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
Nunca estuve en una situación en la que necesitaba ajustar un modelo lineal generalizado con codificación de efectos ( contr.sum
para R
usuarios). ¿Puedo aplicar la misma interpretación que en el caso del modelo lineal? En un modelo lineal normal de la intersección sería la gran media y la s (parámetros para , , y los efectos, es decir, cómo los niveles de factor se desvían de la gran media.period1
period2
period3
period4 = (Intercept) - period1 - period2 - period3
Así es como creo que va la interpretación análoga para los modelos lineales generalizados. (Expondré todos los parámetros y, por lo tanto, transformaré el log-odds (-ratios) en odds (-ratios).) La intercepción sería la probabilidad general de éxito frente a fracaso ( pegando aquí a la terminología binomial clásica) y las s las log-odds-ratios . Y obtenemos las probabilidades para, por ejemplo, agregando y luego : . ¿Es el realmente las probabilidades generales / medias y el speriod1
odds-ratios ?