0

我正在使用 lmfit 来查找适合的置信区间，但它经常给我一个错误：ValueError: f(a) and f(b) must have different signs

这是一个最小的工作示例（运行几次以获取错误）：

import lmfit
import numpy as np

def residual(p, X):
        a1, a2, t1, t2 = [i.value for i in p.values()]
        return a1*np.exp(-x/t1)+a2*np.exp(-x/t2)-y

if __name__ == '__main__':
    x = np.linspace(0.3,10,100)
    y = 3*np.exp(-x/2.)-5*np.exp(-x/10.)+0.2*np.random.randn(x.size)
    p = lmfit.Parameters()
    p.add_many(('a1', 5), ('a2', -5), ('t1', 2), ('t2', 5))
    mi = lmfit.minimize(residual, p, args=(x,))
    lmfit.printfuncs.report_fit(mi.params, show_correl=False)
    ci, trace = lmfit.conf_interval(mi, sigmas=[0.68,0.95], trace=True, verbose=False)
    lmfit.printfuncs.report_ci(ci)

为什么这样做？有解决办法吗？

谢谢

4

1 回答 1

0

我不知道具体的答案或解决方法。

我会将其报告为问题https://github.com/lmfit/lmfit-py/issues。

置信区间是通过找到轮廓似然等于一个值的点来计算的。这需要找到非线性函数的零或根。这种方法通常提供比依赖局部导数更准确的置信区间，例如 scipy 的 curve_fit。

为 scipy 的 brentq 指定适用于所有情况的边界通常很棘手。我的猜测是，对于您的示例，默认界限太紧了。

我对方法只有一个大概的了解，但不知道 lmfit 中的具体代码。

于 2013-12-17T03:26:37.327 回答