Difficulties in summing log-normal distributions for abundance and potential solutions

Abstract

The log-normal distribution, often used to model animal abundance and its uncertainty, is central to ecological modeling and conservation but its statistical properties are less intuitive than those of the normal distribution. The right skew of the log-normal distribution can be considerable for highly uncertain estimates and the median is often chosen as a point estimate. However, the use of the median can become complicated when summing across populations since the median of the sum of log-normal distributions is not the sum of the constituent medians. Such estimates become sensitive to the spatial or taxonomic scale over which abundance is being summarized and the naive estimate (the median of the distribution representing the sum across populations) can become grossly inflated. Here we review the statistical issues involved and some alternative formulations that might be considered by ecologists interested in modeling abundance. Using a recent estimate of global avian abundance as a case study (Callaghan et al. 2021), we investigate the properties of several alternative methods of summing across species' abundance, including the sorted summing used in the original study (Callaghan et al. 2021) and the use of shifted log-normal distributions, truncated normal distributions, and rectified normal distributions. The appropriate method of summing across distributions was intimately tied to the use of the mean or median as the measure of central tendency used as the point estimate. Use of the shifted log-normal distribution, however, generated scale-consistent estimates for global abundance across a spectrum of contexts. Our paper highlights how seemingly inconsequential decisions regarding the estimation of abundance yield radically different estimates of global abundance and its uncertainty, with conservation consequences that are underappreciated and require careful consideration.