source: trunk/test/statistics_test.cc @ 865

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

/*
  Copyright (C) 2004 Jari Häkkinen, Peter Johansson
  Copyright (C) 2005 Peter Johansson
  Copyright (C) 2006 Jari Häkkinen, Markus Ringnér, 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/statistics/utility.h"
#include "yat/utility/vector.h"

#include <vector>
#include <cstdlib>
#include <iostream>
#include <cmath>

int main()
{  
  using namespace theplu::yat;
  utility::vector gsl_vec(10);
  std::vector<double> data;
  for (unsigned int i=0; i<10; i++){
    data.push_back(static_cast<double>(i));
    gsl_vec(i)=i;
  }

  double m=statistics::median(data);
  double m_gsl=statistics::median(gsl_vec);
  if (m!=4.5 || m!=m_gsl)
    return -1;

  double tolerance=1e-10;
  double skewness_gsl=statistics::skewness(gsl_vec);
  if (fabs(skewness_gsl)>tolerance)
    return -1;
  double kurtosis_gsl=statistics::kurtosis(gsl_vec);
  if (fabs(kurtosis_gsl+1.5616363636363637113)>tolerance) 
    return -1;  
  return 0;
}