vokal = read.table(file.path(pfad, "vokal.txt"))
table(vokal$V)

summary(aov(F2 ~ V, data=vokal))

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

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

with(vok, table(Vokal, Gen))
par(mfrow = c(1,2))
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)
summary(vok.aov)



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