if (FALSE) {
library(agridat)
data(khan.brassica.uniformity)
dat <- khan.brassica.uniformity
# Slightly different results than Khan Table 1.
## dat
## mutate(yield=yield/8)
## group_by(field)
## summarize(mn=mean(yield), sd=sd(yield))
libs(desplot)
desplot(dat, yield ~ col*row | field,
flip=TRUE, aspect=1,
main="khan.brassica.uniformity")
}
Run the code above in your browser using DataLab