我正在使用正交距离回归法(scipy.odr)来拟合我的数据,在拟合之后,我在计算95%置信区间时遇到了麻烦,请帮我没办法计算它〜
这里是代码:
#define fit function
def def f(p,x):
return p*x
linear = scipy.odr.Model(f)
mydata = scipy.odr.RealData(x, y, sx=dx, sy=dy)
myodr = scipy.odr.ODR(mydata, linear, beta0=[1.])
myoutput = myodr.run()
myoutput.pprint()
[solpe] = myoutput.beta
现在我得到了解,得到了拟合函数:y = p * x,但我应该如何获得95%置信区间?