regressionWithoutIntercept.hpp

Go to the documentation of this file.

#ifndef JULIAN_REGRESSIONWITHOUTINTERCEPT_HPP
#define JULIAN_REGRESSIONWITHOUTINTERCEPT_HPP

#include <mathematics/regressions/regression.hpp>
#include <iostream>
#include <vector>
#include <string>
#include <gsl/gsl_fit.h>

namespace julian {

  class RegressionWithoutIntercept : public DeeplyCopyableRegression<RegressionWithoutIntercept> {
  public:
    RegressionWithoutIntercept(): chi_sq_(0.0), cov00_(0.0), cov01_(0.0), cov11_(0.0),  c1_(0.0) {};

    void estimate(const std::vector<double>& x,const std::vector<double>& y);
    void estimate(const std::vector<double>& x,const std::vector<double>& y,const std::vector<double>& w);
    
    std::vector<double> getCoefficient() const;
    double operator()(double) const;

    friend std::ostream& operator<<(std::ostream& s, julian::RegressionWithoutIntercept& r);
    
  private:
    double chi_sq_; 
    double cov00_;  
    double cov01_;  
    double cov11_;  
    double c1_;     
  };
}  //  namespace julian

#endif