在 Stata 中,我知道如果使用以下命令,我可以获得因变量 ( thkbins
) 和两个预测变量 ( cc
& tv
) 之间每个可能组合的 logits:
melogit thkbins cc#tv || school:,
有没有办法在 R 中产生类似的输出?我一直在使用包中的glmer
命令lme4
,虽然我可以使用交互项获得输出,但这并不是我在 Stata 中可以产生的。
model1 <- glmer(thkbin ~ cc + tv + cc*tv + (1|school),
data=thkdata, family = binomial, nAGQ = 7)
summary(model1)