work/mixpp/bdm/stat/libBM.h

Go to the documentation of this file.

#ifndef BM_H
#define BM_H

#include <itpp/itbase.h>
//#include <std>

using namespace itpp;

class RV {
        int size;
        int len;
        ivec ids;
        ivec sizes;
        ivec times;
        ivec obs;
        Array<std::string> names;

private:
        void init ( ivec in_ids, Array<std::string> in_names, ivec in_sizes, ivec in_times, ivec in_obs );
public:
        RV ( ivec in_ids, Array<std::string> in_names, ivec in_sizes, ivec in_times, ivec in_obs );
        RV ( ivec ids );
        RV ();
        
        friend std::ostream &operator<< ( std::ostream &os, const RV &rv );

        int count() const {return size;} ;
        //TODO why not inline and later??
        
        ivec find(RV rv2);
        RV add(RV rv2);
        RV subt(RV rv2);
        RV subselect(ivec ind);
        RV operator()(ivec ind);
        void t(int delta);
        ivec indexlist();
};




class fnc {
protected:
        int dimy;
public: 
        virtual vec eval(const vec &cond){}; //Fixme: virtual?
        int _dimy()const{return dimy;}
};

class BM {
public:
        double ll;
        bool evalll;

        BM():ll(0),evall(true){};
        
        virtual void bayes ( const vec &dt) = 0;
        void bayes ( mat Dt );
        epdf* _epdf();
};

class epdf {
        RV rv;
public:
//      virtual vec moment ( const int order = 1 );
        virtual vec sample ()=0;
        virtual double eval(const vec &val){};
};

class mpdf {
        RV rv;
        RV rvc;
public:

//      virtual fnc moment ( const int order = 1 );
        virtual vec samplecond (vec &cond, double lik){};
        virtual void condition (vec &cond){};
};

class DS {
protected:
        RV Drv; 
        RV Urv; //
public:
        void getdata(vec &dt);
        void getdata(vec &dt, ivec &indeces);
        void write(vec &ut);
        void write(vec &ut, ivec &indeces);
        void linkrvs(RV &drv, RV &urv);
        
        void step();
};


#endif // BM_H

---