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.
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.
# NOT RUN {
# See examples in package 'plotMCMC'
# }
Run the code above in your browser using DataLab