R：绘制置信带问题

Question

我正在尝试使用函数绘制回归的置信带，但我的上下预测带的线条看起来很奇怪和模糊。谁能发现这个问题，因为我似乎找不到任何问题？

代码如下：

drawPrediction <- function(x, y) {

  lm.model <- lm(y ~ x)
  ci <- predict(lm.model, interval = "prediction")
  plot(x,y)
  abline(lm.model, col = "red")
  lines(x, ci[,2], lty = 2, col = 'blue')
  lines(x,ci[,3], lty = 2, col = 'blue' )

}

针对一些测试数据（R state.x77 数据，文盲收入回归）给出此结果：

这是我用于上面测试的代码：

> test <- as.data.frame(state.x77)
> drawPrediction(test$Illiteracy, test$Income)

Answer 1

问题是 x 没有排序，所以你在整个地方画线，而不仅仅是在连续的 x 之间。相反，使用

lm.model <- lm(Income ~ Illiteracy, data=as.data.frame(state.x77))
  ci <- predict(lm.model, interval = "prediction")
  plot(x,y)
  abline(lm.model, col = "red")
  lines(sort(x), ci[order(x),2], lty = 2, col = 'blue')
  lines(sort(x), ci[order(x),3], lty = 2, col = 'blue' )