Skip to contents

Computes summed probability distributions (SPD) of radiocarbon dates.

Usage

c14_spd(object, ...)

# S4 method for class 'CalibratedAges'
c14_spd(object, normalize_date = FALSE, normalize_spd = FALSE)

Arguments

object

A CalibratedAges object.

...

Currently not used.

normalize_date

A logical scalar: should the total probability mass of the calibrated dates be normalised (to sum to unity within the time-span of analysis)?

normalize_spd

A logical scalar: should the total probability mass of the SPD be normalised (to sum to unity)?

Details

Summed probability distributions (SPD) are not statistically valid estimators of the calendar age of a potential future sample. They should not be used in any dates-as-data approach to provide a population proxy.

See also

Author

N. Frerebeau

Examples

## Radiocarbon data from Bosch et al. 2015
data("ksarakil")

## Calibrate
cal <- c14_calibrate(
  values = ksarakil$date,
  errors = ksarakil$error,
  names = ksarakil$code,
  curves = "marine13",
  reservoir_offsets = 53,
  reservoir_errors = 43,
  from = 50000, to = 0
)
plot(cal, level = 0.68, flip = TRUE)
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter
#> Warning: "flip" is not a graphical parameter


## SPD
s <- c14_spd(cal)
plot(s)