[SciPy-User] lmfit: confidence intervals issue

josef.pktd at gmail.com josef.pktd at gmail.com
Wed Jun 18 20:09:04 EDT 2014


On Wed, Jun 18, 2014 at 7:25 PM, Sturla Molden <sturla.molden at gmail.com> wrote:
> lmfit is not a part of SciPy, but as long as it's scientific use of Python
> and no other appropriate list, I don't think anyone will mind. I cannot
> answer your question, though.
>
> Sturla
>
> Antonino Ingargiola <tritemio at gmail.com> wrote:
>> Hi,
>>
>> I'm not sure this is the right list for this kind of issues.
>>
>> I'm using lmfit to solve a non-linear least square problem to fit a "noisy"
>> histogram to a model that is a single exponential decay convoluted with a
>> "sharp" peak and a baseline noise floor.
>>
>> The problem is that unless I set constrains for the 'offset' parameter the
>> confidence interval calculation fails in a strange way in scipy brentq
>> function.
>>
>> You can see the full example notebook and the error message in this ipython
>> notebook:
>>
>> <a
>> href="http://nbviewer.ipython.org/urls/gist.githubusercontent.com/tritemio/901c2eb2a43775e81844/raw/755eede8f170cc33b96316601286ae534e901c49/lmfit%20issue">http://nbviewer.ipython.org/urls/gist.githubusercontent.com/tritemio/901c2eb2a43775e81844/raw/755eede8f170cc33b96316601286ae534e901c49/lmfit%20issue</a>
>>
>> What can I do to fix the problem? If this is the normal behavior, a more
>> informative error message would be helpful.

http://stackoverflow.com/questions/20619156/python-brentq-problems/20625362

I have no idea about the details in lmfit, but it's not an infrequent
problem also outside of lmfit.

Josef

>>
>> Thanks in advance for any help,
>> Antonio
>>
>> PS: Let me congratulate with the lmfit authors as this is a really amazing
>> piece of software that should be included in scipy IMHO!
>>
>> _______________________________________________
>> SciPy-User mailing list
>> SciPy-User at scipy.org
>> <a
>> href="http://mail.scipy.org/mailman/listinfo/scipy-user">http://mail.scipy.org/mailman/listinfo/scipy-user</a>
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user



More information about the SciPy-User mailing list