Mixed effects: a unifying framework for statistical modelling in fisheries biology

Fisheries biology encompasses a tremendous diversity of research questions, methods, and models. Many sub-fields use observational or experimental data to make inference about biological characteristics that are not directly observed (called “latent states”), such as heritability of phenotypic traits, habitat suitability, and population densities to name a few. Latent states will generally cause model residuals to be correlated, violating the assumptionof statistical independencemade inmany statisticalmodelling approaches. In this exposition,we argue thatmixed-effectmodelling (i) is an important and generic solution to non-independence caused by latent states; (ii) provides a unifying framework for disparate statistical methods such as time-series, spatial, and individual-based models; and (iii) is increasingly practical to implement and customize for problemspecific models. We proceed by summarizing the distinctions between fixed and random effects, reviewing a generic approach for parameter estimation, and distinguishing general categories of non-linear mixed-effect models. We then provide four worked examples, including state-space, spatial, individual-level variability, and quantitative genetics applications (with working code for each), while providing comparison with conventional fixed-effect implementations. We conclude by summarizing directions for future research in this important framework for modelling and statistical analysis in fisheries biology.

[1]  Giovanni Petris,et al.  Dynamic linear models , 2009 .

[2]  J. Thorson,et al.  The importance of spatial models for estimating the strength of density dependence. , 2015, Ecology.

[3]  R. Beverton,et al.  A Review of the Lifespans and Mortality Rates of Fish in Nature, and Their Relation to Growth and Other Physiological Characteristics , 2008 .

[4]  André E. Punt,et al.  Model performance analysis for Bayesian biomass dynamics models using bias, precision and reliability metrics , 2012 .

[5]  Brian Dennis,et al.  Bayesian methods for hierarchical models: are ecologists making a Faustian bargain? , 2009, Ecological applications : a publication of the Ecological Society of America.

[6]  Stephen R. Carpenter,et al.  ECOLOGICAL FUTURES: BUILDING AN ECOLOGY OF THE LONG NOW1 , 2002 .

[7]  R. Law On the quantitative genetics of correlated characters under directional selection in age-structured populations , 1991 .

[8]  Terrance J. Quinn,et al.  Quantitative Fish Dynamics , 1999 .

[9]  E. Charnov,et al.  Evolutionary assembly rules for fish life histories , 2013 .

[10]  Catherine A Calder,et al.  Accounting for uncertainty in ecological analysis: the strengths and limitations of hierarchical statistical modeling. , 2009, Ecological applications : a publication of the Ecological Society of America.

[11]  L. Madsen,et al.  Models for Estimating Abundance from Repeated Counts of an Open Metapopulation , 2011, Biometrics.

[12]  Randall M. Peterman,et al.  Comparison of parameter estimation methods for detecting climate-induced changes in productivity of Pacific salmon (Oncorhynchus spp.) , 2000 .

[13]  Haran Murali,et al.  Gaussian Random Field Models for Spatial Data , 2011 .

[14]  A. McKane,et al.  Stochastic formulation of ecological models and their applications. , 2012, Trends in ecology & evolution.

[15]  J. Thorson,et al.  Spatial semiparametric models improve estimates of species abundance and distribution , 2014 .

[16]  Russell B. Millar,et al.  Bayesian state-space modeling of age-structured data: fitting a model is just the beginning , 2000 .

[17]  L. Mark Berliner,et al.  Hierarchical Bayesian Time Series Models , 1996 .

[18]  Perry de Valpine Shared challenges and common ground for Bayesian and classical analysis of hierarchical statistical models. , 2009, Ecological applications : a publication of the Ecological Society of America.

[19]  George Sugihara,et al.  Detecting Causality in Complex Ecosystems , 2012, Science.

[20]  Zhenming Su,et al.  Spatial hierarchical Bayesian models for stock-recruitment analysis of pink salmon (Oncorhynchus gorbuscha) , 2004 .

[21]  Hans J. Skaug,et al.  Automatic approximation of the marginal likelihood in non-Gaussian hierarchical models , 2006, Comput. Stat. Data Anal..

[22]  J. Andrew Royle,et al.  Hierarchical Modeling and Inference in Ecology: The Analysis of Data from Populations, Metapopulations and Communities , 2008 .

[23]  J. Andrew Royle N‐Mixture Models for Estimating Population Size from Spatially Replicated Counts , 2004, Biometrics.

[24]  D. P. Swain,et al.  Stocks and chameleons: the use of phenotypic variation in stock identification , 1999 .

[25]  Robert D. Holt,et al.  RESOLVING ECOLOGICAL QUESTIONS THROUGH META‐ANALYSIS: GOALS, METRICS, AND MODELS , 1999 .

