# NOT RUN {
data(federer.tobacco)
dat <- federer.tobacco
# RCB analysis. Treatment factor not signficant.
dat <- transform(dat, dosef=factor(dose), rowf=factor(row),
blockf=factor(block))
m1 <- lm(height ~ blockf + dosef, data=dat)
anova(m1)
# RCB residuals show strong spatial trends
require(lattice)
levelplot(resid(m1) ~ row * block, data=dat,
col.regions=RedGrayBlue(21), main="federer.tobacco")
# Row-column analysis. Treatment now significant
m2 <- lm(height ~ rowf + blockf + dosef, data=dat)
anova(m2)
# }
Run the code above in your browser using DataLab