Description Usage Arguments Details Value Function version How to cite Note Author(s) References See Also Examples

View source: R/analyse_SAR.CWOSL.R

The function performs a SAR CW-OSL analysis on an RLum.Analysis object including growth curve fitting.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 |

`object` |
RLum.Analysis ( |

`signal.integral.min` |
integer ( |

`signal.integral.max` |
integer ( |

`background.integral.min` |
integer ( |

`background.integral.max` |
integer ( |

`rejection.criteria` |
list ( Allowed arguments are |

`dose.points` |
numeric ( |

`mtext.outer` |
character ( |

`plot` |
logical ( |

`plot.single` |
logical ( |

`onlyLxTxTable` |
logical (with default): If |

`...` |
further arguments that will be passed to the function
plot_GrowthCurve or calc_OSLLxTxRatio
(supported: |

The function performs an analysis for a standard SAR protocol measurements
introduced by Murray and Wintle (2000) with CW-OSL curves. For the
calculation of the Lx/Tx value the function calc_OSLLxTxRatio is
used. For **changing the way the Lx/Tx error is calculated** use the argument
`background.count.distribution`

and `sigmab`

, which will be passed to the function
calc_OSLLxTxRatio.

**Argument object is of type list**

If the argument `object`

is of type list containing **only**
RLum.Analysis objects, the function re-calls itself as often as elements
are in the list. This is usefull if an entire measurement wanted to be analysed without
writing separate for-loops. To gain in full control of the parameters (e.g., `dose.points`

) for
every aliquot (corresponding to one RLum.Analysis object in the list), in
this case the arguments can be provided as list. This `list`

should
be of similar length as the `list`

provided with the argument `object`

,
otherwise the function will create an own list of the requested lenght.
Function output will be just one single RLum.Results object.

Please be careful when using this option. It may allow a fast an efficient data analysis, but the function may also break with an unclear error message, due to wrong input data.

**Working with IRSL data**

The function was originally designed to work just for 'OSL' curves, following the principles of the SAR protocol. An IRSL measurement protocol may follow this procedure, e.g., post-IR IRSL protocol (Thomsen et al., 2008). Therefore this functions has been enhanced to work with IRSL data, however, the function is only capable of analysing curves that follow the SAR protocol structure, i.e., to analyse a post-IR IRSL protocol, curve data have to be pre-selected by the user to fit the standards of the SAR protocol, i.e., Lx,Tx,Lx,Tx and so on.

Example: Imagine the measurement contains pIRIR50 and pIRIR225 IRSL curves. Only one curve type can be analysed at the same time: The pIRIR50 curves or the pIRIR225 curves.

**Supported rejection criteria**

`[recycling.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).

`[testdose.error]`

: set the allowed error for the testdose, which per
default should not exceed 10\

`[palaeodose.error]`

: set the allowed error for the De value, which per
default should not exceed 10\

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

`data` |
data.frame containing De-values, De-error 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. |

`Formula` |
formula formula that have been used for the growth curve fitting |

The output should be accessed using the function get_RLum.

0.8.8

Kreutzer, S., 2020. analyse_SAR.CWOSL(): Analyse SAR CW-OSL measurements. Function version 0.8.8. In: Kreutzer, S., Burow, C., Dietze, M., Fuchs, M.C., Schmidt, C., Fischer, M., Friedrich, J., 2020. Luminescence: Comprehensive Luminescence Dating Data Analysis. R package version 0.9.7. https://CRAN.R-project.org/package=Luminescence

This function must not be mixed up with the function Analyse_SAR.OSLdata, which works with Risoe.BINfileData objects.

**The function currently does support only 'OSL', 'IRSL' and 'POSL' data!**

Sebastian Kreutzer, IRAMAT-CRP2A, UMR 5060, CNRS-Université Bordeaux Montaigne (France) , RLum Developer Team

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

Duller, G., 2003. Distinguishing quartz and feldspar in single grain luminescence measurements. Radiation Measurements, 37 (2), 161-165.

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.

Thomsen, K.J., Murray, A.S., Jain, M., Boetter-Jensen, L., 2008. Laboratory fading rates of various luminescence signals from feldspar-rich sediment extracts. Radiation Measurements 43, 1474-1486. doi:10.1016/j.radmeas.2008.06.002

calc_OSLLxTxRatio, plot_GrowthCurve, RLum.Analysis, RLum.Results, get_RLum

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 | ```
##load data
##ExampleData.BINfileData contains two BINfileData objects
##CWOSL.SAR.Data and TL.SAR.Data
data(ExampleData.BINfileData, envir = environment())
##transform the values from the first position in a RLum.Analysis object
object <- Risoe.BINfileData2RLum.Analysis(CWOSL.SAR.Data, pos=1)
##perform SAR analysis and set rejection criteria
results <- analyse_SAR.CWOSL(
object = object,
signal.integral.min = 1,
signal.integral.max = 2,
background.integral.min = 900,
background.integral.max = 1000,
log = "x",
fit.method = "EXP",
rejection.criteria = list(
recycling.ratio = 10,
recuperation.rate = 10,
testdose.error = 10,
palaeodose.error = 10,
exceed.max.regpoint = TRUE)
)
##show De results
get_RLum(results)
##show LnTnLxTx table
get_RLum(results, data.object = "LnLxTnTx.table")
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.