ede*_*esz 2 python curve-fitting lmfit
这是关于从lmfit fit_report()(1)对象提取拟合统计信息的问题
在此 lmfit示例中,返回以下部分输出:
[[Model]]
Model(gaussian)
[[Fit Statistics]]
# function evals = 31
# data points = 101
# variables = 3
chi-square = 3.409
reduced chi-square = 0.035
Akaike info crit = -336.264
Bayesian info crit = -328.418
.
.
.
.
.
.
Run Code Online (Sandbox Code Playgroud)
我试图提取该Fit Statistics部分中的所有数量作为单独的变量。
for key in fit.params:
print(key, "=", fit.params[key].value, "+/-", fit.params[key].stderr)
Run Code Online (Sandbox Code Playgroud)
但是,这仅给出了模型参数。它没有提供拟合统计参数,这也很有用。我似乎在文档中找不到此内容。
有没有类似的方式来提取拟合统计参数(chi-square,reduced chi-square,function evals分别,等等)?
结果包含所有适合的统计信息。您可以获取所需的参数,如下所示
result = gmodel.fit(y, x=x, amp=5, cen=5, wid=1)
# print number of function efvals
print result.nfev
# print number of data points
print result.ndata
# print number of variables
print result.nvarys
# chi-sqr
print result.chisqr
# reduce chi-sqr
print result.redchi
#Akaike info crit
print result.aic
#Bayesian info crit
print result.bic
Run Code Online (Sandbox Code Playgroud)