Mapping phylogenetic endemism in R using georeferenced branch extents

Applications are needed to map biodiversity from large-scale species occurrence datasets whilst seamlessly integrating with existing functions in R. Phylogenetic endemism (PE) is a biodiversity measure based on range-restricted phylogenetic diversity (PD). Current implementations use area of occupancy (AOO) or frequency to estimate the spatial range of branch-length (i.e.  phylogenetic range-rarity), rather than extent of occurrence (EOO; i.e.  georeferenced phylogenetic endemism), which is known to produce different range estimates. We present R functions to map PD or PE weighted by AOO or EOO (new georeferenced implementation), taking as inputs georeferenced species occurrences and a phylogeny. Non-parametric statistics distinguish PD/PE from trivial correlates of species richness and sampling intensity. Keywords: Phylogenetic diversity, Phylogenetic endemism, Area of occupancy, Extent of occurrence, Minimum convex polygon, R function, Non-parametric statistics, Interquartile range