0

我想测试一个 2x3 因子设计并像这样对比变量

    library(lme4)
    library(emmeans)
    my.helmert = matrix(c(2, -1, -1, 0, -1, 1), ncol = 2)

    contrasts(Target3$mask) = my.helmert
    contrasts(Target3$length)

因此，对于掩码，我想将第一组与其他两组的平均值进行比较，然后在第二步中将第二组与第三组进行比较。

这在我的 LMM 中运行良好

Target3.2_TT.lmer = lmer(logTotalTime ~ mask*length+ (1+length|Subject) +(1|Trialnum), data = Target3)

掩码和长度之间存在显着的相互作用，这就是为什么我想看看这种效果并计算这样的事后测试（土耳其）：

emmeans(Target3.2_TT.lmer, pairwise ~ mask : length)

这对一个问题也很有效：现在我的对比消失了。文本计算所有掩码的差异，而不仅仅是 1 与 2 和 3 以及 2 与 3。是否有可能在事后测试中保留我的对比？

这是数据的样子：

> dput(Target3)
structure(list(mask = structure(c(2L, 1L, 2L, 3L, 1L, 2L, 3L, 
2L, 1L, 3L, 3L, 3L, 3L, 2L, 2L, 3L, 3L, 1L, 3L, 1L, 1L, 2L, 3L, 
3L, 2L, 1L, 3L, 2L, 3L, 2L), contrasts = structure(c(2, -1, -1, 
0, -1, 1), .Dim = c(3L, 2L), .Dimnames = list(c("keine Maske", 
"syntaktisch\n korrekt", "syntaktisch \n inkorrekt"), NULL)), .Label = c("keine Maske", 
"syntaktisch\n korrekt", "syntaktisch \n inkorrekt"), class = "factor"), 
    length = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
    2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 
    2L, 2L, 1L, 1L, 2L, 1L), .Label = c("kurzes \n N+1", "langes\n N+1"
    ), class = "factor"), logTotalTime = c(4.969813299576, 5.37989735354046, 
    5.14166355650266, 5.40717177146012, 5.27299955856375, 5.72358510195238, 
    5.4249500174814, 6.18001665365257, 5.67675380226828, 5.44241771052179, 
    5.66988092298052, 5.04985600724954, 5.78996017089725, 5.03043792139244, 
    5.92958914338989, 5.15329159449778, 6.11146733950268, 5.26269018890489, 
    5.17614973257383, 6.18001665365257, 6.03068526026126, 5.68697535633982, 
    5.17614973257383, 5.19849703126583, 5.29330482472449, 5.89989735358249, 
    5.73979291217923, 5.65599181081985, 5.94017125272043, 5.72031177660741
    )), .Names = c("mask", "length", "logTotalTime"), row.names = c(2L, 
4L, 6L, 8L, 9L, 11L, 13L, 15L, 16L, 18L, 20L, 22L, 27L, 29L, 
31L, 33L, 35L, 37L, 39L, 41L, 42L, 44L, 47L, 49L, 51L, 54L, 55L, 
57L, 59L, 61L), class = "data.frame")
4

1 回答 1

0

好吧，如果你要求成对比较，那就是你得到的，赫尔默特对比与成对比较不同。此外，Tukey（不是土耳其）方法仅适用于成对比较，不适用于其他类型的对比。

这里有一些可以尝试的东西，它可能会给你你想要的东西。

emm = emmeans(Target3.2_TT.lmer, 
    ~ mask | length)
contrast(emm, list(
    c1 = c(2, -1, -1)/2,
    c2 = c(0, 1, -1)),
  adjust = “mvt”)

contrasts这将独立于拟合模型时使用的任何参数化（即设置）。模型参数化影响模型矩阵的设置方式和系数的解释，但不影响来自emmeans或其亲属的结果。

于 2019-01-21T14:07:43.590 回答