GeneralizedChiSquared.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2020 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/General/GeneralizedChiSquared.h
 
#include "BinaryStatistic.h"
#include "ndArray.h"
 
namespace BinaryStatistics
{

  class GeneralizedChiSquared : public BinaryStatistic
  {
    bool robust_linear_fit;
    double outlier_threshold;
    double max_zap_fraction;
 
    std::vector<double> residual;
 
  public:
 
    GeneralizedChiSquared ();
 
    double get (const std::vector<double>&, const std::vector<double>&);

    const std::vector<double>& get_residual () const { return residual; }
    
    GeneralizedChiSquared* clone () const { return new GeneralizedChiSquared; }

    unsigned get_neigen () const { return eigenvalues.size(); }
    
    ndArray<2,double> eigenvectors;
    ndArray<1,double> eigenvalues;
    
  };
}