if (FALSE) {
library(agridat)
data(barrantesaguilar.sugarcane.uniformity)
dat <- barrantesaguilar.sugarcane.uniformity
mean(dat$yield) # 19.96 match 19.9 from Barrantes-Aguilar
libs(desplot)
# Match figure 13
desplot(dat, yield ~ col*row,
flip=TRUE, tick=TRUE, aspect=(40*2)/(40*1.5),
main="barrantesaguilar.sugarcane.uniformity")
}
Run the code above in your browser using DataLab