[SciPy-user] curve_fit failure when supplying uncertainties

josef.pktd@gmai... josef.pktd@gmai...
Sat Mar 7 19:20:07 CST 2009


On Sat, Mar 7, 2009 at 7:18 PM, ElMickerino <elmickerino@hotmail.com> wrote:
>
> I think I've identified a problem with curve_fit which occurs when one
> attempts to fit normally distributed data with a gaussian.  From the
> documentation of curve_fit, it appears that 'sigma' should be the
> uncertainties on the y-values of the data; however, the following example
> (see attached code) should make it clear that there's a problem with this.
>
> My best guess is that the sigma are actually weights (=1.0/sigma).  Can
> anyone confirm or deny this?  Seems like from the name it should be
> uncertainties but from the behavior of the code it appears otherwise.
>
> Also, I was wondering if there's a way to supply asymmetric errors to
> curve_fit (or for that matter, to leastsqr or any wrapper thereof).
>
> Thanks very much,
> Michael
>
> http://www.nabble.com/file/p22380378/testFit.py testFit.py
> --
> View this message in context: http://www.nabble.com/curve_fit-failure-when-supplying-uncertainties-tp22380378p22380378.html
> Sent from the Scipy-User mailing list archive at Nabble.com.
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>

sigma is the vector of standard deviations for the observations. The
weights in the least squares objective function are 1/sigma.
Essentially it is a heteroscedastic non-linear generalized least
squares estimation.

In your case err contains zeros, which means the variance is zero and
the weight 1/sigma for these observations would be infinite. curve_fit
cannot handle zero standard deviation. Infinite standard deviation
works because then the weight becomes zero.

If you force your err to be strictly positive, then this works, e.g.

>>> curve_fit(gaus, centers, data, sigma=1e-6+err)
(array([  1.91520894e+03,   1.68957830e-01,   6.19128687e-01]),
array([[  3.03590409e-09,   5.11571694e-18,  -6.69243949e-18],
       [  5.11571694e-18,   4.43479563e-16,   6.18126633e-17],
       [ -6.69243949e-18,   6.18126633e-17,   3.02641036e-17]]))


I haven't looked closely at the rest of your file, so I'm not sure
about the interpretation of your err

Josef


More information about the SciPy-user mailing list