Makes a calibration curve probability array
Generates CalArray containing a 2D probability array of the calibration curve.
makeCalArray(calcurve, calrange, inc=5)
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 |
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.
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 |
# generate a CalArray of the intcal20 curve covering 5500 calBP to 7000 calBP x <- makeCalArray(intcal13, c(5500,7000), inc=1 )
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.