[Numpy-discussion] Numpy fitting

Olivier Delalleau shish@keba...
Thu Mar 1 09:39:13 CST 2012


Sorry I can't help, but I'd just suggest to post this on the scipy mailing
list as you may get more replies there.

-=- Olivier

Le 1 mars 2012 10:24, Pierre Barthelemy <barthpi@gmail.com> a écrit :

> Dear all,
>
> i am writing a program for data analysis. One of the functions of this
> program gives the possibility to fit the functions. I followed the recipe
> described in : h*MailScanner soupçonne le lien suivant d'être une
> tentative de fraude de la part de "www.scipy.org" *ttp://www.scipy.org/Cookbook/FittingData<http://www.scipy.org/Cookbook/FittingData> under
> the section "Simplifying the syntax".
>
> To fit, i use the function:
>
> out=optimize.leastsq(f, p, full_output=1)
>
> where f is my function and p a list of parameters.
>
> One thing that i would like to know is how can i get the error on the
> parameters ? From what i understood from the "Cookbook" page, and from the
> scipy manual (
> http://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.leastsq.html#scipy.optimize.leastsq),
> the second argument returned by the leastsq function gives access to these
> errors.
> std_error=std(y-function(x))
> param_error=sqrt(diagonal(out[1])*std_error)
>
> The param_errors that i get in this case are extremely small. Much smaller
> than what i expected, and much smaller than what i can get fitting the
> function with matlab. So i guess i made an error here.
>
> Can someone tell me how i should do to retrieve the parameter errors ?
>
> Bests,
>
> Pierre
>
> PS: i got the impression something went wrong with my previous message,
> sorry for that.
>
>
> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion@scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/numpy-discussion/attachments/20120301/a5d5c98a/attachment.html 


More information about the NumPy-Discussion mailing list