R 中的置信區間 plot 的線性回歸預測

Question

說示例數據

library("robustbase")
data(education)

我創建回歸 model

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

現在我需要得到 plot 預測值與置信區間。 即，我想要這樣的 plot：

我怎樣才能創建它？

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

我不需要結果。

Answer 1

您可以使用lattice package 和mosaic package ，例如

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))