# NOT RUN {
data(cochran.lattice)
dat <- cochran.lattice
desplot(y~row*col|rep, dat, text=trt, cex=2, main="cochran.lattice")
# Random rep,row,column model often used by Federer
if(require(lme4)){
dat <- transform(dat, rowf=factor(row), colf=factor(col))
m1 <- lmer(y ~ trt + (1|rep) + (1|rep:row) + (1|rep:col), data=dat)
summary(m1)
}
# }
Run the code above in your browser using DataLab