temp = dip.l %in% c("aU", "OY") & dip.spkr == "67"
d = dur(dip[temp,])
lab = dip.l[temp]
boxplot(d ~ lab)

tapply(d, lab, shapiro.test)
var.test(d ~ lab)
t.test(d ~ lab, var.equal=T)


temp = timevow.l %in% c("I", "a")
d = dur(timevow[temp,])
lab = timevow.l[temp]
boxplot(d ~ lab)

tapply(d, lab, shapiro.test)
var.test(d ~ lab)
t.test(d ~ lab)

temp = polhom.l %in% c("c", "x")
d = dur(polhom[temp,])
lab = polhom.l[temp]
boxplot(d ~ lab)

tapply(d, lab, shapiro.test)
wilcox.test(d ~ lab)

##########################################################
temp = vowlax.word %in% c("bin", "ist", "nicht", "Kinder", "richten", "findet", "finden", "binden")
lab = vowlax.word[temp]
lab[lab %in% c("bin", "ist", "nicht")] = "F"
lab[lab != "F"] = "I"
table(lab)
d = dur(vowlax[temp,])
boxplot(d ~ lab)

##########################################################
on = dcut(isol.fdat[,2], 0, prop=T)
fin = dcut(isol.fdat[,2], 1, prop=T)
boxplot(on - fin)
shapiro.test(on)
shapiro.test(fin)
var.test(on, fin)
t.test(on, fin, paired=T)

par(mfrow=c(1,2))
dplot(isol.fdat[,2])
dplot(isol.fdat[,2], offset=1)
##########################################################
temp = vowlax.l == "I" & vowlax.left %in% c("b", "d")
d = dcut(vowlax.fdat[temp,2], 0, prop=T)
lab = vowlax.left[temp]
boxplot(d ~ lab)
tapply(d, lab, shapiro.test)
var.test(d ~ lab)
t.test(d ~ lab, var.equal=T)

##########################################################
temp = vowlax.left %in% c("k", "g")
d = dcut(vowlax.fund[temp,], start(vowlax[temp,])+15)
lab = vowlax.left[temp]
boxplot(d ~ lab)
temp = lab=="g"
hist(d[temp])
hist(d[!temp])
temp = d > 20
d  = d[temp]
lab = lab[temp]
boxplot(d ~ lab)
tapply(d, lab, shapiro.test)
var.test(d ~ lab)
t.test(d ~ lab, var.equal=T)
t.test(d ~ lab, alt="less", var.equal=T)
##########################################################