Estimator of instantaneous total mortality (Z) from a time series of mean length data.
MLCR(MLZ_data, ncp, CPUE.type = c(NA, "WPUE", "NPUE"),
loglikeCPUE = c("lognormal", "normal"), start = NULL,
grid.search = TRUE, parallel = ifelse(ncp > 2, TRUE, FALSE),
min.time = 3, Z.max = 5, figure = TRUE)
An object of class MLZ_model
.
An object of class MLZ_data
containing mean lengths and
life history data of stock.
The number of change points in total mortality in the time series. ncp + 1
total
mortality rates will be estimated.
Indicates whether CPUE time series is abundance or biomass based.
Indicates whether the log-likelihood for the CPUE will be lognormally or normally distributed.
An optional list of starting values. See details.
If TRUE
, a grid search will be performed using the profile_MLCR
function to find the best starting values for the change points (the years when mortality changes).
Ignored if ncp = 0
or if start
is provided.
Whether grid search is performed with parallel processing. Ignored if grid.search = FALSE
.
The minimum number of years between each change point for the grid search, passed
to profile_MLCR
. Not used if grid.search = FALSE
.
The upper boundary for Z estimates.
If TRUE
, a call to plot
of observed and predicted mean lengths will be produced.
For a model with I
change points, the starting values in
start
is a list with the following entries:
Z a vector of length = I+1
.
yearZ a vector of length = I
.
start
can be NULL
, in which case, the supplied starting values depend on
the value of grid.search
. If grid.search = TRUE
, starting values will use the
values for yearZ
which minimize the negative log-likelihood from the grid search.
Otherwise, the starting values for yearZ
evenly divide the time series.
Huynh, Q.C., Gedamke, T., Porch, C.E., Hoenig, J.M., Walter, J.F, Bryan, M., and Brodziak, J. In revision. Estimating Total Mortality Rates from Mean Lengths and Catch Rates in Non-equilibrium Situations. Transactions of the American Fisheries Society.
profile_MLCR
if (FALSE) {
data(MuttonSnapper)
MLCR(MuttonSnapper, ncp = 2, CPUE.type = "WPUE", grid.search = TRUE)
}
Run the code above in your browser using DataLab