powered by
data(goulden.barley)
dat <- goulden.barley desplot(yield ~ col*row, data=dat) # 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