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.