dat <- gomez.stripplot
# Gomez figure 3.7
desplot(gen~x*y, data=dat, out1=rep, num=nitro, cex=1,
main="gomez.stripplot")
# Gomez table 3.12
tapply(dat$yield, dat$rep, sum)
tapply(dat$yield, dat$gen, sum)
tapply(dat$yield, dat$nitro, sum)
# Gomez table 3.15. Anova table for strip-plot
dat <- transform(dat, nf=factor(nitro))
m1 <- aov(yield ~ gen * nf + Error(rep + rep:gen + rep:nf), data=dat)
summary(m1)Run the code above in your browser using DataLab