dat <- goulden.barley.uniformity
desplot(yield ~ col*row, data=dat, main="goulden.barley.uniformity")
# See LeClerg, Leonard, Clark
hist(dat$yield, 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