if (FALSE) {
library(agridat)
data(baker.strawberry.uniformity)
dat <- baker.strawberry.uniformity
# Match mean and cv of Baker p 414.
libs(dplyr)
dat <- group_by(dat, trial)
summarize(dat, mn=mean(yield), cv=sd(yield)/mean(yield))
dat <- mutate(dat, col = ifelse(trial=="T1" & col==4, 7, col),
col = ifelse(trial=="T1" & col==3, 6, col))
libs(desplot)
desplot(dat, yield ~ col*row, subset=trial=="T1",
flip=TRUE, aspect=500/(5*12),
main="baker.strawberry.uniformity - trial T1")
desplot(dat, yield ~ col*row, subset=trial=="T2",
flip=TRUE, aspect=500/300,
main="baker.strawberry.uniformity - trial T2")
}
Run the code above in your browser using DataLab