Statistical methods in spatial genetics

The joint analysis of spatial and genetic data is rapidly becoming the norm in population genetics. More and more studies explicitly describe and quantify the spatial organization of genetic variation and try to relate it to underlying ecological processes. As it has become increasingly difficult to keep abreast with the latest methodological developments, we review the statistical toolbox available to analyse population genetic data in a spatially explicit framework. We mostly focus on statistical concepts but also discuss practical aspects of the analytical methods, highlighting not only the potential of various approaches but also methodological pitfalls.     

When mechanism matters: Bayesian forecasting using models of ecological diffusion

Ecological diffusion is a theory that can be used to understand and forecast spatio‐temporal processes such as dispersal, invasion, and the spread of disease. Hierarchical Bayesian modelling provides a framework to make statistical inference and probabilistic forecasts, using mechanistic ecological models. To illustrate, we show how hierarchical Bayesian models of ecological diffusion can be implemented for large data sets that are distributed densely across space and time. The hierarchical Bayesian approach is used to understand and forecast the growth and geographic spread in the prevalence of chronic wasting disease in white‐tailed deer (Odocoileus virginianus). We compare statistical inference and forecasts from our hierarchical Bayesian model to phenomenological regression‐based methods that are commonly used to analyse spatial occurrence data. The mechanistic statistical model based on ecological diffusion led to important ecological insights, obviated a commonly ignored type of collinearity, and was the most accurate method for forecasting.     

Statistical ecology comes of age

The desire to predict the consequences of global environmental change has been the driver towards more realistic models embracing the variability and uncertainties inherent in ecology. Statistical ecology has gelled over the past decade as a discipline that moves away from describing patterns towards modelling the ecological processes that generate these patterns. Following the fourth International Statistical Ecology Conference (1–4 July 2014) in Montpellier, France, we analyse current trends in statistical ecology. Important advances in the analysis of individual movement, and in the modelling of population dynamics and species distributions, are made possible by the increasing use of hierarchical and hidden process models. Exciting research perspectives include the development of methods to interpret citizen science data and of efficient, flexible computational algorithms for model fitting. Statistical ecology has come of age: it now provides a general and mathematically rigorous framework linking ecological theory and empirical data.     

Applicability of the Poisson distribution to model the data of the German Children's Cancer Registry

Since 1980 the German Children's Cancer Registry has documented all childhood malignancies in the Federal Republic of Germany. Various statistical procedures have been proposed to identify municipalities or other geographic units with increased numbers of malignancies. Usually the Poisson distribution, which requires the malignancies to be distributed homogeneously and uncorrelated, is applied. Other discrete statistical distributions (so-called cluster distributions) like the generalized or compound Poisson distributions are applicable more generally. In this paper we present a first explorative approach to the question of whether it is necessary to use one of these cluster distributions to model the data of the German Children's Cancer Registry. In conclusion, we find no indication that the Poisson approach is insufficient.     

Mapping natural capital: optimising the use of national scale datasets

Understanding the spatial distribution of specific environmental variables and the interdependencies of these variables is crucial for managing the environment in a sustainable way. Here we discuss two methods of mapping – a Geographical Information System classification‐based approach and a statistical model‐based approach. If detailed, spatially comprehensive covariate datasets exist to complement the ecological‐response data, then using a statistical model‐based analysis provides the potential for greater understanding of underlying relationships, as well as the uncertainty in the spatial predictions. Further, the model‐based approach facilitates scenario testing. Although similar methods are already adopted in species distribution modeling, the flexibility of the model framework used is rarely exploited to go beyond modeling occupancy or suitability for a single species, into modeling complex derived metrics such as community composition and indicators of natural capital. As an example, we assess the potential benefits of the statistical model‐based approach to mapping natural capital through the use of two national survey datasets; The Centre for Ecology and Hydrology (CEH) Land Cover Map (LCM) and the British Geological Survey's (BGS) Parent Material Model (PMM), to predict national soil microbial community distributions based on data from a sample of > 1000 soils covering Great Britain. The results are mapped and compared against a more traditional, land classification‐based approach. The comparison shows that, although the maps look broadly similar, the model‐based approach provides better overall spatial prediction, and the contribution of individual model terms (along with their uncertainty) are far easier to understand and interpret, whilst also facilitating any scenario testing. We therefore both recommend the use of spatial statistical modelling techniques to map natural capital and anticipate that they will become more prominent over the forthcoming years.     

Incorporating ecological principles into statistical models for the prediction of species’ distribution and abundance

