# NOT RUN {
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
# require(lattice)
# qqmath( ~ yield, dat) # clearly non-normal, skewed right
require(desplot)
desplot(yield ~ col*row, dat,
flip=TRUE, aspect=(40*5)/(19*10), # true aspect
main="davies.pasture.uniformity")
# }
Run the code above in your browser using DataLab