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

Question

在Stata中，我知道如果使用以下命令，则可以获取因变量（ thkbins ）和两个预测变量（ cc ＆ tv ）之间每种可能组合的对数：

melogit thkbins cc#tv || school:,

有没有办法在R中产生类似的输出？ 我一直在使用lme4软件包中的glmer命令，尽管我可以使用交互项获得输出，但这并不是我可以在Stata中产生的。

model1 <- glmer(thkbin ~ cc + tv + cc*tv + (1|school),
            data=thkdata, family = binomial, nAGQ = 7)
summary(model1)

Answer 1

我会用clmm从包ordinal （教程这里 ）：

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

我希望这有帮助。