run.Proj: Schaefer's projection function

In habeebollah/montiR: Biomass dynamic models for fisheries using time series data fitting methods

Schaefer's projection function

Description

Function to create projection in Schaefer's model based on different reference points (MSY, Emsy) and different Total Allowable Catch setting (TAC). This projection can be used to assist in the development of limit and target reference point.

Usage

run.Proj(inpars, df, nyears, nsims, stoch, TAC = 1, plot = TRUE)

Arguments

inpars	fitted surplus production parameters which consist of K (carrying capacity), B0 (biomass when fishing is started), r (intrinsic growth rate), q (catchability coefficient), and sigma (observation error)
df	dataframe containing three columns; year, catch and unit of effort
nyears	number of years the projection for the fishery
nsims	number of iteration performed
stoch	stochastic value as a sigma
TAC	number to drive the management level
plot	whether the projection plot is produced or not

Value

The TAC is set at 1 in default. It would depend on the fishery to set whether the TAC can be set lower to be more conservative or set higher to increase more catch. For instance, when the TAC is set conservative at 0.8 of reference point, it will return 0.8MSY, and 0.8Emsy. In contrary, the TAC will increase catch when it is set at 1.2 (or other value higher than 1) and return 1.2MSY, and 1.2Emsy.

Examples

fit <- calc.MSY(K=1000, B0=1000, r=0.2, q=0.00025, s.sigma=0.1,
df=df.goodcontrast, OWT=FALSE, currentF = 0.7, weight=0.5, plot=TRUE)

run.Proj(inpars=fit[[1]][,2], df=df.goodcontrast, nyears=30, nsims=100, stoch=0.2, TAC=1, plot=TRUE)

habeebollah/montiR documentation built on Dec. 11, 2022, 7:55 p.m.