Luminescence (version 0.6.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:

summary

data.frame summary of all relevant model results.

data

data.frame original input data

args

list used arguments

call

call the function call

tables

list a list of data frames containing all calculation tables

The output should be accessed using the function get_RLum.

Function version

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

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)