Function reference
Calibration of a Single Radiocarbon Sample
Function to independently calibrate a single radiocarbon determination against a given calibration curve. The rest of the library concerns calibration of multiple (related) radiocarbon samples.
-
CalibrateSingleDetermination() - Calibrate a Single Radiocarbon Determination
Bayesian Non-Parametric Joint Density Estimation
Functions providing a rigorous Bayesian non-parametric alternative to summed probability distributions (SPDs). These functions enable the calibration and calendar age summarisation of multiple related radiocarbon samples (and provide plotting of the results)
-
PolyaUrnBivarDirichlet() - Calibrate and Summarise Multiple Radiocarbon Samples via a Bayesian Non-Parametric DPMM (with Polya Urn Updating)
-
WalkerBivarDirichlet() - Calibrate and Summarise Multiple Radiocarbon Samples via a Bayesian Non-Parametric DPMM (with Walker Updating)
-
FindPredictiveCalendarAgeDensity() - Find Predictive Estimate of Shared Calendar Age Density from Bayesian Non-Parametric DPMM Output
-
PlotPredictiveCalendarAgeDensity() - Plot Predictive Estimate of Shared Calendar Age Density from Bayesian Non-Parametric DPMM Output
-
PlotNumberOfClusters() - Plot Number of Calendar Age Clusters Estimated in Bayesian Non-Parametric DPMM Output
-
PlotCalendarAgeDensityIndividualSample() - Plot Posterior Calendar Age Estimate for an Individual Determination after Joint Calibration
-
PlotConvergenceData() - Plot KL Divergence of Predictive Density to Assess Convergence of Bayesian Non-Parametric DPMM Sampler
Poisson Process Modelling
Functions for modelling the occurrence of radiocarbon samples as a variable-rate (inhomogeneous) Poisson process. This is a further (analogous) approach that enables the rigorous summarisation of calendar age information from multiple related radiocarbon samples.
-
PPcalibrate() - Model Occurrence of Multiple Radiocarbon Samples as a Variable-Rate Poisson Process
-
FindPosteriorMeanRate() - Find Posterior Mean Rate of Sample Occurrence for Poisson Process Model
-
PlotPosteriorMeanRate() - Plot Posterior Mean Rate of Sample Occurrence for Poisson Process Model
-
PlotNumberOfInternalChanges() - Plot Number of Changepoints in Rate of Sample Occurrence for Poisson Process Model
-
PlotPosteriorChangePoints() - Plot Calendar Ages of Changes in Rate of Sample Occurrence for Poisson Process Model
-
PlotPosteriorHeights() - Plot Heights of Segments in Rate of Sample Occurrence for Poisson Process Model
-
PlotCalendarAgeDensityIndividualSample() - Plot Posterior Calendar Age Estimate for an Individual Determination after Joint Calibration
Common convergence plotting functions
Functions that can be used to plot information on convergence of the MCMC for both the Bayesian non-parametric DPMM and the Poisson process approaches to calibration and summarisation.
-
PlotGelmanRubinDiagnosticSingleChain() - Plot Histogram of the Gelman-Rubin Convergence Diagnostic for a Single MCMC Chain
-
PlotGelmanRubinDiagnosticMultiChain() - Plot Histogram of the Gelman-Rubin Convergence Diagnostic for Multiple Independent MCMC Chains
Helper functions
General functions that may be of use. However, please do not use SPDs in your inference (use either the Bayesian non-parametric or Poisson Process approaches described above).
-
FindSummedProbabilityDistribution() - Find the summed probability distribution (SPD) for a set of radiocarbon observations
-
InterpolateCalibrationCurve() - Interpolate a calibration curve at a set of calendar ages
-
GenerateOxcalCode() - Outputs code suitable for running in OxCal from a series of radiocarbon determinations
Example datasets
Sets of radiocarbon determinations, either from real-life examples or artificially-generated, that can be used as examples with the calibration functions.
-
alces - Example real-life data - Alces in Yukon and Alaska
-
armit - Example real-life data - Population Decline in Iron Age Ireland
-
bison - Example real-life data - Bison in Yukon and Alaska
-
buchanan - Example real-life data - Palaeo-Indian demography
-
cervus - Example real-life data - Cervus in Yukon and Alaska
-
equus - Example real-life data - Equus in Yukon and Alaska
-
human - Example real-life data - Humans in Yukon and Alaska
-
kerr - Example real-life data - Irish Rath
-
mammuthus - Example real-life data - Mammuthus in Yukon and Alaska
-
pp_uniform_phase - Example artificial data - Uniform Phase
-
two_normals - Example artificial data - Mixture of Normal Phases