Changeset 78 for trunk/src/ROC.cc

Timestamp:

May 4, 2004, 1:26:48 PM (19 years ago)

Author:

Peter

Message:

modified get_p_approx

File:

• trunk/src/ROC.cc (modified) (1 diff)

trunk/src/ROC.cc

@@ -62 +62 @@
   {
     double x = area - 0.5;
-    double sigma = (pow((value_.size()-nof_pos_)* nof_pos_ *
-                        (value_.size()+1)/12, 0.5) /
+    // Not integrating from the middle of the bin, but from the inner edge.
+    if (x>0)
+      x -= 0.5/nof_pos_/(value_.size()-nof_pos_);
+    else if(x<0)
+      x += 0.5/nof_pos_/(value_.size()-nof_pos_);
+
+    double sigma = (std::sqrt((value_.size()-nof_pos_)* nof_pos_ *
+                        (value_.size()+1)/12) /
                     (value_.size() - nof_pos_ ) / nof_pos_);
     double p = gsl_cdf_gaussian_Q(x, sigma);