2016-07-22 5 views
1

私は線形回帰直線の信頼区間をプロットしようとしますが、lines()関数は機能しません。この問題を解決するにはどうすればよいですか?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") 
+0

:-)はいますがlinkesに感謝。私の問題は、それをプロットするのに成功していないということです。 – Kaja

答えて

1

問題は、あなたが唯一の2016年に一つのデータポイントを予測し、次のとおりです。

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

linesラインを生成するために単一の点を並べることができますか?既定ではlinestype = "l"です。 type = "b"を設定すると、b番目の点と線が表示されます。 1つのデータムしかない場合は、ポイントが表示されます。

もう1つの問題は、ベクトルではなくデータフレームを最初の引数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") 
+0

Thnksたくさんあります。あなたの解決策は私に2点しか与えませんが、私は2本の線が必要です – Kaja

関連する問題