Estimate_mortality_dists {SharkDemography}

Estimate the range of mortality estimates produced by the different estimators

Description

This function performs a Monte Carlo simulation analysis which determines the distributions of natural mortality (M) produced by each estimator.

Usage

Estimate_mortality_dists(n = 1000, data, M.estimators = NULL)

Arguments

n

The number of simulations to be run

data

A multi-level list of the class 'Demography.inputs' produced from the 'create_data_input' function and then manually completed.

M.estimators

Any specific natural mortality estimators to be included in the analysis. Must be a single estimator or a vector of estimators. These can include: "Pet.Wro", "Jensen.mat","Chen.Yuan","Then_hoenig","Then_pauly", "Jensen.mat","Charnov" or "Chen.Want". If none are specified then all applicable estimators could be chosen.

Value

A list which includes a data.frame of the distributions age invariant estimators summarised as the mean and 95% quantiles and dataframes of any age dependent estimators with the mean and 95% quantiles for each age class

Examples

# load Silky shark data produced by create_data_input()
# Type `?create_data_input()` for details
data("Silky_data")

# Run function to get natural mortality distributions from  Monte Carlo
# simulations for all available natural mortality estimators.
# Set n = at least 1000 for full analysis but use n = 100 for testing

Estimate_mortality_dists(n = 100, Silky_data)