[SciPy-user] Numerical Recipes robust fit implementation
Andrew Straw
strawman at astraw.com
Mon Jul 16 02:37:03 EDT 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