code = c("d","s","w")
d = zmean[,c("cvoice","subj","voiced")]
d.t = Anova.prepare(d,code)
d.lm = lm(d.t$d~1)
Anova(d.lm,idata=d.t$w,idesign=~voiced)

code = c("d","s","w","w")
d = zmean.c[,c("cvoice","subj","voiced","c1")]
d.t = Anova.prepare(d,code)
d.lm = lm(d.t$d~1)
Anova(d.lm,idata=d.t$w,idesign=~voiced*c1)

accmean = anova.mean(accdur$dur,accdur$subj,accdur$focus,accdur$age
names(accmean) = names(c("subj","focus","dur","age"))
code = c("d","s","w","b")
d.t = Anova.prepare(accmean,code)
d.lm = lm(d.t$d~d.t$b)
Anova(d.lm,idata=d.t$w,idesign=~focus)

focmean = anova.mean(all$f0max2.hz,all$subj,all$focus,all$sent)


focmean = anova.mean((all$acc1.end-all$acc1.start),all$subj,all$focus,all$sent)
names(focmean) = c("dur","subj","focus","sent")
code = c("d","s","w","w")
d = focmean[,c("dur","subj","focus","sent")]
d.t = Anova.prepare(focmean,code)
d.lm = lm(d.t$d~1)
Anova(d.lm,idata=d.t$w,idesign=~focus*sent)


## no broad focus contained ##
nob = all$focus != "broad"
datfr = all[nob,]
levels(datfr$focus)[1] = NA
focmean = anova.mean((datfr$acc1.end-datfr$acc1.start),datfr$subj,datfr$focus,datfr$sent)
code = c("d","s","w","w")
names(focmean) = c("dur","subj","focus","sent")
d = focmean[,c("dur","subj","focus","sent")]
d.t = Anova.prepare(focmean,code)
d.lm = lm(d.t$d~1)
Anova(d.lm,idata=d.t$w,idesign=~focus*sent)