“…These procedures were performed using the binaryPGLMM function of the R package “rr2” and the phyloglm function of the R package “phylolm” (Ho & Ané, 2014; Ives & Garland, 2010, 2014; Ives, Helmus, & Ves, 2011; Paradis, Claude, & Strimmer, 2004). In the binaryPGLMM function, s 2 is the scaling component of the variance in the PGLMM, where s 2 = 0 suggests no phylogenetic signal and a high s 2 value implies strong phylogenetic signal (Jamrozy et al., 2017). In the phyloglm function, alpha is the phylogenetic correlation parameter (an alpha value close to 0 suggests strong phylogenetic signal, alpha = 1 indicates a phylogenetic signal of trait evolution consistent with the expectation under Brownian motion, and an alpha value close to infinity implies low phylogenetic signal) (Blumstein et al., 2015; Gallien et al., 2016; O'Meara, Graham, Pellis, & Burghardt, 2015).…”