Sensitivity — sensitivity • ArchaeoPhases

Calculates the ranges of summary statistics from the output of two or more runs of the MCMC algorithm.

Usage

sensitivity(...)

# S4 method for class 'EventsMCMC'
sensitivity(..., positions = NULL, level = 0.95)

Arguments

...: Any EventsMCMC object.
positions: A numeric vector specifying the positions of the columns corresponding to the MCMC chains of interest, or a character vector of column names.
level: A length-one numeric vector giving the confidence level.

Value

A data.frame.

Details

This function is useful for estimating the sensitivity of calibration results to different model parameters.

See also

Other statistics: interval_credible(), interval_hdr(), summary()

Author

T. S. Dye, N. Frerebeau

Examples

## Coerce to MCMC
eve <- as_events(mcmc_events, calendar = CE(), iteration = 1)

## Returns 0's
sensitivity(eve, eve)
#>    mad mean sd min q1 median q3 max start end
#> E1   0    0  0   0  0      0  0   0     0   0
#> E2   0    0  0   0  0      0  0   0     0   0
#> E3   0    0  0   0  0      0  0   0     0   0
#> E4   0    0  0   0  0      0  0   0     0   0