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