if (FALSE) {
library(agridat)
data(eden.tea.uniformity)
dat <- eden.tea.uniformity
# sum(dat$yield) # 140050.6 matches total yield in appendix A
# mean(dat$yield) # 972.574 match page 5554
m1 <- aov(yield ~ factor(entry) + factor(row) + factor(col), data=dat)
summary(m1)
libs(desplot)
desplot(dat, yield ~ col*row,
aspect=1,
main="eden.tea.uniformity")
}
Run the code above in your browser using DataLab