Fig. 3 in Divergence in Body Mass, Wing Loading, and Population Structure Reveals Species-Specific and Potentially Adaptive Trait Variation Across Elevations in Montane

Fig. 3. Effect predictions (with 95% CI) of mass (field and empty), thorax size (ITS), forewing area, and transformed wing loading (pw-Field, pw-Empty) against (A) latitude, (B) AMT (Worldclim BIO1 variable), and (C) elevation from univariate linear mixed effects models for B. vancouverensis (blue) and B. vosnesenskii (red). Statistical analyses were conducted on log-transformed mass, ITS, and wing loading (see SuppTable S4 [online only]), but responses were back-transformed to the original scale for plotting (plot_model function in sjPlot). Lines are labeled with significance estimate of the β parameter estimated from lmerTest (*P <0.05; ***P <0.001; unlabeled = not significant, CI encompasses zero).