source: branches/peters_vector/lib/statistics/WilcoxonFoldChange.cc @ 466

// $Id: WilcoxonFoldChange.cc 465 2005-12-16 23:19:24Z peter $

#include <c++_tools/statistics/WilcoxonFoldChange.h>
#include <c++_tools/statistics/utility.h>

#include <cmath>

namespace theplu {
namespace statistics {


  WilcoxonFoldChange::WilcoxonFoldChange(bool absolute)
    : Score(absolute)
  {
  }



  double WilcoxonFoldChange::score(const gslapi::vector& target, 
                                   const gslapi::vector& value,
                                   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;

    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);
  }

  double WilcoxonFoldChange::score(const gslapi::vector& target, 
                                   const gslapi::vector& value,
                                   const gslapi::vector& weight,
                                   const std::vector<size_t>& train_set) 
  {
    return 0;
  }

} // of namespace statistics
} // of namespace theplu