UNCERTAINTY AND VARIABILITY IN DEMOGRAPHY AND POPULATION GROWTH: A HIERARCHICAL APPROACH

Estimates of uncertainty are the basis for inference of population risk. Uncertainty is estimated from models fitted to data that typically include a deterministic model (e.g., population growth) and stochastic elements, which should accommodate errors in sampling and any sources of variability that affect observations. Prediction from fitted models (of, say, demography) to new variables (say, population growth) requires propagation of these stochastic elements. Ecological models ignore most forms of variability, because they make statistical models complex, and they pose computational challenges. Variability associated with space, time, and among individuals that is not accommodated by demographic models can make parameter estimates and growth rate predictions unrealistic. I adapt a hierarchical approach to the problem of estimating population growth rates and their uncertainties when individuals vary and that variability cannot be assigned to specific causes. In contrast to an overfitted model that would assign a different parameter value to each individual, hierarchical models accommodate individual differences, but assume that those differences derive from an underlying distribution—they belong to a “population.” The hierarchical model can be implemented in classical (frequentist) and Bayesian frameworks (I demonstrate both) and analyzed using Markov chain Monte Carlo simulation. Results show that population growth models that rely on standard propagation of estimation error but ignore variability among individuals can misrepresent uncertainties in ways that erode credibility. Corresponding Editor: O. N. Bjornstad.

[1]  S. Schwartz The fallacy of the ecological fallacy: the potential misuse of a concept and the consequences. , 1994, American journal of public health.

[2]  G. Jolly EXPLICIT ESTIMATES FROM CAPTURE-RECAPTURE DATA WITH BOTH DEATH AND IMMIGRATION-STOCHASTIC MODEL. , 1965, Biometrika.

[3]  S. Ellner,et al.  SIZE‐SPECIFIC SENSITIVITY: APPLYING A NEW STRUCTURED POPULATION MODEL , 2000 .

[4]  James S. Clark,et al.  The relationship between growth and mortality for seven co‐occurring tree species in the southern Appalachian Mountains , 2002 .

[5]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[6]  C. Krebs,et al.  Can the Solar Cycle and Climate Synchronize the Snowshoe Hare Cycle in Canada? Evidence from Tree Rings and Ice Cores , 1993, The American Naturalist.

[7]  Alan B. Franklin,et al.  Population Regulation in Northern Spotted Owls: Theoretical Implications for Management , 1992 .

[8]  Jay M. Ver Hoef,et al.  Parametric Empirical Bayes Methods for Ecological Applications , 1996 .

[9]  B. Silverman Density estimation for statistics and data analysis , 1986 .

[10]  N. Stenseth,et al.  Cycles and trends in cod populations. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Stephen R. Carpenter,et al.  MICROCOSM EXPERIMENTS HAVE LIMITED RELEVANCE FOR COMMUNITY AND ECOSYSTEM ECOLOGY: REPLY , 1999 .

[12]  C. Krebs,et al.  What Drives the 10-year Cycle of Snowshoe Hares? , 2001 .

[13]  Christopher K. Wikle,et al.  Hierarchical Bayesian Models for Predicting The Spread of Ecological Processes , 2003 .

[14]  H. G. Andrewartha,et al.  The distribution and abundance of animals. , 1954 .

[15]  Bradley P. Carlin,et al.  BAYES AND EMPIRICAL BAYES METHODS FOR DATA ANALYSIS , 1996, Stat. Comput..

[16]  David R. Anderson,et al.  Modeling Survival and Testing Biological Hypotheses Using Marked Animals: A Unified Approach with Case Studies , 1992 .

[17]  David R. Anderson,et al.  CLIMATE, HABITAT QUALITY, AND FITNESS IN NORTHERN SPOTTED OWL POPULATIONS IN NORTHWESTERN CALIFORNIA , 2000 .

[18]  M. Kendall,et al.  Kendall's advanced theory of statistics , 1995 .

[19]  D. R. Strong,et al.  Density vagueness abiding the variance in the demography of real populations , 1986 .

[20]  A. Gelfand,et al.  Identifiability, Improper Priors, and Gibbs Sampling for Generalized Linear Models , 1999 .

[21]  Michael C. Dietze,et al.  COEXISTENCE: HOW TO IDENTIFY TROPHIC TRADE-OFFS , 2003 .

[22]  Scott L. Zeger,et al.  Generalized linear models with random e ects: a Gibbs sampling approach , 1991 .

[23]  J. Ogden,et al.  Population increase and nesting patterns of the black noddy Anous minutus in Pisonia forest on Heron Island: Observations in 1978, 1979 and 1992 , 1993 .

[24]  A. Diez-Roux,et al.  Bringing context back into epidemiology: variables and fallacies in multilevel analysis. , 1998, American journal of public health.

[25]  R F Costantino,et al.  Lattice Effects Observed in Chaotic Dynamics of Experimental Populations , 2001, Science.

[26]  John Vandermeer,et al.  Choosing category size in a stage projection matrix , 2004, Oecologia.

[27]  C. Morris Parametric Empirical Bayes Inference: Theory and Applications , 1983 .

[28]  Roland H. Lamberson,et al.  Conservation planning for species occupying fragmented landscapes: the case of the northern spotted owl , 1993 .

[29]  R. May,et al.  Stability and Complexity in Model Ecosystems , 1976, IEEE Transactions on Systems, Man, and Cybernetics.

[30]  S. E. Hills,et al.  Illustration of Bayesian Inference in Normal Data Models Using Gibbs Sampling , 1990 .

[31]  P. Holgate,et al.  Matrix Population Models. , 1990 .

[32]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[33]  Craig A. Stow,et al.  Fitting Predator-Prey Models to Time Series with Observation Errors , 1994 .

[34]  J. Clark,et al.  Interpreting recruitment limitation in forests. , 1999, American journal of botany.

[35]  R. Cousins Why isn’t every physicist a Bayesian? , 1995 .

[36]  Ransom A. Myers,et al.  Density-Dependent Juvenile Mortality in Marine Demersal Fish , 1993 .

[37]  Michael Lavine,et al.  Statistical modeling of seedling mortality , 2002 .

[38]  Kirk A. Moloney,et al.  A generalized algorithm for determining category size , 1986, Oecologia.

[39]  Peter Müller,et al.  INCORPORATING MULTIPLE SOURCES OF STOCHASTICITY INTO DYNAMIC POPULATION MODELS , 2003 .

[40]  Jonathan M. Chase,et al.  Spatial scale dictates the productivity–biodiversity relationship , 2002, Nature.

[41]  Peter Turchin,et al.  Population Regulation" Old Arguments and a New Synthesis , 1995 .

[42]  Taylor,et al.  Dynamical role of predators in population cycles of a forest insect: An experimental test , 1999, Science.

[43]  A. F. Smith,et al.  Conjugate likelihood distributions , 1993 .

[44]  R. Lande,et al.  Demographic models of the northern spotted owl (Strix occidentalis caurina) , 1988, Oecologia.

[45]  Thomas W.F. Stroud,et al.  Bayesian analysis of binary survey data , 1994 .