library(agridat)
data(gomez.groupsplit)
dat <- gomez.groupsplit
# Gomez figure 3.10. Obvious fert and group effects
libs(desplot)
desplot(dat, group ~ col*row,
out1=rep, col=fert, text=gen, # aspect unknown
main="gomez.groupsplit")
# Gomez table 3.19 (not partitioned by group)
m1 <- aov(yield ~ fert*group + gen:group + fert:gen:group +
Error(rep/fert/group), data=dat)
summary(m1)
Run the code above in your browser using DataLab