Description Usage Arguments Details Value Required Data Optional Data Note Author(s) References See Also Examples
A surplus production model that uses only a time-series of catches and a relative abundance index
and coded in TMB. The base model, SP
, is conditioned on catch and estimates a predicted index.
Continuous surplus production and fishing is modeled with sub-annual time steps which should approximate
the behavior of ASPIC (Prager 1994). The Fox model, SP_Fox
, fixes BMSY/K = 0.37 (1/e).
The state-space version, SP_SS
estimates annual deviates in biomass. An option allows for setting a
prior for the intrinsic rate of increase.
The function for the spict
model (Pedersen and Berg, 2016) is available in DLMextra.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 | SP(
x = 1,
Data,
AddInd = "B",
rescale = "mean1",
start = NULL,
fix_dep = TRUE,
fix_n = TRUE,
LWT = NULL,
n_seas = 4L,
n_itF = 3L,
use_r_prior = FALSE,
r_reps = 100,
SR_type = c("BH", "Ricker"),
silent = TRUE,
opt_hess = FALSE,
n_restart = ifelse(opt_hess, 0, 1),
control = list(iter.max = 5000, eval.max = 10000),
...
)
SP_SS(
x = 1,
Data,
AddInd = "B",
rescale = "mean1",
start = NULL,
fix_dep = TRUE,
fix_n = TRUE,
fix_sigma = TRUE,
fix_tau = TRUE,
LWT = NULL,
early_dev = c("all", "index"),
n_seas = 4L,
n_itF = 3L,
use_r_prior = FALSE,
r_reps = 100,
SR_type = c("BH", "Ricker"),
integrate = FALSE,
silent = TRUE,
opt_hess = FALSE,
n_restart = ifelse(opt_hess, 0, 1),
control = list(iter.max = 5000, eval.max = 10000),
inner.control = list(),
...
)
SP_Fox(x = 1, Data, ...)
|
x |
An index for the objects in |
Data |
An object of class Data. |
AddInd |
A vector of integers or character strings indicating the indices to be used in the model. Integers assign the index to the corresponding index in Data@AddInd, "B" (or 0) represents total biomass in Data@Ind, "VB" represents vulnerable biomass in Data@VInd, and "SSB" represents spawning stock biomass in Data@SpInd. |
rescale |
A multiplicative factor that rescales the catch in the assessment model, which
can improve convergence. By default, |
start |
Optional list of starting values. Entries can be expressions that are evaluated in the function. See details. |
fix_dep |
Logical, whether to fix the initial depletion (ratio of biomass to carrying capacity in the
first year of the model). If |
fix_n |
Logical, whether to fix the exponent of the production function. If |
LWT |
A vector of likelihood weights for each survey. |
n_seas |
Integer, the number of seasons in the model for calculating continuous surplus production. |
n_itF |
Integer, the number of iterations to solve F conditional on the observed catch given multiple seasons within an annual time step.
Ignored if |
use_r_prior |
Logical, whether a prior for the intrinsic rate of increase will be used in the model. See details. |
r_reps |
If |
SR_type |
If |
silent |
Logical, passed to |
opt_hess |
Logical, whether the hessian function will be passed to |
n_restart |
The number of restarts (calls to |
control |
A named list of parameters regarding optimization to be passed to
|
... |
For |
fix_sigma |
Logical, whether the standard deviation of the index is fixed. If |
fix_tau |
Logical, the standard deviation of the biomass deviations is fixed. If |
early_dev |
Character string describing the years for which biomass deviations are estimated in |
integrate |
Logical, whether the likelihood of the model integrates over the likelihood of the biomass deviations (thus, treating it as a state-space variable). |
inner.control |
A named list of arguments for optimization of the random effects, which
is passed on to newton via |
To provide starting values for the SP
, a named list can be provided for FMSY
,
MSY
, dep
, and n
via the start argument (see example).
For SP_SS
, a start value can also be provided for sigma
and tau
, the standard deviation
of the index and log-biomass deviates, respectively. Default for tau is 0.1. Deviations are estimated beginning in the year when index
data are available.
If use_r_prior = TRUE
, SP
and SP_SS
will use a prior for the intrinsic rate of increase in the objective function.
A vector of length two can be passed in the start
list for the mean and standard deviation of the prior (see example). The normal
distribution is used.
If no values are provided, a prior is created using the Euler-Lotka method (Equation 15a of McAllister et al. 2001).
The Euler-Lotka method is modified to multiply the left-hand side of equation 15a by the alpha parameter of the
stock-recruit relationship (Stanley et al. 2009). Natural mortality and steepness are sampled in order to generate
a prior distribution for r. See vignette("Surplus_production")
for more details.
An object of Assessment
containing objects and output from TMB.
SP
: Cat, Ind
SP_SS
: Cat, Ind
SP_SS
: CV_Ind
The model uses the Fletcher (1978) formulation and is parameterized with FMSY and MSY as leading parameters. The default conditions assume unfished conditions in the first year of the time series and a symmetric production function (n = 2).
Tip: to create the Fox model (Fox 1970), just fix n = 1. See example.
Q. Huynh
Fletcher, R. I. 1978. On the restructuring of the Pella-Tomlinson system. Fishery Bulletin 76:515:521.
Fox, W.W. 1970. An exponential surplus-yield model for optimizing exploited fish populations. Transactions of the American Fisheries Society 99:80-88.
McAllister, M.K., Pikitch, E.K., and Babcock, E.A. 2001. Using demographic methods to construct Bayesian priors for the intrinsic rate of increase in the Schaefer model and implications for stock rebuilding. Can. J. Fish. Aquat. Sci. 58: 1871-1890.
Pedersen, M. W. and Berg, C. W. 2017. A stochastic surplus production model in continuous time. Fish and Fisheries. 18:226-243.
Pella, J. J. and Tomlinson, P. K. 1969. A generalized stock production model. Inter-Am. Trop. Tuna Comm., Bull. 13:419-496.
Prager, M. H. 1994. A suite of extensions to a nonequilibrium surplus-production model. Fishery Bulletin 92:374-389.
Stanley, R.D., M. McAllister, P. Starr and N. Olsen. 2009. Stock assessment for bocaccio (Sebastes paucispinis) in British Columbia waters. DFO Can. Sci. Advis. Sec. Res. Doc. 2009/055. xiv + 200 p.
SP_production plot.Assessment summary.Assessment retrospective profile make_MP
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 | data(swordfish)
#### Observation-error surplus production model
res <- SP(Data = swordfish)
# Provide starting values, assume B/K = 0.875 in first year of model
# and symmetrical production curve (n = 2)
start <- list(dep = 0.875, n = 2)
res <- SP(Data = swordfish, start = start)
## Not run:
plot(res)
## End(Not run)
profile(res, FMSY = seq(0.1, 0.4, 0.01))
retrospective(res)
#### State-space version
res_SS <- SP_SS(Data = swordfish, start = list(dep = 0.875, sigma = 0.1, tau = 0.1))
## Not run:
plot(res_SS)
## End(Not run)
#### Fox model
res_Fox <- SP(Data = swordfish, start = list(n = 1), fix_n = TRUE)
res_Fox2 <- SP_Fox(Data = swordfish)
#### SP with r_prior
res_prior <- SP(Data = SimulatedData, use_r_prior = TRUE)
#### Pass an r_prior to the model with mean = 0.35, sd = 0.10
res_prior2 <- SP(Data = SimulatedData, use_r_prior = TRUE, start = list(r_prior = c(0.35, 0.10)))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.