if (FALSE) {
library(agridat)
data(strickland.grape.uniformity)
dat <- strickland.grape.uniformity
libs(desplot)
desplot(dat, yield ~ col*row,
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