if (FALSE) {
library(agridat)
data(masood.rice.uniformity)
dat <- masood.rice.uniformity
libs(desplot)
desplot(dat, yield ~ col*row,
flip=TRUE, tick=TRUE, aspect=24/12, # true aspect
main="masood.rice.uniformity - yield heatmap")
libs(agricolae)
libs(reshape2)
dmat <- acast(dat, row~col, value.var='yield')
index.smith(dmat,
main="masood.rice.uniformity",
col="red") # CVs match Table 3
}
Run the code above in your browser using DataLab