if (FALSE) {
libs(agridat)
data(beckett.maize.uniformity)
dat <- beckett.maize.uniformity
# QC. germination, earnum match published values.
# stalks published value is 33091, but here 33101. (Data were hand-checked)
colSums(dat)
# Examine correlations. earwt,yield high cor
pairs(dat[ , c("germination","stalks","earnum","earwt","yield")])
libs(desplot)
desplot(dat, germination ~ col*row,
flip=TRUE, aspect=19/6,
main="beckett.maize.uniformity - stalks")
desplot(dat, yield ~ col*row,
flip=TRUE, aspect=19/6,
main="beckett.maize.uniformity - yield")
}
Run the code above in your browser using DataLab