pfad = "pfad wo anovaobjekte gespeichert wurde

attach(paste(pfad, "anovaobjekte", sep="/"))
library(car)

#### Beziehung zwischen gepaartem t-test und ANOVA
ba = c(10, -20, 5, -10, -25, 10, -5, 0)
pa = c(20, -10, 15, 0, -20, 16, 7, 5)
vot = c(ba, pa)
vot.l = factor(c(rep("ba", length(ba)), rep("pa", length(pa))))
t.test(vot ~ vot.l, var.equal=T)

t.test(vot ~ vot.l, var.equal=T, paired=T)

Sprecher = factor(rep(1:8, 2))

# RM-ANOVA
vot.aov = aov(vot ~ vot.l + Error(Sprecher/vot.l))
summary(vot.aov)

# RM-MANOVA
vdaten = data.frame(vot, Sp= factor(Sprecher), Voice = factor(vot.l))
code = c("d", "s", "w")
vdaten.t = Anova.prepare(vdaten, code)
vdaten.lm = lm(vdaten.t$d ~ 1)
vdaten.aov = Anova(vdaten.lm, idata=vdaten.t$w, idesign=~Voice)
vdaten.aov

summary(vdaten.aov, mult=F)

#######################################

dr
names(dr)
attach(dr)

boxplot(D ~ Dialekt * Position)
interaction.plot(Dialekt, Position, D)

# Befehl fuer RM-MANOVA
code = c("d", "b", "s", "w")
dr.t = Anova.prepare(dr, code)
Dialekt = factor(dr.t$b)
dr.lm = lm(dr.t$d ~ Dialekt)
dr.aov = Anova(dr.lm, idata = dr.t$w, idesign= ~Position)
dr.aov


# Post-hoc
phoc = phoc.prepare(dr.t)
labs = paste(phoc$fac[,1], phoc$fac[,2], sep=".")
pairwise.t.test(phoc$d, labs, paired=T, p.adj="bonf")

# Bestaetigen. 
# Wir vergleichen SH-initial mit SH-final
temp = dr.t$b == "SH"
t.test(dr.t$d[temp,1], dr.t$d[temp,2], paired=T)
0.06192 * 6

# Wir vergleichen SH und B in initialer Position
t.test(dr.t$d[,2] ~ factor(dr.t$b), paired=T)
0.0002315 * 6