我不相信有一种方法可以指定要用于计算预测标准误差的协方差矩阵。请注意,预测代码仍在 statsmodels 存储库的“沙盒”文件夹中。我相信 Github 拉取请求会受到欢迎:)
无论如何,这应该很简单。这是您链接到的预测功能的底层代码的链接。本质上,您只需要替换要使用的协方差矩阵而不是covb
变量。
然后,您可以使用您在其他 SO 帖子中看到的相同 matplotlib 花絮。
https://github.com/statsmodels/statsmodels/blob/master/statsmodels/sandbox/regression/predstd.py#L27
predvar = res.mse_resid/weights + (exog * np.dot(covb, exog.T).T).sum(1)
predstd = np.sqrt(predvar)
tppf = stats.t.isf(alpha/2., res.df_resid)
interval_u = predicted + tppf * predstd
interval_l = predicted - tppf * predstd
return predstd, interval_l, interval_u