Changeset 1624 for branches/0.4-stable/yat/statistics/utility.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/yat/statistics/utility.cc (modified) (2 diffs)

branches/0.4-stable/yat/statistics/utility.cc

@@ -31 +31 @@
 #include <gsl/gsl_statistics_double.h>
 
+#include <algorithm>
 #include <cassert>
 
@@ -41 +42 @@
   {
     double p=0;
-    for (size_t i=0; i<=k; i++)
+    size_t top = std::min(k, std::min(n1, t));
+    for (size_t i=std::max(0, static_cast<int>(t-n2)); i<=top; i++)
       p+= gsl_ran_hypergeometric_pdf(i, n1, n2, t);
     return p;