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
levelplot(resid(m1) ~ row * block, data=dat, col.regions=RedGrayBlue(21))
# 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