if (FALSE) {
library(agridat)
library(agridat)
data(goulden.latin)
dat <- goulden.latin
libs(desplot)
desplot(dat, yield ~ col*row,
text=trt, cex=1, # aspect unknown
main="goulden.latin")
# Matches Goulden.
m1 <- lm(yield~ trt + factor(row) + factor(col), data=dat)
anova(m1)
}
Run the code above in your browser using DataLab