library(lme4)
out <- glmer(cbind(incidence, size - incidence)
~ period
+ (1 | herd),
data = cbpp,
family = binomial,
contrasts = list(period = "contr.sum"))
summary(out)
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.32337 0.22129 -10.499 < 2e-16 ***
period1 0.92498 0.18330 5.046 4.51e-07 ***
period2 -0.06698 0.22845 -0.293 0.769
period3 -0.20326 0.24193 -0.840 0.401
私は私が(コーディング効果と一般化線形モデルにフィットするために必要な状況ではなかったcontr.sum
ためR
、ユーザーが)。線形モデルの場合と同じ解釈を適用できますか?通常の線形モデルでは、切片は総平均と s(、、およびのパラメーター)であり、因子レベルが総平均からどのように逸脱するかを示します。period1
period2
period3
period4 = (Intercept) - period1 - period2 - period3
これが、一般化線形モデルの類似の解釈がどのように行われるかを私が考える方法です。(私はすべてのパラメーターをべき乗するので、log-odds(-ratios)をodds(-ratios)に変換します。)切片は、成功と失敗の全体的なオッズになります(ここでは古典的な二項用語に固執します)、はlog-odds-ratiosです。そして、例えばを指数をすることにより、オッズを取得し:。ある本当に全体的な/中オッズと sがperiod1
オッズ比?