sensSAM {numOSL}

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)