ROB_CHECKFIT
Name
ROB_CHECKFIT
Purpose
Used by ROBUST_... routines to determine the quality of a fit and to
return biweights.
Calling Sequence
status = ROB_CHECKFIT( Y, YFIT, EPS, DEL, SIG, FRACDEV, NGOOD, W, B
BISQUARE_LIMIT = )
Input
Y = the data
YFIT = the fit to the data
EPS = the "too small" limit
DEL = the "close enough" for the fractional median abs. deviations
Returns
Integer status. if =1, the fit is considered to have converged
Outputs
SIG = robust standard deviation analog
FRACDEV = the fractional median absolute deviation of the residuals
NGOOD = the number of input point given non-zero weight in the
calculation
W = the bisquare weights of Y
B = residuals scaled by sigma
Optional Input Keyword
BISQUARE_LIMIT = allows changing the bisquare weight limit from
default 6.0
Procedures Used
ROBUST_SIGMA()
Revision History
Written, H.T. Freudenreich, HSTX, 1/94