if (FALSE) {
library(agridat)
data(christidis.wheat.uniformity)
dat <- christidis.wheat.uniformity
# sum(dat$yield) # Matches Christidis
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, aspect=16/90, # true aspect
main="christidis.wheat.uniformity")
}
Run the code above in your browser using DataLab