library(ez) library(lattice) ssb = read.table(file.path(pfadu, "ssb.txt")) # Tabelle: wieviele Stichproben pro Vpn? with(ssb, table(Vpn, interaction(Wort, Alter))) # Über die Wiederholungen mitteln ssbm = aggregate(F2 ~ Wort * Alter * Vpn, mean, data = ssb) head(ssbm); dim(ssbm) # Tabelle: wieviele Stichproben pro Vpn? with(ssbm, table(Vpn, interaction(Wort, Alter))) # Abbildung bwplot(F2 ~ Alter | Wort, data = ssbm) # Statistik ezANOVA(ssbm, .(F2), .(Vpn), .(Wort), between = .(Alter)) # Greenhouse-Geisser-Korrektur 2 * 0.6860511 20 * 0.6860511 round(2 * 0.6860511, 1) round(20 * 0.6860511, 1) 1 - pf(9.8908882, 2 * 0.6860511, 20 * 0.6860511)