为什么 lines() 不绘制置信区间

Question

我尝试绘制线性回归线的置信区间，但 lines() 函数不起作用。我应该怎么做才能解决这个问题？

Jahr<-Jahr[2:26]
Menge<-Menge[1:25]
plot(Jahr,Menge,xlim=c(1991,2016),ylim=c(25000000,550000000))

coef<-lm(Menge~Jahr)
abline(coef,col="blue")

PredictPoint<-predict(coef,data.frame(Jahr=c(2016)),interval="none")
points.default(x=2016,y=PredictPoint,type="p",col="red",pch=1)
conf<-predict(coef, data.frame(Jahr=c(2016)), interval = "confidence")
lines(data.frame(Jahr=c(2016)),conf[,2],col="red")
lines(data.frame(Jahr=c(2016)),conf[,3],col="red")

Answer 1

问题是，您只预测了 2016 年的一个数据点：

conf<-predict(coef, data.frame(Jahr=c(2016)), interval = "confidence")

lines怎么可能把一个点排成一条线呢？默认情况下，lines 有 type = "l"。如果你设置type = "b"，它会显示b个点和线。如果你只有一个数据，它会显示点。

另一个问题是您将数据框而不是向量传递给 lines 的第一个参数。

总而言之，做：

lines(2016, conf[,2], col="red", type = "b")
lines(2016, conf[,3], col="red", type = "b")

如果你只想要这个单点的置信区间，也许这样更好：

segments(2016, conf[,2], 2016, conf[,3], lty = 2, col = "red")

为什么 lines() 不绘制置信区间

Why does lines() not draw confidence interval

r

plot

regression

linear-regression

lm