Estimating ancient biogeographic patterns with statistical model discrimination

The geographic ranges in which species live is a function of many factors underlying ecological and evolutionary contingencies. Observing the geographic range of an individual species provides valuable information about these historical contingencies for a lineage, determining the distribution of many distantly related species in tandem provides information about large-scale constraints on evolutionary and ecological processes generally. We present a linear regression method that allows for the discrimination of various hypothetical biogeographical models for determining which landscape distributional pattern best matches data from the fossil record. The linear regression models used in the discrimination rely on geodesic distances between sampling sites (typically geologic formations) as the independent variable and three possible dependent variables: Dice/Sorensen similarity; Euclidean distance; and phylogenetic community dissimilarity. Both the similarity and distance measures are useful for full-community analyses without evolutionary information, whereas the phylogenetic community dissimilarity requires phylogenetic data. Importantly, the discrimination method uses linear regression residual error to provide relative measures of support for each biogeographical model tested, not absolute answers or p-values. When applied to a recently published dataset of Campanian pollen, we find evidence that supports two plant communities separated by a transitional zone of unknown size. A similar case study of ceratopsid dinosaurs using phylogenetic community dissimilarity provided no evidence of a biogeographical pattern, but this case study suffers from a lack of data to accurately discriminate and/or too much temporal mixing. Future research aiming to reconstruct the distribution of organisms across a landscape has a statistical-based method for determining what biogeographic distributional model best matches the available data.