Changeset 1624 for branches/0.4-stable/test/fisher_test.cc

Timestamp:

Nov 12, 2008, 11:10:52 PM (15 years ago)

Author:

Peter

Message:

fixes #461. Also modified implementation of cdf_hypergeometric_P, which may cause conflict with modifications done in trunk (refs #87). If so, go with the trunk version (which uses GSL 1.8).

File:

• branches/0.4-stable/test/fisher_test.cc (modified) (2 diffs)

branches/0.4-stable/test/fisher_test.cc

@@ -24 +24 @@
 #include "yat/statistics/Fisher.h"
 
+using namespace theplu::yat;
+void test_p_value(test::Suite&);
+void test_p_value_approximative(test::Suite&);
+void test_p_value_exact(test::Suite&);
+
 int main(int argc, char* argv[])
 {  
-  using namespace theplu::yat;
   test::Suite suite(argc, argv);
 
@@ -60 +64 @@
     suite.err() << "minimum_size failed\n";
   }
+  test_p_value(suite);
   return suite.return_value();
 }
+
+
+void test_p_value(test::Suite& suite)
+{
+  test_p_value_exact(suite);
+  test_p_value_approximative(suite);
+}
+
+
+void test_p_value_approximative(test::Suite& suite)
+{
+  suite.err() << "testing p_value_approximative\n";
+  statistics::Fisher f;
+  f.minimum_size() = 0;
+  f.oddsratio(10,20,20,50);
+  suite.add(suite.equal(f.p_value(), 2*f.p_value_one_sided()));
+  f.oddsratio(10,20,10,20);
+  suite.add(suite.equal(f.p_value(), 1.0));
+  suite.add(suite.equal(f.p_value_one_sided(), 0.5));
+}
+
+
+void test_p_value_exact(test::Suite& suite)
+{
+  suite.err() << "test p_value_exact\n";
+  statistics::Fisher f;
+  f.minimum_size() = 1000;
+  f.oddsratio(10,20,10,20);
+  suite.add(suite.equal(f.p_value(), 1.0));
+
+  f.oddsratio(10, 20, 20, 200);
+  suite.add(suite.equal(f.p_value(), 0.0008119060627676223));
+  suite.add(suite.equal(f.p_value_one_sided(), 0.0008119060627676223));
+
+  // testing symmetry
+  statistics::Fisher f2;
+  f2.minimum_size() = 1000;
+  f2.oddsratio(20, 200, 10, 20);
+  suite.add(suite.equal(f2.p_value(), f.p_value()));
+
+  f.oddsratio(1, 1, 1, 2);
+  suite.add(suite.equal(f.p_value(), 1.0));
+  suite.add(suite.equal(f.p_value_one_sided(), 0.7));
+
+  f.oddsratio(1, 1, 2, 1);
+  suite.add(suite.equal(f.p_value(), 1.0));
+  suite.add(suite.equal(f.p_value_one_sided(), 0.9));
+
+}