Changeset 103 for trunk/src/ROC.cc

Timestamp:

Jun 15, 2004, 4:29:15 PM (19 years ago)

Author:

Peter

Message:

output operator added

File:

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

trunk/src/ROC.cc

@@ -2 +2 @@
 
 // System includes
+#include <iostream>
 //#include <algorithm>
 //#include <utility>
@@ -19 +20 @@
     
   {
-    for (unsigned int i=0; i<target.size(); i++){
-      int targ=static_cast<int>(target(i));
-      std::pair<int, double> tmp(targ, value(i));
-      value_.push_back(tmp);
-      if (targ==1)
-        nof_pos_++;
-    }
-    sort(value_.begin(),value_.end(),
-         pair_value_compare<int,double>());
-    
+    sort(target, value);
   }
 
@@ -36 +28 @@
   {
   }
-
-
-  double ROC::score(const gslapi::vector& target, const gslapi::vector& value)
-  {
-    for (unsigned int i=0; i<target.size(); i++){
-      int targ=static_cast<int>(target(i));
-      std::pair<int, double> tmp(targ, value(i));
-      value_.push_back(tmp);
-      if (targ==1)
-        nof_pos_++;
-    }
-    sort(value_.begin(),value_.end(),
-         pair_value_compare<int,double>());
-    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::score()
-  {
-    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::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_), 
-                          nof_pos_, value_.size()-nof_pos_);
-    else
-    p = get_p_approx(area_);
-    return p;
-  }
-
 
   double ROC::get_p_approx(const double area) const
@@ -104 +46 @@
   }
 
-
   double ROC::get_p_exact(const double block, const double nof_pos, 
                             const double nof_neg)
@@ -121 +62 @@
   }
 
+  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_), 
+                          nof_pos_, value_.size()-nof_pos_);
+    else
+    p = get_p_approx(area_);
+    return p;
+  }
+
+  double ROC::score()
+  {
+    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::score(const gslapi::vector& target, const gslapi::vector& value)
+  {
+    sort(target, value);
+    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_;
+  }
+
+  void ROC::sort(const gslapi::vector& target, const gslapi::vector& value)
+  {
+    for (unsigned int i=0; i<target.size(); i++){
+      int targ=static_cast<int>(target(i));
+      std::pair<int, double> tmp(targ, value(i));
+      value_.push_back(tmp);
+      if (targ==1)
+        nof_pos_++;
+    }
+    std::sort(value_.begin(),value_.end(),
+         pair_value_compare<int,double>());
+  }
+
+  gslapi::vector ROC::target(void) const
+  {
+    gslapi::vector target(value_.size());
+    for(size_t i=0; i<target.size(); ++i) 
+      target(i) = value_[i].first;
+    return target;
+  }
+
+  std::ostream& operator<<(std::ostream& s, const ROC& r)
+  {
+    s.setf( std::ios::dec );
+    s.precision(12);
+    double sens = 1;
+    double spec = 0;
+    gslapi::vector target = r.target();
+    size_t n = target.size();
+    double nof_pos = (target.sum()+n)/2;
+    for(size_t i=0; i<n-1; ++i) {
+      s << sens << "\t";
+      s << spec << std::endl;
+      if (target(i)==1)
+        spec -= 1/(n-nof_pos);
+      else 
+        sens -= 1/nof_pos;
+    }
+    s << sens << "\t";
+    s << spec << std::endl;
+    return s;
+  }
+
+
 }} // of namespace cpptools and namespace theplu