sensitivity_Lc | R Documentation |
The function re-calculates mean lengths for each alternative value of Lc and re-estimates mortality. Currently supports only the ML estimator.
sensitivity_Lc(MLZ_data, MLZ_model, Lc.vec, grid.search = FALSE, figure = TRUE)
MLZ_data |
An object of class |
MLZ_model |
An object of class |
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 |
figure |
Whether a figure will be produced, similar to Figure 6 of Huynh et al. (2017). |
A matrix of mortality and change point estimates with each value Lc.
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.
ML
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.