Processing math: 100%

An adaptive MP that uses trajectory in inferred suplus production and fishing mortality rate to update a TAC

Fadapt(x, Data, reps = 100, plot = FALSE, yrsmth = 7, gg = 1)

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

Years over which to smooth recent estimates of surplus production

gg

A gain parameter controlling the speed in update in TAC.

Value

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

A numeric vector of quota recommendations

Details

Fishing rate is modified each year according to the gradient of surplus production with biomass (aims for zero). F is bounded by FMSY/2 and 2FMSY and walks in the logit space according to dSP/dB. This is derived from the theory of Maunder 2014.

The TAC is calculated as: TACy=FyBy1 where By1 is the most recent biomass, estimated with a loess smoother of the most recent yrsmth years from the index of abundance (Data@Ind) and estimate of current abundance (Data@Abun), and

Fy=Flim1+(expFmod21+expFmod2Flim3)

where Flim1=0.5FMSYMM, Flim2=2FMSYMM, Flim3 is Flim2Flim1, Fmod2 is Fmod1+gG where g is gain parameter gg, G is the predicted surplus production given current abundance, and: Fmod1={2if Fold<Flim12if Fold>Flim2logFfrac1Ffracif Flim1FoldFlim2 where Ffrac=FoldFlim1Flim3, Fold=ChistBhist/n where Chist and Bhist are smooth catch and biomass over last yrsmth, and n is yrsmth.

Tested in Carruthers et al. 2015.

Required Data

See Data for information on the Data object

Fadapt: Abun, Cat, FMSY_M, Ind, Mort, Year

Rendered Equations

See Online Documentation for correctly rendered equations

References

Carruthers et al. 2015. Performance evaluation of simple management procedures. ICES J. Mar Sci. 73, 464-482.

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

See also

Other Fmsy/M methods: DynF(), Fratio()

Other Surplus production MPs: Rcontrol(), SPMSY(), SPSRA(), SPmod(), SPslope()

Examples

Fadapt(1, Data=DLMtool::Atlantic_mackerel, plot=TRUE)
#> TAC (median) #> 8.117012