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')