Question

我有一个表格 av.rt ，包含3个因素：主题（N = 28），奖励（奖励或未奖励），一致（全等或不一致）和数值 - 反应时间。我的表有108个观察所有主题/奖励/一致性组合。

aov_CSRA <- with(av.rt, aov(TTime.Shifted ~ Reward*Congruency + Error(Subject / (Reward*Congruency))))
summary(aov_CSRA)
library(lsr)
etaSquared( aov_CSRA )

我运行了一个rANOVA和一个函数来计算eta平方：

Error in etaSquared(aov_CSRA) : "x" must be a linear model object

然而，我收到错误：

ImageUploader

我该如何解决？

Answer 1

在快速谷歌搜索后，我发现还有另一个包含计算eta squared，包DescTools，函数EtaSq的函数的包。

library(DescTools)

aov_CSRA <- aov(TTime.Shifted ~ Reward*Congruency + Error(Subject / (Reward*Congruency)), data = av.rt)
summary(aov_CSRA)

EtaSq(aov_CSRA, type = 1, anova = FALSE)
#                       eta.sq eta.sq.part eta.sq.gen
#Reward            0.281762511   0.6430538 0.31363362
#Congruency        0.094438684   0.7177132 0.13281446
#Reward:Congruency 0.007180224   0.1280018 0.01151048

EtaSq(aov_CSRA, type = 1, anova = TRUE)
#                       eta.sq eta.sq.part eta.sq.gen         SS df         MS
#Reward            0.281762511   0.6430538 0.31363362 193165.658  1 193165.658
#Congruency        0.094438684   0.7177132 0.13281446  64743.569  1  64743.569
#Reward:Congruency 0.007180224   0.1280018 0.01151048   4922.489  1   4922.489
#                        SSE dfE         F            p
#Reward            107222.34  26 46.840119 2.895631e-07
#Congruency         25464.56  26 66.104912 1.304265e-08
#Reward:Congruency  33533.93  26  3.816574 6.158667e-02

它在第一次尝试时起作用请注意，我已将呼叫更改为aov一点。而不是with我使用aov参数data。我发现它使代码更具可读性。

部分eta平方误差

1 个答案: