# NOT RUN {
library(agridat)
data(lehmann.millet.uniformity)
dat <- lehmann.millet.uniformity
libs(desplot)
dat$year = factor(dat$year)
desplot(dat, yield ~ range*plot|year,
aspect=(22*50)/(6*200),
main="lehmann.millet.uniformity",
flip=TRUE, tick=TRUE)
desplot(dat, total ~ range*plot|year,
aspect=(22*50)/(6*200),
main="lehmann.millet.uniformity",
flip=TRUE, tick=TRUE)
# libs(dplyr)
# group_by(dat, year) <!-- %>% summarize(yield=mean(yield, na.rm=TRUE)) -->
# }
Run the code above in your browser using DataLab