Changeset 176 for bdm/stat/libEF.h

Timestamp:

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

Author:

smidl

Message:

Corrections to mixtures & merger

Files:

• bdm/stat/libEF.h (modified) (5 diffs)

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
 };