if (FALSE) {
library(agridat)
data("sawyer.multi.uniformity")
dat <- sawyer.multi.uniformity
libs(desplot)
# The field plan shows square plots
desplot(dat, grain~col*row,
subset= year==1925,
main="sawyer.multi.uniformity - 1925 wheat grain yield",
aspect=(6)/(8)) # true aspect
desplot(dat, rootwt~col*row,
subset= year==1926,
main="sawyer.multi.uniformity - 1926 root weight of swedes",
aspect=(6)/(8))
desplot(dat, grain~col*row, subset= year==1927,
main="sawyer.multi.uniformity - 1927 oats grain yield",
aspect=(6)/(8))
# This plot shows the "outlier" in the wheat data reported by Mackenzie.
libs(lattice)
xyplot(grain ~ straw, data=subset(dat, year==1925))
round(cor(dat[,7:9], use="pair"),2) # Matches McCullagh p 2121
## leafwt rootwt rootct
## leafwt 1.00 0.66 0.47
## rootwt 0.66 1.00 0.43
## rootct 0.47 0.43 1.00
## pairs(dat[,7:9],
## main="sawyer.multi.uniformity")
}
Run the code above in your browser using DataLab