Apply the model after Fuchs & Lang (2001) to a given De distribution.

Description

This function applies the method according to Fuchs & Lang (2001) for heterogeneously bleached samples with a given coefficient of variation threshold.

Usage

calc_FuchsLang2001(data, cvThreshold = 5, startDeValue = 1, plot = TRUE,
  ...)

Arguments

Details

Used values
If the coefficient of variation (c[v]) of the first two values is larger than the threshold c[v_threshold], the first value is skipped. Use the startDeValue argument to define a start value for calculation (e.g. 2nd or 3rd value).

Basic steps of the approach

(1) Estimate natural relative variation of the sample using a dose recovery test
(2) Sort the input values ascendingly
(3) Calculate a running mean, starting with the lowermost two values and add values iteratively.
(4) Stop if the calculated c[v] exceeds the specified cvThreshold

Value

Returns a plot (optional) and terminal output. In addition an RLum.Results object is returned containing the following elements:

Function version

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

Note

Please consider the requirements and the constraints of this method (see Fuchs & Lang, 2001)

Author(s)

Sebastian Kreutzer, IRAMAT-CRP2A, Universite Bordeaux Montaigne (France) Christoph Burow, University of Cologne (Germany)
R Luminescence Package Team

References

Fuchs, M. & Lang, A., 2001. OSL dating of coarse-grain fluvial quartz using single-aliqout protocols on sediments from NE Peloponnese, Greece. In: Quaternary Science Reviews 20, 783-787.

Fuchs, M. & Wagner, G.A., 2003. Recognition of insufficient bleaching by small aliquots of quartz for reconstructing soil erosion in Greece. Quaternary Science Reviews 22, 1161-1167.

See Also

plot, calc_MinDose, calc_FiniteMixture, calc_CentralDose, calc_CommonDose, RLum.Results

Examples

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

##calculate De according to Fuchs & Lang (2001)
temp<- calc_FuchsLang2001(ExampleData.DeValues$BT998, cvThreshold = 5)

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 new pIRIR derivative: '20000 miles below the sea.'
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Luminescence/calc_FuchsLang2001.Rd_%03d_medium.png", width=480, height=480)
> ### Name: calc_FuchsLang2001
> ### Title: Apply the model after Fuchs & Lang (2001) to a given De
> ###   distribution.
> ### Aliases: calc_FuchsLang2001
> ### Keywords: dplot
> 
> ### ** Examples
> 
> 
> 
> ##load example data
> data(ExampleData.DeValues, envir = environment())
> 
> ##calculate De according to Fuchs & Lang (2001)
> temp<- calc_FuchsLang2001(ExampleData.DeValues$BT998, cvThreshold = 5)

 [calc_FuchsLang2001]

----------- meta data --------------
 cvThreshold:             5 %
 used values:             22
----------- dose estimate ----------
 mean:                    2866.11
 sd:                      157.35
 weighted mean:           2846.66
 weighted sd:             20.58
------------------------------------

> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>