sensitivity_Lc {MLZ}

Sensitivity to Lc

Description

The function re-calculates mean lengths for each alternative value of Lc and re-estimates mortality. Currently supports only the ML estimator.

Usage

sensitivity_Lc(MLZ_data, MLZ_model, Lc.vec, grid.search = FALSE,
  figure = TRUE)

Arguments

MLZ_data

An object of class MLZ_data containing mean lengths and life history data of stock. Must contain length composition data.

MLZ_model

An object of class MLZ_model with base value of Lc.

Lc.vec

A vector of alternative Lc values.

grid.search

Whether a grid search is performed or not. By default, the starting values in the sensitivity analysis are the estimates from object MLZ_model.

figure

Whether a figure will be produced, similar to Figure 6 of Huynh et al. (2017).

Value

A matrix of mortality and change point estimates with each value Lc.

References

Huynh, Q.C, Gedamke, T., Hoenig, J.M, and Porch C. 2017. Multispecies Extensions to a Nonequilibrium Length-Based Mortality Estimator. Marine and Coastal Fisheries 9:68-78.

See Also

ML

Examples

## Not run: 
data(SilkSnapper)
new.dataset <- new("MLZ_data", Year = 1983:2013, Len_df = SilkSnapper, length.units = "mm",
vbLinf = 794, vbK = 0.1)

new.dataset@Lc <- 310
new.dataset <- calc_ML(new.dataset)

first.MLZmodel <- ML(new.dataset, 1)
Lc.vec <- seq(240, 340, 5)

sensitivity_Lc(new.dataset, first.MLZmodel, Lc.vec)

## End(Not run)