sensSAM: Investigate of the sensitivity of a statistical age model to...

In numOSL: Numeric Routines for Optically Stimulated Luminescence Dating

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)

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

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

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.

See Also

RadialPlotter; EDdata; optimSAM

Examples

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

numOSL documentation built on Sept. 18, 2023, 9:07 a.m.