x = 1:6
n = length(x)
v = var(x)
v * (n-1)

m = mean(x)
ssx = sum((x-m)^2)

vokal = read.table(file.path(pfad, "vokal.txt"))
attach(vokal)
SSY = var(F2) * (length(F2) -1)
SSE = sum(tapply(F2, V, var) * 9)
SSR = SSY - SSE

MSE = SSE/18
MSR = SSR/1
Fratio = MSR/MSE
detach(vokal)

summary(aov(F2 ~ V, data=vokal))
# oder
anova(aov(F2 ~ V, data=vokal))
# das gleiche
reg = lm(F2 ~ V, data = vokal)
anova(reg)

t.test(F2 ~ V, var.equal=T, data = vokal)

vok = read.table(file.path(pfad, "vok.txt"))
names(vok)

with(vok, table(Vokal, Gen))
boxplot(F2 ~ Gen * Vokal, data = vok)
with(vok, interaction.plot(Vokal, Gen, F2))

vok.aov = aov(F2 ~ Vokal * Gen, data = vok)
# das gleiche
vok.aov = aov(F2 ~ Vokal + Gen + Vokal:Gen, data = vok)


anova(vok.aov)
# das gleiche
reg = lm(F2 ~ Vokal * Gen, data = vok)
anova(reg)

tk = TukeyHSD(vok.aov)
tk
source(file.path(pfad, "tksel.txt"))
tksel(tk[[3]], 1)
tksel(tk[[3]], 2)