与emmeans的成对比较

Question

我对emmeans的输出有疑问。我进行了MCMCglmm，后跟emmeans进行了成对比较。输出包含“ estimate”，“ lower.HPD”和“ upper.HPD”。这里的估计是多少？是日志吗？

谢谢！ 仁

Answer 1

MCMCglmm 使用贝叶斯方法，因此， emmeans 默认情况下会产生贝叶斯摘要，显示每个EMM的后中位数和最高后密度间隔。如果您想要（可能更熟悉的）常客摘要，请使用类似

emm = emmeans(...)
summary(emm, freq = TRUE)
summary(pairs(emm), freq = TRUE)

如果您为模型配备了链接功能，则汇总将处于链接功能的范围内（例如，对数），除非您还向这些type = “response”调用中添加了summary。

希望有帮助。