F2 = c(2300, 2212, 2005, 2010, 2440, 1920, 1855, 1761, 1880, 2010, 2401, 2415, 2308, 2100, 2520, 2420, 2332, 2505, 2210, 2325)
Alter = rep(rep(c("j", "a"), each=5), 2)
Dialekt = rep(c("M", "W"), each=10)
ac = data.frame(F2, Alter=factor(Alter), Dialekt=factor(Dialekt))

par(mfrow=c(1,2))
boxplot(F2 ~ Alter * Dialekt, data = ac)
with(ac, interaction.plot(Alter, Dialekt, F2))

ac.aov = aov(F2 ~ Alter * Dialekt, data=ac)
summary(ac.aov)
tk = TukeyHSD(ac.aov)
phsel(tk[[3]])
phsel(tk[[3]], 2)

Alter (F[1, 16] = 5.5, p < 0.05) und Dialekt (F[1,16] = 24.4, p < 0.001)
hatten einen signifikanten Einfluss auf F2 und es gab eine
signifikante (F[1,16] = 6.23, p < 0.05) Interaktion zwischen diesen Faktoren.
Post-hoc Tukey Tests zeigten signifikante F2-Unterschiede zwischen
den Altersgruppen nur für München (p < 0.05) jedoch nicht für Wien. Es
gab signifikante F2-Unterschiede zwischen den Dialekten nur für die alten
(p < 0.001) jedoch nicht für die jungen Sprecher.

drei = read.table(file.path(pfad, "drei.txt"))
boxplot(Response ~ Region * Alter, data=drei)
boxplot(Response ~ Region * Gender, data=drei)
drei.aov = aov(Response ~ Region * Alter * Gender, data = drei)
summary(drei.aov)

F2 wurde von Alter (F[2,48] = 28.2, p < 0.001) jedoch weder von Geschlecht
noch von der Region signifikant beeinflusst. Es gab auch signifikante Interaktionen zwischen Region und Alter (F[6, 48] = 4.9) sowie zwischen Alter und Gender (F[2,48] = 9.9, p < 0.001). 

tk = TukeyHSD(drei.aov)
names(tk)
phsel(tk[[4]])
phsel(tk[[4]], 2)

Post-Hoc Tukey-Tests zeigten signifikante Unterschiede  zwischen jung und alt für Region A (p < 0.01) und Region B (p < 0.05);  zwischen mittelalt und alt (p < 0.001) für Region D (p < 0.001); und zwischen mittelalt und jung für Region C (p < 0.05).

phsel(tk[[6]])
phsel(tk[[6]], 2)

Post-hoc Tukey-Tests zeigten ebenfalls signifikante Unterschiede zwischen alt und jung sowohl in Männern (p < 0.01) als auch in Frauen (p < 0.05); sie zeigten aber signifikante Unterschiede zwischen alt und mittalt  sowie zwischen mittelalt und jung nur in Frauen (p < 0.05; p < 0.01) aber nicht in Männern.