library(agridat)
data(lechadolopez.sorghum.uniformity)
dat <- lechadolopez.sorghum.uniformity
mean(dat$yield) # 5382 lechado p 14, 4384 here (because we rounded to 0 dec).
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, tick=TRUE, aspect=(32*.75)/(16*3), # true aspect
main="lechadolopez.sorghum.uniformity")
Run the code above in your browser using DataLab