Flower and host specialization and interacting species’ distribution range
In general, the contribution of the interacting species to the modelled species increased as the geographic distribution of the interacting species decreased and the flower specializationincreased . Both the GLMs of the flower visiting bees and the cleptoparasitic bees resulted in a set of candidate models with the body size, distribution of the interacting species and diversity of the interaction network as explanatory variables with the coefficient (β) and the contribution of the interacting species as the response variable (Appendix S7). The inverse link function was used for the GLM and the model for the flower visiting bees with the lowest AICc (AICc = 960.8) resulted in the formula with flower specialization (β1) and distribution of the interacting species (β2): y = (0.063 – 0.0099 * β1 + 0.013 * β2)-1 (R2 = 0.39; fig. 6A & 6B). The selected model of the cleptoparasitic bees contained the distribution of the host bee (β) and it had the formula: y = (0.13 + 0.0099*β)-1 (R2: 0.47) (fig. 6C; Appendix S7).