Analysis of data
Analyses were carried out using the statistical software “R” (version 3.2.3. R Development Core Team 2016), with methods following Thomas et al. (2017). The MASS package (Venables and Ripley 2002) and the mgcv package (Wood 2011) were used to implement for each species a negative binomial Generalised Additive Model (GAM) analysis of the number of echo location calls. Model selection was tailored to identify the detail of seasonal patterns (week-to-week variation), by choosing K-values (degree of non-linearity) that was higher than needed to minimise AIC. These models examined non-linear temporal (seasonal and/or overnight) variation in calling activity, and tested for the linear effects of both temperature and rainfall on activity. We used the contrast package (Kuhn 2016) to examine pairwise differences between specific weeks in the number of calls recorded. Capture data were analysed using a Chi-squared approach to test for sex-ratio differences for each Myotisspecies within and between years, and using negative binomial GAM models to examine seasonal and between-year variation in capture rate. Outputs of all GAM models are presented in the Appendix. Data files and R script files for running the analyses are archived in the Supplementary Materials.