[SciPy-User] Least median of squares for regression in scipy?

Jason Rennie jrennie@gmail....
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

> Hi,
> I have to perform a linear regression on noisy data. On the last paper I
> read
> least median of squares was suggested for dealing with outliers. I have
> searched
> the scipy docs but it seems nothing is readily available. Searching the web
> for
> "(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
> to
> wrap (I have zero knowledge about fortran or wrapping it into python,
> except
> 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
> into
> one of the (probably more time demanding) possibilities I mentioned above.
> Thanks in advance,
> Jorge
> _______________________________________________
> SciPy-User mailing list
> SciPy-User@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user

Jason Rennie
Research Scientist, ITA Software
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-user/attachments/20100824/bc8e23cc/attachment.html 

More information about the SciPy-User mailing list