my apologies. anova() doesn't work on class "aovlist". It is necessary to use summary(). They are similar for "aov", but not for "aovlist". wing <- read.table("tmp.dat", header=TRUE) wing.aov <- aov( wing ~ cage + Error(cage:female), data=wing ) summary(wing.aov)