Changeset 878 for library/bdm/stat/exp_family.cpp
- Timestamp:
- 03/24/10 14:30:43 (14 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/stat/exp_family.cpp
r850 r878 21 21 dimx = dimx0; 22 22 nPsi = V0.rows() - dimx; 23 dim = dimx * ( dimx + nPsi ); // size(R) + size(Theta)24 23 25 24 V = V0; … … 53 52 eWishartCh Omega; //inverse Wishart, result is R, 54 53 Omega.set_parameters ( iChLam, nu - 2*nPsi - dimx ); // 2*nPsi is there to match numercial simulations - check if analytically correct 54 Omega.validate(); 55 55 56 56 mat OmChi; … … 87 87 eWishartCh Omega; //inverse Wishart, result is R, 88 88 Omega.set_parameters ( iCh, nu - 2*nPsi - dimx ); // 2*nPsi is there to match numercial simulations - check if analytically correct 89 Omega.validate(); 89 90 90 91 chmat Omi; … … 420 421 k = k0; 421 422 iepdf.set_parameters ( k * ones ( beta0.length() ), beta0 ); 422 dimc = iepdf.dimension();423 dim = iepdf.dimension();424 423 } 425 424 … … 532 531 UI::get ( ref, set, "ref" , UI::compulsory ); 533 532 set_parameters ( set["k"], ref, set["l"] ); 533 validate(); 534 534 } 535 535 … … 538 538 UI::get ( mu0, set, "mu0", UI::compulsory ); 539 539 set_parameters ( mu0.length(), set["k"] ); 540 validate(); 540 541 condition ( mu0 ); 541 542 }