CLMM 事后多重比较，包括响应变量的水平

Posted 2023-02-18

【中文标题】CLMM 事后多重比较，包括响应变量的水平

【英文标题】：CLMM post hoc multiple comparisons including the levels of the response variable

【发布时间】：2021-11-19 17:28:05

【问题描述】：

我正在为具有 3 个序数水平的响应变量运行一些 clmm 模型

磁阻面 没有防水布 没有毯子 一些防水布 一些防水布 很多毯子 很多毯子

fmm1

emmeans(fmm1b,pairwise~Surface,adjust="tukey")

在摘要和 emmeans 中，我只得到了毯子与防水布的比较，但它从不包括我的可变响应水平之间的比较。

它为我提供了值，告诉我每个表面内发生的事情是否因表面而异。我得到一个 p 值，告诉我毯子与防水布不同，但这不是我感兴趣的。 enter image description here

但我想比较这里展示的每个表面的相同水平

enter image description here

我想得到一个 p 值，告诉我“无”不情愿不会因表面而异，其他级别也一样。

我知道如果我改变我的模型并将不情愿作为因变量并将表面作为响应变量，我可以获得这些结果。 fmm1

有没有办法使用我的初始模型进行这些比较？

【参考方案1】：

见vignette("models", "emmeans")。默认是潜在模式，但是你需要mode = "prob")

【讨论】：

感谢您的回答，但更改模式并不会更改我在 emmeans 摘要中得到的比较。它仍然不包括变量响应的水平。

另请参阅vignette("sophisticated") 在序数模型部分。那里的例子会告诉你你应该期待什么。

谢谢！ emmeans(fmm1, ~ reluctance | Surface, mode = "prob") 我设法找到了接近我想要的东西。

更新：使用 emmeans(fmm1,pairwise ~ Surface | reluctance2, mode = "prob")我设法得到了我正在寻找的东西！