set.seed(18)
x <- rnorm(15)
y <- x + rnorm(15)
new <- seq(-3, 3, 0.5)
pred.w.plim <- predict(lm(y ~ x), newdata=data.frame(x=new), interval="prediction")
pred.w.clim <- predict(lm(y ~ x), newdata=data.frame(x=new), interval="confidence")
plot(y ~ x)
DrawBand(y = c(pred.w.plim[,2], rev(pred.w.plim[,3])),
x=c(new, rev(new)), col= SetAlpha("grey90", 0.5))
# passing y as matrix interface allows more intuitive arguments
DrawBand(y = pred.w.clim[, 2:3],
x = new, col= SetAlpha("grey80", 0.5))
abline(lm(y~x), col="brown")
Run the code above in your browser using DataLab