# NOT RUN {
data(masood.rice.uniformity)
dat <- masood.rice.uniformity
require(desplot)
desplot(yield ~ col*row, data=dat,
flip=TRUE, tick=TRUE, aspect=24/12, # true aspect
main="masood.rice.uniformity - yield heatmap")
require(agricolae)
require(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