if (FALSE) {
library(agridat)
data(lyon.potato.uniformity)
dat <- lyon.potato.uniformity
# Yield per square foot, assuming 72 foot rows
sum(dat$yield)/(72*96) # 0.67 # seems about right
# Yield per square foot, assuming 72 foot plots
sum(dat$yield)/(6*72*96) # 0.11
libs(desplot)
desplot(dat, yield ~ col*row,
tick=TRUE, flip=TRUE, aspect=96/72, # true aspect
main="lyon.potato.uniformity")
}
Run the code above in your browser using DataLab