R/calc_IEU.R
calc_IEU.Rd
Function to calculate the IEU De for a De data set.
calc_IEU(data, a, b, interval, decimal.point = 2, plot = TRUE, ...)
data | RLum.Results or data.frame (required):
for data.frame: two columns with De |
---|---|
a | numeric (required): slope |
b | numeric (required): intercept |
interval | numeric (required):
fixed interval (e.g. 5 Gy) used for iteration of |
decimal.point | numeric (with default): number of decimal points for rounding calculations (e.g. 2) |
plot | logical (with default): plot output |
... | further arguments ( |
Returns a plot (optional) and terminal output. In addition an RLum.Results object is returned containing the following elements:
data.frame summary of all relevant model results.
data.frame original input data
list used arguments
call the function call
list a list of data frames containing all calculation tables
This function uses the equations of Thomsen et al. (2007). The parameters a and b are estimated from dose-recovery experiments.
0.1.1
Smedley, R.K., 2021. calc_IEU(): Apply the internal-external-uncertainty (IEU) model after Thomsen et al. (2007) to a given De distribution. Function version 0.1.1. In: Kreutzer, S., Burow, C., Dietze, M., Fuchs, M.C., Schmidt, C., Fischer, M., Friedrich, J., Mercier, N., Riedesel, S., Autzen, M., Mittelstrass, D., Gray, H.J., 2021. Luminescence: Comprehensive Luminescence Dating Data Analysis. R package version 0.9.11. https://CRAN.R-project.org/package=Luminescence
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.
Rachel Smedley, Geography & Earth Sciences, Aberystwyth University (United Kingdom)
Based on an excel spreadsheet and accompanying macro written by Kristina Thomsen.
, RLum Developer Team
## 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)#> #>#> #> #> #> #> #>