# NOT RUN {
data(strickland.grape.uniformity)
dat <- strickland.grape.uniformity
require(desplot)
desplot(yield ~ col*row, dat,
main="strickland.grape.uniformity",
flip=TRUE, aspect=(31*8)/(5*10) )
# CV 43.4<!-- % -->
sd(dat$yield, na.rm=TRUE)/mean(dat$yield, na.rm=TRUE)
# anova like Strickland, appendix 1
anova(aov(yield ~ factor(row) + factor(col), data=dat))
# numbers ending in .5 much more common than .0
# table(substring(format(na.omit(dat$yield)),4,4))
# 0 5
# 25 100
# }
Run the code above in your browser using DataLab