Learn R Programming

MLZ (version 0.1.5)

MLCR: Mean length with catch rate mortality estimator

Description

Estimator of instantaneous total mortality (Z) from a time series of mean length data.

Usage

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)

Value

An object of class MLZ_model.

Arguments

MLZ_data

An object of class MLZ_data containing mean lengths and life history data of stock.

ncp

The number of change points in total mortality in the time series. ncp + 1 total mortality rates will be estimated.

CPUE.type

Indicates whether CPUE time series is abundance or biomass based.

loglikeCPUE

Indicates whether the log-likelihood for the CPUE will be lognormally or normally distributed.

start

An optional list of starting values. See details.

grid.search

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.

parallel

Whether grid search is performed with parallel processing. Ignored if grid.search = FALSE.

min.time

The minimum number of years between each change point for the grid search, passed to profile_MLCR. Not used if grid.search = FALSE.

Z.max

The upper boundary for Z estimates.

figure

If TRUE, a call to plot of observed and predicted mean lengths will be produced.

Details

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.

References

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.

See Also

profile_MLCR

Examples

Run this code
if (FALSE) {
data(MuttonSnapper)
MLCR(MuttonSnapper, ncp = 2, CPUE.type = "WPUE", grid.search = TRUE)
}

Run the code above in your browser using DataLab