如何用can方法对lme4模型进行正交对比？

Question

我有一个多变量的lme4混合模型。我对模型中的一个变量的正交对比感兴趣，它有三个组(比如A，B，C)。我要把A和B和C的平均值进行比较，我怎么能用R中的do来做这件事呢？

Answer 1

以下是我想象中的一种可行的方法：

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 -.5x2-.5x3

这样做的结果如下：

contrast estimate   SE df t.ratio p.value
 cyl          5.16 1.37 28 3.781   0.0008 

这符合一个模型，使用emmeans开发边缘，然后在contrast函数中，我们查看了4柱体与6和8的平均值之间的差异。