source: branches/peters_vector/lib/classifier/InputRanker.cc @ 469

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

#include <c++_tools/classifier/InputRanker.h>

#include <c++_tools/classifier/DataView.h>
#include <c++_tools/classifier/VectorView.h>
#include <c++_tools/classifier/Target.h>
#include <c++_tools/statistics/ROC.h>
#include <c++_tools/utility/stl_utility.h>

#include <vector>
#include <utility>

namespace theplu {
namespace classifier {  


  InputRanker::InputRanker(const DataView& data, 
                           const Target& target, 
                           statistics::Score& score_object)
  {
    size_t nof_genes = data.rows();
    size_t nof_samples = data.columns();
    if (!train_set_.size()){
      train_set_.resize(nof_samples);
      for (size_t i=0; i<nof_samples; i++)
        train_set_[i]=i;  
    }

    //scoring each input
    std::vector<std::pair<size_t, double> > score;
    for (size_t i=0; i<nof_genes; i++) {
      double area = score_object.score(target,VectorView(data,i,false));
      std::pair<size_t, double> tmp(i,area);
      score.push_back(tmp);
    }

    //sort the scores and assign id_ and rank_
    sort(score.begin(), score.end(), 
         utility::pair_value_compare<size_t, double>());
    
    for (size_t i=0; i<nof_genes; i++){
      id_[i]=score[nof_genes-i-1].first;
      rank_[id_[i]]=i;            
    }
  }



  InputRanker::InputRanker(const DataView& data, 
                           const Target& target, 
                           statistics::Score& score_object,
                           const DataView& weight)
  {
    size_t nof_genes = data.rows();
    size_t nof_samples = data.columns();
    if (!train_set_.size()){
      train_set_.resize(nof_samples);
      for (size_t i=0; i<nof_samples; i++)
        train_set_[i]=i;  
    }
    //scoring each input
    std::vector<std::pair<size_t, double> > score;
    for (size_t i=0; i<nof_genes; i++) {
      double area = score_object.score(target, VectorView(data,i,false),
                                       VectorView(weight,i,false));
      std::pair<size_t, double> tmp(i,area);
      score.push_back(tmp);
    }
    //sort the scores and assign id_ and rank_
    sort(score.begin(), score.end(), 
         utility::pair_value_compare<size_t, double>());

    for (size_t i=0; i<nof_genes; i++){
      id_[i]=score[nof_genes-i-1].first;
      rank_[id_[i]]=i;            
    }
  }


}} // of namespace classifier and namespace theplu