[SciPy-user] Numerical Recipes robust fit implementation

Andrew Straw strawman@astraw....
Mon Jul 16 01:37:03 CDT 2007


Hi Angus,

I don't know if the MEDFIT algorithm is something you really want, but I
wrote a simple implementation of the RANSAC algorithm:
http://www.scipy.org/Cookbook/RANSAC

This is not to say I recommend the use of such algorithms as RANSAC! :)

-Andrew

PS There seems to be some criticism of the MEDFIT algorithm in general
and as implemented in NR in particular online.

Angus McMorland wrote:
> Hi all,
> 
> I'm in need of a robust linear-fit algorithm- Numerical Recipes (in C)
> medfit routine will do nicely. I've had a quick look and don't see any
> similar routines in scipy. Is there one I'm missing?
> 
> If not, I'd lke to implement the NR routine for use with my
> numpy-arrayed data. What are my options? Obviously I could translate
> the routine into python, but I realise a potentially easier
> alternative is to interface the existing NR C code. Having never done
> that before, can someone give me some pointers/tips on how to start -
> there's quite a range of linking tools listed on
> http://www.scipy.org/Topical_Software.
> 
> Thanks,
> 
> Angus.



More information about the SciPy-user mailing list