Changeset 956 for library/bdm/stat/exp_family.cpp

Timestamp:

05/19/10 03:56:03 (14 years ago)

Author:

sarka

Message:

to_setting

Files:

• library/bdm/stat/exp_family.cpp (modified) (5 diffs)

library/bdm/stat/exp_family.cpp

@@ -326 +326 @@
 }
 
+void egiw::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set );
+                UI::get ( dimx, set, "dimx", UI::compulsory );
+                if ( !UI::get ( nu, set, "nu", UI::optional ) ) {
+                        nu = -1;
+                }
+                mat V;
+                if ( !UI::get ( V, set, "V", UI::optional ) ) {
+                        vec dV;
+                        UI::get ( dV, set, "dV", UI::compulsory );
+                        set_parameters ( dimx, ldmat ( dV ), nu );
+
+                } else {
+                        set_parameters ( dimx, V, nu );
+                }
+        }
+
+void egiw::to_setting ( Setting& set ) const {
+                epdf::to_setting ( set );
+                UI::save ( dimx, set, "dimx" );
+                UI::save ( V.to_mat(), set, "V" );
+                UI::save ( nu, set, "nu" );
+        };
+
+void egiw::validate() {
+                eEF::validate();
+                dim = dimx * ( dimx + nPsi );   
+        
+            // check sizes, rvs etc.
+                // also check if RV are meaningful!!!
+                // meaningful =  rv for theta  and rv for r are split!
+        }
 void multiBM::bayes ( const vec &yt, const vec &cond ) {
         if ( frg < 1.0 ) {
@@ -421 +453 @@
 }
 
+void egamma::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set ); // reads rv
+                UI::get ( alpha, set, "alpha", UI::compulsory );
+                UI::get ( beta, set, "beta", UI::compulsory );
+        }
+        
+void egamma::to_setting ( Setting &set ) const
+        {                       
+                epdf::to_setting( set );
+                UI::save( alpha, set, "alpha" );
+                UI::save( beta, set, "beta" );
+        } 
+        
+        
+void egamma::validate() {
+                eEF::validate();
+                bdm_assert ( alpha.length() == beta.length(), "parameters do not match" );
+                dim = alpha.length();
+        }
+
 void mgamma::set_parameters ( double k0, const vec &beta0 ) {
         k = k0;
         iepdf.set_parameters ( k * ones ( beta0.length() ), beta0 );
 }
