Luminescence (version 0.7.1)

calc_IEU: Apply the internal-external-uncertainty (IEU) model after Thomsen et al. (2007) to a given De distribution

Description

Function to calculate the IEU De for a De data set.

Usage

calc_IEU(data, a, b, interval, decimal.point = 2, plot = TRUE, ...)

Arguments

data

RLum.Results or data.frame (required): for data.frame: two columns with De (data[,1]) and De error (values[,2])

a

numeric: slope

b

numeric: intercept

interval

numeric: fixed interval (e.g. 5 Gy) used for iteration of Dbar, from the mean to Lowest.De used to create Graph.IEU [Dbar.Fixed vs Z]

decimal.point

numeric (with default): number of decimal points for rounding calculations (e.g. 2)

plot

logical (with default): plot output

...

further arguments (trace, verbose).

Value

Returns a plot (optional) and terminal output. In addition an RLum.Results object is returned containing the following element:The output should be accessed using the function get_RLum.

Function version

0.1.0 (2016-05-02 09:36:06)

How to cite

Smedley, R.K. (2017). calc_IEU(): Apply the internal-external-uncertainty (IEU) model after Thomsen et al. (2007) to a given De distribution. Function version 0.1.0. In: Kreutzer, S., Dietze, M., Burow, C., Fuchs, M.C., Schmidt, C., Fischer, M., Friedrich, J. (2017). Luminescence: Comprehensive Luminescence Dating Data Analysis. R package version 0.7.1. https://CRAN.R-project.org/package=Luminescence

Details

This function uses the equations of Thomsen et al. (2007). The parameters a and b are estimated from dose-recovery experiments.

References

Smedley, R.K., 2015. A new R function for the Internal External Uncertainty (IEU) model. Ancient TL 33, 16-21.

Thomsen, K.J., Murray, A.S., Boetter-Jensen, L. & Kinahan, J., 2007. Determination of burial dose in incompletely bleached fluvial samples using single grains of quartz. Radiation Measurements 42, 370-379.

See Also

plot, calc_CommonDose, calc_CentralDose, calc_FiniteMixture, calc_FuchsLang2001, calc_MinDose

Examples

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

## apply the IEU model
ieu <- calc_IEU(ExampleData.DeValues$CA1, a = 0.2, b = 1.9, interval = 1)