Try example.aov <- aov(y~year*site*treatment, data = df) example.aov
Then if significant, follow up with TukeyHSD(example.aov) Hope this helps. Jim Milks Sent from my iPhone > 1. Syntax of R for RCBD (Dinberu Million)
> _______________________________________________ > R-sig-ecology mailing list > R-sig-ecology@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
_______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology