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).
Details
This function uses the equations of Thomsen et al. (2007). The parameters a
and b are estimated from dose-recovery experiments.
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)
Author(s)
Rachel Smedley, Geography & Earth Sciences, Aberystwyth University
(United Kingdom) Based on an excel spreadsheet and accompanying macro
written by Kristina Thomsen.
R Luminescence Package Team
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.
## 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)
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]
A motivated R-Team member: 'We are doing this not just for statistical reasons, there is real science behind it!'
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Luminescence/calc_IEU.Rd_%03d_medium.png", width=480, height=480)
> ### Name: calc_IEU
> ### Title: Apply the internal-external-uncertainty (IEU) model after
> ### Thomsen et al. (2007) to a given De distribution
> ### Aliases: calc_IEU
>
> ### ** 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)
[calc_IEU]
Dbar: 46.67
IEU.De (Gy): 46.67
IEU.Error (Gy): 2.55 Number of De: 24
a: 0.2000
b: 1.9000
>
>
>
>
>
>
> dev.off()
null device
1
>