[26]  John Hinde,et al.  Statistical Modelling in R , 2009 .

[27]  Ray Hilborn,et al.  THE DARK SIDE OF REFERENCE POINTS , 2002 .

[28]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[29]  C. Walters,et al.  Quantitative fisheries stock assessment: Choice, dynamics and uncertainty , 2004, Reviews in Fish Biology and Fisheries.

[30]  I. MacKenzie Occupancy estimation and modeling , 2013 .

[31]  Christopher K. Wikle,et al.  Hierarchical Bayesian Spatio-Temporal Conway–Maxwell Poisson Models with Dynamic Dispersion , 2013 .

[32]  Devin S Johnson,et al.  Continuous-time correlated random walk model for animal telemetry data. , 2008, Ecology.

[33]  Ian D. Jonsen,et al.  ROBUST STATE-SPACE MODELING OF ANIMAL MOVEMENT DATA , 2005 .

[34]  Noel A Cressie,et al.  Statistics for Spatio-Temporal Data , 2011 .

[35]  Anders Nielsen,et al.  Estimation of time-varying selectivity in stock assessments using state-space models , 2014 .

[36]  A. Punt Refocusing Stock Assessment in Support of Policy Evaluation , 2008 .

[37]  André E. Punt,et al.  Extending production models to include process error in the population dynamics , 2003 .

[38]  B. Worm,et al.  Productivity dynamics of Atlantic cod , 2014 .

[39]  Michael Schaub,et al.  Bayesian Population Analysis using WinBUGS: A Hierarchical Perspective , 2011 .

[40]  Richard D. Methot,et al.  Hiding or dead? A computationally efficient model of selective fisheries mortality , 2013 .

[41]  D. Bates,et al.  Mixed-Effects Models in S and S-PLUS , 2001 .

[42]  Stephan B. Munch,et al.  Sustaining Fisheries Yields Over Evolutionary Time Scales , 2002, Science.

[43]  M. Mangel,et al.  Separating Intrinsic and Environmental Contributions to Growth and Their Population Consequences , 2013, The American Naturalist.

[44]  Henrik Madsen,et al.  Estimation methods for nonlinear state-space models in ecology , 2011 .

[45]  K. Andersen,et al.  Asymptotic Size Determines Species Abundance in the Marine Size Spectrum , 2006, The American Naturalist.

[46]  D. Fraser,et al.  Multigenerational hybridisation and its consequences for maternal effects in Atlantic salmon , 2013, Heredity.

[47]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[48]  Tyler Wagner,et al.  Accounting for Multilevel Data Structures in Fisheries Data using Mixed Models FEATURE: FISHERIES RESEARCH , 2006 .

[49]  G. Gudmundsson,et al.  Selection and estimation of sequential catch-at-age models , 2012 .

[50]  S. Levin,et al.  Theories of Simplification and Scaling of Spatially Distributed Processes , 2011 .

[51]  Mathew A. Leibold,et al.  Metacommunities: Spatial Dynamics and Ecological Communities , 2005 .

[52]  Russell B. Millar,et al.  Remedies for pseudoreplication , 2004 .

[53]  H. Rue,et al.  An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach , 2011 .

[54]  Casper W. Berg Theta logistic population model writeup , 2012 .

[55]  Martin I. Taylor,et al.  Experimental harvesting of fish populations drives genetically based shifts in body size and maturation , 2013 .

[56]  M. Morgan,et al.  Does increased information about reproductive potential result in better prediction of recruitment , 2011 .

[57]  Ransom A. Myers,et al.  When Do Environment–recruitment Correlations Work? , 1998, Reviews in Fish Biology and Fisheries.

[58]  Jon T. Schnute,et al.  The importance of noise in fish population models , 1991 .

[59]  Martyn Plummer,et al.  JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling , 2003 .

[60]  Elise F. Zipkin,et al.  How variable is recruitment for exploited marine fishes? A hierarchical model for testing life history theory , 2014 .

[61]  A. Zuur,et al.  Dynamic factor analysis to estimate common trends in fisheries time series , 2003 .

[62]  David Peel,et al.  A Model-Based Approach to Designing a Fishery-Independent Survey , 2013 .

[63]  A. Gelman Analysis of variance: Why it is more important than ever? , 2005, math/0504499.

[64]  W. Tobler A Computer Movie Simulating Urban Growth in the Detroit Region , 1970 .

[65]  James S. Clark,et al.  A future for models and data in environmental science. , 2006, Trends in ecology & evolution.

[66]  J. Andrew Royle,et al.  Efficient statistical mapping of avian count data , 2005, Environmental and Ecological Statistics.

[67]  H. Akaike A new look at the statistical model identification , 1974 .

