Question

我一直在使用emmeans软件包为我的数据创建估计的边缘平均值，置信水平为0.95％。尽管我似乎无法将其更改为.99％的置信度。任何帮助将非常感激。通常我会使用“ levels =”函数，但对于emmeans来说似乎不存在。

library(emmeans)

emmeans(AcidLevels1, specs=~MAScore,)

最好， -内森（Nathan）

Answer 1

emmeans提供了方法confint.emmGrid重新计算置信区间，并且（可能更重要的是）还针对多种假设检验进行了调整。

由于您没有提供示例数据，因此以下是使用warpbreaks数据的示例。

library(emmeans)
lm <- lm(breaks ~ wool * tension, data = warpbreaks)
emm <- emmeans(lm, ~ wool | tension);

要重新计算99％的置信区间（无需多次校正）

confint(emm, adjust = "none", level = 0.99)
#tension = L:
# wool   emmean       SE df lower.CL upper.CL
# A    44.55556 3.646761 48 34.77420 54.33691
# B    28.22222 3.646761 48 18.44086 38.00358
#
#tension = M:
# wool   emmean       SE df lower.CL upper.CL
# A    24.00000 3.646761 48 14.21864 33.78136
# B    28.77778 3.646761 48 18.99642 38.55914
#
#tension = H:
# wool   emmean       SE df lower.CL upper.CL
# A    24.55556 3.646761 48 14.77420 34.33691
# B    18.77778 3.646761 48  8.99642 28.55914
#
#Confidence level used: 0.99

要使用Bonferroni校正在99％的水平上重新计算CI，并针对多种假设检验进行校正，

confint(emm, adjust = "bonferroni", level = 0.99)
#tension = L:
# wool   emmean       SE df lower.CL upper.CL
# A    44.55556 3.646761 48 33.82454 55.28657
# B    28.22222 3.646761 48 17.49120 38.95324
#
#tension = M:
# wool   emmean       SE df lower.CL upper.CL
# A    24.00000 3.646761 48 13.26898 34.73102
# B    28.77778 3.646761 48 18.04676 39.50880
#
#tension = H:
# wool   emmean       SE df lower.CL upper.CL
# A    24.55556 3.646761 48 13.82454 35.28657
# B    18.77778 3.646761 48  8.04676 29.50880
#
#Confidence level used: 0.99
#Conf-level adjustment: bonferroni method for 2 estimates

如何设置“平均数”中的置信区间

1 个答案: