library(ez)


#### Beziehung zwischen gepaartem t-test und ANOVA
voice = read.table(file.path(pfad, "voice.txt"))
t.test(vot ~ Stimm, var.equal=T, data = voice)


t.test(vot ~ Stimm, var.equal=T, paired=T, data = voice)

# t-test (hier falsch)
par(mfrow=c(1,2))
boxplot(vot ~ Stimm, data = voice, main="falsch")
# gepaarter t-test
with(voice, boxplot(vot[9:16] - vot[1:8], main = "/ba-pa/ Unterschiede im selben Sprecher"))

ezANOVA(voice, .(vot), .(Vpn), .(Stimm))

dr = read.table(file.path(pfad, "dr.txt"))
par(mfrow=c(1,2))
boxplot(D ~ Dialekt * Position, data=dr)
with(dr, interaction.plot(Dialekt, Position, D))

# besser
temp = with(dr, Dialekt=="B" & Position == "final")
a = dr[temp,]
temp = with(dr, Dialekt=="B" & Position == "initial")
b = dr[temp,]
# Reihenfolge der Vpn prüfen
# alles OK, sonst b = b[m,]
m = match(a$Vpn, b$Vpn)
boxplot(a$D - b$D)



dr.ez = ezANOVA(dr, .(D), .(Vpn), .(Position), .(Dialekt))
dr.ez

# Lädt die Funktionen phoc() phsel()
source(file.path(pfad, "phoc.txt"))
dr.p = phoc(dr, .(D), .(Vpn), .(Dialekt, Position))
dr.p

phsel(dr.p[[1]])
phsel(dr.p[[1]], 2)

# Wiederholung in derselben Zelle
ssb = read.table(file.path(pfad, "ssb.txt"))

with(ssb, table(Wort, Vpn))

ssbm = with(ssb, aggregate(F2, list(Alter, Wort, Vpn), mean))
dim(ssbm)
names(ssbm) = c("Alter", "Wort", "Vpn", "F2")
ezANOVA(ssbm, .(F2), .(Vpn), .(Wort), .(Alter))

# Sphericity
1 - pf(9.8908882, 2 * 0.6860511, 20 * 0.6860511)