Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/81.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 如何使用emmeans限制估计边际均值之间的比较次数_R_Nlme_Emmeans - Fatal编程技术网

R 如何使用emmeans限制估计边际均值之间的比较次数

R 如何使用emmeans限制估计边际均值之间的比较次数,r,nlme,emmeans,R,Nlme,Emmeans,我正在使用lme(作为分类协变量输入的时间)分析纵向激素数据,并希望测试研究组(对照组和两个实验组)之间的差异,以及使用emmeans可视化差异。我对组内或组间连续时间点之间的比较一点也不感兴趣。因此,我希望改变对比度,以便对多重比较的校正只考虑我所看到的那些比较。然而,即使在查阅了emmeans的小插曲之后,我似乎也不明白如何做到这一点 这是一个可重复的例子,它确实产生了组间的比较,但可能过度调整了p值。当然,这只是一个很小的问题,我不介意在这方面过于保守 ortho <- nlme::

我正在使用lme(作为分类协变量输入的时间)分析纵向激素数据,并希望测试研究组(对照组和两个实验组)之间的差异,以及使用emmeans可视化差异。我对组内或组间连续时间点之间的比较一点也不感兴趣。因此,我希望改变对比度,以便对多重比较的校正只考虑我所看到的那些比较。然而,即使在查阅了emmeans的小插曲之后,我似乎也不明白如何做到这一点

这是一个可重复的例子,它确实产生了组间的比较,但可能过度调整了p值。当然,这只是一个很小的问题,我不介意在这方面过于保守

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