This is a single-sex model with 10 age classes, the catch data starting in 1971 and ending in 2003. The model was fitted to three data components: survey abundance index, commercial catch at age, and survey catch at age.
x.cod
scape
containing:
N |
predicted numbers at age |
B |
predicted biomass, recruitment, and observed landings (year things) |
Sel |
predicted selectivity and observed maturity (age things) |
Dev |
predicted recruitment deviates from the stock-recruitment curve |
Survey |
survey abundance index and fit |
CAc |
commercial C@A (catch at age) and fit |
CAs |
survey C@A and fit |
A maturity vector of zeros and ones was used to predict the biomass of age 4 and older, the quantity of main interest for the management of this stock.
Estimated parameters: R0, Rinit, uinit, Sleft[commercial], Sfull[c], Sleft[survey], Sfull[s], q, and 41 recruitment deviates.
Magnusson, A. and Hilborn, R. (2007) What makes fisheries data informative? Fish and Fisheries 8, 337--358.
importCol
was used to import the fitted model. x.cod
, x.ling
, x.oreo
, and
x.sbw
are fitted scape
models to explore.
xmcmc
and xproj
are MCMC results and
projections for the x.cod
model run.
scape-package
gives an overview of the package.
plotB(x.cod)
plotCA(x.cod, "c")
plotCA(x.cod, "s")
plotIndex(x.cod, "s")
plotN(x.cod)
plotSel(x.cod)
Run the code above in your browser using DataLab