求拟合函数的根和估计误差

Question

我对Python有点陌生，我想找出线性拟合函数包含错误的根源。

我用

import numpy as np
import scipy.optimize as op

def lin_fit(x, a, b):
    return a * x + b

...

popt, pconv = op.curve_fit(lin_fit, U, sqrt_I)
x = np.linspace(np.min(U), np.max(U), 100)
y = lin_fit(x, *popt)

然后我找到了根

^{pr2}$

这没问题，但我现在不知道如何找出这些根的错误。

有人能帮忙吗？

Answer 1

误差包含在pconv中，其中对角线包含每个估计系数的方差，在本例中是a和{}。在

根据MathWorld，测量拟合整体质量的相关系数r可计算为：

sxy = pconv[0,1]*x.shape[0]
sxx = pconv[0,0]
syy = pconv[1,1]
r = (sxy**2/(sxx*syy))**0.5