Changeset 310 for mpdm

Timestamp:

04/09/09 14:33:08 (15 years ago)

Author:

smidl

Message:

Published example 2d

Location:

mpdm/TR2244

Files:

• merger_iter_cond.cpp (modified) (4 diffs)
• merger_iter_debug.m (added)

mpdm/TR2244/merger_iter_cond.cpp

@@ -25 +25 @@
 
         //Differneces in constant term are essential
-        f1.set_parameters ( "0.4","1",mat ( "0.04" ) );
-        f2.set_parameters ( "0.2", mat ( "0.08" ) );
+        f1.set_parameters ( "1","1",mat ( "0.01" ) );
+        f2.set_parameters ( "2", mat ( "0.02" ) );
 
         mepdf mf2(&f2);
@@ -40 +40 @@
 /*      Cov(1,2)=0.29;
         Cov(2,1)=0.29;*/
-        g0.set_parameters ( vec ( "0 0" ),Cov);// +1*ones ( 3,3 ) );
+        g0.set_parameters ( vec ( "2 1" ),Cov);// +1*ones ( 3,3 ) );
 
 //      enorm<ldmat>* teste=g0.marginal(concat(u1,u2));
 //      mlnorm<ldmat>* testm=(mlnorm<ldmat>*)teste->condition(u2);
         
-        M.set_parameters ( 1.2,1000,1);
+        M.set_parameters ( 1.2,1000,1,0.5);
         
         Array<vec> YUU(2); 
-        YUU(0)=linspace(-2.9,3.1,20);
-        YUU(1)=linspace(-5,5,20);
-        //YUU(2)=linspace(-3,3,20);
+        YUU(0)=linspace(1.5,2.5,20);
+        YUU(1)=linspace(0.5,1.5,20);
+        
+        YUU(0)=linspace(0,4,20);
+        YUU(1)=linspace(-1,3,20);
+//YUU(2)=linspace(-3,3,20);
 
-        M.set_grid(YUU);
+        //M.set_grid(YUU);
         
         int Ntrials=1;
@@ -68 +71 @@
                 epdf* MP = MM._Coms ( 0 )->epredictor ( );
                 MP->set_rv(all);
+                
+                //Missing piece
+                enorm<ldmat>* miss=     (enorm<ldmat>*) MP->marginal(u1);
+                cout << "Missing: " << miss->_mu() << " +- " << miss->_R().to_mat() << endl;
 
                 RV yu1 = y; yu1.add ( u1 );
@@ -102 +109 @@
         cout << "R1: " << R1mean << " +- " << 2*sqrt(sum_sqr(R1s)/Ntrials-R1mean*R1mean) <<endl;
         cout << "R2: " << R2mean << " +- " << 2*sqrt(sum_sqr(R2s)/Ntrials-R2mean*R2mean) <<endl;
+        
+        
 }