предсказания линейной регрессии с доверительным интервалом в 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))