# NOT RUN {
library(agridat)
data(li.millet.uniformity)
dat <- li.millet.uniformity
mean(dat$yield) # matches Li et al.
libs(desplot)
desplot(dat, yield~col*row,
aspect=100/100, # true aspect
main="li.millet.uniformity")
# }
Run the code above in your browser using DataLab