An official website of the United States government
Abstract Structured demographic models are among the most common and useful tools in population biology. However, the introduction of integral projection models (IPMs) has caused a profound shift in the way many demographic models are conceptualized. Some researchers have argued that IPMs, by explicitly representing demographic processes as continuous functions of state variables such as size, are more statistically efficient, biologically realistic, and accurate than classic matrix projection models, calling into question the usefulness of the many studies based on matrix models. Here, we evaluate how IPMs and matrix models differ, as well as the extent to which these differences matter for estimation of key model outputs, including population growth rates, sensitivity patterns, and life spans. First, we detail the steps in constructing and using each type of model. Second, we present a review of published demographic models, concentrating on size‐based studies, which shows significant overlap in the way IPMs and matrix models are constructed and analyzed. Third, to assess the impact of various modeling decisions on demographic predictions, we ran a series of simulations based on size‐based demographic data sets for five biologically diverse species. We found little evidence that discrete vital rate estimation is less accurate than continuous functions across a wide range of sample sizes or size classes (equivalently bin numbers or mesh points). Most model outputs quickly converged with modest class numbers (≥10), regardless of most other modeling decisions. Another surprising result was that the most commonly used method to discretize growth rates for IPM analyses can introduce substantial error into model outputs. Finally, we show that empirical sample sizes generally matter more than modeling approach for the accuracy of demographic outputs. Based on these results, we provide specific recommendations to those constructing and evaluating structured population models. Both our literature review and simulations question the treatment of IPMs as a clearly distinct modeling approach or one that is inherently more accurate than classic matrix models. Importantly, this suggests that matrix models, representing the vast majority of past demographic analyses available for comparative and conservation work, continue to be useful and important sources of demographic information.
more »
« less
This content will become publicly available on May 7, 2026
My, how you've grown: A practical guide to modeling size transitions for integral projection model ( IPM ) applications
Abstract Integral projection models (IPMs) are widely used for studying continuously size‐structured populations. IPMs require a growth sub‐model that describes the probability of future size conditional on current size and any covariates. Most IPM studies assume that this distribution is Gaussian, despite calls for non‐Gaussian models that accommodate skewness and excess kurtosis. We provide a general workflow for accommodating non‐Gaussian growth patterns while retaining important covariates and random effects. Our approach emphasizes visual diagnostics from pilot Gaussian models and quantile‐based metrics of skewness and kurtosis that guide selection of a non‐Gaussian alternative, if necessary. Across six case studies, skewness and excess kurtosis were common features of growth data, and non‐Gaussian models consistently generated simulated data that were more consistent with real data than pilot Gaussian models. However, effects of “improved” growth modeling on IPM results were moderate to weak and differed in direction or magnitude between different outputs from the same model. Using tools not available when IPMs were first developed, it is now possible to fit non‐Gaussian models to growth data without sacrificing ecological complexity. Doing so, as guided by careful interrogation of the data, will result in models that better represent the populations for which they are intended.
more »
« less
- PAR ID:
- 10594964
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Ecology
- Volume:
- 106
- Issue:
- 5
- ISSN:
- 0012-9658
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Structured population models are among the most widely used tools in ecology and evolution. Integral projection models (IPMs) use continuous representations of how survival, reproduction and growth change as functions of state variables such as size, requiring fewer parameters to be estimated than projection matrix models (PPMs). Yet, almost all published IPMs make an important assumption that size‐dependent growth transitions are or can be transformed to be normally distributed. In fact, many organisms exhibit highly skewed size transitions. Small individuals can grow more than they can shrink, and large individuals may often shrink more dramatically than they can grow. Yet, the implications of such skew for inference from IPMs has not been explored, nor have general methods been developed to incorporate skewed size transitions into IPMs, or deal with other aspects of real growth rates, including bounds on possible growth or shrinkage.Here, we develop a flexible approach to modelling skewed growth data using a modified beta regression model. We propose that sizes first be converted to a (0,1) interval by estimating size‐dependent minimum and maximum sizes through quantile regression. Transformed data can then be modelled using beta regression with widely available statistical tools. We demonstrate the utility of this approach using demographic data for a long‐lived plant, gorgonians and an epiphytic lichen. Specifically, we compare inferences of population parameters from discrete PPMs to those from IPMs that either assume normality or incorporate skew using beta regression or, alternatively, a skewed normal model.The beta and skewed normal distributions accurately capture the mean, variance and skew of real growth distributions. Incorporating skewed growth into IPMs decreases population growth and estimated life span relative to IPMs that assume normally distributed growth, and more closely approximate the parameters of PPMs that do not assume a particular growth distribution. A bounded distribution, such as the beta, also avoids the eviction problem caused by predicting some growth outside the modelled size range.Incorporating biologically relevant skew in growth data has important consequences for inference from IPMs. The approaches we outline here are flexible and easy to implement with existing statistical tools.more » « less
-
Abstract In many applications of data assimilation, especially when the size of the problem is large, a substantial assumption is made: all variables are well‐described by Gaussian error statistics. This assumption has the advantage of making calculations considerably simpler, but it is often not valid, leading to biases in forecasts or, even worse, unphysical predictions. We propose a simple, but effective, way of replacing this assumption, by making use of transforming functions, while remaining consistent with Bayes' theorem. This method allows the errors to have any value of the skewness and kurtosis, and permits physical bounds for the variables. As such, the error distribution can conform better to the underlying statistics, reducing biases introduced by the Gaussian assumption. We apply this framework to a 3D variational data assimilation method, and find improved performance in a simple atmospheric toy model (Lorenz‐63), compared to an all‐Gaussian technique.more » « less
-
The isothermal compressibility (i.e., related to the asymptotic number variance) of equilibrium liquid water as a function of temperature is minimal under near-ambient conditions. This anomalous non-monotonic temperature dependence is due to a balance between thermal fluctuations and the formation of tetrahedral hydrogen-bond networks. Since tetrahedrality is a many-body property, it will also influence the higher-order moments of density fluctuations, including the skewness and kurtosis. To gain a more complete picture, we examine these higher-order moments that encapsulate many-body correlations using a recently developed, advanced platform for local density fluctuations. We study an extensive set of simulated phases of water across a range of temperatures (80–1600 K) with various degrees of tetrahedrality, including ice phases, equilibrium liquid water, supercritical water, and disordered nonequilibrium quenches. We find clear signatures of tetrahedrality in the higher-order moments, including the skewness and excess kurtosis, which scale for all cases with the degree of tetrahedrality. More importantly, this scaling behavior leads to non-monotonic temperature dependencies in the higher-order moments for both equilibrium and non-equilibrium phases. Specifically, under near-ambient conditions, the higher-order moments vanish most rapidly for large length scales, and the distribution quickly converges to a Gaussian in our metric. However, under non-ambient conditions, higher-order moments vanish more slowly and hence become more relevant, especially for improving information-theoretic approximations of hydrophobic solubility. The temperature non-monotonicity that we observe in the full distribution across length scales could shed light on water’s nested anomalies, i.e., reveal new links between structural, dynamic, and thermodynamic anomalies.more » « less
-
Browman, Howard (Ed.)Ocean acidification may impact the fitness of marine fish, however, studies reporting neutral to moderate effects have mostly performed short-term exposures to elevated CO2, whereas longer-term studies across life stages are still scarce. We performed a CO2 exposure experiment, in which a large number (n > 2200) of Atlantic silverside Menidia menidia offspring from wild spawners were reared for 135 days through their embryonic, larval, and juvenile stages under control (500 µatm) and high CO2 conditions (2300 µatm). Although survival was high across treatments, subtle but significant differences in length, weight, condition factor and fatty acid (FA) composition were observed. On average, fish from the acidified treatment were 4% shorter and weighed 6% less, but expressed a higher condition factor than control juveniles. In addition, the metrics of length and weight distributions differed significantly, with juveniles from the high CO2 treatment occupying more extreme size classes and the length distribution shifting to a positive kurtosis. Six of twenty-seven FAs differed significantly between treatments. Our results suggest that high CO2 conditions alter long-term growth in M. menidia, particularly in the absence of excess food. It remains to be shown whether and how these differences will impact fish populations in the wild facing size-selective predation and seasonally varying prey abundance.more » « less
