if (FALSE) {
library(agridat)
data(summerby.multi.uniformity)
dat <- summerby.multi.uniformity
libs(desplot)
dat <- mutate(dat, env=paste(range, year, crop))
desplot(dat, yield ~ col*row|env, aspect=(5*20)/(35*20),
main="summerby.multi.uniformity")
# Show all ranges for a single year.
# dat
# Compare the variance for each dataset in Summerby, page 18, column (a)
# with what we calculate. Very slight differences.
# libs(dplyr)
# dat
## range year var summerby
## 1 R2 1922 82404 82404
## 2 R2 1923 254780. 254780
## 3 R2 1924 111978. 111978
## 4 R2 1925 84515. 84515
## 5 R2 1926 101008. 100960
## 6 R3 1922 185031. 185031
## 7 R3 1923 154777. 154784
## 8 R3 1924 252451. 252451
## 9 R3 1926 472087. 472088
## 10 R4 1924 19.3 19.341
## 11 R4 1925 14.2 14.234
## 12 R4 1926 14.2 14.236
## 13 R5 1924 134472. 134472
## 14 R5 1925 289001. 289026
## 15 R5 1926 131714. 131714
## 16 R5 1927 8.62 8.622
}
Run the code above in your browser using DataLab