track.aov <- aov(mice ~ coon + block*veget*fruit*time - block:veget:fruit:time + Error(block/plot), data = track) anova(track.aov)
I think this is what you are looking for. This model in words says, What is the effect of the four-way crossing after adjusting for the covariate coon? We don't have enough degrees of freedom for the full crossing, so remove the four-way interaction. anova() gives the sequential sums of squares of (multi-degree of freedom) effects. Be sure that the five variables (block, veget, fruit, time, plot) are factors. I am assuming time is a factor here. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.