Changeset 98 for trunk/src/ROC.cc

Timestamp:

Jun 10, 2004, 5:24:05 PM (19 years ago)

Author:

Peter

Message:

tScore added

File:

• trunk/src/ROC.cc (modified) (3 diffs)

trunk/src/ROC.cc

@@ -2 +2 @@
 
 // System includes
-#include <algorithm>
-#include <utility>
-#include <vector>
+//#include <algorithm>
+//#include <utility>
+//#include <vector>
 
 // Thep C++ Tools
@@ -16 +16 @@
   ROC::ROC(const gslapi::vector& target, const gslapi::vector& value) 
   
-    : value_(), nof_pos_(0), minimum_size_(10)
+    : Score(), value_(), nof_pos_(0), minimum_size_(10), area_(-1)
 
        
@@ -33 +33 @@
 
 
-  double ROC::area()
+  double ROC::score()
   {
     using namespace std;
-    double area=0;
-    for (unsigned int i=0; i<value_.size(); i++)
-      if (value_[i].first==1)
-        area+=i;
-    // Normalizing the area to 0-1
-    area = (area/nof_pos_ - (nof_pos_ - 1)/2 )/(value_.size() - nof_pos_);
-
-    return area;
+    if (area_==-1){
+      double area_=0;
+      for (unsigned int i=0; i<value_.size(); i++)
+        if (value_[i].first==1)
+          area_+=i;
+      // Normalizing the area to 0-1
+      area_ = (area_/nof_pos_ - (nof_pos_ - 1)/2 )/(value_.size() - nof_pos_);
+    }
+    return area_;
   }
 
 
-  double ROC::get_p(const double area)
+  double ROC::p_value(void)
   {
+    if (area_==-1)
+      area_ = score();
     double p;
     if (nof_pos_ < minimum_size_ & value_.size()-nof_pos_ < minimum_size_)
-      p = get_p_exact(area*nof_pos_*(value_.size()-nof_pos_), 
+      p = get_p_exact(area_*nof_pos_*(value_.size()-nof_pos_), 
                           nof_pos_, value_.size()-nof_pos_);
     else
-    p = get_p_approx(area);
+    p = get_p_approx(area_);
     return p;
   }