if (FALSE) {
library(agridat)
data(panse.cotton.uniformity)
dat <- panse.cotton.uniformity
# match the CV of Panse 1954
# sd(dat$yield)/mean(dat$yield) * 100
# 32.1
# match the fertility map of Hutchinson, fig 1
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, aspect=207/150, # true aspect
main="panse.cotton.uniformity")
}
Run the code above in your browser using DataLab