if (FALSE) {
library(agridat)
data(holtsmark.timothy.uniformity)
dat <- holtsmark.timothy.uniformity
# Define diagonal 'check' plots like Holtsmark does
dat <- transform(dat,
check = ifelse(floor((row+col)/3)==(row+col)/3, "C", ""))
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, text=check, show.key=FALSE,
aspect=30/200, # true aspect
main="holtsmark.timothy.uniformity")
# sd(dat$yield) # 2.92 matches Holtsmark p. 348
}
Run the code above in your browser using DataLab