if (FALSE) {
library(agridat)
data(lehmann.rice.uniformity)
dat <- lehmann.rice.uniformity
libs(desplot)
dat$year = factor(dat$year)
desplot(dat, yield ~ range*plot|year,
aspect=(17*50)/(2*200),
main="lehmann.rice.uniformity",
flip=TRUE, tick=TRUE)
desplot(dat, total ~ range*plot|year,
aspect=(17*50)/(2*200),
main="lehmann.rice.uniformity",
flip=TRUE, tick=TRUE)
# libs(dplyr)
# group_by(dat, year)
}
Run the code above in your browser using DataLab