R/calc_HomogeneityTest.R
calc_HomogeneityTest.Rd
A simple homogeneity test for De estimates
calc_HomogeneityTest(data, log = TRUE, ...)
data | RLum.Results or data.frame (required):
for data.frame: two columns with De |
---|---|
log | logical (with default): perform the homogeneity test with (un-)logged data |
... | further arguments (for internal compatibility only). |
Returns a 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
For details see Galbraith (2003).
0.3.0
Burow, C., Kreutzer, S., 2021. calc_HomogeneityTest(): Apply a simple homogeneity test after Galbraith (2003). Function version 0.3.0. 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
Galbraith, R.F., 2003. A simple homogeneity test for estimates of dose obtained using OSL. Ancient TL 21, 75-77.
Christoph Burow, University of Cologne (Germany), Sebastian Kreutzer, IRAMAT-CRP2A, Université Bordeaux Montaigne (France) , RLum Developer Team
## load example data data(ExampleData.DeValues, envir = environment()) ## apply the homogeneity test calc_HomogeneityTest(ExampleData.DeValues$BT998)#> #> [calc_HomogeneityTest()] #> #> --------------------------------- #> n: 25 #> --------------------------------- #> mu: 7.9812 #> G-value: 155.5127 #> Degrees of freedom: 24 #> P-value: 0 #> --------------------------------- #>#> #> [RLum.Results-class] #> originator: calc_HomogeneityTest() #> data: 3 #> .. $summary : data.frame #> .. $data : data.frame #> .. $args : list #> additional info elements: 1## using the data presented by Galbraith (2003) df <- data.frame( x = c(30.1, 53.8, 54.3, 29.0, 47.6, 44.2, 43.1), y = c(4.8, 7.1, 6.8, 4.3, 5.2, 5.9, 3.0)) calc_HomogeneityTest(df)#> #> [calc_HomogeneityTest()] #> #> --------------------------------- #> n: 7 #> --------------------------------- #> mu: 3.7727 #> G-value: 19.2505 #> Degrees of freedom: 6 #> P-value: 0.0038 #> --------------------------------- #>#> #> [RLum.Results-class] #> originator: calc_HomogeneityTest() #> data: 3 #> .. $summary : data.frame #> .. $data : data.frame #> .. $args : list #> additional info elements: 1