dat <- cochran.latin
desplot(diff~col*row, dat, text=operator, cex=1)
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