我有一个带有多个变量的lme4混合模型。我对模型中一个变量的正交对比感兴趣,该模型有三组(比如A、B、C(。我想将A与B和C的平均值进行比较。如何使用R中的emmeans来实现这一点?
这里有一种我认为可行的方法:
mtcars$cyl <-as.factor(mtcars$cyl)
fit <- lm(mpg ~ wt + cyl, mtcars)
library(emmeans)
fit
emm_fit <- emmeans(fit, specs = "cyl")
contrast(emm_fit, list(cyl = c(1, -.5, -.5)))
这相当于x_1-.5x_2-.5x_3
结果如下:
contrast estimate SE df t.ratio p.value
cyl 5.16 1.37 28 3.781 0.0008
这符合一个模型,使用emmeans来开发边缘,然后在contrast
函数中,我们从6和8的平均值来看4个圆柱体之间的差异。有关更多详细信息,请参阅https://cran.r-project.org/web/packages/emmeans/vignettes/comparisons.html#contrasts