numOSL (version 2.0)

calSARED: SAR equivalent doses calculation

Description

Calculating a series of equivalent doses using the single aliquot regenerative-dose (SAR) method (Murray and Wintle, 2000) using a wrapper of function calED.

Usage

calSARED(Data, model = "gok", origin = FALSE,  ErrorMethod = "mc", nsim = 1000, weight = TRUE,  trace = TRUE, outpdf = NULL, outfile = NULL)

Arguments

Data

data.frame(required): data used for SAR equivalent dose calculation, it should contain five columns (i.e., Grain.NO, SAR.Cycle, Dose, Signal, and Signal.Err), see SARdata for details

model

character(with default): model used for growth curve fitting, see fitGrowth for available models

origin

logical(with default): force the fitted growth curve to pass the origin or not

ErrorMethod

character(with default): method used for assessing the standard error of an equivalent dose, see function calED for details

nsim

integer(with default): number of Monte Carlo simulations used for equivalent dose error estimation if ErrorMethod="mc", see function calED for details

weight

logical(with default): perform growth curve fitting using a weighted or an un-weighted procedure, see function fitGrowth for details

trace

logical(with default): logical value indicating if a trace of the calculating progress should be printed

outpdf

character(optional): if specified, results of SAR equivalent dose calculation will be drawn to a PDF file named "outpdf" and saved to the current work directory

outfile

character(optional): if specified, calculated equivalent doses and relevant quantities will be written to a CSV file named "outfile" and saved to the current work directory

Value

Return an invisible list of S3 class "calSARED" that contains the following elements:

LMpars

parameters of fitted growth curves for each Grain.NO

N

number of analyzed Grain.NO

failFit.NO

Grain.NO failed in growth curve fitting

saturate.NO

Grain.NO with saturated standardised natural dose signal

failED.NO

Grain.NO failed in equivalent dose (error) calculation

extrapolate.NO

Grain.NO from which equivalent dose is calculated by extrapolation

largeRcy.NO

Grain.NO with recycling ratio falls outside the range 0.9-1.1

largeRcp5.NO

Grain.NO with recuperation (i.e., ratio of the standardised zero regenerative dose signal to natural dose signal) larger than 5 percent

largeRcp10.NO

Grain.NO with recuperation larger than 10 percent

tab

calculated SAR equivalent doses and relevant quantities, such as recycling ratio, recuperation, etc

References

Murray AS, Wintle AG, 2000. Luminescence dating of quartz using improved single-aliquot regenerative-dose protocol. Radiation Measurements, 32(1): 57-73.

See Also

fitGrowth; calED; calSGCED

Examples

  data(SARdata)
  # Use only the first two aliquots of SARdata.
  Data <- SARdata[1:16,]
  res <- calSARED(Data, model="gok", trace=FALSE)
  summary(res)
  print(res$LMpars)