Question

我有一个心理学实验的数据，其中有一个对照组和一个实验组。两组均在反应时间（RTs）上进行了测试，以定位点在屏幕上出现的位置（左/右）。在出现该点之前，会为参与者展示两张图像，一幅在左侧，一幅在右侧。图像可以是快乐，悲伤或中性的。在三个不同的时间点（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评分为自变量，以三个时间点作为组内变量进行重复测量ANOVA。

>

我对应该如何在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个单元格...

这是怎么了？

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

0 个答案: