Changeset 178 for bdm/stat/libBM.h

Timestamp:

10/15/08 19:04:09 (16 years ago)

Author:

smidl

Message:

Changes in basic structures! new methods

Files:

• bdm/stat/libBM.h (modified) (6 diffs)

bdm/stat/libBM.h

@@ -140 +140 @@
 };
 
+class mpdf;
 
 //! Probability density function with numerical statistics, e.g. posterior density.
@@ -173 +174 @@
                 return x;
         }
+        //! Return conditional density on the given RV, the remaining rvs will be in conditioning
+        mpdf* condition(const RV &rv){it_warning("Not implemented"); return NULL;}
+        //! Return marginal density on the given RV, the remainig rvs are intergrated out
+        epdf* marginal(const RV &rv){it_warning("Not implemented"); return NULL;}
 
         //! return expected value
@@ -183 +188 @@
         //! modifier function - useful when copying epdfs
         void _renewrv(const RV &in_rv){rv=in_rv;}
+        //!
 };
 
@@ -213 +219 @@
         };
         //! Update \c ep so that it represents this mpdf conditioned on \c rvc = cond
-        virtual void condition ( const vec &cond ) {};
+        virtual void condition ( const vec &cond ) {it_error("Not implemented");};
 
         //! Shortcut for conditioning and evaluation of the internal epdf. In some cases,  this operation can be implemented efficiently.
@@ -252 +258 @@
                 //! Indeces of rvc in common rv
                 Array<ivec> rvcsinrv;
+                //! Indeces of common rv in rvc
+                Array<ivec> rvinrvcs;
         public:
-                compositepdf(Array<mpdf*> A0): n(A0.length()), mpdfs(A0), rvsinrv(n), rvcsinrv(n){};
+                compositepdf(Array<mpdf*> A0): n(A0.length()), mpdfs(A0), rvsinrv(n), rvcsinrv(n),rvinrvcs(n){};
                 RV getrv(bool checkoverlap=false);
                 void setrvc(const RV &rv, RV &rvc);
@@ -329 +337 @@
         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;}
         
+        //!Constructs a predictive density (marginal density on data)
+        virtual epdf* predictor(const RV &rv){it_error("Not implemented");return NULL;};
+        
         //! Destructor for future use;
         virtual ~BM() {};