if (FALSE) {
library(agridat)
data(robinson.peanut.uniformity)
dat <- robinson.peanut.uniformity
# Mean yield per year. Robinson has 703.9, 787.3
# tapply(dat$yield, dat$year, mean)
# 1939 1940
# 703.7847 787.8125
libs(desplot)
desplot(dat, yield ~ col*row|year,
flip=TRUE, tick=TRUE, aspect=200/48,
main="robinson.peanut.uniformity")
}
Run the code above in your browser using DataLab