Skip to contents

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

Usage

sensitivity(...)

# S4 method for 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

summary()

Other statistics: credible(), hpdi(), summary()

Author

T. S. Dye, N. Frerebeau

Examples

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

## Returns 0's
sensitivity(eve, eve)
#>    mad mean sd min q1 median q3 max lower upper
#> 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