# NOT RUN {
data(goulden.barley.uniformity)
dat <- goulden.barley.uniformity
desplot(yield ~ col*row, data=dat, main="goulden.barley.uniformity")
# See LeClerg, Leonard, Clark
hist(dat$yield, main="goulden.barley.uniformity",
breaks=c(21,40,59,78,97,116,135,154,173,192,211,230,249)+.5)
# }
Run the code above in your browser using DataLab