Testing species abundance models: a new bootstrap approach applied to Indo-Pacific coral reefs

Patterns in the commonness and rarity of species are a fundamental characteristic of ecological assemblages; however, testing between alternative models for such patterns remains an important challenge. Conventional approaches to fitting or testing species abundance models often assume that species, not individuals, are the units that are sampled and that species' abundances are independent of one another. Here we test three different models (the Poisson lognormal, the negative binomial, and the neutral, “zero-sum multinomial” [ZSM]) against species abundance distributions of Indo-Pacific corals and reef fishes. We derive and apply several alternative bootstrap analyses of model fit, each of which makes different assumptions about how species abundance data are sampled, and we assess the extent to which tests of model fit are sensitive to such assumptions. For all models, goodness of fit is remarkably consistent, regardless of whether one assumes that species or individuals are the units that are sampled or whether or not one assumes that species' abundances are statistically independent of one another. However, goodness-of-fit estimates are approximately twice as precise and detect lack of model fit more frequently, when based on sampling of individuals, rather than species. Bootstrap analyses indicate that the Poisson lognormal distribution exhibits substantially better fit to species abundance patterns, consistent with model selection analyses. In particular, heterogeneity in species abundances (many rare and few highly abundant species) is too great to be captured by the ZSM model or the negative binomial model and is best explained by models that predict species abundance patterns that are much closer, but not identical, to the lognormal distribution. More broadly, our bootstrap analyses suggest that estimates of model fit are likely to be robust to assumptions about the statistical interdependence of species abundances, but that tests of model fit are more powerful when they assume sampling of individuals, rather than species. Such individual-based tests therefore may be able to identify lack of model fit where previous tests have been inconclusive