以下代碼將爲您生成好看的迴歸圖。我在代碼中的評論應該清楚地解釋一切。該代碼將使用value
,model
在您的問題。
## all date you are interested in, 4 years with observations, 10 years for prediction
all_date <- seq(as.Date("2012-12-31"), by="1 year", length.out = 14)
## compute confidence bands (for all data)
pred.c <- predict(model, data.frame(date=all_date), interval="confidence")
## compute prediction bands (for new data only)
pred.p <- predict(model, data.frame(date=all_date[5:14]), interval="prediction")
## set up regression plot (plot nothing here; only set up range, axis)
ylim <- range(range(pred.c[,-1]), range(pred.p[,-1]))
plot(1:nrow(pred.c), numeric(nrow(pred.c)), col = "white", ylim = ylim,
xaxt = "n", xlab = "Date", ylab = "prediction",
main = "Regression Plot")
axis(1, at = 1:nrow(pred.c), labels = all_date)
## shade 95%-level confidence region
polygon(c(1:nrow(pred.c),nrow(pred.c):1), c(pred.c[, 2], rev(pred.c[, 3])),
col = "grey", border = NA)
## plot fitted values/lines
lines(1:nrow(pred.c), pred.c[, 1], lwd = 2, col = 4)
## add 95%-level confidence bands
lines(1:nrow(pred.c), pred.c[, 2], col = 2, lty = 2, lwd = 2)
lines(1:nrow(pred.c), pred.c[, 3], col = 2, lty = 2, lwd = 2)
## add 95%-level prediction bands
lines(4 + 1:nrow(pred.p), pred.p[, 2], col = 3, lty = 3, lwd = 2)
lines(4 + 1:nrow(pred.p), pred.p[, 3], col = 3, lty = 3, lwd = 2)
## add original observations on the plot
points(1:4, rev(value), pch = 20)
## finally, we add legend
legend(x = "topleft", legend = c("Obs", "Fitted", "95%-CI", "95%-PI"),
pch = c(20, NA, NA, NA), lty = c(NA, 1, 2, 3), col = c(1, 4, 2, 3),
text.col = c(1, 4, 2, 3), bty = "n")
的JPEG由代碼生成:
jpeg("regression.jpeg", height = 500, width = 600, quality = 100)
## the above code
dev.off()
## check your working directory for this JPEG
## use code getwd() to see this director if you don't know
正如你可以從圖中看到,
- 置信帶變寬了,當你試圖做預測遠離你觀察的數據;
- 預測間隔寬於置信區間。
如果您想了解更多關於predict.lm()
如何計算內部置信度/預測時間間隔,請閱讀How does predict.lm() compute confidence interval and prediction interval?,並在那裏回答我的答案。
感謝Alex的演示,簡單的使用了visreg
包;但我仍然喜歡使用R基地。
的可能的複製[如何計算斜率的95%置信區間在R中的線性迴歸模型中(http://stackoverflow.com/questions/15180008/how-to-calculate-the-95-confidence-interval-for-the-slope-in-a-linear- regressio) – majom
它不是上面提到的重複。 OP是要求一個置信區間而不是置信區間。 – Alex
@Zheyuan Li我也不會認爲它是重複的。首先,術語被使用得很奇怪。如果有人正在尋找它並使用正確的術語,他可能不會找到它。其次,我們在這裏有一個線性迴歸的特例,在其他答案中沒有提到。令人驚訝的是(對我來說),如果你尋找'[r]置信區間線性迴歸',你將不會得到好的結果。所以我認爲這個問題很好。但我也確信在SO上有很多解決方案。但是如果你找不到它們,他們又能怎麼幫忙? – Alex