Generates a summed probability distribution (SPD) after phasing dates
Generates a Summed Probability Distribution (SPD) after phasing dates.
summedPhaseCalibrator(data, calcurve, calrange, inc=5, width=200)
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. |
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.
Returns a single-column data frame of SPD probabilities. Row names are the calendar years.
data <- subset(SAAD, site %in% c('Carrizal','Pacopampa'))
SPD <- summedPhaseCalibrator(data, shcal20, c(2000,6000))
plotPD(SPD)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.