work/mixpp/bdm/stat/emix.h

Go to the documentation of this file.

#ifndef MX_H
#define MX_H

#include "libBM.h"
#include "libEF.h"
//#include <std>

using namespace itpp;

class emix : public epdf {
protected:
        vec w;
        Array<epdf*> Coms;
public:
        emix ( RV &rv ): epdf(rv){};
        void set_parameters ( const vec &w, const Array<epdf*> &Coms );

        vec sample() const;
        vec mean() const {
                int i; vec mu=zeros ( rv.count() );
                for ( i=0;i<w.length();i++ ) {mu+=w ( i ) *Coms ( i )->mean(); }
                return mu;
        }
        double evalpdflog ( const vec &val ) const {int i; double sum=0.0; for ( i=0;i<w.length();i++ ) {sum+=w ( i ) *Coms ( i )->evalpdflog ( val );} return log ( sum );};

//Access methods
        vec& _w() {return w;}
};


#endif //MX_H

---