if (FALSE) {
library(agridat)
data(palenciaortiz.sugarcane.uniformity)
dat <- palenciaortiz.sugarcane.uniformity
sum(dat$yield) # matches source document 8878.03
var(dat$yield) # matches source document 39.85
# Gradient shows higher yield in lower right
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, tick=TRUE, aspect=(18*2)/(26*2),
main="palenciaortiz.sugarcane.uniformity")
}
Run the code above in your browser using DataLab