Changeset 3236 for trunk/test/statistics.cc

Timestamp:

May 23, 2014, 3:42:51 PM (9 years ago)

Author:

Peter

Message:

BH correction on unsorted range. closes #796

File:

• trunk/test/statistics.cc (modified) (3 diffs)

trunk/test/statistics.cc

@@ -48 +48 @@
 using namespace theplu::yat;
 void test_benjamini_hochberg(test::Suite&);
+void test_benjamini_hochberg_unsorted(test::Suite&);
 void test_entropy(test::Suite&);
 void test_mad(test::Suite&);
@@ -135 +136 @@
   }
   test_benjamini_hochberg(suite);
+  test_benjamini_hochberg_unsorted(suite);
   test_entropy(suite);
   test_median_empty(suite);
@@ -168 +170 @@
 
 
+void test_benjamini_hochberg_unsorted(test::Suite& suite)
+{
+  std::vector<double> p;
+  p.push_back(0.015);
+  p.push_back(0.0001);
+  p.push_back(0.01);
+  p.push_back(0.5);
+  p.push_back(0.99);
+  std::vector<double> q(p.size());
+  statistics::benjamini_hochberg_unsorted(p.begin(), p.end(), q.begin());
+  suite.add(suite.equal(q[1], p[1]*5));
+  suite.add(suite.equal(q[2], p[2]*2.5));
+  suite.add(suite.equal(q[0], 0.025));
+  suite.add(suite.equal(q[3], p[3]*1.25));
+  suite.add(suite.equal(q[4], 0.99));
+
+  // do nut run compiler test
+  if (false) {
+    using statistics::benjamini_hochberg_unsorted;
+    boost::random_access_iterator_archetype<double> input;
+    boost::mutable_random_access_iterator_archetype<double> result;
+    benjamini_hochberg_unsorted(input, input, result);
+  }
+}
+
+
 void test_entropy(test::Suite& suite)
 {