Analysis of deviance table of hierarchical generalized mixed-effect models applied to the data.

Maximal model (with all possible main effects and interactions of fixed factors) was simplified and less complex models were created. First, a three-way interaction between treatment, time and year was removed, followed by two-way interactions and then main effects of treatment, time and year. The table shows the growth of residual deviance associated with omitting each term or interaction, which was tested using Chi-square tests.