# NOT RUN {
data(kadam.millet.uniformity)
dat <- kadam.millet.uniformity
# similar to Kadam fig 1
require(desplot)
desplot(yield ~ col*row, dat,
subset=year==1933,
flip=TRUE, aspect=(10*33)/(8*16.5), # true aspect
main="kadam.millet.uniformity 1933")
desplot(yield ~ col*row, dat,
subset=year==1934,
flip=TRUE, aspect=(20*16.5)/(8*16.5), # true aspect
main="kadam.millet.uniformity 1934")
# }
Run the code above in your browser using DataLab