python拟合显示置信区间_python – StatsModels的置信度和预测间隔

iv_l,iv_u给出每个点的预测间隔的限制.

预测间隔是观测的置信区间,包括误差估计.

我认为,平均预测的置信区间在statsmodels中还不可用.

(实际上,拟合值的置信区间隐藏在effects_outlier的summary_table内,但我需要验证这一点.)

统计模型的正确预测方法在TODO列表中.

加成

有信心的间隔是OLS,但访问有点笨拙.

在运行你的脚本后包括在内：

from statsmodels.stats.outliers_influence import summary_table

st, data, ss2 = summary_table(re, alpha=0.05)

fittedvalues = data[:,2]

predict_mean_se = data[:,3]

predict_mean_ci_low, predict_mean_ci_upp = data[:,4:6].T

predict_ci_low, predict_ci_upp = data[:,6:8].T

# check we got the right things

print np.max(np.abs(re.fittedvalues - fittedvalues))

print np.max(np.abs(iv_l - predict_ci_low))

print np.max(np.abs(iv_u - predict_ci_upp))

plt.plot(x, y, 'o')

plt.plot(x, fittedvalues, '-', lw=2)

plt.plot(x, predict_ci_low, 'r--', lw=2)

plt.plot(x, predict_ci_upp, 'r--', lw=2)

plt.plot(x, predict_mean_ci_low, 'r--', lw=2)

plt.plot(x, predict_mean_ci_upp, 'r--', lw=2)

plt.show()