if (FALSE) {
library(agridat)
data(andres.rice.uniformity)
dat <- andres.rice.uniformity
var(dat$yield) # 5611.803 match andres
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, aspect=(30*1)/(33*0.9),
main="andres.rice.uniformity")
}
Run the code above in your browser using DataLab