if (FALSE) {
library(agridat)
data(buesocampos.tomato.uniformity)
dat <- buesocampos.tomato.uniformity
mean(dat$yield) # 5.1185 match page 43
var(dat$yield) # 2.7347 matches
require(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, tick=TRUE, aspect=(26*0.9)/(10*4),
main="buesocampos.tomato.uniformity")
}
Run the code above in your browser using DataLab