if (FALSE) {
library(agridat)
data(christidis.cotton.uniformity)
dat <- christidis.cotton.uniformity
# Match the mean yields in table 2. Not sure why '16' is needed
# sapply(split(dat$yield, dat$block), mean)*16
libs(desplot)
dat$yld <- dat$yield/4*1000 # re-scale to match Christidis fig 1
desplot(dat, yld ~ col*row|block,
flip=TRUE, aspect=(16)/(64),
main="christidis.cotton.uniformity")
}
Run the code above in your browser using DataLab