Changeset 200 for bdm/stat/libEF.h

Timestamp:

11/10/08 15:40:24 (16 years ago)

Author:

smidl

Message:

BM has now function _e() returning posterior of correct type

Files:

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

bdm/stat/libEF.h

@@ -167 +167 @@
         int nPsi;
 public:
-        //!Default constructor, assuming
-        egiw ( RV rv, mat V0, double nu0 ) : eEF ( rv ), V ( V0 ), nu ( nu0 ) {
+        //!Default constructor, if nu0<0 a minimal nu0 will be computed
+        egiw ( RV rv, mat V0, double nu0=-1.0 ) : eEF ( rv ), V ( V0 ), nu ( nu0 ) {
                 xdim = rv.count() /V.rows();
                 it_assert_debug ( rv.count() ==xdim*V.rows(),"Incompatible V0." );
                 nPsi = V.rows()-xdim;
+                //set mu to have proper normalization and 
+                if (nu0<0){
+                        nu = 0.1 +nPsi +2*xdim +2; // +2 assures finite expected value of R
+                        // terms before that are sufficient for finite normalization
+                }
         }
         //!Full constructor for V in ldmat form
-        egiw ( RV rv, ldmat V0, double nu0 ) : eEF ( rv ), V ( V0 ), nu ( nu0 ) {
+        egiw ( RV rv, ldmat V0, double nu0=-1.0 ) : eEF ( rv ), V ( V0 ), nu ( nu0 ) {
                 xdim = rv.count() /V.rows();
                 it_assert_debug ( rv.count() ==xdim*V.rows(),"Incompatible V0." );
                 nPsi = V.rows()-xdim;
+                if (nu0<0){
+                        nu = 0.1 +nPsi +2*xdim +2; // +2 assures finite expected value of R
+                        // terms before that are sufficient for finite normalization
+                }
         }
 
@@ -275 +284 @@
         }
         const epdf& _epdf() const {return est;};
+        const eDirich* _e() const {return &est;};
         void set_parameters ( const vec &beta0 ) {
                 est.set_parameters ( beta0 );