Fitting population dynamic models to time-series data by gradient matching

We describe and test a method for fitting noisy differential equation models to a time series of population counts, motivated by stage-structured models of insect and zooplankton populations. We consider semimechanistic models, in which the model struc- ture is derived from knowledge of the life cycle, but the rate equations are estimated nonparametrically from the time-series data. The method involves smoothing the population time series x(t) in order to estimate the gradient dx/dt, and then fitting rate equations using penalized regression splines. Computer-intensive methods are used to estimate and remove the biases that result from the data being discrete time samples with sampling errors from a continuous time process. Semimechanistic modeling makes it possible to test assumptions about the mechanisms behind population fluctuations without the results being confounded by possibly arbitrary choices of parametric forms for process-rate equations. To illustrate this application, we analyze time-series data on laboratory populations of blowflies Lucilia cuprina and Lucilia sericata. The models assume that the populations are limited by com- petition among adults affecting their current birth and death rates. The results correspond to the actual experimental conditions. For L. cuprina (where the model's structure is ap- propriate) a good fit can be obtained, while for L. sericata (where the model is inappropriate), the fitted model does not reproduce some major features of the observed cycles. A docu- mented set of R functions for all steps in the model-fitting process is provided as a sup- plement to this article.

[1]  A. Nicholson.  An outline of the dynamics of animal populations. , 1954 .

[2]  A. Nicholson,et al.  The Self-Adjustment of Populations to Change , 1957 .

[3]  J. Readshaw,et al.  A Model of Nicholson's Blowfly Cycles and its Relevance to Predation Theory , 1980 .

[4]  S. P. Blythe,et al.  Nicholson's blowflies revisited , 1980, Nature.

[5]  J. Readshaw,et al.  Age-specific survival, fecundity and fertility of the adult blowfly, Lucilia cuprina, in relation to crowding, protein food and population cycles , 1983 .

[6]  William Gurney,et al.  Fluctuation periodicity, generation separation, and the expression of larval competition , 1985 .

[7]  W. Gurney,et al.  Parameter evolution in a laboratory insect population , 1988 .

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

[9]  L. Olsen,et al.  Chaos versus noisy periodicity: alternative hypotheses for childhood epidemics. , 1990, Science.

[10]  Larry L. Schumaker,et al.  Spline Models for Observational Data (Grace Wahba) , 1991, SIAM Review.

[11]  A. Gallant,et al.  Finding Chaos in Noisy Systems , 1992 .

[12]  W M Schaffer,et al.  The case for chaos in childhood epidemics. II. Predicting historical epidemics from mathematical models , 1993, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[13]  S. Wood,et al.  Obtaining birth and mortality patterns from structured population trajectories. , 1994, Ecological monographs.

[14]  Brian Dennis,et al.  DENSITY DEPENDENCE IN TIME SERIES OBSERVATIONS OF NATURAL POPULATIONS: ESTIMATION AND TESTING' , 1994 .

[15]  William W. Murdoch,et al.  POPULATION REGULATION IN THEORY AND PRACTICE , 1994 .

[16]  J. R. Cook,et al.  Simulation-Extrapolation Estimation in Parametric Measurement Error Models , 1994 .

[17]  J. R. Cook,et al.  Simulation-Extrapolation: The Measurement Error Jackknife , 1995 .

[18]  C. Pfister.  Estimating Competition Coefficients from Census Data: A Test with Field Manipulations of Tidepool Fishes , 1995, The American Naturalist.

[19]  D. Ruppert,et al.  Measurement Error in Nonlinear Models , 1995 .

[20]  S. Wood,et al.  Mortality estimation for planktonic copepods: Pseudocalanus newmani in a temperate fjord , 1996 .

[21]  William Gurney,et al.  Structured Population Models of Herbivorous Zooplankton , 1996 .

[22]  A. Hastings.  Population Biology: Concepts and Models , 1996 .

[23]  Paul H. C. Eilers,et al.  Flexible smoothing with B-splines and penalties , 1996 .

[24]  W. Murdoch,et al.  Theory for Biological Control: Recent Developments , 1996 .

[25]  Ross Ihaka,et al.  Gentleman R: R: A language for data analysis and graphics , 1996 .

[26]  Michael Bruce. Forrest.  Toxins and blowfly population dynamics. , 1996 .

[27]  R. Hilborn,et al.  The Ecological Detective: Confronting Models with Data , 1997 .

[28]  M. J. Hatcher,et al.  Modeling Biological Systems: Principles and Applications , 1997 .

[29]  Simon N. Wood,et al.  Inverse Problems and Structured-Population Dynamics , 1997 .

[30]  Stephen P. Ellner,et al.  Inferring mechanism from time-series data: delay-differential equations , 1997 .

[31]  Shripad Tuljapurkar,et al.  Structured-Population Models in Marine, Terrestrial, and Freshwater Systems , 1997, Population and Community Biology Series.

[32]  Erkki Korpimäki,et al.  Experimental reduction of predators reverses the crash phase of small-rodent cycles , 1998 .

[33]  A. Dobson,et al.  Prevention of population cycles by parasite removal. , 1998, Science.

[34]  A. R. Gallant,et al.  Noise and Nonlinearity in Measles Epidemics: Combining Mechanistic and Statistical Approaches to Population Modeling , 1998, The American Naturalist.

[35]  J. Timothy Wootton,et al.  THEORETICAL CONCEPTS AND EMPIRICAL APPROACHES TO MEASURING INTERACTION STRENGTH , 1998 .

[36]  H. Müller,et al.  Local Polynomial Modeling and Its Applications , 1998 .

[37]  Theoretical Approaches to Biological Control: Recent developments in theory for biological control of insect pests by parasitoids , 1999 .

[38]  Simon N. Wood,et al.  Super–sensitivity to structure in biological models , 1999, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[39]  B. Kendall,et al.  WHY DO POPULATIONS CYCLE? A SYNTHESIS OF STATISTICAL AND MECHANISTIC MODELING APPROACHES , 1999 .

[40]  D. Ruppert,et al.  Nonparametric regression in the presence of measurement error , 1999 .

[41]  David Ruppert,et al.  Penalized regression splines , 1999 .

[42]  Anthony R. Ives,et al.  COMMUNITY INTERACTION WEBS AND ZOOPLANKTON RESPONSES TO PLANKTIVORY MANIPULATIONS , 1999 .

[43]  Stephen P. Ellner,et al.  Living on the edge of chaos: population dynamics of fennoscandian voles , 2000 .

[44]  Susan Daniels,et al.  Blowflies as a Case Study in Non-Linear Population Dynamics , 2000 .

[45]  David Ruppert,et al.  Theory & Methods: Spatially‐adaptive Penalties for Spline Fitting , 2000 .

[46]  Stephen P. Ellner,et al.  Modelling Time-Series Data , 2000 .

[47]  J. Perry,et al.  Chaos in real data : analysis of non-linear dynamics from short ecological time-series , 2000 .

[48]  N. Stenseth,et al.  Exploring the density-dependent structure of blowfly populations by nonparametric additive modeling , 2001 .

[49]  Simon N. Wood,et al.  PARTIALLY SPECIFIED ECOLOGICAL MODELS , 2001 .

[50]  Semi-parametric population models , 2001 .

[51]  Alan Hastings,et al.  FITTING POPULATION MODELS INCORPORATING PROCESS NOISE AND OBSERVATION ERROR , 2002 .