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).