if (FALSE) {
library(agridat)
data(vishnaadevi.rice.uniformity)
dat <-vishnaadevi.rice.uniformity
# CV in Table 2 for 1x1 is reported as 2.8
# sd(dat$yield)/mean(dat$yield) = .0277
# If we remove the outlier yield 685, then we calculate .0256
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, aspect=24/12,
main="vishnaadevi.rice.uniformity")
}
Run the code above in your browser using DataLab