By fitting a glm model, I obtained the following result:
anova(m1.glm, test="F")
Analysis of Deviance Table
Model: Negative Binomial(4.0016), link: log
Response: KDE95
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev F Pr(>F)
NULL 54 63.952
Season 3 1.2058 51 62.746 0.4661 0.70735
Year 1 5.2026 50 57.544 6.0327 0.01779 *
Season:Year 3 10.5437 47 47.000 4.0753 0.01184 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
I would like to plot the interaction among the factor "Season" and the discrete variable "Year" on the dependent variable KDE 95. In essence, I would like to plot the difference in term of KDE 95 among seasons within each year. Is this possible?
Many thanks!