x = read.table(file.path(pfad, "int.txt"))
attach(x)
f0ab = abs(G1f - Tf)
dab = abs(Tt - G1t)
f0an = abs(G2f - Tf)
dan = abs(Tt - G2t)
detach(x)

mab.f = with(x, tapply(f0ab, Satz, mean))
mab.t = with(x, tapply(dab, Satz, mean))

man.f = with(x, tapply(f0an, Satz, mean))
man.t = with(x, tapply(dan, Satz, mean))

par(mar=rep(2, 4))
par(mfrow=c(1,3))
plot(mab.t, mab.f, main="Abstieg", xlab="Dauer (s)", ylab="f0 (Hz)", type="n", ylim=c(50, 90), xlim = c(0.15, 0.4))
text(mab.t, mab.f, c("A", "B", "C", "D"), cex=3)
plot(man.t, man.f, main="Anstieg", ylab="", xlab="Dauer (s)", type="n", ylim=c(45, 100), xlim = c(.1, .35))
text(man.t, man.f, c("A", "B", "C", "D"), cex=3)
# gesamt Dauer vs. gesamt f0-Änderung
plot(mab.t + man.t, mab.f + man.f, ylab="", main="Abstieg+Anstieg", xlab="Dauer (s)",type="n", ylim=c(90, 170), xlim=c(.3, .7))
text(man.t+mab.t, man.f+mab.f, c("A", "B", "C", "D"), cex=3)



# Geschwindigkeit
gan = f0an/dan
gab = f0ab/dab

g = data.frame(Gan = gan, Gab = gab, Satz = with(x, Satz))
par(mfrow=c(1,2))
boxplot(gan ~ Satz, data = g, ylim = c(0, 800))
boxplot(gab ~ Satz, data = g, ylim= c(50, 600))