Social animal models for quantifying plasticity, assortment, and selection on interacting phenotypes.

Both assortment and plasticity can facilitate social evolution, as each may generate heritable associations between the phenotypes and fitness of individuals and their social partners. However, it currently remains difficult to empirically disentangle these distinct mechanisms in the wild, particularly for complex and environmentally responsive phenotypes subject to measurement error. To address this challenge, we extend the widely used animal model to facilitate unbiased estimation of plasticity, assortment and selection on social traits, for both phenotypic and quantitative genetic (QG) analysis. Our social animal models (SAMs) estimate key evolutionary parameters for the latent reaction norms underlying repeatable patterns of phenotypic interaction across social environments. As a consequence of this approach, SAMs avoid inferential biases caused by various forms of measurement error in the raw phenotypic associations between social partners. We conducted a simulation study to demonstrate the application of SAMs and investigate their performance for both phenotypic and QG analyses. With sufficient repeated measurements, we found desirably high power, low bias and low uncertainty across model parameters using modest sample and effect sizes, leading to robust predictions of selection and adaptation. Our results suggest that SAMs will readily enhance social evolutionary research on a variety of phenotypes in the wild. We provide detailed coding tutorials and worked examples for implementing SAMs in the Stan statistical programming language.