# Example 1:
# create dataframe
dat <- createDF(nVP = 50, nTrl = 1,
design = list("Comp" = c("comp", "incomp")))
dat <- addDataDF(dat, RT = list("Comp comp" = c(500, 100, 100),
"Comp incomp" = c(520, 100, 100)))
aovRT <- aov(RT ~ Comp + Error(VP/(Comp)), dat)
aovDispMeans(aovRT)
# or with ezANOVA
library(ez)
aovRT <- ezANOVA(dat, dv=.(RT), wid = .(VP), within = .(Comp),
return_aov = TRUE, detailed = TRUE)
aovRT <- aovTable(aovRT)
aovDispMeans(aovRT)
Run the code above in your browser using DataLab