surplus_production_model_sigma: surplus_production_model_sigma

In Craig44/stockassessmenthelper: Package to help plot stock assessment models and functionality

surplus_production_model_sigma

Description

simulates a surplus production model can be used as an operating model, Different from the above uses a standard deviation for process error instead of CV

Usage

surplus_production_model_sigma(
  seed = 123,
  n_t,
  B1,
  r,
  k,
  catches,
  obs_std,
  catchability,
  observation_likelihood = 1,
  process_error = NULL,
  process_std = NULL
)

Arguments

seed	integer random number seed
n_t	<integer> number of time-steps to run the model
B1	<real> starting exploited biomass
r	<real> intrinsic growth parameter.
k	<real> carrying capacity parameter.
catches	<vector> observed catches, defines number of years, so if you have a period of no catches set values of 0, number of years = length(catches+ 1)
obs_std	<real> Standard deviation which is used to generate simulated observations.
catchability	<vector> a scalar to generate the relative index observation
observation_likelihood	<int> likelihood type, 1 = lognormal, 2 = normal, 3 = ...
process_error	(optional) time specific time deviations (epsilon_t,p) for each time step
process_std	(optional) standard deviation for process deviations corrections

Value

Returns a list of model quantities and observations that can be used as an operating model

Craig44/stockassessmenthelper documentation built on April 14, 2023, 10:57 a.m.