計算R中的Eta Squared用於重復測量

Question

我正在使用aov()函數使用R執行重復測量ANOVA，例如在http://personality-project.org/r/r.anova.html網站上。

datafilename <- "http://personality-project.org/r/datasets/R.appendix4.data"
data.ex4 <- read.table(datafilename, header=TRUE)   
data.ex4                                      

計算標准ANOVA允許使用etaSquared()函數從lsr包或者可替換地etasq()從函數heplots ：

aov.standard=aov(Recall~(Task*Valence),data.ex4 ) 
etaSquared(aov.standard)
etasq(aov.standard)

但是，在實現主體內模型時，兩個功能都無法識別線性對象模型。

aov.within=aov(Recall~(Task*Valence)+Error(Subject/(Task*Valence)),data.ex4 )
etaSquared(aov.within)
etasq(aov.within)

有人有一個簡單的解決方案嗎？ 提前致謝。

Answer 1

當您使用ez包中的ezANOVA函數時，將自動計算eta平方：

library(ez)

m1 <- ezANOVA(data = data.ex4, dv = Recall, wid = Subject, within = .(Task,Valence), type = 3, detailed = TRUE)
m1$ANOVA

這給出了以下輸出：

        Effect DFn DFd    SSn       SSd          F           p p<.05         ges
1  (Intercept)   1   4 4177.2 349.13333 47.8579339 0.002291098     * 0.910303347
2         Task   1   4   30.0  16.33333  7.3469388 0.053508297       0.067934783
3      Valence   2   8    9.8  26.86667  1.4590571 0.288250073       0.023255814
4 Task:Valence   2   8    1.4  19.26667  0.2906574 0.755343673       0.003389831

ges列具有您需要的廣義eta平方值。