DLMtool (version 5.2)

SPslope: Slope in surplus production MP

Description

A management procedure that makes incremental adjustments to TAC recommendations based on the apparent trend in recent surplus production. Based on the theory of Mark Maunder (IATTC)

Usage

SPslope(x, Data, reps = 100, plot = FALSE, yrsmth = 4, alp = c(0.9,
  1.1), bet = c(1.5, 0.9))

Arguments

x

A position in a data-limited methods data object

Data

A data-limited methods data object

reps

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

plot

Logical. Show the plot?

yrsmth

Years over which to smooth recent estimates of surplus production

alp

Condition for modifying the Data (bounds on change in abundance)

bet

Limits for how much the Data can change among years

Value

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

Required Data

See '>Data for information on the Data object

SPslope: Abun, Cat, Ind, Year

Rendered Equations

See Online Documentation for correctly rendered equations

Details

Note that this isn't exactly what Mark has previously suggested and is stochastic in this implementation.

The TAC is calculated as: $${\text{TAC}}_y=\{\begin{array}{ll}M\overline{C}& \text{if }r<{\alpha }_1\\ \overline{C}& \text{if }{\alpha }_1<r<{\alpha }_2\\ {\text{bet}}_2\text{SP}& \text{if }r>{\alpha }_2\end{array}$$ where $r$ is the ratio of predicted biomass in next year to biomass in current year $\overline{C}$ is the mean catch over the last yrmsth years, ${\alpha }_1$ and ${\alpha }_2$ are specified in alp, ${\text{bet}}_1$ and ${\text{bet}}_2$ are specified in bet, $\text{SP}$ is estimated surplus production in most recent year, and: $$M=1-{\text{bet}}_1\frac{B_y-{\tilde{B}}_y}{B_y}$$ where $B_y$ is the most recent estimate of biomass and $\tilde{B}$ is the predicted biomass in the next year.

References

http://www.iattc.org/Meetings/Meetings2014/MAYSAC/PDFs/SAC-05-10b-Management-Strategy-Evaluation.pdf

See Also

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

Examples

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