Skip to contents

Bayesian Credible Interval

Usage

# S4 method for class 'CalibratedAges'
interval_credible(x, level = 0.954, n = 100, ...)

Arguments

x

A CalibratedAges object.

level

A length-one numeric vector giving the confidence level.

n

An integer specifying the number of random samples.

...

Currently not used.

Value

A CalibratedIntervals object.

See also

Author

N. Frerebeau

Examples

## Calibrate multiple dates
cal <- c14_calibrate(
  values = c(5000, 4500),
  errors = c(45, 35),
  names = c("X", "Y")
)

## Credible intervals
crd68 <- interval_credible(cal, level = 0.683)
crd95 <- interval_credible(cal, level = 0.954)
crd99 <- interval_credible(cal, level = 0.997)

## Coerce to data.frame
as.data.frame(crd95, calendar = BC())
#>   label start  end     p
#> 1     X  3917 3663 0.954
#> 2     Y  3340 3094 0.954

## Plot
plot(cal, interval = "credible")