r - 从多级序数 logit (clmm) 估计和绘制预测概率

Question

我正在使用 clmm 估计具有一个交互项的多级序数 logit。我得到的系数很有趣且具有统计意义，我现在想计算和绘制预测概率。在这样做时，考虑到我所有的自变量都是分类的可能是相关的。如何计算和绘制使用 clmm 估计的模型的预测概率？

我在网上查看了相对广泛的内容，但找不到计算和绘制 clmm 模型的预测概率的方法。“效果”包似乎也不适用于 clmm 。

这是我用来估计模型的代码：

mm3c = clmm(var1 ~ var3 + var4 + var6 + var2*var5 + var7 + var6
            +  (1|var8), data=my.data)

我的目标是估计和绘制 var2 和 var5 相互作用的概率。任何帮助将不胜感激 - 感谢您的考虑！

标签： rmixed-modelsmulti-levelordinalmultilevel-analysis