Accounting for uncertainty due to data processing in virtual population analysis using Bayesian multiple imputation

Virtual population analysis (VPA) is used in many stock assessment settings and requires a total catch-at-age data set where an age is assigned to each fish that has been caught. These data sets are typically constructed using ad hoc methods that rely on numerous assumptions. Although approaches are available to account for observation error in these data, no statistically rigorous methods have been developed to account for uncertainty from data processing. To address this, we investigated a Bayesian multiple imputation approach to filling missing size data. Using Atlantic yellowfin tuna (Thunnus albacares) and bigeye tuna (Thunnus obesus) as case studies, we evaluated the hypothesis that data processing is as important in determining management reference points in stock assessments as conventional sources of uncertainty. Size imputation models accounting for location, season, and year provided good predictive capacity. Uncertainty from data processing could be large; however, the circumstances for this w...

[1]  André E. Punt,et al.  The performance of VPA-based management , 1997 .

[2]  John Sibert,et al.  AD Model Builder: using automatic differentiation for statistical inference of highly parameterized complex nonlinear models , 2012, Optim. Methods Softw..

[3]  J. Rice,et al.  Coexistence in North Sea fish communities: Implications for growth and natural mortality , 2008 .

[4]  Laurence T. Kell,et al.  The spectre of uncertainty in management of exploited fish stocks: The illustrative case of Atlantic bluefin tuna , 2014 .

[5]  Murdoch K. McAllister,et al.  Atlantic Bluefin Tuna: A Novel Multistock Spatial Model for Assessing Population Biomass , 2011, PloS one.

[6]  Roderick J A Little,et al.  A Review of Hot Deck Imputation for Survey Non‐response , 2010, International statistical review = Revue internationale de statistique.

[7]  O. Maury,et al.  GROWTH OF BIGEYE TUNA (THUNNUS OBESUS) IN THE EASTERN ATLANTIC OCEAN FROM TAGGING-RECAPTURE DATA AND OTOLITH READINGS , 2005 .

[8]  T. Matsuishi,et al.  Effects of sampling errors on abundance estimates from virtual population analysis for walleye pollock in northern waters of Sea of Japan , 2007, Fisheries Science.

[9]  J. Ianelli,et al.  A scientific alternative to moratoria for rebuilding depleted international tuna stocks , 2016 .

[10]  D. Rubin,et al.  Statistical Analysis with Missing Data , 1988 .

[11]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[12]  André E. Punt,et al.  Experiences in the evaluation and implementation of management procedures , 1999 .

[13]  Emmanuel Chassot,et al.  Global spatio-temporal patterns in tropical tuna purse seine fisheries on drifting fish aggregating devices (DFADs): Taking a historical perspective to inform current challenges , 2013 .

[14]  M. Prager.  Comparison of logistic and generalized surplus-production models applied to swordfish, Xiphias gladius, in the north Atlantic Ocean , 2002 .

[15]  ASSESSING COMPONENTS OF UNCERTAINTY IN VPA ABUNDANCE AND MORTALITY ESTIMATES USING AN ALTERNATIVE EXPLOITATION RATE‐BASED ALGORITHM , 2012 .

[16]  André E. Punt,et al.  Management strategy evaluation: best practices , 2016 .

[17]  Daniel K. Kimura,et al.  Statistical Assessment of the Age–Length Key , 1977 .

[18]  A. Dobson,et al.  Multiple imputation for body mass index: lessons from the Australian Longitudinal Study on Women's Health , 2004, Statistics in medicine.

[19]  C. Davies,et al.  Spatial and Sex-Specific Variation in Growth of Albacore Tuna (Thunnus alalunga) across the South Pacific Ocean , 2012, PloS one.

[20]  S. Richardson,et al.  Strategy for Modelling Nonrandom Missing Data Mechanisms in Observational Studies Using Bayesian Methods , 2012 .

[21]  David A. Fournier,et al.  MULTIFAN a Likelihood-Based Method for Estimating Growth Parameters and Age Composition from Multiple Length Frequency Data Sets Illustrated using Data for Southern Bluefin Tuna (Thunnus maccoyii) , 1990 .

[22]  John Shepherd,et al.  An Alternative Perspective on Recruitment Overfishing and Biological Reference Points , 1987 .

[23]  J. Fromentin,et al.  An evaluation of management strategies for Atlantic tuna stocks , 2003 .

[24]  Cody S. Szuwalski,et al.  Examining common assumptions about recruitment: a meta‐analysis of recruitment dynamics for worldwide marine fisheries , 2015 .

[25]  Douglas M. Bates,et al.  Unconstrained parametrizations for variance-covariance matrices , 1996, Stat. Comput..

[26]  D. Rubin,et al.  Multiple Imputation for Nonresponse in Surveys , 1989 .

[27]  Richard D. Methot,et al.  Stock synthesis: A biological and statistical framework for fish stock assessment and fishery management , 2013 .

[28]  D. Rubin.  INFERENCE AND MISSING DATA , 1975 .

[29]  A. Fonteneau,et al.  Modélisation d'une croissance en deux stances chez l'albacore (Thunnus albacares) de l'Atlantique Est , 1992 .

[30]  J. Rooker,et al.  Origin of Atlantic bluefin tuna (Thunnus thynnus) in the Bay of Biscay , 2015 .

[31]  F. Abascal,et al.  Stereological assessment of the reproductive status of female Atlantic northern bluefin tuna during migration to Mediterranean spawning grounds through the Strait of Gibraltar , 2002 .

[32]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[33]  Thomas R. Carruthers,et al.  Imputing recreational angling effort from time-lapse cameras using an hierarchical Bayesian model , 2015 .

[34]  John Shepherd,et al.  Extended survivors analysis: An improved method for the analysis of catch-at-age data and abundance indices , 1999 .

[35]  T. Pitcher,et al.  Estimating the Worldwide Extent of Illegal Fishing , 2009, PloS one.