if (FALSE) {
library(agridat)
data(mcclelland.corn.uniformity)
dat <- mcclelland.corn.uniformity
# McClelland table 3, first row, gives 11.2
# Probable error = 0.67449 * sd(). Relative to mean.
# 0.67449 * sd(dat$yield)/mean(dat$yield) # 11.2
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE,
aspect=(219*44/12)/132, # true aspect, 219 rows * 44 inches x 132 feet
main="mcclelland.corn.uniformity")
}
Run the code above in your browser using DataLab