Changeset 1553 for trunk/yat/statistics

Timestamp:

Oct 7, 2008, 6:06:50 PM (13 years ago)

Author:

Peter

Message:

fixes #450

Location:

trunk/yat/statistics

Files:

• utility.cc (modified) (1 diff)
• utility.h (modified) (1 diff)

trunk/yat/statistics/utility.cc

@@ -47 +47 @@
     if (n<=2)
       return std::numeric_limits<double>::quiet_NaN();
+    if (r>=1.0)
+      return 0.0;
+    if (r<=-1.0)
+      return 1.0;
     return gsl_cdf_tdist_Q(r*sqrt((n-2)/(1-r*r)), n-2);
   }

trunk/yat/statistics/utility.h

@@ -100 +100 @@
 
      \return Probability that correlation is larger than \a r by
-     chance when having \a n samples.
+     chance when having \a n samples. For \a r larger or equal to 1.0,
+     0.0 is returned. For \a r smaller or equal to -1.0, 1.0 is
+     returned.
    */
   double pearson_p_value(double r, unsigned int n);