要计算我想使用的线性回归(使用scipy)的斜率误差
slope, intercept, r_value, p_value, std_err = stats.linregress(x,y)
sd_slope = (len(x)*std_err**2)/(len(x)*(x**2).sum()-x.sum()**2)
。该方程式取自维基百科。为什么这么错?
答案 0 :(得分:0)
最后要计算的是估计斜率的方差,即斜率的std误差的平方。事实证明,std_err
返回的.linregress
已经是斜率的标准误差。
但是,如果您仍然想手动计算,则需要用std_err
(s
(其中{{ 1}},样本大小。因此,在代码中:
sqrt(SSE / (n-2))