PUBLIC

WeightedLeastSqr2dC::WeightedLeastSqr2dC(const DListC &) WeightedLeastSqr2dC::WeightedLeastSqr2dC(const Array1dC &) WeightedLeastSqr2dC::Line(void) const WeightedLeastSqr2dC::LackOfFit(RealT) const WeightedLeastSqr2dC::InLimits(const LineABC2dC &) const WeightedLeastSqr2dC::JointConfidence(const LineABC2dC &) const WeightedLeastSqr2dC::JointConfidenceLucia(const LineABC2dC &) const WeightedLeastSqr2dC::SigmaIntercept(void) const WeightedLeastSqr2dC::SigmaSlope(void) const WeightedLeastSqr2dC::Chi2(void) const WeightedLeastSqr2dC::Q(void) const WeightedLeastSqr2dC::myfit(void)

WeightedLeastSqr2dC
Performs weighted least squares fit on a set of data points
include	"amma/Num/WeightedLeastSqr2d.hh"
User Level:	Default
Library:	NumRBF
Example:	doRbfValidate.cc
Section:	default.Kieron J Messer
In Scope:	std

Comments:
Algorithm taken from Bevington P.R. 1969, Data Reduction and Error Analysis for the Physical Sciences (New York: McGraw-Hill), Chapter 6.

RealArray1dC x;

x-values

RealArray1dC y;

y-values

RealArray1dC stdev;

standard deviation

LineABC2dC line;

estimated line

RealT lof;

the quality of the fit

RealT chi2;

the chi2 statistic

RealT siginter;

intercept sigma

RealT sigslope;

slope sigma

RealT covar;

the covar of slope/intercept

WeightedLeastSqr2dC(const DListC<Point3dC> & pts)

Perform weighted least squares on data

WeightedLeastSqr2dC(const Array1dC<Point3dC> & pts)

Perform weighted least squares on data

LineABC2dC Line() const

Access to the fitted line

BooleanT LackOfFit(RealT prob = 0.01) const

test for lack of fit at desired probability level

BooleanT InLimits(const LineABC2dC & test) const

test line is in confidence limits

BooleanT JointConfidence(const LineABC2dC & test) const

test line is in confidence limits

RealT JointConfidenceLucia(const LineABC2dC & test) const

this routine returns the Fratio on joint confidence level do not use does not work properly, use InLimits()

RealT SigmaIntercept() const

return sigma of intercept

RealT SigmaSlope() const

return sigma of intercept

RealT Chi2() const

return the Chi2 value for LOF

RealT Q() const

return probability for LOF

void myfit()

perform alg