Changeset 176 for bdm/stat

Timestamp:

10/09/08 11:32:12 (16 years ago)

Author:

smidl

Message:

Corrections to mixtures & merger

Location:

bdm/stat

Files:

• emix.h (modified) (1 diff)
• libBM.cpp (modified) (4 diffs)
• libBM.h (modified) (5 diffs)
• libEF.h (modified) (5 diffs)
• merger.cpp (deleted)
• merger.h (deleted)

bdm/stat/emix.h

@@ -78 +78 @@
         Additional parameter overlap is left for future use. Do not set to true for mprod.
         */
-        mprod ( Array<mpdf*> mFacs): compositepdf(mFacs), mpdf(getrv(true),RV()), rvcinds(n), epdfs(n)
+        mprod ( Array<mpdf*> mFacs): compositepdf(mFacs), mpdf(getrv(true),RV()), epdfs(n), rvcinds(n)
         {       setrvc(rv,rvc);
                 setrvcinrv(rvc,rvcinds);

bdm/stat/libBM.cpp

@@ -91 +91 @@
 // }
 
-RV RV::subselect ( ivec ind ) const {
+RV RV::subselect (const ivec &ind ) const {
         RV ret;
         ret.init ( ids ( ind ), names ( to_Arr ( ind ) ), sizes ( ind ), times ( ind ) );
@@ -99 +99 @@
 void RV::t ( int delta ) { times += delta;}
 
-RV RV::operator() ( ivec ind ) const {
+RV RV::operator() (const ivec &ind ) const {
         RV ret;
         if ( ind.length() >0 ) {
@@ -140 +140 @@
 }
 
-ivec RV::dataind ( RV rv2 ) const {
+ivec RV::dataind (const RV &rv2 ) const {
         ivec res ( 0 );
         if ( rv2.count()>0 ) {
@@ -155 +155 @@
 }
 
-RV RV::subt ( const RV rv2 ) const {
+RV RV::subt ( const RV &rv2 ) const {
         ivec res = this->findself ( rv2 ); // nonzeros
         ivec valid = itpp::find ( res == -1 ); //-1 => value not found => it remains

bdm/stat/libBM.h

@@ -80 +80 @@
         bool add ( const RV &rv2 );
         //! Subtract  another variable from the current one
-        RV subt ( const RV rv2 ) const;
+        RV subt ( const RV &rv2 ) const;
         //! Select only variables at indeces ind
-        RV subselect ( ivec ind ) const;
+        RV subselect ( const ivec &ind ) const;
         //! Select only variables at indeces ind
-        RV operator() ( ivec ind ) const;
+        RV operator() ( const ivec &ind ) const;
         //! Shift \c time shifted by delta.
         void t ( int delta );
@@ -90 +90 @@
         str tostr() const;
         //! generate indeces into \param crv data vector that form data vector of self.
-        ivec dataind(RV crv) const;
+        ivec dataind(const RV &crv) const;
 
         //!access function
@@ -103 +103 @@
         //!access function
         std::string name ( int at ) {return names ( at );};
+        
+        //!access function
+        void set_id ( int at, int id0 ) {ids ( at )=id0;};
+        //!access function
+        void set_size ( int at, int size0 ) {sizes ( at )=size0; tsize=sum(sizes);};
+        //!access function
+        void set_time ( int at, int time0 ) {times ( at )=time0;};
+
         //!Assign unused ids to this rv 
         void newids();
@@ -160 +168 @@
 
         //! Compute log-probability of multiple values argument \c val
-        virtual vec evalpdflog ( const mat &Val ) const {
+        virtual vec evalpdflog_m ( const mat &Val ) const {
                 vec x ( Val.cols() );
                 for ( int i=0;i<Val.cols();i++ ) {x ( i ) =evalpdflog( Val.get_col(i) ) ;}
@@ -318 +326 @@
         //! I.e. marginal likelihood of the data with the posterior integrated out.
         virtual double logpred(const vec &dt)const{it_error("Not implemented");return 0.0;}
+        //! Matrix version of logpred
+        vec logpred_m(const mat &dt)const{vec tmp(dt.cols());for(int i=0;i<dt.cols();i++){tmp(i)=logpred(dt.get_col(i));}return tmp;}
         
         //! Destructor for future use;

bdm/stat/libEF.h

@@ -91 +91 @@
         void bayes ( const vec &dt );
         //!Flatten the posterior
-        virtual void flatten ( BMEF * B) {it_error ( "Not implemented" );}
+        virtual void flatten ( BMEF * B ) {it_error ( "Not implemented" );}
 };
 
@@ -206 +206 @@
         vec sample() const {it_error ( "Not implemented" );return vec_1 ( 0.0 );};
         vec mean() const {return beta/sum ( beta );};
-        //! In this instance, val= [theta, r]. For multivariate instances, it is stored columnwise val = [theta_1 theta_2 ... r_1 r_2 ]
+        //! In this instance, val is ...
         double evalpdflog_nn ( const vec &val ) const {return ( beta-1 ) *log ( val );};
         double lognc () const {
@@ -216 +216 @@
         //!access function
         vec& _beta() {return beta;}
+        //!Set internal parameters
+        void set_parameters(const vec &beta0){
+                if(beta0.length()!=beta.length()){
+                        it_assert_debug(rv.length()==1,"Undefined");
+                        rv.set_size(0,beta0.length());
+                }
+                beta= beta0;
+        }
 };
 
@@ -239 +247 @@
                 eDirich pred ( est );
                 vec &beta = pred._beta();
-                
+
                 double lll;
                 if ( frg<1.0 )
@@ -250 +258 @@
                 return pred.lognc()-lll;
         }
-        void flatten (BMEF* B ) {
-                eDirich* E=dynamic_cast<eDirich*>(B);
+        void flatten ( BMEF* B ) {
+                eDirich* E=dynamic_cast<eDirich*> ( B );
                 // sum(beta) should be equal to sum(B.beta)
                 const vec &Eb=E->_beta();
-                est.pow ( sum(beta)/sum(Eb) );
+                est.pow ( sum ( beta ) /sum ( Eb ) );
+                if ( evalll ) {last_lognc=est.lognc();}
+        }
+        const epdf& _epdf() const {return est;};
+        void set_parameters ( const vec &beta0 ) {
+                est.set_parameters(beta0);
+                rv = est._rv();
                 if(evalll){last_lognc=est.lognc();}
         }
-        const epdf& _epdf() const {return est;};
-        //!access funct
 };