dat <- gomez.groupsplit
# Gomez figure 3.10. Obvious fert and group effects
desplot(group~x*y, dat, out1=rep, col=fert, text=gen,
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