r - 重复测量方差分析 - ezANOVA 失败

Question

我有一个心理学实验的数据，其中有一个控制组和一个实验组。两组都在反应时间 (RT) 上进行测试，以确定屏幕上出现点的位置（左/右）。在点出现之前，参与者会看到两张图像，一张在左边，一张在右边。这些图像可以是快乐的、悲伤的或中性的。在三个不同的时间点（T0、T1 和 T2）重复测量：

     Emotion    Location     Time Group    RT     Score.D
1     Sad       Incongruent   T0    XP     625   22.83693
2     Happy     Congruent     T1    CTRL   783   22.83693
3     Sad       Congruent     T2    XP     578   22.83693
4     Happy     Incongruent   T0    XP     780   22.83693
5     Sad       Incongruent   T2    CTRL   566   22.83693
6     Happy     Congruent     T1    XP     535   22.83693

这是我可以直接在 R 中导入的数据的子集：

structure(list(Emotion = structure(c(2L, 2L, 2L, 2L, 2L, 2L), .Label = c("Happy", 
"Sad"), class = "factor"), Location = structure(c(2L, 2L, 2L, 
2L, 2L, 2L), .Label = c("Congruent", "Incongruent"), class = "factor"), 
    Time = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = c("T0", 
    "T1", "T2"), class = "factor"), Group = structure(c(2L, 2L, 
    2L, 2L, 2L, 2L), .Label = c("CTRL", "XP"), class = "factor"), 
    RT = c(625L, 783L, 578L, 780L, 566L, 535L), Score.D = c(22.83693, 
    22.83693, 22.83693, 22.83693, 22.83693, 22.83693)), row.names = c(NA, 
6L), class = "data.frame")

我的主管告诉我进行重复测量方差分析，其中治疗类型（对照/实验）作为组间因素，D 分数作为自变量，三个时间点作为组内变量。

我对如何在 R 中实现这一点感到困惑。这是我的尝试：

Incong.Cong.XP.CTRL.T0T1T2.formatted$Group <- factor(Incong.Cong.XP.CTRL.T0T1T2.formatted$Group)
Incong.Cong.XP.CTRL.T0T1T2.formatted$Time <- factor(Incong.Cong.XP.CTRL.T0T1T2.formatted$Time)
Incong.Cong.XP.CTRL.T0T1T2.formatted$Emotion <- factor(Incong.Cong.XP.CTRL.T0T1T2.formatted$Emotion)
rANOVA.BA.ech1 <- with(data = Incong.Cong.XP.CTRL.T0T1T2.formatted, aov(Score.D ~ Time*Emotion + Error(Group/(Time*Emotion))))
summary(rANOVA.BA.ech1)

这行得通，但我无法计算 Mauchly 的球形度测试，这是我非常想做的。该软件包ez允许 rANOVA 与 Mauchly 的测试一起使用。这是我尝试过的代码：

ezANOVA(Incong.Cong.XP.CTRL.T0T1T2.formatted, 
        dv = .(Score.D), 
        wid = .(generated_uid),
        within = .(Time,Emotion), 
        between = .(Group), type = 1)

但我收到此错误消息：

一个或多个单元格缺少数据。尝试使用 ezDesign() 检查您的数据。

然而数据中没有NA单元格......

这里有什么问题？

标签： ranova