3

回帰直線上の 1 点の CI を取得するにはどうすればよいですか? そのためにconfint()を使用する必要があると確信していますが、これを試してみると

confint(model,param=value)

入力した場合と同じ番号が表示されるだけです

confint(model)

値なしで試しても、まったく値が得られません。

私は何を間違っていますか?

4

1 に答える 1

21

predict()の代わりに必要ですconfint()。また、Joranが指摘したように、特定のxの信頼区間と予測区間のどちらが必要かを明確にする必要があります。(信頼区間は、特定のxでのy値の期待値に関する不確実性を表します。予測区間は、xの値を持つ単一のサンプリングされたポイントの予測y値を取り巻く不確実性を表します。)

これをRで行う方法の簡単な例を次に示します。

df <- data.frame(x=1:10, y=1:10 + rnorm(10))

f <- lm(y~x, data=df)

predict(f, newdata=data.frame(x=c(0, 5.5, 10)), interval="confidence")
#         fit       lwr       upr
# 1 0.5500246 -1.649235  2.749284
# 2 5.7292889  4.711230  6.747348
# 3 9.9668688  8.074662 11.859075

predict(f, newdata=data.frame(x=c(0, 5.5, 10)), interval="prediction")
#         fit       lwr       upr
# 1 0.5500246 -3.348845  4.448895
# 2 5.7292889  2.352769  9.105809
# 3 9.9668688  6.232583 13.701155
于 2012-10-04T17:37:58.560 に答える