R中带有置信区间图的线性回归预测

Question

说示例数据

library("robustbase")
data(education)

我创建回归模型

model=lm(Y~X1+X2+X3,data=education)

现在我需要绘制带有置信区间的预测值。 即，我想要这样的情节：

如何创建它？

model=lm(Y~X1+X2+X3,data=education)
plot(model, which = 1)

我不需要这个结果。

Answer 1

您可以使用lattice 包和mosaic 包，比如

library("lattice")
library(mosaic)
library(robustbase)

data(education)
mylm=lm(Y~X1+X2+X3,data=education)
pred <- predict(mylm, data=education)
df <- data.frame(Observed=education$Y, Predicted=pred)

xyplot(Predicted ~ Observed, data = df, pch = 19,  panel=panel.lmbands,
                  band.lty = c(conf =2, pred = 1))