[68]  Haavard Rue,et al.  A toolbox for fitting complex spatial point process models using integrated nested Laplace approximation (INLA) , 2012, 1301.1817.

[69]  Giovanni Petris,et al.  Dynamic Linear Models with R , 2009 .

[70]  John C. Nash,et al.  Strategies for fitting nonlinear ecological models in R, AD Model Builder, and BUGS , 2013 .

[71]  Erik Postma,et al.  An ecologist's guide to the animal model. , 2010, The Journal of animal ecology.

[72]  S. R. Searle,et al.  The estimation of environmental and genetic trends from records subject to culling. , 1959 .

[73]  Andrew Gelman,et al.  Handbook of Markov Chain Monte Carlo , 2011 .

[74]  Ian J. Stewart,et al.  Development and application of an agent-based model to evaluate methods for estimating relative abundance indices for shoaling fish such as Pacific rockfish (Sebastes spp.) , 2012, ICES Journal of Marine Science.

[75]  D. MacKenzie Occupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species Occurrence , 2005 .

[76]  John Skilling,et al.  Maximum Entropy and Bayesian Methods , 1989 .

[77]  Jameal F. Samhouri,et al.  Giants' shoulders 15 years later: lessons, challenges and guidelines in fisheries meta‐analysis , 2015 .

[78]  R. Frankham Introduction to quantitative genetics (4th edn): by Douglas S. Falconer and Trudy F.C. Mackay Longman, 1996. £24.99 pbk (xv and 464 pages) ISBN 0582 24302 5 , 1996 .

[79]  Andrew Gelman,et al.  Data Analysis Using Regression and Multilevel/Hierarchical Models , 2006 .

[80]  U. Dieckmann,et al.  Managing Evolving Fish Stocks , 2007 .

[81]  Tim Coulson,et al.  Towards a general, population-level understanding of eco-evolutionary change. , 2013, Trends in ecology & evolution.

[82]  Jonas Knape,et al.  ESTIMABILITY OF DENSITY DEPENDENCE IN MODELS OF TIME SERIES DATA. , 2008, Ecology.

[83]  K. Radway Allen,et al.  Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations , 1996 .

[84]  C. Clark,et al.  Dynamic State Variable Models in Ecology: Methods and Applications , 2001 .

[85]  Russell B. Millar,et al.  BUGS in Bayesian stock assessments , 1999 .

[86]  U. Dieckmann,et al.  Ecology: Managing Evolving Fish Stocks , 2007, Science.

[87]  Ian G. Taylor,et al.  A comparison of stock assessment uncertainty estimates using maximum likelihood and Bayesian methods implemented with the same model framework , 2013 .

[88]  André E. Punt,et al.  A stepwise‐selected spline approximation to time‐varying parameters, with application to occupancy modelling , 2013 .

[89]  HighWire Press Philosophical Transactions of the Royal Society of London , 1781, The London Medical Journal.

[90]  Peter Kareiva,et al.  Spatial ecology : the role of space in population dynamics and interspecific interactions , 1998 .

[91]  Päivi Elisabet Haapasaari,et al.  Integrated, age-structured, length-based stock assessment model with uncertain process variances, structural uncertainty, and environmental covariates: case of Central Baltic herring , 2013 .

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

[93]  Colin Fox,et al.  Posterior Exploration for Computationally Intensive Forward Models , 2011 .

[94]  Ray Hilborn,et al.  Standing on the Shoulders of Giants: Learning from Experience in Fisheries , 1998, Reviews in Fish Biology and Fisheries.

[95]  J. Thorson,et al.  Spatial variation buffers temporal fluctuations in early juvenile survival for an endangered Pacific salmon. , 2014, The Journal of animal ecology.

[96]  Jarrod Had MCMC Methods for Multi-Response Generalized Linear Mixed Models: The MCMCglmm R Package , 2010 .

[97]  Milner B. Schaefer Some aspects of the dynamics of populations important to the management of the commercial marine fisheries , 1991 .

[98]  N. Breslow,et al.  Approximate inference in generalized linear mixed models , 1993 .

[99]  GudmundssonGudmundur,et al.  Selection and estimation of sequential catch-at-age models , 2012 .

[100]  Karl J. Friston,et al.  Variance Components , 2003 .

[101]  Daniel Pauly,et al.  Gasping Fish and Panting Squids: Oxygen, Temperature and the Growth of Water-Breathing Animals , 2010 .

[102]  M. Mangel,et al.  Maternal age, fecundity, egg quality, and recruitment: linking stock structure to recruitment using an age-structured Ricker model , 2012 .

[103]  Murray Hill More AD of Nonlinear AMPL Models: Computing Hessian Information and Exploiting Partial Separability † , 1996 .