# Ist die proportionale /i: - u:/ Entfernung kleiner
# für gespanntes /y:/ als ungespanntes /Y/?

efun <- function(a, b) 
{ 
sqrt(sum((a - b)^2)) 
}


# 1. Segmentliste
vok.s = emu.query("florian", "*", "phonetic != n") 

# 2. Label-Vektor
vok.l = label(vok.s)

# 3. Formanten-Trackdatei
vok.fm = emu.track(vok.s, "fm")

# 4. Formanten zum zeitlichen Mittelpunkt
vok.m = dcut(vok.fm, 0.5, prop=T)
eplot(vok.m[,1:2], vok.l, dopoints=T, doellipse=F, form=T)
vok.fmb = bark(vok.fm)

# F1-F2 Werte zum zeitlichen Mittelpunkt in Bark
vok.b = dcut(vok.fmb[,1:2], 0.5, prop=T)
# F1 x F2 Abbildung in Bark
eplot(vok.b[,1:2], vok.l, dopoints=T, doellipse=F, form=T)

# zentroid
zen.i = apply(vok.b[vok.l=="i:",], 2, mean)
zen.u = apply(vok.b[vok.l=="u:",], 2, mean)
zen.I = apply(vok.b[vok.l=="i",], 2, mean)
zen.U = apply(vok.b[vok.l=="u",], 2, mean)

mat = NULL
lab = c("i:", "u:", "i", "u")
for (j in lab){
print(j)
zen = apply(vok.b[vok.l==j,], 2, mean)
mat = rbind(mat, zen)
}
rownames(mat) = lab


# Entfernungen zu i:
di = apply(vok.b, 1, efun, zen.i)

# Entfernungen zu u:
du = apply(vok.b, 1, efun, zen.u)

# Entfernungen zu I
dI = apply(vok.b, 1, efun, zen.I)

# Entfernungen zu U
dU = apply(vok.b, 1, efun, zen.U)

# log Euk. Entfernungen gespannt
ent.g = log(di/du)

# log Euk. Enfernungen, ungespannt
ent.u = log(dI/dU)

# Daten für /y:/
ent.y = ent.g[vok.l=="y:"]

# Daten für /Y/
ent.Y = ent.u[vok.l=="y"]

mean(ent.y)
mean(ent.Y)