source: trunk/test/test_regression_linear.cc @ 238

// $Id: test_regression_linear.cc 238 2005-02-21 16:42:00Z peter $

// C++ tools include
////////////////////
#include "RegressionLinear.h"
#include "vector.h"


// Standard includes
////////////////////
#include <gsl/gsl_fit.h>
#include <iostream>


using namespace std;

int main()

{  
  bool ok=true;

  theplu::gslapi::vector x(4);
  x(0)=1970;
  x(1)=1980;
  x(2)=1990;
  x(3)=2000;
  theplu::gslapi::vector y(4);
  y(0)=12;
  y(1)=11;
  y(2)=14;
  y(3)=13;
  theplu::gslapi::vector w(4);
  w(0)=0.1;
  w(1)=0.2;
  w(2)=0.3;
  w(3)=0.4;

  theplu::statistics::RegressionLinear r;
  r.fit(x,w,y);
  double y_predicted=0;
  double y_predicted_err=0;
  r.predict(0,y_predicted,y_predicted_err);

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