if (FALSE) {
library(agridat)
data(davies.pasture.uniformity)
dat <- davies.pasture.uniformity
# range(dat$yield) # match Davies
# mean(dat$yield) # 227.77, Davies has 221.7
# sd(dat$yield)/mean(dat$yield) # 33.9, Davies has 32.5
# libs(lattice)
# qqmath( ~ yield, dat) # clearly non-normal, skewed right
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, aspect=(40*5)/(19*10), # true aspect
main="davies.pasture.uniformity")
}
Run the code above in your browser using DataLab