Understanding the determinants of species’ distributions and abundances is a central theme in ecology. The development of statistical models to achieve this has a long history and the notion that the model should closely reflect underlying scientific understanding has encouraged ecologists to adopt complex statistical methods as they arise. In this paper we describe a Bayesian hierarchical model that reflects a conceptual ecological model of multi‐scaled environmental determinants of riverine fish species’ distributions and abundances. We illustrate this with distribution and abundance data of a small‐bodied fish species, the Empire gudgeon Hypseleotris galii, in the Mary and Albert Rivers, Queensland, Australia. Specifically, the model sought to address; 1) the extent that landscape‐scale abiotic variables can explain the species’ distribution compared to local‐scale variables, 2) how local‐scale abiotic variables can explain species’ abundances, and 3) how are these local‐scale relationships mediated by landscape‐scale variables. Overall, the model accounted for around 60% of variation in the distribution and abundance of H. galii. The findings show that the landscape‐scale variables explain much of the distribution of the species; however, there was considerable improvement in estimating the species’ distribution with the addition of local‐scale variables. There were many strong relationships between abundance and local‐scale abiotic variables; however, several of these relationships were mediated by some of the landscape‐scale variables. The extent of spatial autocorrelation in the data was relatively low compared to the distances among sampling reaches. Our findings exemplify that Bayesian statistical modelling provides a robust framework for statistical modelling that reflects our ecological understanding. This allows ecologists to address a range of ecological questions with a single unified probability model rather than a series of disconnected analyses.     

Intrinsic scaling complexity in animal dispersion and abundance

Ecological theory related to animal distribution and abundance is at present incomplete and to some extent naive. We suggest that this may partly be due to a long tradition in the field of model development for choosing mathematical and statistical tools for convenience rather than applicability. Real population dynamics are influenced by nonlinear interactions, nonequilibrium conditions, and scaling complexity from system openness. Thus, a coherent theory for individual-, population-, and community-level processes should rest on mathematical and statistical methods that explicitly confront these issues in a manner that satisfies principles from statistical mechanics for complex systems. Instead, ecological theory is traditionally based on premises from simpler statistical mechanical theory for memory-free, scale-specific, random-walk, and diffusion processes, while animals from many taxa generally express strategic homing, site fidelity, and conspecific attraction in direct violation of primary model assumptions. Thus, the main challenge is to generalize the theory for memory-free physical, many-body systems to include a more realistic memory-influenced framework that better satisfies ecological realism. We describe, simulate, and discuss three testable aspects of a model for multiscaled habitat use at the individual level: (1) scale-free distribution of movement steps under influence of self-reinforcing site fidelity, (2) fractal spatial dispersion of intra-home range relocations, and (3) nonasymptotic expansion of observed intra-home range patch use with increasing set of relocations. Examples of literature data apparently supporting the conjecture that multiscaled, strategic space use is widespread among many animal taxa are also described. We suggest that the present approach, which provides a protocol to test for influence from scale-free, memory-dependent habitat use at the individual level, may also point toward a guideline for development of a generalized theoretical framework for complex population kinetics and spatiotemporal population dynamics.     

A new multivariate zero‐adjusted Poisson model with applications to biomedicine

Recently, although advances were made on modeling multivariate count data, existing models really has several limitations: (i) The multivariate Poisson log‐normal model (Aitchison and Ho, 1989) cannot be used to fit multivariate count data with excess zero‐vectors; (ii) The multivariate zero‐inflated Poisson (ZIP) distribution (Li et al., 1999) cannot be used to model zero‐truncated/deflated count data and it is difficult to apply to high‐dimensional cases; (iii) The Type I multivariate zero‐adjusted Poisson (ZAP) distribution (Tian et al., 2017) could only model multivariate count data with a special correlation structure for random components that are all positive or negative. In this paper, we first introduce a new multivariate ZAP distribution, based on a multivariate Poisson distribution, which allows the correlations between components with a more flexible dependency structure, that is some of the correlation coefficients could be positive while others could be negative. We then develop its important distributional properties, and provide efficient statistical inference methods for multivariate ZAP model with or without covariates. Two real data examples in biomedicine are used to illustrate the proposed methods.     

Accounting for dispersal and biotic interactions to disentangle the drivers of species distributions and their abundances

