Changeset 529 for trunk/lib/statistics/tScore.cc

Timestamp:

Mar 1, 2006, 3:03:27 PM (16 years ago)

Author:

Markus Ringnér

Message:

Added score for signal-to-noise ratio (Golub score). Fixed documentation. Improved weighted version of tScore

File:

• trunk/lib/statistics/tScore.cc (modified) (2 diffs)

trunk/lib/statistics/tScore.cc

@@ -50 +50 @@
     statistics::AveragerWeighted positive;
     statistics::AveragerWeighted negative;
-    dof_=target.size()-2;
     for(size_t i=0; i<target.size(); i++){
       if (target.binary(i))
@@ -60 +59 @@
     double s2=(positive.sum_xx_centered()+negative.sum_xx_centered())/
       (positive.n()+negative.n()-2);
-    t_=diff/sqrt(s2*(1.0/positive.sum_w()+1.0/negative.sum_w()));
+    t_=diff/sqrt(s2*(1.0/positive.sum_w()+1.0/negative.sum_w()));   
     if (t_<0 && absolute_)
       t_=-t_;
+
+    if(positive.n()==0 || negative.n()==0 ||
+       positive.sum_w()==0 || positive.sum_w()==0)
+      t_=0;
+    dof_=target.size()-2;
      
     return t_;