if (FALSE) {
library(agridat)
data(magistad.pineapple.uniformity)
dat <- magistad.pineapple.uniformity
# match table page 641
## dat
## summarize(number=mean(number),
## weight=mean(weight))
## field number weight
## 1 1 596.4062 2499.922
## 2 19 171.1667 2100.250
## 3 21 171.1600 2056.800
## 4 82 220.7500 1264.500
libs(desplot)
desplot(dat, weight ~ col*row,
subset=field==19,
aspect=300/39,
main="magistad.pineapple.uniformity - field 19")
desplot(dat, weight ~ col*row,
subset=field==82,
aspect=228/176,
main="magistad.pineapple.uniformity - field 82")
desplot(dat, weight ~ col*row,
subset=field==21,
aspect=300/97.5,
main="magistad.pineapple.uniformity - field 21")
desplot(dat, weight ~ col*row,
subset=field==1,
aspect=1200/96,
main="magistad.pineapple.uniformity - field 1")
}
Run the code above in your browser using DataLab