sensitivity_Lc: Sensitivity to Lc

In quang-huynh/MLZ: Mean Length-Based Estimators of Mortality using TMB

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)

quang-huynh/MLZ documentation built on April 10, 2022, 7:39 p.m.