rice (version 2.3.0)

younger: Find the probability of a calibrated distribution being of a certain age or younger than it

Description

Find the probability that a date or distribution is of a certain calendar age x or younger than it, by calculating the proportion of the calibrated distribution up to and including x (i.e., summing the calibrated distribution up to year x).

Usage

younger(
  x,
  y,
  er = c(),
  cc = 1,
  postbomb = FALSE,
  glue = 0,
  bombalert = TRUE,
  deltaR = 0,
  deltaSTD = 0,
  normal = TRUE,
  as.F = FALSE,
  is.F = FALSE,
  t.a = 3,
  t.b = 4,
  BCAD = FALSE,
  zero = FALSE,
  threshold = 0
)

Value

The probability of a date being of a certain calendar age or younger than it. Returns a vector if x contains multiple ages.

Arguments

x

The year or years of interest, in cal BP by default.

y

The radiocarbon date's mean. Alternatively, can be an age distribution (two columns: ages and their probabilities) or a vector of ages.

er

The radiocarbon date's lab error. Leave empty if y is an age distribution.

cc

calibration curve for the radiocarbon date(s) (see the rintcal package).

postbomb

Whether or not to use a postbomb curve (see caldist()).

glue

Glue postbomb and prebomb curves together. Defaults to 0 (none), can be 1 (IntCal20 + NH1), 2 (IntCal20 + NH2), 3 (IntCal20 + NH3), 4 (SHCal20 + SH1-2) or 5 (SHCal20 + SH3). Note that this will override the value of cc.

bombalert

Warn if a date is close to the lower limit of the calibration curve. Defaults to postbomb=TRUE.

deltaR

Age offset (e.g. for marine samples).

deltaSTD

Uncertainty of the age offset (1 standard deviation).

normal

Use the normal distribution to calibrate dates (default TRUE). The alternative is to use the t model (Christen and Perez 2016).

as.F

Whether or not to calculate ages in the F14C timescale. Defaults to as.F=FALSE, which uses the C14 timescale.

is.F

Use this if the provided date is in the F14C timescale.

t.a

Value a of the t distribution (defaults to 3).

t.b

Value b of the t distribution (defaults to 4).

BCAD

Which calendar scale to use. Defaults to cal BP, BCAD=FALSE.

zero

Whether or not zero BC/AD should be included if using BCAD=TRUE. Defaults to zero=FALSE.

threshold

Report only values above a threshold. Defaults to threshold=0.

Author

Maarten Blaauw

Details

Instead of providing a to-be-calibrated age (y +- er), one can also provide an age distribution (as two columns: cal BP ages and their probabilities) or a vector of ages.

Examples

younger(2800, 2450, 20)
younger(2400, 2450, 20)
calibrate(160, 20, BCAD=TRUE)
younger(1750, 160, 20, BCAD=TRUE)