R 如何使用emmeans限制估计边际均值之间的比较次数
我正在使用lme(作为分类协变量输入的时间)分析纵向激素数据,并希望测试研究组(对照组和两个实验组)之间的差异,以及使用emmeans可视化差异。我对组内或组间连续时间点之间的比较一点也不感兴趣。因此,我希望改变对比度,以便对多重比较的校正只考虑我所看到的那些比较。然而,即使在查阅了emmeans的小插曲之后,我似乎也不明白如何做到这一点 这是一个可重复的例子,它确实产生了组间的比较,但可能过度调整了p值。当然,这只是一个很小的问题,我不介意在这方面过于保守R 如何使用emmeans限制估计边际均值之间的比较次数,r,nlme,emmeans,R,Nlme,Emmeans,我正在使用lme(作为分类协变量输入的时间)分析纵向激素数据,并希望测试研究组(对照组和两个实验组)之间的差异,以及使用emmeans可视化差异。我对组内或组间连续时间点之间的比较一点也不感兴趣。因此,我希望改变对比度,以便对多重比较的校正只考虑我所看到的那些比较。然而,即使在查阅了emmeans的小插曲之后,我似乎也不明白如何做到这一点 这是一个可重复的例子,它确实产生了组间的比较,但可能过度调整了p值。当然,这只是一个很小的问题,我不介意在这方面过于保守 ortho <- nlme::
ortho <- nlme::Orthodont
ortho <-ortho %>%
mutate(Sex = as.character(Sex)) %>%
mutate(Sex = replace(Sex, Subject %in%
c("M10", "M11", "M12", "M13", "M14", "M15", "M16", "F01", "F02"),
"Other"))
fit <- lme(fixed = distance ~ factor(age) * Sex,
random = ~ 1 | Subject, data = ortho)
fit %>%
emmeans(trt.vs.ctrl ~ factor(age) * Sex)
ortho%
变异(性别=替换(性别,受试者百分比百分比)
c(“M10”、“M11”、“M12”、“M13”、“M14”、“M15”、“M16”、“F01”、“F02”),
“其他”))
适合%
艾美指数(trt.vs.ctrl~因子(年龄)*性别)
似乎我没有足够长的时间阅读@aosmith16的教程。这是由
fit %>%
emmeans(trt.vs.ctrl ~ Sex|factor(age)) %>%
confint()