sensSAM: Investigate of the sensitivity of a statistical age model to the additional uncertainty (sigmab)

Description

Estimate of the parameters of a statistical age model using a number of sigmab values.

Usage

sensSAM(EDdata, model, sigmaVEC = NULL, iflog = TRUE, 
        maxcomp = 8, plot = TRUE, outfile = NULL)

Value

Return an invisible list that contains the following elements:

pars: a list that contains the optimized parameters for each sigmab value
mat: a matrix that contains the optimized parameters, the maximum logged likelihood value, and the corresponding Bayesian Information Criterion (BIC) value

Arguments

EDdata: matrix(required): a two-column matrix (i.e., equivalent dose values and
associated standard errors)
model: character(with default): the fitting model, one of "com", "cam", "mam3", "mam4", "mxam3", "mxam4", "fmm0", "fmm1", "fmm2", ..., "fmm9"
sigmaVEC: vector(with default): a series of sigmab values that will be used as inputs for the model. For example, sigmaVEC=seq(from=0,to=0.3,by=0.01)
iflog: logical(with default): transform equivalent dose values to log-scale or not
maxcomp: integer(with default): the maximum number of components in the FMM
plot: logical(with default): logical value indicating if the results should be plotted
outfile: character(optional): if specified, the results will be written to a CSV file named "outfile" and saved to the current work directory

References

Peng J, Li B, Jacobs Z, 2020. Modelling heterogeneously bleached single-grain equivalent dose distributions: Implications for the reliability of burial dose determination. Quaternary Geochronology, 60: 101108.

Peng J, Li B, Jacobs Z, Gliganic LA, 2023. Optical dating of sediments affected by post-depositional mixing: Modelling, synthesizing and implications. Catena, 232: 107383.

Examples

Run this code


  # Not run.
  # data(EDdata)
  # sensSAM(EDdata$al3, model="mam4", iflog=TRUE)