Introducing untb, an R package for simulating ecological drift under the unified neutral theory of biodiversity

The distribution of abundance amongst species with similar ways of life is a classical problem in ecology. The unified neutral theory of biodiversity, due to Hubbell, states that observed popula-tion dynamics may be explained on the assumption of per capita equivalence amongst in-dividuals. One can thus dispense with differences between species, and differences between abundant and rare species: all individuals behave alike in respect of their probabilities of reproducing and death. It is a striking fact that such a parsimonious theory results in a non-trivial dominance-diversity curve (that is, the simultaneous existence of both abundant and rare species) and even more striking that the theory predicts abundance curves that match observations across a wide range of ecologies. Here I introduce the untb package of R routines, for numerical simulation of ecological drift under the unified neutral theory. A range of visualization, analytical, and simulation tools are provided in the package and these are presented with examples in the paper. This vignette is based on Hankin (2007a). For reasons of performance, some of the more computationally expensive results are pre-loaded. To calculate them from scratch, change “calc_from_scratch <- TRUE ” to “calc_from_scratch <- FALSE ” in chunk time_saver