## tmp.dat is your data in a file.
wing <- read.table("tmp.dat", header=TRUE)
wing.aov <- aov( wing ~ cage + Error(cage:female), data=wing )
anova(wing.aov)
## Notice that at full precision, the R values of the mean squares
## give the number you are looking for
332.84 / 1.3017
______________________________________________
[email protected] 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.
