R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Analyse SAR TL measurements

analyse_SAR.TL

R Documentation

Analyse SAR TL measurements

Description

The function performs a SAR TL analysis on a RLum.Analysis object including growth curve fitting.

Usage

analyse_SAR.TL(object, object.background, signal.integral.min,
  signal.integral.max, sequence.structure = c("PREHEAT", "SIGNAL",
  "BACKGROUND"), rejection.criteria = list(recycling.ratio = 10,
  recuperation.rate = 10), dose.points, log = "", ...)

Arguments

`object`	`RLum.Analysis`(required): input object containing data for analysis
`object.background`	currently not used
`signal.integral.min`	integer (required): requires the channel number for the lower signal integral bound (e.g. `signal.integral.min = 100`)
`signal.integral.max`	integer (required): requires the channel number for the upper signal integral bound (e.g. `signal.integral.max = 200`)
`sequence.structure`	vector character (with default): specifies the general sequence structure. Three steps are allowed ( `"PREHEAT"`, `"SIGNAL"`, `"BACKGROUND"`), in addition a parameter `"EXCLUDE"`. This allows excluding TL curves which are not relevant for the protocol analysis. (Note: None TL are removed by default)
`rejection.criteria`	list (with default): list containing rejection criteria in percentage for the calculation.
`dose.points`	`numeric` (optional): option set dose points manually
`log`	character (with default): a character string which contains "x" if the x axis is to be logarithmic, "y" if the y axis is to be logarithmic and "xy" or "yx" if both axes are to be logarithmic. See plot.default).
`...`	further arguments that will be passed to the function `plot_GrowthCurve`

Details

This function performs a SAR TL analysis on a set of curves. The SAR procedure in general is given by Murray and Wintle (2000). For the calculation of the Lx/Tx value the function calc_TLLxTxRatio is used.

Provided rejection criteria

‘recyling.ratio’: calculated for every repeated regeneration dose point.
‘recuperation.rate’: recuperation rate calculated by comparing the Lx/Tx values of the zero regeneration point with the Ln/Tn value (the Lx/Tx ratio of the natural signal). For methodological background see Aitken and Smith (1988)

Value

A plot (optional) and an RLum.Results object is returned containing the following elements:

`De.values`	data.frame containing De-values and further parameters
`LnLxTnTx.values`	data.frame of all calculated Lx/Tx values including signal, background counts and the dose points.
`rejection.criteria`	data.frame with values that might by used as rejection criteria. NA is produced if no R0 dose point exists.

note: the output should be accessed using the function get_RLum

Function version

0.1.4 (2015-11-29 17:27:48)

Note

THIS IS A BETA VERSION

None TL curves will be removed from the input object without further warning.

Author(s)

Sebastian Kreutzer, IRAMAT-CRP2A, Universite Bordeaux Montaigne (France)
R Luminescence Package Team

References

Aitken, M.J. and Smith, B.W., 1988. Optical dating: recuperation after bleaching. Quaternary Science Reviews 7, 387-393.

Murray, A.S. and Wintle, A.G., 2000. Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol. Radiation Measurements 32, 57-73.

Examples



##load data
data(ExampleData.BINfileData, envir = environment())

##transform the values from the first position in a RLum.Analysis object
object <- Risoe.BINfileData2RLum.Analysis(TL.SAR.Data, pos=3)

##perform analysis
analyse_SAR.TL(object,
               signal.integral.min = 210,
               signal.integral.max = 220,
               log = "y",
               fit.method = "EXP OR LIN",
               sequence.structure = c("SIGNAL", "BACKGROUND"))

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(Luminescence)
Welcome to the R package Luminescence version 0.6.0 [Built: 2016-05-30 16:47:30 UTC]
Anonymous student hotel employee: 'Let me double check this.'
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Luminescence/analyse_SAR.TL.Rd_%03d_medium.png", width=480, height=480)
> ### Name: analyse_SAR.TL
> ### Title: Analyse SAR TL measurements
> ### Aliases: analyse_SAR.TL
> ### Keywords: datagen plot
> 
> ### ** Examples
> 
> 
> 
> ##load data
> data(ExampleData.BINfileData, envir = environment())
> 
> ##transform the values from the first position in a RLum.Analysis object
> object <- Risoe.BINfileData2RLum.Analysis(TL.SAR.Data, pos=3)
> 
> ##perform analysis
> analyse_SAR.TL(object,
+                signal.integral.min = 210,
+                signal.integral.max = 220,
+                log = "y",
+                fit.method = "EXP OR LIN",
+                sequence.structure = c("SIGNAL", "BACKGROUND"))
[plot_GrowthCurve()] Fit: EXP OR LIN | De = 406.86 | D01 = 615183610.2

 [RLum.Results]
	 originator: analyse_SAR.TL()
	 data: 3
 	 .. $De.values : data.frame
	 .. $LnLxTnTx.table : data.frame
	 .. $rejection.criteria : data.frame
	 additional info elements:  0There were 50 or more warnings (use warnings() to see the first 50)
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>