if (FALSE) {
library(agridat)
data("sawyer.multi.uniformity")
dat <- sawyer.multi.uniformity
libs(desplot)
# The field plan shows square plots
desplot(dat, yield~col*row,
subset= (year==1925 & crop=="wheat_grain"),
main="sawyer.multi.uniformity - 1925 wheat grain yield",
aspect=(6)/(8)) # true aspect
desplot(dat, yield~col*row,
subset= (year==1926 & crop=="swedes_rootwt"),
main="sawyer.multi.uniformity - 1926 root weight of swedes",
aspect=(6)/(8))
desplot(dat, yield~col*row,
subset= (year==1927 & crop=="oats_grain"),
main="sawyer.multi.uniformity - 1927 oats grain yield",
aspect=(6)/(8))
# Match Mackenzie, p. 279. This plot shows a possible "outlier".
d1 <- subset(dat, year==1925 & crop=="wheat_grain")
d2 <- subset(dat, year==1925 & crop=="wheat_straw")
d12 <- merge(d1[ , c("row","col","yield")],
d2[ , c("row","col","yield")], by=c("row","col"))
plot(d12$yield.y, d12$yield.x, xlab="1925 straw", ylab="1925 grain")
# Swede correlations, 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
}
Run the code above in your browser using DataLab