[SciPy-user] Curve fitting and LaTeX output

[Warren Weckesser]

Hi David,

On Thu, Aug 28, 2008 at 8:56 AM, David Lonie <loniedavid at gmail.com> wrote:

> Thanks for the fast reply! I plugged in my data to that script and
> added some output, but I get an invalid fit out of it. The fit I got
> from oocalc is approx
>
> y = 32 - 0.85^x, and fmin is returning
> y = 12 - 0.65^x, which for some reason returns an array of non-numbers
> when I try to generate plot data?


If c[0]=12 and c[1]=-0.65, then your function is not 12-0.65**x,
but rather 12 + (-0.65)**x, which is not defined when x is not
an integer.   As Ryan suggested, try a better starting guess.
You know the function must be decreasing, so c[1] must be
between 0 and 1.  If I use c_initial_guess = [10, 0.5], I get
c_final =  [12.20289742   0.65819867]
This is still a bad fit; I think this sum is not a good model for
your data.  You'll get a better fit with a simple scaled
exponential:  c[0]*(c[1]**x)


Cheers,

Warren
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20080828/75f4e946/attachment.html>