+
+
+void mgamma::from_setting ( const Setting &set ) {
+                pdf::from_setting ( set ); // reads rv and rvc
+                vec betatmp; // ugly but necessary
+                UI::get ( betatmp, set, "beta", UI::compulsory );
+                UI::get ( k, set, "k", UI::compulsory );
+                set_parameters ( k, betatmp );
+        }
+
+void mgamma::to_setting  (Setting  &set) const {
+                pdf::to_setting(set);
+                UI::save( _beta, set, "beta");
+                UI::save( k, set, "k");
+
+        }
+
+void mgamma::validate() {
+                pdf_internal<egamma>::validate();
+
+                dim = _beta.length();
+                dimc = _beta.length();
+        }
 
 void eEmp::resample ( RESAMPLING_METHOD method ) {
@@ -536 +611 @@
 
 void migamma_ref::from_setting ( const Setting &set ) {
+        migamma::from_setting(set);
         vec ref;
         UI::get ( ref, set, "ref" , UI::compulsory );
@@ -541 +617 @@
         validate();
 }
-
+void  migamma_ref::to_setting  (Setting  &set) const {
+        migamma::to_setting(set);
+        UI::save ( pow ( refl, 1/(1.0 - l) ), set, "ref");      
+        UI::save(l,set,"l");
+        UI::save(k,set,"k");
+}
 void mlognorm::from_setting ( const Setting &set ) {
+        pdf_internal<elognorm>::from_setting(set);
         vec mu0;
         UI::get ( mu0, set, "mu0", UI::compulsory );
         set_parameters ( mu0.length(), set["k"] );
-        validate();
         condition ( mu0 );
 }
+
+
+void mlognorm::to_setting  (Setting  &set) const {
+        pdf_internal<elognorm>::to_setting(set);
+        UI::save ( exp(mu + sig2), set, "mu0");
+}
+
 
 void mlstudent::condition ( const vec &cond ) {
@@ -586 +674 @@
 }
 
+void eEmp::to_setting ( Setting &set ) const {
+                epdf::to_setting( set );
+                UI::save ( samples, set, "samples" );
+                UI::save ( w, set, "w" );
+        }
+
+void eEmp::from_setting ( const Setting &set ) {
+                epdf::from_setting( set );
+                
+                UI::get( samples, set, "samples", UI::compulsory );
+                UI::get ( w, set, "w", UI::compulsory );
+        }
+
+void    eEmp::validate (){
+          epdf::validate();
+          bdm_assert (samples.length()==w.length(),"samples and weigths are of different lengths");
+          n = w.length();
+          if (n>0)
+                pdf::dim = samples ( 0 ).length();
+        }
+        
+        void eDirich::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set );
+                UI::get ( beta, set, "beta", UI::compulsory );
+        }
+void eDirich::validate() {
+                //check rv
+                eEF::validate();
+                dim = beta.length();
+        }
+
+void eDirich::to_setting ( Setting &set ) const
+        {                       
+                eEF::to_setting( set );
+                UI::save( beta, set, "beta" );
+        } 
+
+void euni::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set ); // reads rv and rvc
+
+                UI::get ( high, set, "high", UI::compulsory );
+                UI::get ( low, set, "low", UI::compulsory );
+                set_parameters ( low, high );
+                
+        }
+        
+void    euni::to_setting  (Setting  &set) const {
+                epdf::to_setting ( set );
+                UI::save ( high, set, "high" );
+                UI::save ( low, set, "low" );
+        }
+        
+void euni::validate() {
+                epdf::validate();
+                bdm_assert ( high.length() == low.length(), "Incompatible high and low vectors" );
+                dim = high.length();
+                bdm_assert ( min ( distance ) > 0.0, "bad support" );
+        }
+
+void mgdirac::from_setting(const Setting& set){
+                        pdf::from_setting(set);
+                        g=UI::build<fnc>(set,"g",UI::compulsory);
+                        validate();
+                }
+void mgdirac::to_setting(Setting &set) const{
+                        pdf::to_setting(set);
+                        UI::save(g.get(), set, "g");
+                }
+void mgdirac::validate() {
+                        pdf::validate();
+                        dim = g->dimension();
+                        dimc = g->dimensionc();
+                }
+
+void mDirich::from_setting ( const Setting &set ) {
+                pdf::from_setting ( set ); // reads rv and rvc
+                if ( _rv()._dsize() > 0 ) {
+                        rvc = _rv().copy_t ( -1 );
+                }
+                vec beta0;
+                if ( !UI::get ( beta0, set, "beta0", UI::optional ) ) {
+                        beta0 = ones ( _rv()._dsize() );
+                }
+                if ( !UI::get ( betac, set, "betac", UI::optional ) ) {
+                        betac = 0.1 * ones ( _rv()._dsize() );
+                }
+                _beta = beta0;
+
+                UI::get ( k, set, "k", UI::compulsory );
+        }
+
+void mDirich::to_setting  (Setting  &set) const {
+                pdf::to_setting(set);
+                UI::save( _beta, set, "beta0");
+                UI::save( betac, set, "betac");
+                UI::save ( k, set, "k" );
+        }
+
+
+void mDirich::validate() {
+                pdf_internal<eDirich>::validate();
+                bdm_assert ( _beta.length() == betac.length(), "beta0 and betac are not compatible" );
+                if ( _rv()._dsize() > 0 ) {
+                        bdm_assert ( ( _rv()._dsize() == dimension() ) , "Size of rv does not match with beta" );
+                }
+                dimc = _beta.length();
+        }
 
 };