if (FALSE) {
library(agridat)
data(beaven.barley)
dat <- beaven.barley
# Match the means shown in Richey table IV
tapply(dat$yield, dat$gen, mean)
## a b c d e f g h
## 298.080 300.710 318.685 295.260 306.410 276.475 304.605 271.820
# Compare to Student 1923, diagram I,II
libs(desplot)
desplot(dat, yield ~ col*row,
aspect=15/96, # true aspect
main="beaven.barley - variety trial", text=gen)
}
Run the code above in your browser using DataLab