source: trunk/test/consensus_inputranker_test.cc @ 2119

// $Id: consensus_inputranker_test.cc 2119 2009-12-12 23:11:43Z peter $

/*
  Copyright (C) 2004, 2005 Peter Johansson
  Copyright (C) 2006, 2007, 2008 Jari Häkkinen, Peter Johansson

  This file is part of the yat library, http://dev.thep.lu.se/yat

  The yat library is free software; you can redistribute it and/or
  modify it under the terms of the GNU General Public License as
  published by the Free Software Foundation; either version 3 of the
  License, or (at your option) any later version.

  The yat library is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with yat. If not, see <http://www.gnu.org/licenses/>.
*/

#include "Suite.h"

#include "yat/classifier/ConsensusInputRanker.h"
#include "yat/statistics/AUC.h"
#include "yat/utility/Matrix.h"
#include "yat/classifier/MatrixLookup.h"
#include "yat/classifier/CrossValidationSampler.h"
#include "yat/classifier/IRRank.h"
#include "yat/statistics/VectorFunction.h"

#include <cstdlib>
#include <fstream>
#include <iostream>

using namespace std;

int main(int argc, char* argv[])
{  
  using namespace theplu::yat;
  theplu::yat::test::Suite suite(argc, argv);
  suite.err() << "testing consensus_inputranker" << std::endl;

  ifstream is(test::filename("data/rank_data.txt").c_str());
  theplu::yat::utility::Matrix data_tmp(is);
  theplu::yat::classifier::MatrixLookup data(data_tmp);
  is.close();

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

  
  theplu::yat::statistics::AUC roc;
  theplu::yat::classifier::CrossValidationSampler sampler(target,30,3);
  suite.err() << "Building Consensus_Inputranker" << std::endl;
  theplu::yat::classifier::IRRank retrieve;
  theplu::yat::statistics::Median median;
  theplu::yat::classifier::ConsensusInputRanker cir(retrieve,median);
  cir.add(sampler,data,roc);

  suite.err() << "test ids... ";
  if (cir.id(0)!=2 || cir.id(1)!=0 || cir.id(2)!=1){
    suite.err() << "\nincorrect id for weighted" << endl;
    suite.add(false);
  }
  else 
    suite.err() << "ok." << std::endl;

  suite.err() << "test ranks... ";
  if (cir.rank(0)!=1 || cir.rank(1)!=2 || cir.rank(2)!=0){
    suite.err() << "\nincorrect rank for weighted" << endl;
    suite.add(false);
  }
  else 
    suite.err() << "ok." << std::endl;

  theplu::yat::utility::Matrix flag(data.rows(),data.columns(),1);
  // Peter, fix weighted version instead
  theplu::yat::classifier::ConsensusInputRanker cir2(retrieve,median);
  cir2.add(sampler,data,roc);

  suite.err() << "test ids... ";
  if (cir2.id(0)!=2 || cir2.id(1)!=0 || cir2.id(2)!=1){
    suite.err() << "\nincorrect id for weighted" << endl;
    suite.add(false);
  }
  else 
    suite.err() << "ok." << std::endl;
  
  suite.err() << "test ranks... ";
  if (cir2.rank(0)!=1 || cir2.rank(1)!=2 || cir2.rank(2)!=0){
    suite.err() << "\nincorrect rank for weighted" << endl;
    suite.add(false);
  }
  else 
    suite.err() << "ok." << std::endl;

  return suite.return_value();
}