※上記の広告は60日以上更新のないWIKIに表示されています。更新することで広告が下部へ移動します。

#predict
.glmについて
#作成した予測式の応答変数や説明変数に
$を入れるとpredict.glmが上手く動かない

#ダメな例
glm.res <- glm(dm$adductor~dm$temp4+I(log(dm$chl4))+dm$soya4+dm$adductor4,family=Gamma(link="log"), data=dm)

#良い例
ls <- glm(adductor~temp4+I(log(chl4))+soya4+adductor4,family=Gamma(link="log"), data=dm)

#予測の仕方
new.d <- data.frame(temp4=2.6, chl4=2.48, soya4=181, adductor4=12.50)
pred <- predict(ls, newdata=new.d, type="response")
|