Comparing Bayesian and frequentist methods of fisheries models: Hierarchical catch curves

Bayesian inference is an emerging statistical paradigm and is becoming an increasingly used alternative to frequentist inference. Unfortunately, little is known about the efficacy of Bayesian inference and how it relates to the historical methodology of evaluating fisheries related models. Mortality information is routinely used in fisheries management to describe fish population abundance over time and has been historically estimated using catch curves and frequentist inference (i.e., maximum likelihood estimation). The objective of this study was to compare frequentist and Bayesian inference approaches to estimate instantaneous mortality (Z) from a hierarchical catch curve model. The data used in the comparison were from a long term monitoring program of yellow perch Perca flavescens from southern Lake Michigan in addition to a simulated dataset where parameter estimates were compared to known values. Point estimates of Z were similar among both methods. Similarly, Bayesian inference 95% credible intervals were concordant with frequentist 95% confidence intervals. However, the root mean squared error of frequentist inference increased at a higher rate than Bayesian inference with increasing variability in the simulated dataset. Our study builds on the literature that seeks to compare results between these two paradigms to assist managers to make the best decision possible when deciding what statistical paradigm to employ.

[1]  J. Kruschke Bayesian estimation supersedes the t test. , 2013, Journal of experimental psychology. General.

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

[3]  Bradley P. Carlin,et al.  Bayesian Methods for Data Analysis , 2008 .

[4]  James R. Bence,et al.  Detecting Temporal Trends in Freshwater Fisheries Surveys: Statistical Power and the Important Linkages between Management Questions and Monitoring Objectives , 2013 .

[5]  E. Rivot,et al.  Hierarchical Bayesian analysis of capture-mark-recapture data , 2002 .

[6]  Andrew Gelman,et al.  General methods for monitoring convergence of iterative simulations , 1998 .

[7]  J. Berger The case for objective Bayesian analysis , 2006 .

[8]  Zhenming Su,et al.  Performance of a Bayesian state-space model of semelparous species for stock-recruitment data subject to measurement error , 2012 .

[9]  I. Potter,et al.  A Bayesian approach for overcoming inconsistencies in mortality estimates using, as an example, data for Acanthopagrus latus , 2004 .

[10]  R. O’Hara,et al.  Do not log‐transform count data , 2010 .

[11]  S. Jacquemin,et al.  Long-term fish assemblages respond to habitat and niche breadth in the West Fork White River, Indiana , 2013 .

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

[13]  D. G. Chapman,et al.  Catch Curves and Mortality Rates , 1961 .

[14]  T. Lauer,et al.  Abiotic and biotic correlates of yellow perch recruitment to age-2 in southern Lake Michigan, 1984–2007 , 2012 .

[15]  Yong Liu,et al.  Predicting the Hypoxic‐Volume in Chesapeake Bay with the Streeter–Phelps Model: A Bayesian Approach 1 , 2011 .

[16]  H. Buckley,et al.  Local to Continental-Scale Variation in the Richness and Composition of an Aquatic Food Web , 2010 .

[17]  Y. Benjamini Discovering the false discovery rate , 2010 .

[18]  F. Lutzoni,et al.  Bayes or bootstrap? A simulation study comparing the performance of Bayesian Markov chain Monte Carlo sampling and bootstrapping in assessing phylogenetic confidence. , 2003, Molecular biology and evolution.

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

[20]  Mark E. Borsuk,et al.  A Bayesian network of eutrophication models for synthesis, prediction, and uncertainty analysis , 2004 .

[21]  C. P. Madenjian,et al.  Exploring Mechanisms Underlying Sex-Specific Differences in Mortality of Lake Michigan Bloaters , 2012 .

[22]  T. Lauer,et al.  Bayesian Estimation of Age and Length at 50% Maturity , 2012 .

[23]  Daniel Stegmueller,et al.  How Many Countries for Multilevel Modeling? A Comparison of Frequentist and Bayesian Approaches , 2013 .

[24]  Angelo J. Canty,et al.  Comparison of Bayesian and frequentist approaches in modelling risk of preterm birth near the Sydney Tar Ponds, Nova Scotia, Canada , 2007, BMC medical research methodology.

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

[26]  Aaron M. Ellison,et al.  AN INTRODUCTION TO BAYESIAN INFERENCE FOR ECOLOGICAL RESEARCH AND ENVIRONMENTAL , 1996 .

