if (FALSE) {
library(agridat)
data(bancroft.peanut.uniformity)
dat <- bancroft.peanut.uniformity
# match means Bancroft page 3
## dat
## # A tibble: 2 x 2
## block mn
##
## 1 B1 2.46
## 2 B2 2.05
libs(desplot)
desplot(dat, yield ~ col*row|block,
flip=TRUE, aspect=(18*3)/(6*16.66), # true aspect
main="bancroft.peanut.uniformity")
}
Run the code above in your browser using DataLab