RNGkind("Mersenne-Twister") set.seed(1) X <- rnorm(100,10,1) Y <- X + runif(100,0,1) summary(res <- lm(Y~X)) XX <- data.frame(X=seq(min(X),max(X),length=20)) plim <- predict(res, XX, interval="prediction") clim <- predict(res, XX, interval="confidence") par(family="sans") plot(X,Y) matlines(XX,plim,col=1,lty=c(1,2,2)) matlines(XX,clim,col=1,lty=c(1,3,3))