optimSAM {numOSL}

Optimization of statistical age models

Description

Estimating the parameters of statistical age models, including the common age model (COM), the central age model (CAM), the minimum age model (MAM), the maximum age model (MXAM), and the finite mixture age model (FMM), using the Maximum Likelihood Estimation method.

Usage

optimSAM(EDdata, model = "cam", addsigma = 0, iflog = TRUE, maxcomp = 6)

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"

addsigma

numeric(with default): additional uncertainty, i.e., the sigmab value

iflog

logical(with default): transform equivalent dose values to log-scale or not

maxcomp

integer(with default): the maximum number of components in FMM

Value

Return a list that contains the following elements:

pars

optimized parameters, the name of the parameter of COM is "COM.De", that of CAM are c("CAM.OD","CAM.De"), that of MAM3 are c("Prop","MAM3.De","Sigma"), that of MXAM3 are c("Prop","MXAM3.De","Sigma"), that of MAM4 are
c("Prop","MAM4.De","Mu","Sigma"), that of MXAM4 are
c("Prop","MXAM4.De","Mu","Sigma"), and that of FMM are c("Prop","FMM.De")

maxlik

optimized maximum logged likelihood value

bic

calculated Bayesian Information Criterion (BIC) value

References

Galbraith RF, 1988. Graphical display of estimates having differing standard errors. Technometrics, 30(3): 271-281.

Galbraith RF, 1990. The radial plot: Graphical assessment of spread in ages. International Journal of Radiation Applications and Instrumentation. Part D. Nuclear Tracks and Radiation Measurements, 17(3): 207-214.

Galbraith RF, Green P, 1990. Estimating the component ages in a finite mixture. International Journal of Radiation Applications and Instrumentation. Part D. Nuclear Tracks and Radiation Measurements, 17: 197-206.

Galbraith RF, Laslett GM, 1993. Statistical models for mixed fission track ages. Nuclear Tracks And Radiation Measurements, 21(4): 459-470.

Galbraith RF, 1994. Some applications of radial plots. Journal of the American Statistical Association, 89(428): 1232-1242.

Galbraith RF, Roberts RG, Laslett GM, Yoshida H, Olley JM, 1999. Optical dating of single grains of quartz from Jinmium rock shelter, northern Australia. Part I: experimental design and statistical models. Archaeometry, 41(2): 339-364.

Galbraith RF, 2005. Statistics for fission track analysis. Chapman & Hall/CRC Press.

Galbraith RF, 2010. On plotting OSL equivalent doses. Ancient TL, 28(1): 1-10.

Galbraith RF, Roberts RG, 2012. Statistical aspects of equivalent dose and error calculation and display in OSL dating: an overview and some recommendations. Quaternary Geochronology, 11: 1-27.

See Also

mcMAM; mcFMM; dbED; psRadialPlot; RadialPlotter; EDdata

Examples

  data(EDdata)
  
  ### Fitting a 3-component FMM.
  optimSAM(EDdata$al3, model="fmm3", addsigma=0)

  ### Fitting a 4-parameter MXAM.
  optimSAM(EDdata$al3, model="mxam4", addsigma=0.1)