if (FALSE) {
library(agridat)
data(heath.cabbage.uniformity)
dat <- heath.cabbage.uniformity
# Heath Fig 3-1, p. 50
libs(desplot)
desplot(dat, yield ~ col*row,
aspect=(8*12)/(6*40),
main="heath.cabbage.uniformity")
}
Run the code above in your browser using DataLab