source: trunk/lib/statistics/Naive.cc @ 430

// $Id: Naive.cc 430 2005-12-08 22:53:08Z peter $

#include <c++_tools/statistics/Naive.h>

#include <c++_tools/statistics/Averager.h>
#include <c++_tools/statistics/AveragerWeighted.h>
#include <c++_tools/statistics/OneDimensional.h>
#include <c++_tools/gslapi/vector.h>

#include <iostream>


namespace theplu {
namespace statistics {
namespace regression {


  void Naive::fit(const gslapi::vector& x, const gslapi::vector& y)
  {
    ap_.reset();
    for (size_t i=0; i<y.size(); i++)
      ap_.add(x(i),y(i));
  }

  void Naive::predict(const double x, double& y, double& y_err) const
  { 
    y = ap_.x_averager().mean(); 
    y_err = ap_.x_averager().std()*(1.0+1.0/ap_.n());
  }
  
  std::ostream& Naive::print_header(std::ostream& s) const
  {
    s << "# column 1: x\n"
      << "# column 2: y\n"
      << "# column 3: y_err";
    return s;
  }


}}} // of namespaces regression, statisitcs and thep