if (FALSE) {
library(agridat)
data(hernandezdavila.potato.uniformity)
dat <- hernandezdavila.potato.uniformity
# Match the mean and sd in the source
libs(dplyr)
dat |> group_by(expt) |> summarize(mn=mean(yield), sd=sd(yield))
libs(desplot)
desplot(dat, yield ~ col*row, subset=expt=="E1",
tick=TRUE, flip=TRUE, aspect=(42*0.9)/(12*0.9),
main="hernandezdavila.potato.uniformity - Expt E1")
desplot(dat, yield ~ col*row, subset=expt=="E2",
tick=TRUE, flip=TRUE, aspect=(42*0.9)/(12*0.9),
main="hernandezdavila.potato.uniformity - Expt E2")
}
Run the code above in your browser using DataLab