l.calib: Find the calibrated probability of a calendar age for a 14C...

In rice: Radiocarbon Equations

Find the calibrated probability of a calendar age for a 14C date.

Description

Find the calibrated probability of a cal BP age for a radiocarbon date. Can handle either multiple calendar ages for a single radiocarbon date, or a single calendar age for multiple radiocarbon dates.

Usage

l.calib(
  x,
  y,
  er,
  cc = 1,
  postbomb = FALSE,
  deltaR = 0,
  deltaSTD = 0,
  thiscurve = c(),
  cc.dir = c(),
  normal = TRUE,
  as.F = FALSE,
  is.F = FALSE,
  t.a = 3,
  t.b = 4
)

Arguments

x	The cal BP year.
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. Required for negative radiocarbon ages.
deltaR	Age offset (e.g. for marine samples).
deltaSTD	Uncertainty of the age offset (1 standard deviation).
thiscurve	As an alternative to providing cc and/or postbomb, the data of a specific curve can be provided (3 columns: cal BP, C14 age, error).
cc.dir	Directory of the calibration curves. Defaults to where the package's files are stored (system.file), but can be set to, e.g., cc.dir="curves".
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 realm. Defaults to as.F=FALSE, which uses the C14 realm.
is.F	Use this if the provided date is in the F14C realm.
t.a	Value a of the t distribution (defaults to 3).
t.b	Value b of the t distribution (defaults to 4).

Details

The function cannot deal with multiple calibration curves if multiple calendar years or radiocarbon dates are entered.

Value

The calibrated probability of a calendar age for a 14C age

Author(s)

Maarten Blaauw

Examples

  l.calib(100, 130, 20)
  l.calib(100:110, 130, 20) # multiple calendar ages of a single date
  l.calib(100, c(130,150), c(15,20)) # multiple radiocarbon ages and a single calendar age
  plot(0:300, l.calib(0:300, 130, 20), type='l')

rice documentation built on April 11, 2025, 5:53 p.m.