在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)

有关更多详细信息，请查看文档