# NOT RUN {
data(kulkarni.sorghum.uniformity)
dat <- kulkarni.sorghum.uniformity
# match means on page 462
# tapply(dat$yield, dat$year, mean)
# 1930 1931 1932
# 116.2875 67.2250 126.3688
require(reshape2)
require(lattice)
dmat <- acast(dat, row+col ~ year, value.var="yield")
splom(dmat, main="kulkarni.sorghum.uniformity")
cor(dmat)
require(desplot)
desplot(yield ~ col*row|year, dat,
flip=TRUE, aspect=150/290,
main="kulkarni.sorghum.uniformity")
# }
Run the code above in your browser using DataCamp Workspace