if (FALSE) {
library(agridat)
data(cochran.latin)
dat <- cochran.latin
libs(desplot)
desplot(dat, diff~col*row,
text=operator, cex=1, # aspect unknown
main="cochran.latin")
dat <- transform(dat, rf=factor(row), cf=factor(col))
aov.dat <- aov(diff ~ operator + Error(rf*cf), dat)
summary(aov.dat)
model.tables(aov.dat, type="means")
}
Run the code above in your browser using DataLab