## Abstract

Statistical models are helping palaeontologists to elucidate the history of biodiversity. Sampling standardization has been extensively applied to remedy the effects of uneven sampling in large datasets of fossil invertebrates. However, many vertebrate datasets are smaller, and the issue of uneven sampling has commonly been ignored, or approached using pairwise comparisons with a numerical proxy for sampling effort. Although most authors find a strong correlation between palaeodiversity and sampling proxies, weak correlation is recorded in some datasets. This has led several authors to conclude that uneven sampling does not influence our view of vertebrate macroevolution. We demonstrate that multi-variate regression models incorporating a model of underlying biological diversification, as well as a sampling proxy, fit observed sauropodomorph dinosaur palaeodiversity best. This bivariate model is a better fit than separate univariate models, and illustrates that observed palaeodiversity is a composite pattern, representing a biological signal overprinted by variation in sampling effort. Multi-variate models and other approaches that consider sampling as an essential component of palaeodiversity are central to gaining a more complete understanding of deep time vertebrate diversification.

## 1. Introduction

Understanding how biodiversity has fluctuated over extended intervals of deep time, and how it responds to extinction events, is central to understanding the significance of modern biodiversity loss and responses to climate change. Fossils provide the only data on these issues, and attempts to understand the diversification of life using fossil data have a long pedigree (e.g. [1–3]). Initial taxonomic compilations (e.g. [4–6]) formed the basis of ‘palaeodiversity curves’, which were commonly interpreted as a literal reading of ancient biodiversity [2,3]. However, it is possible that much of the variation in observed ‘global’ palaeodiversity is driven by temporal and spatial variation in the amount of available fossiliferous rock, or by disparities in collection effort by palaeontologists (e.g. [7,8]). To remedy this, recent palaeodiversity databases, exemplified by *The Palaeobiology Database* (PBDB; http://www.paleodb.org/), incorporate subsampling mechanisms that allow for the simulation of even fossil sampling through time. This has resulted in substantial revisions to palaeodiversity curves, and to our understanding of the history of life on Earth [9,10].

Large-scale patterns across the entire Phanerozoic (542 million years ago (Mya)–present) have been the focus of palaeodiversity research (e.g. [2,3,9,10]). The Phanerozoic record primarily comprises shallow marine invertebrates inhabiting the submerged continental shelves, where most rock deposition occurs. Terrestrial and open ocean palaeodiversity have received comparatively little attention, and studies of vertebrate palaeodiversity were relatively rare until recently. Vertebrate fossils are less abundant. For some taxonomic groups, facies and time intervals, many collecting localities have yielded only a single specimen or taxon, instead of a more complete faunal sample (PBDB, accessed on 27 May 2011). These factors restrict the utility of subsampling approaches. Thus, many vertebrate palaeodiversity studies have relied on modelling approaches to ‘correct’ data for uneven sampling (e.g. [11–15]). However, perspectives on the importance of sampling to our understanding of vertebrate macroevolution are polarized. Some workers have stated that macroevolutionary patterns contain a strong biological signal, and can be interpreted at face value [16–18], whereas others suggest that these signals are distorted by sampling biases, and emphasize caution (e.g. [12–14,19,20]).

Because of this lack of consensus, some high-profile macroevolutionary studies of vertebrates recently have ignored sampling biases altogether [18,21]. This runs counter to scientific intuition, and the fact that several studies detected a significant correlation between observed vertebrate palaeodiversity and proxies representing sampling effort (e.g. [12–14,20]). Several authors have argued that the apparent absence of this correlation for some datasets indicates that vertebrate palaeodiversity signals are not influenced by sampling bias [17,18,22]. However, this may not be correct: only if biodiversity was constant through time would we expect to see a perfect correlation between sampling and observed palaeodiversity. In cases where underlying biodiversity exhibits high levels of variation, this correlation should become weaker. However, this would not be a special case of an ‘unbiased’ fossil record, because the strong, ‘genuine’ biodiversity signal should still be overprinted by variation in sampling effort.

We illustrate this principle using multi-variate regression models for sauropodomorph dinosaur palaeodiversity. These show that including even a simple numerical representation of biological diversification in a multi-variate model can improve the fit of a sampling proxy beyond that obtained by a univariate, sampling-only, model. Furthermore, some representations of biodiversity, though justified by palaeontological observations, may show poor fit to observed palaeodiversity unless sampling is explicitly considered by following a multi-variate approach.

## 2. Material and methods

Sauropodomorpha is a dinosaurian clade of primarily long-necked, large-bodied, herbivorous taxa. It includes more than 200 known taxa (e.g. *Diplodocus*, *Brachiosaurus*), first appearing in the Late Triassic (228 Mya) and surviving until the end of the Cretaceous (65.5 Mya) [14,23]. Sauropodomorphs have been the focus of several palaeodiversity studies [14,24]. Compared with other dinosaurian clades, observed sauropodomorph palaeodiversity shows only a weak correlation with sampling proxies [12]. The significance of this is disputed; it may indicate that sampling does not influence sauropodomorph palaeodiversity (e.g. [18], see the electronic supplementary material) or that sauropodomorph biodiversity showed genuinely greater fluctuations, obscuring the relationship between sampling and palaeodiversity in univariate comparisons [14]. Notably, sauropodomorphs suffered a substantial extinction event at the end of the Jurassic (approx. 145 Mya), resulting in the disappearance of broad-toothed forms and their replacement by ornithischians in Cretaceous ecosystems [14,24,25]. If this is correct then sauropodomorph biodiversity may be approximated by a simple model (‘TJK’ model) in which Triassic–Jurassic (TJ) biodiversity is assigned one value (‘0’) and Cretaceous (K) biodiversity is assigned a different value (‘1’).

To test between interpretations of sauropodomorph palaeodiversity, we analysed a comprehensive species-level dataset of Late Triassic–Cretaceous sauropodomorphs [14,23] parsed into 26 stage-level time bins (figure 1*a*; electronic supplementary material, appendix S1). We compared this with our simple model of sauropodomorph biodiversity, and to two sampling proxies (counts of geological formations (dinosaur-bearing formations, DBFs) and collections (dinosaur-bearing collections, DBCs) yielding dinosaur body fossils; figure 1*a*) [14,23], using generalized least-squares regression. This approach incorporates an autoregressive model to account for potential non-independence of successive data points in a time series. Variables were log-e transformed prior to analysis. Residuals from the regression models were normally distributed and homoskedastic. Three types of regression model were compared: univariate models comprising (i) a sampling proxy, (ii) our simple biodiversity model, and (iii) a bivariate model comprising a sampling proxy and our simple biodiversity model. All analyses were implemented in R v. 2.10.1 [26], following the approach of Hunt *et al*. [27] and Marx & Uhen [28] (electronic supplementary material, appendix S2).

## 3. Results

All models including a sampling proxy fit significantly better than the null model in which palaeodiversity is constant and variation is subsumed by an error term. However, substantially higher *R*^{2} (proportion of variance explained) and lower Akaike information criterion (AICc) scores [29,30] indicate that the bivariate model is the best explanation of sauropodomorph palaeodiversity (table 1). The univariate model including only our simple biodiversity model is a poorer fit than the null model (figure 1*b* and table 1), indicating that it is significantly different from the impression of sauropodomorph diversity obtained by inspecting ‘raw’ palaeodiversity data. Importantly, both sampling proxies show a greater *t*-value and stronger statistical significance within the bivariate models than they do in the univariate models (table 1). Thus, including an estimate of sauropodomorph biodiversity in the regression model actually improves the fit of the sampling proxy. This is not consistent with the suggestion that sauropodomorph palaeodiversity is independent of sampling bias [18]. These results are independent of the effects of the variable duration of geological stages (electronic supplementary material, appendix S2).

## 4. Discussion

Our results illustrate a general principle relevant to all palaeodiversity studies. Poor fit of univariate sampling models (e.g. pairwise statistical tests of correlation) may not always result from an ‘unbiased’ palaeodiversity curve. Instead, all palaeodiversity curves represent a composite signal, comprising both sampling *and* genuine biological diversity. Our sauropodomorph biodiversity model is coarse, postulating two phases of standing biodiversity, a Late Triassic–Jurassic interval of relatively high diversity followed by a Cretaceous interval of lower diversity. Nonetheless, when considered alongside uneven sampling, this model fits the observed palaeodiversity data extremely well. In principle, more complex biodiversity models could be applied. For example, logistic dynamics, which incorporate an initial phase of increasing diversity followed by the attainment of ‘carrying capacity’, describe sampling-standardized palaeodiversity curves for most invertebrate clades [10,31].

Multi-variate models (e.g. [32,33]), and other approaches, which isolate biological and sampling signals (e.g. subsampling [9,10,14,19,31] and the ‘residuals’ method of Smith & McGowan [11–15]), are key to understanding palaeodiversity. This knowledge is central to correctly interpreting macroevolutionary patterns, and we strongly urge that all palaeodiversity studies employ appropriate methods to account for sampling biases to elucidate the biological significance of the fossil record.

## Footnotes

One contribution to a Special Feature on ‘Methods in Palaeontology’.

- Received April 29, 2011.
- Accepted May 27, 2011.

- This Journal is © 2011 The Royal Society