Changeset 198 for tests

Timestamp:

11/04/08 14:54:34 (16 years ago)

Author:

smidl

Message:

opravy + zavedeni studenta + zakomentovani debug v mergeru

Location:

tests

Files:

• arx_test.cpp (modified) (4 diffs)
• merger_iter_test.cpp (modified) (3 diffs)
• merger_test.cpp (modified) (3 diffs)

tests/arx_test.cpp

@@ -14 +14 @@
         
         //Test constructor
-        mat V0 = 0.00001*eye(ord+1); V0(0.0)*= 10; //
-        double nu0 = ord+1;
+        mat V0 = 0.00001*eye(ord+1); V0(0.0)*= 100; //
+        double nu0 = ord+4;
         
         RV thr("{theta_and_r }",vec_1(ord+1));
@@ -22 +22 @@
                                 
         //Test estimation
-        int ndat = 10000;
+        int ndat = 100;
         int t,j;
         vec Yt(ndat);
@@ -28 +28 @@
         Yt.set_subvector(0,randn(ord)); //initial values
         vec rgr(ord);
+        
         
         cout << Ar_ep.mean()<<endl;
@@ -37 +38 @@
                 Ar.bayes(Psi);
                 LL(t) = Ar._ll();
+                
+                cout << "y: " << Yt(t) << endl;
+                mlstudent*      Pr = Ar.predictor_student(RV("{y }"),RV("{y1 y2 y3 y4 }"));
+                cout << Ar._ll() <<" , " << log(Pr->evalcond(vec_1(Yt(t)),rgr)) <<endl;
+                delete Pr;
         }
         cout << Ar_ep.mean()<<endl;

tests/merger_iter_test.cpp

@@ -11 +11 @@
 int main() {
 
-        //RNG_randomize();
+        RNG_randomize();
         
         RV x ( "{x }","1" );
@@ -19 +19 @@
         
         enorm<fsqmat> f1 ( xy );
-        enorm<fsqmat> f2 ( x );
+        enorm<fsqmat> f2 ( xy );
+        enorm<fsqmat> f3(y);
                 
-        f1.set_parameters ( "3 2",mat ( "0.5 0; 0 0.5" ) );
-        f2.set_parameters ( "3",mat ( "0.5" ) );
+        f1.set_parameters ( "4 3",mat ( "0.4 0.3; 0.3 0.4" ) );
+        f2.set_parameters ( "1 3",mat ( "0.3 -0.2; -0.2 0.3" ) );
+        f3.set_parameters ( "2",mat("0.4") );
         
-        Array<mpdf* > A ( 2 );
+        Array<mpdf* > A ( 3 );
         mepdf A1(f1);
         mepdf A2(f2);
+        mepdf A3(f3);
         A ( 0 ) =&A1;
         A ( 1 ) =&A2;
+        A ( 2 ) =&A3;
         
         int Npoints=100;
@@ -42 +46 @@
         merger M ( A );
         enorm<fsqmat> g0(xy);
-        g0.set_parameters(vec("1 1"),mat("1 0; 0 1"));
+        g0.set_parameters(vec("4 4"),mat("1 0; 0 1"));
         
-        M.set_parameters(1.2,100,1);
+        M.set_parameters(1.2,400,5);
         M.merge(&g0);
         

tests/merger_test.cpp

@@ -11 +11 @@
 int main() {
 
+        RNG_randomize();
+        
         RV x ( "{x }","1" );
 
@@ -42 +44 @@
         g0.set_parameters(vec("0.0"),mat("100.0"));
         
-        M.set_parameters(1.2,1000,3);
+        M.set_parameters(1.2,200,3);
         M.merge(&g0);
         
@@ -49 +51 @@
         
         it_file it("merger_test.it");
+        it << Name("x_grid") << x_grid;
         it << Name("lf1") << l_f1;
         it << Name("lf2") << l_f2;