Description Usage Arguments Details Value Functions Required Data Rendered Equations Author(s) References See Also Examples
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
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 |
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 |
The TAC is calculated as:
\textrm{TAC} = \textrm{SP} (1-gG)
where g is a gain parameter, \textrm{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 - 2bB_\textrm{hist} where B_\textrm{hist}
is the smoothed biomass overlast yrsmth
years and:
b = ∑{\frac{\textrm{SP}}{B_\textrm{hist}} - r} \frac{∑{B_\textrm{hist}}}{∑{B_\textrm{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% of average catch and maximum increase
is 2 x average catch.
An object of class Rec
with the TAC
slot populated with a numeric vector of length reps
Rcontrol
: Base version Rcontrol
Rcontrol2
: This is different from
Rcontrol
because it includes a quadratic approximation of recent trend in
surplus production given biomass
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
See Online Documentation for correctly rendered equations
C. Walters and T. Carruthers
Made-up for this package.
Other Surplus production MPs:
Fadapt()
,
SPMSY()
,
SPSRA()
,
SPmod()
,
SPslope()
1 2 | Rcontrol(1, Data=DLMtool::Atlantic_mackerel, plot=TRUE)
Rcontrol2(1, Data=DLMtool::Atlantic_mackerel, plot=TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.