makeCalArray: Makes a calibration curve probability array

In UCL/ADMUR: Ancient Demographic Modelling Using Radiocarbon

Makes a calibration curve probability array

Description

Generates CalArray containing a 2D probability array of the calibration curve

Usage

makeCalArray(calcurve, calrange, inc=5)

Arguments

calcurve	A calibration curve object. Choose from intcal20 (default), shcal20, intcal13 or shcal13.
calrange	A vector of two calendar dates BP, giving the calendar range of CalArray. Can be in either order.
inc	Increments to interpolate calendar years. Default = 5

Details

Generates an array of probabilities mapping the calibration curve and its error ribbon.

Each column represents a Gaussian PDF constructed from the C14 date and error of the calibration curve (typically a column every 5 cal years).

Row names of CalArray are 14C dates, column names are cal dates.

This function is memory and time costly if used to construct the entire 50,000 year range of the calibration curve at a resolution of 1 cal years, therefore typically used only for a constrained date range, by specifying calrange.

This array only needs constructing once to generate a Summed Probability Distribution of any number of calibrated dates, allowing efficient downstream calibration.

Value

Returns a list including:

probs	a 2D probability array of the calibration curve
calcurve	the calibration curve as provided as an argument
cal	a numeric vector of calendar years
inc	the resolution of the array in calendar years

Examples

	# generate a CalArray of the intcal20 curve covering 5500 calBP to 7000 calBP
	x <- makeCalArray(intcal13, c(5500,7000), inc=1 )
	

UCL/ADMUR documentation built on Sept. 14, 2023, 11:41 a.m.