if (FALSE) {
library(agridat)
data(lord.rice.uniformity)
dat <- lord.rice.uniformity
# match table on page 180
## libs(dplyr)
## dat
## field grain straw
##
## 1 10 590 732
## 2 11 502 600
## 3 12 315 488
## 4 13 291 538
## 5 14 489 670
## 6 26 441 560
## 7 27 451 629
## 8 28 530 718
# There are consistently high yields along all edges of the field
# libs(lattice)
# bwplot(grain ~ factor(col)|field,dat)
# bwplot(grain ~ factor(col)|field,dat)
# Heatmaps
libs(desplot)
desplot(dat, grain ~ col*row|field,
flip=TRUE, aspect=140/50,
main="lord.rice.uniformity")
# bivariate scatterplots
# xyplot(grain ~ straw|field, dat)
}
Run the code above in your browser using DataLab