work/git/mixpp/bdm/estim/mixef.h

Go to the documentation of this file.

#ifndef MEF_H
#define MEF_H

#include <itpp/itbase.h>
#include "../stat/libFN.h"
#include "../stat/libEF.h"
#include "../stat/emix.h"

using namespace itpp;

class MixEF: public BM {
protected:
        int n;
        Array<BMEF*> Coms;
        multiBM weights;
        eprod* est;
        
        void build_est() {
                Array<const epdf*> epdfs ( n+1 );
                for ( int i=0;i<Coms.length();i++ ) {
//                      it_assert_debug(!x,"MixEF::MixEF : Incompatible components");
                        epdfs ( i ) =& ( Coms ( i )->_epdf() );
                }
                // last in the product is the weight
                epdfs ( n ) =& ( weights._epdf() );
                est = new eprod ( epdfs );
        }

public:
        MixEF ( const Array<BMEF*> &Coms0, const vec &alpha0 ) :
                        BM ( RV() ), n ( Coms0.length() ), Coms ( n ),
                        weights ( RV ( "{w }", vec_1 ( n ) ),alpha0 ) {
        //      it_assert_debug ( n>0,"MixEF::MixEF : Empty Component list" );

                for ( int i=0;i<n;i++ ) {Coms ( i ) = ( BMEF* ) Coms0 ( i )->_copy_();}
                build_est();
        };
        MixEF () :
                        BM ( RV() ), n ( 0 ), Coms ( 0 ),
                        weights ( RV ( "{w }", vec_1 ( 0 ) ),vec ( 0 ) ) {build_est();}
        void init ( BMEF* Com0, const mat &Data, int c=5 );
        //Destructor
        ~MixEF() {
                delete est;
                for ( int i=0;i<n;i++ ) {delete Coms ( i );}
        }
        void bayes ( const vec &dt );
        void bayes ( const mat &dt );
        void bayesB ( const mat &dt );
        double logpred ( const vec &dt ) const;
        const epdf& _epdf() const {return *est;}
        emix* predictor(const RV &rv);
        void flatten(double sumw=1.0);
};


#endif // MEF_H

---