schnute_orig: Prepare list for Schnute Original Process Error Model...

In lbatts/sbar: Stage-based Assessments in R

Description

Create an object with TMB framework, including data, gradients and NLL function for a Schnute Original Process Error Model assessment that can be optimised.

Usage

schnute_orig(
  version = 2,
  catch_b,
  indices_b,
  ts,
  mwts,
  tsp = 0,
  mu = 0.5,
  rho,
  W,
  ind_l_wt = 1,
  start_q = 1e-08,
  start_indexsigma = 0.1,
  start_sigma = exp(-0.2),
  start_rec_a,
  start_rec_b,
  spawn_prop = 1,
  fix_sigma = TRUE,
  fix_indexsigma = FALSE,
  adrep = FALSE
)

Arguments

version	numeric, either 1, 2 or 3. This controls what deterministic equations in the model are used to derive population biomass. 1 and 2 use the fraction of of total biomass in a given year due to newly recruited fish. This fraction is derived from mean weights and detailed in the schnute vignette. version = 3 is the more classical population dynamics. 1 whole biomass derived from recruit biomass 2 whole biomass derived from previously exploited biomass 3 whole biomass is a combination of recruit biomass and previously exploited biomass
catch_b	numeric vector of catch biomass over time period of assessment
indices_b	matrix of biomass surveys (CPUE) of dimensions: no. of surveys x no.years
ts	numeric. Survey timing parameters
mwts	matrix of mean weights from sampling with dimensions: 3 x no. years. recruit mean weights \bar{Y} (first row), previously exploited biomass mean weights \bar{Z} (second row) and entire assessed biomass mean weight \bar{X} (third row).
tsp	numeric. Timing of spawning. Default to 0 (start of year).
mu	numeric. proportion of catch taken before natural mortality.
rho	numeric. Growth parameter, slope of linear growth model.
W	numeric. Growth parameter, intercept of linear growth model.
ind_l_wt	numeric. Survey weighting in the likelihood. Defaults to 1 fro each survey, ie.e. equal weighting
start_q	Starting values for survey catchability parameters. Default is 1e-6
start_indexsigma	Starting values for survey sd parameters. Default is 0.1
start_sigma	Starting parameter value fraction of population that survives natural moratlity. Default is e^{0.2}
start_rec_a	Starting parameter value for the 'a' parameter of the Beverton-Holt stock-recruit function. The asymptotic biomass of recruits. Default is 1/5*max(catch_b).
start_rec_b	Starting parameter value for the 'b' parameter of the Beverton-Holt stock-recruit function. The spawning stock biomass needed to produce a/2 on average. Default is 4*max(catch_b).
spawn_prop	proportion of biomass that is mature. Defaults to 1 for each year.
fix_sigma	logical. logical. Should survival be fixed in the model
fix_indexsigma	logical. Should survey standard deviation be fixed in the model
adrep	logical. Whether the user would like the ADreport variables (and their derivatives) reported for starting parameters.

Details

schnute_orig is simply a wrapper function that gives the output from MakeADFun from TMB, i.e. an objective function with derivatives, hessian etc. Otimisation and extraction of values from the assessment models has been kept separate to allow flexibility with optimisation methods as well as easy access to MakeADFun outputs such as the hessian, gradients etc. See http://kaskr.github.io/adcomp/_book/Introduction.html and TMB documentation for details. Users should refer to the sbar vignette by running vignette("intro_to_sbar", "sbar") for for details on the schnute_orig function.

The table below gives the outputs and description of the values that can be extracted from this assessment model after optimisation and summary(TMB::sdreport(x)) of the schnute_orig object.

#'

Output	Description
logitq	logit transformed survey catchability
logitsigma	logit transformed survival of natural mortality
logindex_sigma	log transformed survey standard deviation
lnb	log transformed total biomass
lnpr	log transformed previously-exploited biomass
lnr	log transformed recruit biomass
lnN	log transformed total numbers
lnPR	log transformed previously-exploited or post-recruit numbers
lnR	log transformed recruit numbers
lnC	log transformed observed catch numbers
biomass	total biomass
N	total numbers
ssb	spawning biomass
post_rec	previously-exploited or post-recruit biomass
PR	previously-exploited or post-recruit numbers
rec_bio	recruit biomass
rec_no	recruit numbers
C	observed catch numbers
logpred_survey	log transformed predicted survey indices
omega	fraction of total biomass in a given year due to newly recruited fish
index_sigma	Survey indices standard deviation
sigma	survival of natural mortality
logrec_param	log transformed recruit parameters from a Beverton-Holt SR function if estimated in versions 1 or 3
rec_param	the two recruit parameters fro a Beverton-Holt SR functionif estimated in versions 1 or 3
qhat	survey catchability

Value

List with components for optimiser in R. This output is that of the function MakeADFun from TMB

Access (transformed) starting values of parameters to be estimated with x$par (where x is the schnute_orig object) to see what parameters are to be estimated.

Examples

## Not run:  obj <- schnute_orig(catch_b = catch_biomass, indices_b = obs_fill,
ts = timing, mwts = mwts, rho = rho1, W = W1, start_sigma = sigma)
## End(Not run)

lbatts/sbar documentation built on Feb. 25, 2022, 9:02 a.m.