if (FALSE) {
library(agridat)
data(igue.sugarcane.uniformity)
dat <- igue.sugarcane.uniformity
# match Igue CV top row of page 171
sd(dat$yield)/mean(dat$yield) # 16.4
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, tick=TRUE, aspect=(42*2)/(36*1.5),
main="igue.sugarcane.uniformity")
}
Run the code above in your browser using DataLab