segs = emu.query("ema5", "*", "Segment=aI")
segs.fm = emu.track(segs, "fm")[,1]
fsmooth = trapply(segs.fm, dct, 2, T, returntrack=T)
dplot(fsmooth, ylab="F1 (Hz)", xlab="Dauer (ms)", leg=F, bty="n", lwd=2)

peakfun <- function(fr, maxtime=T)
{
if(maxtime) 
num = which.max(fr)
else num = which.min(fr)
tracktimes(fr)[num]
}

times = trapply(fsmooth, peakfun, simplify=T)
dplot(fsmooth, offset=times, prop=F, xlab="Dauer (ms)", leg=F, bty="n", lwd=2)

jaw = emu.track(segs, "jw_posz")
jmax = dcut(jaw, times)
fmax = dcut(fsmooth, times)
cor.test(fmax, jmax)
plot(jmax, fmax)


llip = emu.track(segs, "ll_posz")
ulip = emu.track(segs, "ul_posz")

llipmax = dcut(llip, times)
ulipmax = dcut(ulip, times)
plot(llipmax, fmax,  lwd=2, xlab="Unterlippenposition (mm)", ylab="F1 (Hz)", bty="n")
plot(jmax, fmax,  lwd=2, xlab="Kieferposition (mm)", ylab="F1 (Hz)", bty="n")

cor.test(llipmax, fmax)

-----------
regm = lm(fmax ~ jmax+llipmax)