Changeset 1626 for trunk

Timestamp:

Nov 15, 2008, 5:01:46 PM (12 years ago)

Author:

Peter

Message:

closes #457 - kolmogorov_smirnov_test fails on Mac OS 10.4

Location:

trunk

Files:

• test/kolmogorov_smirnov_test.cc (modified) (2 diffs)
• yat/statistics/KolmogorovSmirnov.cc (modified) (3 diffs)

trunk/test/kolmogorov_smirnov_test.cc

@@ -42 +42 @@
   test_one_sample(suite);
   test_two_sample(suite);
-  //test_p_value(suite);
+  test_p_value(suite);
   test_reset(suite);
 
@@ -116 +116 @@
     ks.add(i+14.5, false);
   } 
+  suite.add(suite.equal(ks.score(), 0.15, 10));
+
   statistics::Averager a;
   for (size_t n=0; n<100; ++n) {

trunk/yat/statistics/KolmogorovSmirnov.cc

@@ -29 +29 @@
 #include <deque>
 #include <functional>
+#include <limits>
 
 namespace theplu {
@@ -84 +85 @@
       targets.push_back(i->second);
 
+    double score_threshold = score()-10*std::numeric_limits<double>().epsilon();
+
     for (size_t i=0; i<perm; ++i){
       random::random_shuffle(targets.begin(), targets.end()); 
@@ -92 +95 @@
         ks.add(iter->first.first, *target_i, iter->first.second);
       
-      if (ks.score()>=score())
+      if (ks.score()>=score_threshold)
         ++count;
     }