Appendix 5. R code for fitting phenomenological models in Chapter 5 of PhD thesis "Linking phenology to population dynamics and distribution change in a changing climate"

Appendix 5 is a .R file that contains the R code (R Core Team, 2017) required to fit the phenomenological models shown in Chapter 5. This code has been altered from that presented in Dennis et al. (2016) to allow spatial variation alongside year-to-year variation in rho (productivity) for both uni- and bi-voltine species. All alterations are marked in-line in the script with explanatory notes, using the notation "### Start addition" and "### End addition" to highlight where code has been added. Literature Cited Dennis, E. B., Morgan, B. J. T., Freeman, S. N., Roy, D. B. and Brereton, T. (2016). Dynamic models for longitudinal butterfly data. Journal of Agricultural, Biological, and Environmental Statistics, 21(1): 1–21. R Core Team (2017). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Available at: https://www.R-project.org/.