# 5. 
w = read.emusegs("H:/segd/sb.txt")

# 6
# (a) 
dur = mudur(dip)

# (b)
dur[c(10, 11, 20)]

# (c)
sd(dur)

# (d)
n = length(dur)
dur[(n-2):n]

# 7
dip25 = dcut(dip.fdat, .25, prop=T)

# (a)
dip25[11:20,3]

# (b)
log(dip25[,3], base=10) - log(dip25[,2], base=10)

# (c)
sum(dip[,3] - dip[,2] > 200)

# 8
# (a)
sum(dip25[dip.l=="OY",1] > 380)

# (b)
mean(dip25[dip.l=="aI",2])

# (c)
sd(dip25[dip.l %in% c("aI", "aU"),1])

# 9
# (a)
mid = dcut(vowlax.fund, .5, prop=T)
mid[vowlax.spkr=="67" & vowlax.l=="a"]

#(b)
vowlax.fund[vowlax.l=="E" & vowlax.spkr=="68",]

# (c)
on = dcut(vowlax.fdat[,1:2], 0, prop=T)
on[vowlax.l %in% c("E", "a") & vowlax.spkr=="67" & vowlax.left=="d",]


# 10.
temp = vowlax.l %in% c("I", "E") & vowlax.spkr=="67"
dat = dcut(vowlax.fdat[temp,1:2], .5, prop=T)
labs = vowlax.l[temp]
eplot(dat, labs, dopoints=T, form=T)

# 11.
which = dat[,1] < 400 & labs=="E"
sum(which)

# (a)
segmente = vowlax[temp,]
mudur(segmente[which,])

# (b)
utt(segmente[which,])

# 12. 
# (a)
fs = 2 * max(trackfreq(keng.dft.5))
fs

# (b)
N = 2 * (ncol(keng.dft.5) - 1)
N/(fs/1000)

# 13.
plot(keng.dft.5[,500:3800], keng.l, fun="mean", power=T)