source: trunk/src/Regression.h @ 221

// $Id: Regression.h 221 2004-12-30 22:36:25Z peter $

#ifndef _theplu_statistics_regression_
#define _theplu_statistics_regression_

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

// Standard C++ includes
////////////////////////


namespace theplu {
namespace statistics {  

  
  ///
  /// Virtual Class for Regression.   
  ///
 
  class Regression 
  {
  
  public:
    ///
    /// Default Constructor.
    ///
    Regression();

    ///
    /// Destructor 
    ///
    virtual ~Regression(void) {};
         
    ///
    /// This function computes the best-fit given a model (see
    /// specific class for details) by minimizing \f$
    /// \sum{(\hat{y_i}-y_i)^2} \f$, where \f$ \hat{y} \f$ is the fitted value.
    ///
    virtual void fit(const gslapi::vector& x, const gslapi::vector& y)=0; 
    
    ///
    /// This function computes the best-fit given a model (see
    /// specific class for details) by minimizing \f$
    /// \sum{w_i(\hat{y_i}-y_i)^2} \f$, where \f$ \hat{y} \f$ is the
    /// fitted value. The weight \f$ w_i \f$ is should be proportional
    /// to the inverse of the variance for \f$ y_i \f$
    ///
    virtual void fit(const gslapi::vector& x, const gslapi::vector& y, 
                    const gslapi::vector& w)=0;
    ///
    ///
    ///
    virtual void predict(const double x, double& y, double& y_err, 
                         const double w=1) const=0;
             
  protected:
    
  };

}} // of namespace statistics and namespace theplu

#endif