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

// $Id: WilcoxonFoldChange.cc 469 2005-12-19 14:58:29Z peter $

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

#include <cmath>
#include <vector>

namespace theplu {
namespace statistics {


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



  double WilcoxonFoldChange::score(const classifier::Target& target, 
                                   const classifier::VectorAbstract& value)
  {
    std::vector<double> distance;
    //Peter, should reserve the vector to avoid reallocations   
    weighted_=false;
    for (size_t i=0; i<target.size(); i++) {
      for (size_t j=0; i<j; j++) {
        if (target(i)==target(j)) continue;
        distance.push_back(value(i)-value(j));
      }
    }
    if (absolute_)
      return fabs(median(distance));
    return median(distance);
  }

  double WilcoxonFoldChange::score(const classifier::Target& target, 
                                   const classifier::VectorAbstract& value,
                                   const classifier::VectorAbstract& weight)
  {
    return 0;
  }

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