if (FALSE) {
library(agridat)
data(kiesselbach.oats.uniformity)
dat <- kiesselbach.oats.uniformity
range(dat$yield) # 56.7 92.8 match Kiesselbach p 64.
libs(desplot)
desplot(dat, yield ~ col*row,
tick=TRUE, flip=TRUE, aspect=792/469, # true aspect
main="kiesselbach.oats.uniformity")
}
Run the code above in your browser using DataLab