if (FALSE) {
library(agridat)
data(devries.pine)
dat <- devries.pine
libs(desplot)
desplot(dat, volume ~ col*row,
main="devries.pine - expt design and tree volume",
col=spacing, num=thinning, cex=1, out1=block, aspect=200/720)
libs(HH)
HH::interaction2wt(volume ~ spacing+thinning, dat,
main="devries.pine")
# ANOVA matches appendix 5 of DeVries
m1 <- aov(volume ~ block + spacing + thinning + block:factor(row) +
block:factor(col), data=dat)
anova(m1)
}
Run the code above in your browser using DataLab