yatsm.regression.robust_fit module¶
Perform an iteratively re-weighted least squares ‘robust regression’. Basically a clone of statsmodels.robust.robust_linear_model.RLM without all the lovely, but costly, creature comforts.
- Reference:
- http://statsmodels.sourceforge.net/stable/rlm.html http://cran.r-project.org/web/packages/robustreg/index.html http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-robust-regression.pdf
Run this file to test performance gains. Implementation is ~3x faster than statsmodels and can reach ~4x faster if Numba is available to accelerate.
-
yatsm.regression.robust_fit.
bisquare
(resid, c=4.685)[source]¶ Returns weighting for each residual using bisquare weight function
Parameters: Returns: weights for residuals
Return type: weight (ndarray)