模型预测的置信区间
问题内容:
我正在关注statsmodels教程
装有OLS模型
formula = 'S ~ C(E) + C(M) + X'
lm = ols(formula, salary_table).fit()
print lm.summary()
通过以下方式提供预测值:
lm.predict({'X' : [12], 'M' : [1], 'E' : [2]})
结果以单值数组形式返回。
是否有一种方法还可以返回统计模型中预测值的置信区间(预测区间)?
谢谢。
问题答案:
我们一直想让它更容易实现。您应该可以使用
from statsmodels.sandbox.regression.predstd import wls_prediction_std
prstd, iv_l, iv_u = wls_prediction_std(results)
如果您有任何问题,请在github上提交问题。