Changeset 1287 for trunk/yat

Timestamp:

Apr 24, 2008, 11:24:12 PM (15 years ago)

Author:

Peter

Message:

Fixes #87 - replace cumulative hypergeometric with GSL implementation
NOTE: This means that yat 0.5 will require GSL 1.8 or
later. configure.ac is updated to check for that.

File:

• trunk/yat/statistics/utility.cc (modified) (1 diff)

trunk/yat/statistics/utility.cc

@@ -40 +40 @@
                               unsigned int n2, unsigned int t)
   {
-    double p=0;
-    for (size_t i=0; i<=k; i++)
-      p+= gsl_ran_hypergeometric_pdf(i, n1, n2, t);
-    return p;
+    return gsl_cdf_hypergeometric_P(k, n1, n2, t);
   }