if (FALSE) {
library(agridat)
data(obsi.potato.uniformity)
dat <- obsi.potato.uniformity
# libs(dplyr)
# dat <- group_by(dat, loc)
# summarize(dat, yield=mean(yield))
## loc yield
##
## 1 L1 2.54 # Obsi says 2.55
## 2 L2 5.31 # Obsi says 5.36
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