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

Question

我正在使用lsmeans/emmeansR 中的包在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)

查看文档以获取更多详细信息