Although abiotic factors, together with dispersal and biotic interactions, are often suggested to explain the distribution of species and their abundances, species distribution models usually focus on abiotic factors only. We propose an integrative framework linking ecological theory, empirical data and statistical models to understand the distribution of species and their abundances together with the underlying community assembly dynamics. We illustrate our approach with 21 plant species in the French Alps. We show that a spatially nested modelling framework significantly improves the model's performance and that the spatial variations of species presence-absence and abundances are predominantly explained by different factors. We also show that incorporating abiotic, dispersal and biotic factors into the same model bring new insights to our understanding of community assembly. This approach, at the crossroads between community ecology and biogeography, is a promising avenue for a better understanding of species co-existence and biodiversity distribution.     

Prospects and challenges for parametric models in historical biogeographical inference

In historical biogeography, phylogenetic trees have long been used as tools for addressing a wide range of inference problems, from explaining common distribution patterns of species to reconstructing ancestral geographic ranges on branches of the tree of life. However, the potential utility of phylogenies for this purpose has yet to be fully realized, due in part to a lack of explicit conceptual links between processes underlying the evolution of geographic ranges and processes of phylogenetic tree growth. We suggest that statistical approaches that use parametric models to forge such links will stimulate integration and propel hypothesis-driven biogeographical inquiry in new directions. We highlight here two such approaches and describe how they represent early steps towards a more general framework for model-based historical biogeography that is based on likelihood as an optimality criterion, rather than having the traditional reliance on parsimony. The development of this framework will not be without significant challenges, particularly in balancing model complexity with statistical power, and these will be most apparent in studies of regions with many component areas and complex geological histories, such as the Mediterranean Basin.     

The gambin model provides a superior fit to species abundance distributions with a single free parameter: evidence,implementation and interpretation

The species abundance distribution (SAD) has been a central focus of community ecology for over fifty years, and is currently the subject of widespread renewed interest. The gambin model has recently been proposed as a model that provides a superior fit to commonly preferred SAD models. It has also been argued that the model's single parameter (α) presents a potentially informative ecological diversity metric, because it summarises the shape of the SAD in a single number. Despite this potential, few empirical tests of the model have been undertaken, perhaps because the necessary methods and software for fitting the model have not existed. Here, we derive a maximum likelihood method to fit the model, and use it to undertake a comprehensive comparative analysis of the fit of the gambin model. The functions and computational code to fit the model are incorporated in a newly developed free‐to‐download R package (gambin). We test the gambin model using a variety of datasets and compare the fit of the gambin model to fits obtained using the Poisson lognormal, logseries and zero‐sum multinomial distributions. We found that gambin almost universally provided a better fit to the data and that the fit was consistent for a variety of sample grain sizes. We demonstrate how α can be used to differentiate intelligibly between community structures of Azorean arthropods sampled in different land use types. We conclude that gambin presents a flexible model capable of fitting a wide variety of observed SAD data, while providing a useful index of SAD form in its single fitted parameter. As such, gambin has wide potential applicability in the study of SADs, and ecology more generally.     

How to make more out of community data? A conceptual framework and its implementation as models and software

Community ecology aims to understand what factors determine the assembly and dynamics of species assemblages at different spatiotemporal scales. To facilitate the integration between conceptual and statistical approaches in community ecology, we propose Hierarchical Modelling of Species Communities (HMSC) as a general, flexible framework for modern analysis of community data. While non‐manipulative data allow for only correlative and not causal inference, this framework facilitates the formulation of data‐driven hypotheses regarding the processes that structure communities. We model environmental filtering by variation and covariation in the responses of individual species to the characteristics of their environment, with potential contingencies on species traits and phylogenetic relationships. We capture biotic assembly rules by species‐to‐species association matrices, which may be estimated at multiple spatial or temporal scales. We operationalise the HMSC framework as a hierarchical Bayesian joint species distribution model, and implement it as R‐ and Matlab‐packages which enable computationally efficient analyses of large data sets. Armed with this tool, community ecologists can make sense of many types of data, including spatially explicit data and time‐series data. We illustrate the use of this framework through a series of diverse ecological examples.     

Ecological traps and species distribution models: a challenge for prioritizing areas of conservation importance

