x <- rnorm(15)
y <- x + rnorm(15)
fm=lm(y ~ x)
new <- data.frame(x = seq(-3, 3, 0.5))
predict(fm, new, interval="prediction")

多谢小王子！