Akaike weights of all the models (Table 1) were used for calculating the multi-model coefficient estimates, the sign of which are given in the table. The sums of Akaike weights of models containing each of the variables are given (wi). Vegetation productivity (GVI, Figure 1a) of the habitats received the most support. Other abbreviations used: precipitation variation (PSea; Figure 1b), two binary regional variables (Region binary 1/2) and the interactions of the two environmental predictors with the former categorical variable for region (Region*GVI, Region*PSea; Figure 1a & b) of each of the three social species. The coefficient estimates for the main regional and interaction effects are not given, as they are difficult to interpret and only weakly supported. βA denotes the model-averaged regression coefficients; ORsoc is the odds ratio of being social; and Psoc is the probability of finding a social species with one unit increase of the respective environmental variable.