if (FALSE) {
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)
}
Run the code above in your browser using DataLab