Fadapt.Rd
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)
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. |
An object of class Rec
with the TAC
slot populated with a numeric vector of length reps
A numeric vector of quota recommendations
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=FyBy−1
where By−1 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 Flim2−Flim1,
Fmod2 is
Fmod1+g−G
where g is gain parameter gg
, G
is the predicted surplus production given current abundance,
and:
Fmod1={−2if Fold<Flim12if Fold>Flim2logFfrac1−Ffracif Flim1≤Fold≤Flim2
where −Ffrac=Fold−Flim1Flim3,
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.
See Data for information on the Data
object
Fadapt
: Abun, Cat, FMSY_M, Ind, Mort, Year
See Online Documentation for correctly rendered equations
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
Other Fmsy/M methods:
DynF()
,
Fratio()
Other Surplus production MPs:
Rcontrol()
,
SPMSY()
,
SPSRA()
,
SPmod()
,
SPslope()
#> TAC (median) #> 8.117012