# NOT RUN {
data(cochran.lattice)
dat <- cochran.lattice
if(require(desplot)){
desplot(y~row*col|rep, dat,
text=trt, # aspect unknown, should be 2 or .5
main="cochran.lattice")
}
# }
# NOT RUN {
# Random rep,row,column model often used by Federer
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)
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace