Changeset 465 for trunk/lib/statistics/WilcoxonFoldChange.cc

Timestamp:

Dec 17, 2005, 12:19:24 AM (17 years ago)

Author:

Peter

Message:

implemented score for WilcoxonFoldChange?

File:

• trunk/lib/statistics/WilcoxonFoldChange.cc (modified) (3 diffs)

trunk/lib/statistics/WilcoxonFoldChange.cc

@@ -1 +1 @@
 // $Id$
 
-#include "WilcoxonFoldChange.h"
+#include <c++_tools/statistics/WilcoxonFoldChange.h>
+#include <c++_tools/statistics/utility.h>
 
-#include "Score.h"
-#include "Averager.h"
-#include "AveragerWeighted.h"
+#include <cmath>
 
 namespace theplu {
@@ -28 +27 @@
       train_set_=train_set;
 
+    std::vector<double> distance;
+    //Peter, should reserve the vector to avoid reallocations   
     weighted_=false;
-
+    for (size_t i=0; i<train_set_.size(); i++) {
+      for (size_t j=0; i<j; j++) {
+        if (target(train_set_[i])==target(train_set_[j])) continue;
+        distance.push_back(value(train_set_[i])-value(train_set_[j]));
+      }
+    }
+    if (absolute_)
+      return fabs(median(distance));
+    return median(distance);
   }
 
@@ -37 +46 @@
                                    const std::vector<size_t>& train_set) 
   {
-    if (!train_set_.size())
-      for (size_t i=0; i<target_.size(); i++)
-        train_set_.push_back(i);  
-    else
-      train_set_=train_set;
-
-    weighted_=true;
+    return 0;
   }