Species distribution models analyse how species use different types of habitats. Their spatial predictions are often used to prioritize areas for conservation. Individuals may, however, prefer settling in habitat types of low quality compared to other available habitats. This ecological trap phenomenon is usually studied in a small number of habitat patches and consequences at the landscape level are largely unknown. It is therefore often unclear whether the spatial pattern of habitat use is aligned with the behavioural decisions made by the individuals during habitat selection or reflects actual variation in the quality of different habitat types. As species distribution models analyse the pattern of occurrence in different habitats, there is a conservation interest in examining what their predictions mean in terms of habitat quality when ecological traps are operating. Previous work in Belgium showed that red-backed shrikes Lanius collurio are more attracted to newly available clear-cut habitat in plantation forests than to the traditionally used farmland habitat. We developed models with shrike distribution data and compared their predictions with spatial variation in shrike reproductive performance used as a proxy for habitat quality. Models accurately predicted shrike distribution and identified the preferred clear-cut patches as the most frequently used habitat, but reproductive performance was lower in clear-cut areas than in farmland. With human-induced rapid environmental changes, organisms may indeed be attracted to low-quality habitats and occupy them at high densities. Consequently, the predictions of statistical models based on occurrence records may not align with variation in significant population parameters for the maintenance of the species. When species expand their range to novel habitats, such models are useful to document the spatial distribution of the organisms, but data on population growth rates are worth collecting before using model predictions to guide the spatial prioritization of conservation actions.     

基于不同零模型的点格局分析

在种群空间格局研究中,定量分析格局及其形成过程已成为生态学家的主要目标。在量化分析的众多方法中,点格局分析是最常用的方法,而在选择零模型时,完全空间随机模型以外的复杂零模型很少使用,实际上,这些零模型可能有助于认识格局的内在特征。为此,我们在研究实例中,选择完全空间随机模型(complete spatial randomness)、泊松聚块模型(Poisson cluster process)和嵌套双聚块模型(nested double-cluster process)对典型草原处于不同恢复演替阶段的羊草(Leymus chinensis)种群空间格局进行了分析。结果发现:完全空间随机模型仅能检测种群在不同尺度下的格局类型;而通过泊松聚块模型和嵌套双聚块模型检验表明,在恢复演替的初期阶段,羊草种群在小尺度范围内偏离泊松聚块模型,而在整个取样范围内完全符合嵌套双聚块模型;随着恢复演替时间的推移,在恢复演替的后期,在整个取样尺度上,羊草种群与泊松聚块模型相吻合。这是很有意义的生态学现象。这一实例表明在应用点格局分析种群空间格局时,仅通过完全空间随机模型的检验来分析格局特征,或许很难论证复杂的生态过程,而选择一些完全空间随机模型以外的较复杂的零模型,可能发现一些有价值的生态学现象,对揭示格局掩盖下的内在机制有所裨益。     

Chebanov law and Vakar formula in mathematical models of complex systems

Ecological models written in a mathematical language L(M) or model language, with a given style or methodology can be considered as a text. It is possible to apply statistical linguistic laws and the experimental results demonstrate that the behaviour of a mathematical model is the same of any literary text of any natural language. A text has the following characteristics: (a) the variables, its transformed functions and parameters are the lexic units or LUN of ecological models; (b) the syllables are constituted by a LUN, or a chain of them, separated by operating or ordering LUNs; (c) the flow equations are words; and (d) the distribution of words (LUM and CLUN) according to their lengths is based on a Poisson distribution, the Chebanov's law. It is founded on Vakar's formula, that is calculated likewise the linguistic entropy for L(M). We will apply these ideas over practical examples using MARIOLA model. In this paper it will be studied the problem of the lengths of the simple lexic units composed lexic units and words of text models, expressing these lengths in number of the primitive symbols, and syllables. The use of these linguistic laws renders it possible to indicate the degree of information given by an ecological model.     

Embracing asymmetry in nature: How to account for skewness in ecological data

Asymmetric regression is an alternative to conventional linear regression that allows us to model the relationship between predictor variables and the response variable while accommodating skewness. Advantages of asymmetric regression include incorporating realistic ecological patterns observed in data, robustness to model misspecification and less sensitivity to outliers. Bayesian asymmetric regression relies on asymmetric distributions such as the asymmetric Laplace (ALD) or asymmetric normal (AND) in place of the normal distribution used in classic linear regression models. Asymmetric regression concepts can be used for process and parameter components of hierarchical Bayesian models and have a wide range of applications in data analyses. In particular, asymmetric regression allows us to fit more realistic statistical models to skewed data and pairs well with Bayesian inference. We first describe asymmetric regression using the ALD and AND. Second, we show how the ALD and AND can be used for Bayesian quantile and expectile regression for continuous response data. Third, we consider an extension to generalize Bayesian asymmetric regression to survey data consisting of counts of objects. Fourth, we describe a regression model using the ALD, and show that it can be applied to add needed flexibility, resulting in better predictive models compared to Poisson or negative binomial regression. We demonstrate concepts by analyzing a data set consisting of counts of Henslow’s sparrows following prescribed fire and provide annotated computer code to facilitate implementation. Our results suggest Bayesian asymmetric regression is an essential component of a scientist’s statistical toolbox.