[SciPy-User] Least median of squares for regression in scipy?
Tue Aug 24 08:27:59 CDT 2010
"least median of squares" doesn't mean anything to me. But, I know that
minimizing sum of absolute differences will provide a good estimate of the
median and is a good technique for dealing with outliers:
Note that you'll need an LP solver. Another option is a hybrid between the
squared and absolute value loss functions, such as the one that Peter Huber
This loss provides the outlier-insensitivity of L1 while being easy to solve
using gradient descent & line search.
On Tue, Aug 24, 2010 at 9:17 AM, Jorge Scandaliaris
> I have to perform a linear regression on noisy data. On the last paper I
> least median of squares was suggested for dealing with outliers. I have
> the scipy docs but it seems nothing is readily available. Searching the web
> "(python OR scipy OR numpy) least median square" doesn't yield meaningfull
> results. The best I found were fortran and matlab code, which I would need
> wrap (I have zero knowledge about fortran or wrapping it into python,
> that there's a tool called f2py, but I would have to learn that as well) or
> rewrite (I used matlab in the past, so this should be feasible).
> I am asking here in the hope there's something I overlooked before I jump
> one of the (probably more time demanding) possibilities I mentioned above.
> Thanks in advance,
> SciPy-User mailing list
Research Scientist, ITA Software
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-User