# Example 1:
# create dataframe with 2(Comp: comp vs. incomp) and 2(Side: left vs. right) factors/levels
dat <- createDF(nVP = 20, nTrl = 1,
design = list("Comp" = c("comp", "incomp", "neutral"),
"Side" = c("left", "right")))
dat <- addDataDF(dat,
RT = list("Comp:Side comp:left" = c(500, 150, 150),
"Comp:Side comp:right" = c(500, 150, 150),
"Comp:Side incomp:left" = c(550, 150, 150),
"Comp:Side incomp:right" = c(550, 150, 150),
"Comp:Side neutral:left" = c(525, 150, 150),
"Comp:Side neutral:right" = c(525, 150, 150)))
aovRT <- aov(RT ~ Comp * Side + Error(VP/(Comp*Side)), dat)
aovDispMeans(aovRT)
aovRT <- aovEffectSize(aovRT)
aovRT <- aovDispTable(aovRT)
# or with ezANOVA
library(ez)
aovRT <- ezANOVA(dat, dv=.(RT), wid = .(VP), within = .(Comp, Side),
return_aov = TRUE, detailed = TRUE)
aovRT <- aovEffectSize(aovRT)
aovDispTable(aovRT)
Run the code above in your browser using DataLab