optimSAM: Optimization of statistical age models

In numOSL: Numeric Routines for Optically Stimulated Luminescence Dating

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)

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