phyr: An r package for phylogenetic species-distribution modelling in ecological communities

Model-based approaches are increasingly popular in ecological studies. A good example of this trend is the use of joint species distribution models to ask questions about ecological communities. However, most current applications of model-based methods do not include phylogenies despite the well-known importance of phylogenetic relationships in shaping species distributions and community composition. In part, this is due to a lack of accessible tools allowing ecologists to fit phylogenetic species distribution models easily. To fill this gap, the r package phyr (pronounced fire) implements a suite of metrics, comparative methods and mixed models that use phylogenies to understand and predict community composition and other ecological and evolutionary phenomena. The phyr workhorse functions are implemented in C++ making all calculations and model estimations fast. phyr can fit a variety of models such as phylogenetic joint-species distribution models, spatiotemporal-phylogenetic autocorrelation models, and phylogenetic trait-based bipartite network models. phyr also estimates phylogenetically independent trait correlations with measurement error to test for adaptive syndromes and performs fast calculations of common alpha and beta phylogenetic diversity metrics. All phyr methods are united under Brownian motion or Ornstein–Uhlenbeck models of evolution, and phylogenetic terms are modelled as phylogenetic covariance matrices. The functions and model formula syntax we propose in phyr provide an easy-to-use collection of tools that we hope will ignite the use of phylogenies to address a variety of ecological questions.