# NOT RUN {
library(agridat)
data(obsi.potato.uniformity)
dat <- obsi.potato.uniformity
# libs(dplyr)
# dat <- group_by(dat, loc)
# summarize(dat, yield=mean(yield))
## loc yield
## <fct> <dbl>
## 1 L1 2.54 # Obsi says 2.55
## 2 L2 5.31 # Obsi says 5.36
# }
# NOT RUN {
libs(desplot)
desplot(dat, yield ~ col*row, subset=loc=="L1",
main="obsi.potato.uniformity - loc L1",
flip=TRUE, tick=TRUE)
desplot(dat, yield ~ col*row, subset=loc=="L2",
main="obsi.potato.uniformity - loc L2",
flip=TRUE, tick=TRUE)
# }
Run the code above in your browser using DataLab