Changeset 182 for bdm/stat/emix.cpp
- Timestamp:
- 10/17/08 10:28:15 (16 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
bdm/stat/emix.cpp
r181 r182 31 31 } 32 32 33 // mprod::mprod ( Array<mpdf*> mFacs, bool overlap) : mpdf ( RV(), RV() ), n ( mFacs.length() ), epdfs ( n ), mpdfs ( mFacs ), rvinds ( n ), rvcinrv ( n ), rvcsinrvc ( n ) { 33 emix* emix::marginal(const RV &rv) const{ 34 Array<epdf*> Cn(Coms.length()); 35 for(int i=0;i<Coms.length();i++){Cn(i)=Coms(i)->marginal(rv);} 36 emix* tmp = new emix(rv); 37 tmp->set_parameters(w,Coms,false); 38 tmp->ownComs(); 39 return tmp; 40 } 41 42 mratio* emix::condition(const RV &rv) const{ 43 return new mratio(this,rv); 44 }; 45 46 // mprod::mprod ( Array<mpdf*> mFacs, bool overlap) : mpdf ( RV(), RV() ), n ( mFacs.length() ), epdfs ( n ), mpdfs ( mFacs ), rvinds ( n ), rvcinrv ( n ), irvcs_rvc ( n ) { 34 47 // int i; 35 48 // bool rvaddok; … … 54 67 // rvcinrv ( i ) = mpdfs ( i )->_rvc().dataind ( rv ); 55 68 // // find ith rvc in common rv 56 // rvcsinrvc ( i ) = mpdfs ( i )->_rvc().dataind ( rvc );69 // irvcs_rvc ( i ) = mpdfs ( i )->_rvc().dataind ( rvc ); 57 70 // // 58 71 // /* if ( rvcinrv ( i ).length() >0 ) {independent = false;} 59 // if ( rvcsinrvc ( i ).length() >0 ) {independent = false;}*/72 // if ( irvcs_rvc ( i ).length() >0 ) {independent = false;}*/ 60 73 // } 61 74 // };