# NOT RUN {
data(williams.barley.uniformity)
dat <- williams.barley.uniformity
require("lattice")
asp = (15*5)/(48*.75) # true aspect
desplot(yield ~ col*row, dat, aspect=asp,
main="williams.barley.uniformity")
# Smoothed contour/persp plot like Williams Fig 1b, 2b
dat$fit <- fitted(loess(yield~col*row, dat, span=.1))
contourplot(fit~col*row, data=dat, main="williams.barley.uniformity")
wireframe(fit~col*row, data=dat, zlim=c(100, 350),
main="williams.barley.uniformity")
# Williams table 1
anova(aov(yield ~ factor(row) + factor(col), dat))
# }
Run the code above in your browser using DataLab