source: branches/peters_vector/test/consensus_inputranker_test.cc @ 470

// $Id: consensus_inputranker_test.cc 470 2005-12-19 19:46:50Z peter $

// C++ tools include
////////////////////
#include <c++_tools/classifier/ConsensusInputRanker.h>
#include <c++_tools/statistics/ROC.h>
#include <c++_tools/gslapi/matrix.h>
#include <c++_tools/classifier/MatrixLookup.h>
#include <c++_tools/classifier/CrossSplitting.h>

#include <cstdlib>
#include <fstream>
#include <iostream>
//#include <vector>

using namespace std;

int main()
{  
  ifstream is("data/rank_data.txt");
  theplu::gslapi::matrix data_tmp(is);
  theplu::classifier::MatrixLookup data(data_tmp);
  is.close();

  is.open("data/rank_target.txt");
  theplu::classifier::Target target(is);
  is.close();

  
  theplu::statistics::ROC roc;
  theplu::classifier::CrossSplitting sampler(3);
  theplu::classifier::ConsensusInputRanker cir(data,target,roc,sampler,30);

  if (cir.id(0)!=2 || cir.id(1)!=0 || cir.id(2)!=1){
    cerr << "wrong id" << endl;
    return -1;
  }
  
  if (cir.rank(0)!=1 || cir.rank(1)!=2 || cir.rank(2)!=0){
    cerr << "wrong rank" << endl;
    return -1;
  }

  theplu::gslapi::matrix flag(data.rows(),data.columns(),1);
  theplu::classifier::ConsensusInputRanker 
    cir2(data,target,flag,roc,sampler,30);

  if (cir2.id(0)!=2 || cir2.id(1)!=0 || cir2.id(2)!=1){
    cerr << "wrong id" << endl;
    return -1;
  }
  
  if (cir2.rank(0)!=1 || cir2.rank(1)!=2 || cir2.rank(2)!=0){
    cerr << "wrong rank" << endl;
    return -1;
  }

  
  return 0;
}