source: branches/0.4-stable/test/fisher_test.cc @ 1621

// $Id: fisher_test.cc 1621 2008-11-10 21:39:48Z peter $

/*
  Copyright (C) 2008 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/statistics/Fisher.h"

int main(int argc, char* argv[])
{  
  using namespace theplu::yat;
  test::Suite suite(argc, argv);

  statistics::Fisher f;
  if (!suite.equal(f.oddsratio(1,4,5,1), 0.05)) {
    suite.add(false);
    suite.err() << "oddsratio failed\n";
  }
  double a, b, c, d;
  f.expected(a,b,c,d);
  if (!suite.equal(a, 30.0/11.0)) {
    suite.add(false);
    suite.err() << "expected a failed\n";
  }
  if (!suite.equal(b, 25.0/11.0)) {
    suite.add(false);
    suite.err() << "expected b failed\n";
  }
  if (!suite.equal(c, 36.0/11.0)) {
    suite.add(false);
    suite.err() << "expected c failed\n";
  }
  if (!suite.equal(d, 30.0/11.0)) {
    suite.add(false);
    suite.err() << "expected d failed\n";
  }
  if (!suite.equal(f.Chi2(), 4.4122222222222222222222)) { 
    suite.add(false);
    suite.err() << "Chi2 failed\n";
  }
  if (!suite.equal(f.minimum_size(),10)) {
    suite.add(false);
    suite.err() << "minimum_size failed\n";
  }
  return suite.return_value();
}