if (FALSE) {
library(agridat)
data(strickland.peach.uniformity)
dat <- strickland.peach.uniformity
mean(dat$yield) # 131.3, Strickland has 131.3
sd(dat$yield)/mean(dat$yield) # 31.1, Strickland has 34.4
libs(desplot)
desplot(dat, yield ~ col*row,
main="strickland.peach.uniformity",
flip=TRUE, aspect=1)
}
Run the code above in your browser using DataLab