if (FALSE) {
library(agridat)
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
libs(reshape2)
libs(lattice)
dmat <- acast(dat, row+col ~ year, value.var="yield")
splom(dmat, main="kulkarni.sorghum.uniformity")
cor(dmat)
libs(desplot)
desplot(dat, yield ~ col*row|year,
flip=TRUE, aspect=150/290,
main="kulkarni.sorghum.uniformity")
}
Run the code above in your browser using DataLab