Markov chain Monte Carlo results from stock assessment of cod (Gadus morhua) in Icelandic waters.
xmcmc
List containing four data frames:
L | likelihood components: CAc (commercial catch at
age), CAs (survey catch at age), Survey (survey
abundance index), Prior (penalty for priors), Total
(sum of all likelihood components; the objective function). |
P | estimated parameters: R0 (average virgin
recruitment), Rinit (initial recruitment scaler),
uinit (initial harvest rate), cSleft (left-side slope
of commercial selectivity curve), cSfull (age at full
commercial selectivity), sSleft (left-side slope of survey
selectivity curve), sSfull (age at full survey selectivity),
logq (log-transformed survey catchability). |
B | predicted biomass (age 4+) by year. |
R | predicted recruitment by year. |
See the x.cod
help page for details about the data and
model.
Magnusson, A., Punt, A.E. and Hilborn, R. (2013). Measuring uncertainty in fisheries stock assessment: the delta method, bootstrap, and MCMC. Fish and Fisheries, 14, 325--342. tools:::Rd_expr_doi("10.1111/j.1467-2979.2012.00473.x")
importMCMC
was used to import the MCMC results.
xmcmc
and xproj
are MCMC results and projections
for the x.cod
model run.
scape-package
gives an overview of the package.
# See examples in package 'plotMCMC'
Run the code above in your browser using DataLab