Changeset 198 for tests/arx_test.cpp

Show
Ignore:
Timestamp:
11/04/08 14:54:34 (16 years ago)
Author:
smidl
Message:

opravy + zavedeni studenta + zakomentovani debug v mergeru

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • tests/arx_test.cpp

    r170 r198  
    1414         
    1515        //Test constructor 
    16         mat V0 = 0.00001*eye(ord+1); V0(0.0)*= 10; // 
    17         double nu0 = ord+1; 
     16        mat V0 = 0.00001*eye(ord+1); V0(0.0)*= 100; // 
     17        double nu0 = ord+4; 
    1818         
    1919        RV thr("{theta_and_r }",vec_1(ord+1)); 
    … …  
    2222                                 
    2323        //Test estimation 
    24         int ndat = 10000; 
     24        int ndat = 100; 
    2525        int t,j; 
    2626        vec Yt(ndat); 
    … …  
    2828        Yt.set_subvector(0,randn(ord)); //initial values 
    2929        vec rgr(ord); 
     30         
    3031         
    3132        cout << Ar_ep.mean()<<endl; 
    … …  
    3738                Ar.bayes(Psi); 
    3839                LL(t) = Ar._ll(); 
     40                 
     41                cout << "y: " << Yt(t) << endl; 
     42                mlstudent*      Pr = Ar.predictor_student(RV("{y }"),RV("{y1 y2 y3 y4 }")); 
     43                cout << Ar._ll() <<" , " << log(Pr->evalcond(vec_1(Yt(t)),rgr)) <<endl; 
     44                delete Pr; 
    3945        } 
    4046        cout << Ar_ep.mean()<<endl;