如何在python中计算正交距离回归线拟合的置信区间

Owe*_*wen 5 python linear-regression confidence-interval

我正在使用正交距离回归方法(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
Run Code Online (Sandbox Code Playgroud)

现在我得到了 solpe,得到拟合函数:y=p*x,但是我应该如何得到 95% 置信区间?