A vector of 1-sigma values for the ages given above
calCurves
A vector of values containing either 'intcal13', 'shcal13', 'marine13', or 'normal'. Should be the same length the number of ages supplied. Non-standard calibration curves can be used provided they are supplied in the same format as those previously mentioned and are placed in the same directory. Normal indicates a normally-distributed (non-14C) age.
pathToCalCurves
(optional) File path to where the calibration curves are located. Defaults to the system directory where the 3 standard calibration curves are stored.
dfs
(optional) Degrees-of-freedom values for the t-distribution associated with the calibration calculation. A large value indicates Gaussian distributions assumed for the 14C ages
samples
(optional) Number of samples of calibrated dates required
G
(optional) Number of Gaussian mixture components
Details
This is a faster approximate version of BchronDensity that uses the densityMclust function to compute the Gaussian mixtures for a set of calibrated ages. The method is an approximation as it does not fit a fully Bayesian model as BchronDensity does. It is designed to be a probabilistic version of the Oxcal SUM command which takes calibrated ages and sums the probability distributions with the aim of estimating activity through age as a proxy.
Value
An object of class BchronDensityRunFast with the following components:
out
The output from the run of densityMclust with the given number of mixture components
calAges
The calibrated ages from the BchronDensity function
Author(s)
Andrew Parnell <andrew.parnell@ucd.ie>
References
See the mclust package for references about this method.
See Also
Bchronology, BchronCalibrate, BchronRSL, BchronDensity for a slower exact version of this function
Examples
## Not run: # Read in some data from Sluggan Moss
data(Sluggan)
# Run the model
SlugDensFast = BchronDensityFast(ages=Sluggan$ages,ageSds=Sluggan$ageSds,
calCurves=Sluggan$calCurves)
# plot it
plot(SlugDensFast)
## End(Not run)