surplus_production_millar: surplus_production_millar

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

surplus_production_millar

Description

Returns a list of model quantities and observations that can be used as an operating model when looking at different estimation models. A main idea here is changing the lognormal distribution to be more traditional no adjustments.

Usage

surplus_production_millar(K, r, q, P1, obs_std, pro_std, catch, seed = 123)

Arguments

K	<real> carrying capacity parameter.
r	<real> intrinsic growth parameter.
q	<real> relative index catchability.
P1	<real> initial state relative to K (1 = B_t = K)
obs_std	<real> lognormal standard deviation for relative index
pro_std	<real> lognormal standard deviation for process equation
catch	<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)
seed	<int> set.seed for reproducibility

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