webbpsf_ext.robust

Small collection of robust statistical estimators based on functions from Henry Freudenriech (Hughes STX) statistics library (called ROBLIB) that have been incorporated into the AstroIDL User’s Library. Function included are:

  • medabsdev - median absolute deviation

  • biweightMean - biweighted mean estimator

  • mean - robust estimator of the mean of a data set

  • mode - robust estimate of the mode of a data set using the half-sample method

  • std - robust estimator of the standard deviation of a data set

  • checkfit - return the standard deviation and biweights for a fit in order to determine its quality

  • linefit - outlier resistant fit of a line to data

  • polyfit - outlier resistant fit of a polynomial to data

For the fitting routines, the coefficients are returned in the same order as np.polyfit, i.e., with the coefficient of the highest power listed first.

For additional information about the original IDL routines, see:

http://idlastro.gsfc.nasa.gov/contents.html#C17

Functions

biweightMean(inputData[, axis, dtype, iterMax])

Biweight Mean

checkfit(inputData, inputFit, epsilon, delta)

Determine the quality of a fit and biweights.

linefit(inputX, inputY[, iterMax, Bisector, ...])

Outlier resistance two-variable linear regression function.

mean(inputData[, Cut, axis, dtype, ...])

Robust Mean

medabsdev(data[, axis, keepdims, nan])

Median Absolute Deviation

mode(inputData[, axis, dtype])

Robust estimator of the mode of a data set using the half-sample mode.

polyfit(inputX, inputY, order[, iterMax])

Outlier resistance two-variable polynomial function fitter.

std(inputData[, Zero, axis, dtype, ...])

Robust Sigma