# NOT RUN {
library(agridat)
data(bancroft.peanut.uniformity)
dat <- bancroft.peanut.uniformity
# match means Bancroft page 3
## dat <!-- %>% group_by(block) %>% summarize(mn=mean(yield)) -->
## # A tibble: 2 x 2
## block mn
## <chr> <dbl>
## 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