if (FALSE) {
library(agridat)
data(iyer.wheat.uniformity)
dat <- iyer.wheat.uniformity
libs(desplot)
desplot(dat, yield ~ col*row,
main="iyer.wheat.uniformity", tick=TRUE,
aspect=(25*5)/(80*5)) # true aspect
# not exactly the same as Iyer table 1, p. 241
var(subset(dat, col <= 20)$yield)
var(subset(dat, col > 20 & col <= 40)$yield)
var(subset(dat, col > 40 & col <= 60)$yield)
var(subset(dat, col > 60)$yield)
# cv for 1x1 whole-field
# sd(dat$yield)/mean(dat$yield)
# 18.3
}
Run the code above in your browser using DataLab