[R] ANOVA in Randomized-complete blocks design

[Richard M. Heiberger]

## 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