Changeset 78

Timestamp:

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

Author:

Peter

Message:

modified get_p_approx

Location:

trunk/src

Files:

• ROC.cc (modified) (1 diff)
• ROC.h (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);

trunk/src/ROC.h

@@ -61 +61 @@
     ///
     ///
-    /// Implemented as in MatLab 12.1
+    /// Implemented as in MatLab 13.1
     /// @return the p-value
     ///