if (FALSE) {
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