R +多级Logistic回归+组固定效应

Question

在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)

Answer 1

我会使用clmm包中的ordinal（教程here）：

model<-clmm(DepVar~IndVar+(1|WithinVar),data=df)

我希望这有帮助。