如何用emmeans创建手动对比？ -R

Question

假设我有这些数据

library(MASS)
m<-lmer(Y~N*V + (1|B),data=oats)

如何在emmeans中创建手动对比度？例如

Victoria_0.2cwt    1
Victoria_0.4cwt   -1
Marvellous_0.2cwt -1
Marvellous_0.4cwt  1

Answer 1

emm = emmeans(m, ~ V * N)
emm
contrast(emm, list(con = c(0,0,0,0,-1,1,0,0,-1,0,0,0)))

但是，这实际上是线性函数，而不是对比度，因为系数的总和不为零。

注意：我可能会误记了因子水平，如果这样，可能需要重新排列系数。如果第二行，它们应该与您在结果中看到的组合相对应