DLMtool (version 5.4.1)

Rcontrol: Intrinsic rate of Increase MP

Description

An MP proposed by Carl Walters that modifies the TAC according to trends in apparent surplus production that includes information from a demographically derived prior for intrinsic rate of increase

Usage

Rcontrol(x, Data, reps = 100, plot = FALSE, yrsmth = 10, gg = 2,
  glim = c(0.5, 2))
Rcontrol2(x, Data, reps = 100, plot = FALSE, yrsmth = 10, gg = 2,
  glim = c(0.5, 2))

Arguments

x

A position in the data object

Data

A data object

reps

The number of stochastic samples of the MP recommendation(s)

plot

Logical. Show the plot?

yrsmth

The number of years for smoothing catch and biomass data

gg

A gain parameters

glim

Limits for the change in TAC among years

Value

An object of class Rec with the TAC slot populated with a numeric vector of length reps

Functions

• Rcontrol: Base version Rcontrol

• Rcontrol2: This is different from Rcontrol because it includes a quadratic approximation of recent trend in surplus production given biomass

Required Data

See '>Data for information on the Data object

Rcontrol: Abun, Cat, Dep, FMSY_M, Ind, L50, MaxAge, Mort, Year, steep, vbK, vbLinf, vbt0, wla, wlb

Rendered Equations

See Online Documentation for correctly rendered equations

Details

The TAC is calculated as: $$\text{TAC}=\text{SP}(1-gG)$$ where $g$ is a gain parameter, $\text{SP}$ is estimated surplus production, and $G$ is: For Rcontrol: $G=r(1-2D)$ where $r$ is the estimated intrinsic rate of increase, and $D$ is assumed depletion.

For Rcontrol2: $G=r-2b{B}_{\text{hist}}$ where $B_{\text{hist}}$ is the smoothed biomass overlast yrsmth years and: $$b=\sum \frac{\text{SP}}{B_{\text{hist}}}-r\frac{\sum B_{\text{hist}}}{\sum B_{\text{hist}}^2}$$.

The TAC is subject to conditions limit the maximum change from the smoothed catch over the last yrsmth years by the glim argument, e.g, default values of glim = c(0.5, 2) means that maximum decrease in TAC is 50 is 2 x average catch.

References

Made-up for this package.

See Also

Other Surplus production MPs: Fadapt, SPMSY, SPSRA, SPmod, SPslope

Examples

# NOT RUN {
Rcontrol(1, Data=DLMtool::Atlantic_mackerel, plot=TRUE)
Rcontrol2(1, Data=DLMtool::Atlantic_mackerel, plot=TRUE)
# }