root/library/tests/testsuite/egiw_test.cpp @ 1383

Revision 1064, 1.3 kB (checked in by mido, 14 years ago)

astyle applied all over the library

  • Property svn:eol-style set to native
Line 
1#include "base/bdmbase.h"
2#include "base/user_info.h"
3#include "stat/exp_family.h"
4#include "stat/discrete.h"
5#include "itpp_ext.h"
6#include "../epdf_harness.h"
7#include "../mat_checks.h"
8#include "UnitTest++.h"
9
10const double epsilon = 0.00001;
11
12using namespace bdm;
13
14TEST ( egiw_test ) {
15    epdf_harness::test_config ( "egiw.cfg" );
16}
17
18TEST ( egiw_1_2_test ) {
19    // Setup model
20    double mu = 1.1; //unit step parametr
21    double b = 3.0; // sequence of <1 -1 1 -1...>
22    double s = 0.1;
23
24
25    // TEST 1x1 EGIW
26    mat V ( 3, 3 );
27    V ( 0, 0 ) = pow ( mu, 2 ) + pow ( b, 2 ) + s;
28    V ( 1, 0 ) = mu;
29    V ( 2, 0 ) = b;
30
31    V ( 0, 1 ) = V ( 1, 0 );
32    V ( 1, 1 ) = 1.0;
33    V ( 2, 1 ) = 0.0;
34
35    V ( 0, 2 ) = V ( 2, 0 );
36    V ( 1, 2 ) = V ( 2, 1 );
37    V ( 2, 2 ) = 1.0;
38
39    double nu = 20;
40
41    egiw E ( 1, nu * V, nu );
42    CHECK_CLOSE ( vec ( "1.1 3.0 0.142857" ), E.mean(), epsilon );
43    CHECK_CLOSE ( 7.36731, E.lognc(), epsilon );
44
45    int n = 100;
46    rectangular_support Sup;
47    Sup.set_parameters ( "{[-2.,4.], [1.,5.], [0.,2.]} ", n*ones_i ( E.dimension() ) );
48
49    double summ = exp ( E.evallog ( Sup.first_vec() ) );
50    // sum all likelihood at all points of support
51    for ( int k = 1; k < Sup.points(); k++ ) { // ALL b
52        summ += exp ( E.evallog ( Sup.next_vec() ) );
53    }
54
55    CHECK_CLOSE ( 1.0, summ*prod ( Sup._steps() ), 0.01 );
56}
Note: See TracBrowser for help on using the browser.