scape (version 2.3-2)

getSigmaI: Extract Abundance Index Sigma

Description

Extract the sigma (magnitude of observation noise) that was used in a model, from survey or commercial abundance index data.

Usage

getSigmaI(model, what="s", series=NULL, digits=NULL)

Arguments

model

fitted scape model containing element CPUE and/or Survey.

what

which sigma to extract: "c"[ommercial] or "s"[urvey] abundance index.

series

vector of strings indicating which gears or surveys to analyze (all by default).

digits

number of decimal places to use when rounding, or NULL to suppress rounding.

Value

Numeric vector of year-specific sigmas, or a list of such vectors when analyzing multiple series.

See Also

getN, getSigmaI, getSigmaR, estN, estSigmaI, and estSigmaR extract and estimate sample sizes and sigmas.

scape-package gives an overview of the package.

Examples

# NOT RUN {
## Exploring candidate sigmas:

getSigmaI(x.cod)  # sigma used in assessment 0.20
estSigmaI(x.cod)  # model fit implies 0.17

## Multiple series:

getSigmaI(x.oreo, "c")                 # sigma used in assessment
getSigmaI(x.oreo, "c", digits=2)       # rounded
estSigmaI(x.oreo, "c")                 # model fit implies smaller sigma

getSigmaI(x.oreo, "c", series="Series 2-1")  # get one series
estSigmaI(x.oreo, "c", series="Series 2-1")  # estimate one series
# }

Run the code above in your browser using DataCamp Workspace