if (FALSE) {
library(agridat)
cv <- function(x) sd(x)/mean(x)
libs(desplot)
# Pole Bean
data(moore.polebean.uniformity)
cv(moore.polebean.uniformity$yield) # 8.00. Moore says 6.73.
desplot(moore.polebean.uniformity, yield~col*row,
flip=TRUE, tick=TRUE, aspect=120/60, # true aspect
main="moore.polebean.uniformity - yield")
# Bush bean
data(moore.bushbean.uniformity)
cv(moore.bushbean.uniformity$yield) # 12.1. Moore says 10.8
desplot(moore.bushbean.uniformity, yield~col*row,
flip=TRUE, tick=TRUE, aspect=120/72, # true aspect
main="moore.bushbean.uniformity - yield")
# Sweet corn
data(moore.sweetcorn.uniformity)
cv(moore.sweetcorn.uniformity$yield) # 17.5. Moore says 13.6
desplot(moore.sweetcorn.uniformity, yield~col*row,
flip=TRUE, tick=TRUE, aspect=120/72, # true aspect
main="moore.sweetcorn.uniformity - yield")
## desplot(moore.sweetcorn.uniformity, ears~col*row,
## flip=TRUE, tick=TRUE, aspect=120/72, # true aspect
## main="moore.sweetcorn.uniformity - ears")
## libs(lattice)
## xyplot(yield ~ ears, moore.sweetcorn.uniformity)
libs(desplot)
# Carrot
data(moore.carrot.uniformity)
cv(moore.carrot.uniformity$yield) # 33.4. Moore says 27.6
desplot(moore.carrot.uniformity, yield~col*row,
flip=TRUE, tick=TRUE, aspect=60/36, # true aspect
main="moore.carrot.uniformity - yield")
libs(desplot)
# Spring cauliflower
data(moore.springcauliflower.uniformity)
cv(moore.springcauliflower.uniformity$yield) # 21. Moore says 19.5
desplot(moore.springcauliflower.uniformity, yield~col*row,
flip=TRUE, tick=TRUE, aspect=300/36, # true aspect
main="moore.springcauliflower.uniformity - yield")
## desplot(moore.springcauliflower.uniformity, heads~col*row,
## flip=TRUE, tick=TRUE, aspect=300/36, # true aspect
## main="moore.springcauliflower.uniformity - heads")
## libs(lattice)
## xyplot(yield ~ heads, moore.springcauliflower.uniformity)
libs(desplot)
# Fall cauliflower
data(moore.fallcauliflower.uniformity)
cv(moore.fallcauliflower.uniformity$yield) # 17.7. Moore says 17.0
desplot(moore.fallcauliflower.uniformity, yield~col*row,
flip=TRUE, tick=TRUE, aspect=300/36, # true aspect
main="moore.fallcauliflower.uniformity - yield")
## desplot(moore.fallcauliflower.uniformity, heads~col*row,
## flip=TRUE, tick=TRUE, aspect=300/36, # true aspect
## main="moore.fallcauliflower.uniformity - heads")
## libs(lattice)
## xyplot(yield ~ heads, moore.fallcauliflower.uniformity)
}
Run the code above in your browser using DataLab