numOSL (version 2.0)

calSGCED: Equivalent dose calculation using the standardised growth curve (SGC) method

Description

Calculating equivalent doses according to the standardised growth curve (SGC) method using known parameters.

Usage

calSGCED(Data1, pars, model, origin,  ErrorMethod = "sp", outpdf = NULL)

Arguments

Data1

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

pars

vector(required): optimized parameters of the SGC obtained using function lsNORM

model

character(required): fitting model used for obtaining pars

origin

logical(required): if established SGC passes the origin or not

ErrorMethod

character(with default): method used for assessing the standard error of an equivalent dose, only ErrorMethod="sp" is available

outpdf

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

Value

Return a list that contains the following elements:

scaleLtx

scaled standardised natural dose signals and associated standard errors

sgcED

calculated SGC equivalent doses

saturate.NO

Saturated Grain.NO from which SGC equivalent doses cannot be calculated

References

Li B, Roberts RG, Jacobs Z, Li SH, 2015. Potential of establishing a "global standardised growth curve" (gSGC) for optical dating of quartz from sediments. Quaternary Geochronology, 27: 94-104.

Li B, Jacobs Z, Roberts RG, 2016. Investigation of the applicability of standardised growth curves for OSL dating of quartz from Haua Fteah cave, Libya. Quaternary Geochronology, 35: 1-15.

Roberts HM, Duller GAT, 2004. Standardised growth curves for optical dating of sediment using multiple-grain aliquots. Radiation Measurements, 38(2): 241-252.

See Also

lsNORM; fitGrowth; calED; calSARED

Examples

 data(SARdata)
 Data1 <- SARdata[is.element(SARdata[,2], c("N","R3")),]
 pars <- c(137.440874251, 0.007997863, 2.462035263, -0.321536177)
 calSGCED(Data1, pars, model="gok", origin=FALSE)