source: trunk/test/consensus_inputranker_test.cc @ 865

// $Id: consensus_inputranker_test.cc 865 2007-09-10 19:41:04Z peter $

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

  This file is part of the yat library, http://trac.thep.lu.se/trac/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 2 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 this program; if not, write to the Free Software
  Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA
  02111-1307, USA.
*/

#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(const int argc,const char* argv[])
{  
  std::ostream* error;
  if (argc>1 && argv[1]==std::string("-v"))
    error = &std::cerr;
  else {
    error = new std::ofstream("/dev/null");
    if (argc>1)
      std::cout << "consensus_inputranker_test -v : for printing extra " 
                << "information\n";
  }
  *error << "testing consensus_inputranker" << std::endl;
  bool ok = true;

  ifstream is("data/rank_data.txt");
  theplu::yat::utility::matrix data_tmp(is);
  theplu::yat::classifier::MatrixLookup data(data_tmp);
  is.close();

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

  
  theplu::yat::statistics::AUC roc;
  theplu::yat::classifier::CrossValidationSampler sampler(target,30,3);
  *error << "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);

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

  *error << "test ranks... ";
  if (cir.rank(0)!=1 || cir.rank(1)!=2 || cir.rank(2)!=0){
    *error << "\nincorrect rank for weighted" << endl;
    ok=false;
  }
  else 
    *error << "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);

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

  if (error!=&std::cerr)
    delete error;

  if(ok)
    return 0;
  return -1;
}