if (FALSE) {
library(agridat)
data(goulden.barley.uniformity)
dat <- goulden.barley.uniformity
libs(desplot)
desplot(dat, yield ~ col*row,
aspect=48/48, # true aspect
main="goulden.barley.uniformity")
# Left skewed distribution. 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,268,287)+.5)
}
Run the code above in your browser using DataLab