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)