# Aufgabe aus anova1.ppt
form = c(2300, 2212, 2005, 2010, 2440, 1920, 1855, 1761, 1880, 2010,
2401, 2415, 2308, 2100, 2520, 2420, 2332, 2505, 2210, 2325)
dialekt = c(rep("M", 10), rep("W", 10))
dialekt = factor(dialekt)
alter = c(rep("j", 5), rep("a", 5), rep("j", 5), rep("a", 5))
alter = factor(alter)
form.aov = aov(form ~ dialekt * alter)
summary(form.aov)

TukeyHSD(form.aov)


# Aufgabe aus anova2.ppt  2
rt = c(121, 192, 110, 130, 180, 95, 88, 54, 78, 62, 92, 57, 75, 71, 70, 91, 72, 61, 69, 58)
sprache = factor(rep(c(rep("F", 5), rep("E", 5)), 2))
aufent = factor(c(rep(0, 10), rep(6, 10)))
vpn = factor(rep(1:10, 2))


rt.aov = aov(rt ~ sprache * aufent + Error(vpn/aufent))
summary(rt.aov)
# Aufenthaltsdauer hat eine signifikante Wirkung
# auf die Reaktionszeiten (F(1, 8) = 13.69, p < 0.01)
# und es gab eine signifikante Interaktion
# zwischen Aufenthaltsdauer und Sprache (F(1, 8) = 10.32,
p < 0.05).

Tukey.rm(rt, 566.9, 8, sprache, aufent)
# Post-hoc Tukey Tests
# zeigten: Aufenthaltsdauer hat einen signifikanten Einfluss
# auf die Reaktionszeiten der Franzosen, nicht der Englaender.