[27]  Bradley P Carlin,et al.  Ecological boundary detection using Bayesian areal wombling. , 2010, Ecology.

[28]  G. Vivó-Truyols,et al.  Bayesian approach for peak detection in two-dimensional chromatography. , 2012, Analytical chemistry.

[29]  D. Robson THE ANALYSIS OF A CATCH CURVE , 1960 .

[30]  T. Lauer,et al.  Year Three Progress Report: Dynamics and Models of the Yellow Perch in Indiana Waters of Lake Michigan and Near-Shore Fish Community Characteristics , 2001 .

[31]  J. Thorson,et al.  Better Catch Curves: Incorporating Age-Specific Natural Mortality and Logistic Selectivity , 2011 .

[32]  George B. Arhonditsis,et al.  A revaluation of lake-phosphorus loading models using a Bayesian hierarchical framework , 2009, Ecological Research.

[33]  S. Jacquemin,et al.  Body Size and Geographic Range Do Not Explain Long Term Variation in Fish Populations: A Bayesian Phylogenetic Approach to Testing Assembly Processes in Stream Fish Assemblages , 2014, PloS one.

[34]  N. Stenseth,et al.  Life-history variation among local populations of Atlantic cod from the Norwegian Skagerrak coast , 2004 .

[35]  Paul G. Ambrose,et al.  Frequentist and Bayesian Pharmacometric-Based Approaches To Facilitate Critically Needed New Antibiotic Development: Overcoming Lies, Damn Lies, and Statistics , 2011, Antimicrobial Agents and Chemotherapy.

[36]  Etienne Rivot,et al.  Estimation of gillnet efficiency and selectivity across multiple sampling units: A hierarchical Bayesian analysis using mark-recapture data , 2007 .

[37]  Kenneth H. Reckhow,et al.  Bayesian inference in non-replicated ecological studies , 1990 .

[38]  J. Kruschke Doing Bayesian Data Analysis: A Tutorial with R and BUGS , 2010 .

[39]  J. Kruschke What to believe: Bayesian methods for data analysis , 2010, Trends in Cognitive Sciences.

[40]  T. S. McComish,et al.  Relationship between Alewife Abundance and Yellow Perch Recruitment in Southern Lake Michigan , 2000 .

[41]  W. E. Pine,et al.  FISH RECRUITMENT IS INFLUENCED BY RIVER FLOWS AND FLOODPLAIN INUNDATION AT APALACHICOLA RIVER, FLORIDA , 2013 .

[42]  G. Hilley,et al.  Deducing Paleoearthquake Timing and Recurrence from Paleoseismic Data, Part I: Evaluation of New Bayesian Markov-Chain Monte Carlo Simulation Methods Applied to Excavations with Continuous Peat Growth , 2008 .

[43]  M. Cappo,et al.  Age, growth, mortality rates and corresponding yield estimates using otoliths of the tropical red snappers, Lutjanus erythropterus, L. malabaricus and L. sebae, from the central Great Barrier Reef $ , 2000 .

[44]  W. Chaplin,et al.  A comparison of frequentist and Bayesian inference: searching for low-frequency p modes and g modes in Sun-as-a-star data , 2010, 1004.4505.

[45]  Anders Nielsen,et al.  Comparison of the frequentist properties of Bayes and the maximum likelihood estimators in an age-structured fish stock assessment model , 2002 .

[46]  M. McCarthy,et al.  The influence of abundance on detectability , 2013 .

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

[48]  E. Baker,et al.  Precision of Ages Determined from Scales and Opercles for Yellow Perch Perca flavescens , 1998 .

[49]  M. Lee Three case studies in the Bayesian analysis of cognitive models , 2008, Psychonomic bulletin & review.

[50]  C D Naylor,et al.  A comparison of a Bayesian vs. a frequentist method for profiling hospital performance. , 2001, Journal of evaluation in clinical practice.

[51]  J. Jackson,et al.  Relationships between Annual Variations in Reservoir Conditions and Age-0 Largemouth Bass Year-Class Strength , 2000 .

[52]  Marc Kery,et al.  Introduction to WinBUGS for Ecologists: Bayesian approach to regression, ANOVA, mixed models and related analyses , 2010 .

[53]  B. Mesnil,et al.  Estimating stock parameters from trawl cpue-at-age series using year-class curves , 2007 .