r - 予測()の使い方

Question

値を予測したいのですが、これは明らかに解決策ではありません。複数選択テストを行っていますが、0.304... は答えではありません。predict() を正しく使用するには?

library(glm2)
data(crabs)
fit= glm(Satellites~Width,data=crabs, family="poisson")
plot(Satellites~Width,data=crabs)
abline(fit)
predict(fit, newdata=data.frame(Width=c(22)))
1 
0.3042347 

Answer 1

ポアソン回帰の関数predict()（一般にGLMの場合）は、デフォルトで線形予測子のスケール、つまりこの場合は対数スケールで値を計算します（のヘルプファイルを参照predict.glm）。

predict(fit, newdata=data.frame(Width=c(22)))
        1 
0.3042347 

応答変数のスケールで予測値を取得するには、type="response"関数に引数を追加する必要がありますpredict()。

predict(fit, newdata=data.frame(Width=c(22)),type="response")
       1 
1.355587