使用R中的lsmeans / emmeans计算成对比较的置信区间

Question

我正在R中使用lsmeans/emmeans包来创建treatA级别（二进制/因子变量）之间响应的成对比较图。我可以使用lsmeans（对比度）获得差异估计，但它仅提供估计的SE，而不提供置信度限制。置信限仅针对单个效果，不针对对比。任何人都可以帮助生成均值差（对比度）的CL吗？

model = glmer.nb(response ~ treatA * treatB + (1|random), data, family=nbinom1)

lsm <- lsmeans(model, pairwise ~ treatA*treatB)

contrast(lsm, method = "pairwise")

我最终希望制作一个这样的图，但是用CL而不是SE。

Answer 1

confint(lsm[[2]])

plot(lsm[[2]], CIs = TRUE)

查看文档以了解更多详细信息