rice (version 1.6.2)

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

Description

Find the probability that a sample 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,
  cc = 1,
  postbomb = FALSE,
  bombalert = TRUE,
  deltaR = 0,
  deltaSTD = 0,
  normal = TRUE,
  as.F = FALSE,
  is.F = FALSE,
  t.a = 3,
  t.b = 4,
  BCAD = FALSE,
  threshold = 0
)

Value

The probability of a date being of a certain calendar age or younger than it.

Arguments

x

The year of interest, in cal BP by default.

y

The radiocarbon date's mean.

er

The radiocarbon date's lab error.

cc

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

postbomb

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

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.

threshold

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

Author

Maarten Blaauw

Details

The function can only deal with one date at a time.

Examples

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