summedPhaseCalibrator: Generates a summed probability distribution (SPD) after...

In AdrianTimpson/ADMUR: Ancient Demographic Modelling Using Radiocarbon

Generates a summed probability distribution (SPD) after phasing dates

Description

Generates a Summed Probability Distribution (SPD) after phasing dates

Usage

summedPhaseCalibrator(data, calcurve, calrange, inc=5, width=200)

Arguments

data	A dataframe of 14C dates. Requires 'age' and 'sd', and at least one of 'site' and 'phase'. Optional 'datingType' comprising '14C' and/or anything else.
calcurve	A calibration curve object. Choose from intcal20 (default), shcal20, intcal13 or shcal13.
calrange	A vector of two cal dates, giving the calendar range of CalArray. Can be in either order.
inc	Increments to interpolate calendar years. Default = 5.
width	A timespan in 14C years used to automatically bin dates if they have not been phased, i.e., 'phase' is missing from the data. Default = 200.

Details

Wrapper function to generate an overall SPD for phased dates. Internally this first generates an SPD for each phase. Data may be phased already, alternatively if 'phase' is not provided, this function will automatically bin dates into phases, see phaseCalibrator. Each phase's distribution is then summed, and the final SPD is normalised to unity.

Optionally 'datingType' can be provided in the data. Only '14C' will be calibrated in the usual way, anything else is assumed to be provided in calendar time. If 'datingType' is not provided, all dates are assumed to be 14C.

Value

Returns a single-column data frame of SPD probabilities. Row names are the calendar years.

Examples

	data <- subset(SAAD, site %in% c('Carrizal','Pacopampa'))
	SPD <- summedPhaseCalibrator(data, shcal20, c(2000,6000))
	plotPD(SPD)
	

AdrianTimpson/ADMUR documentation built on July 2, 2024, 8:39 p.m.