This function performs the weight transformation of the data matrix after
Klovan & Imbrie (1971) and performs EMMA() with different weight limits to
check if valied results are yielded. It returns the maximum value for which
the transformation remains stable.
Usage
test.l(X, l, ...)
Arguments
X
Numeric matrix with m samples (rows) and n variables (columns).
l
Numeric vector specifying the weight transformation limit, i.e.
quantile; default is 0.
...
Further arguments passed to the function.
Value
A list with objects
step
Numeric scalar with position of last
valid value.
l.max
Numeric scalar with last valid value of l.
Author(s)
Michael Dietze, Elisabeth Dietze
References
Dietze E, Hartmann K, Diekmann B, IJmker J, Lehmkuhl F, Opitz S,
Stauch G, Wuennemann B, Borchers A. 2012. An end-member algorithm for
deciphering modern detrital processes from lake sediments of Lake Donggi
Cona, NE Tibetan Plateau, China. Sedimentary Geology 243-244: 169-180.
Klovan JE, Imbrie J. 1971. An Algorithm and FORTRAN-IV Program for
Large-Scale Q-Mode Factor Analysis and Calculation of Factor Scores.
Mathematical Geology 3: 61-77.
See Also
EMMA, check.data,
test.parameters
Examples
## load example data set
data(X, envir = environment())
test <- test.l(X = X, l = seq(from = 0, to = 0.6, by = 0.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(EMMAgeo)
Loading required package: GPArotation
Loading required package: limSolve
Loading required package: shape
Loading required package: shiny
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/EMMAgeo/test.l.Rd_%03d_medium.png", width=480, height=480)
> ### Name: test.l
> ### Title: Function to test maximum valid l value.
> ### Aliases: test.l
> ### Keywords: EMMA
>
> ### ** Examples
>
> ## load example data set
> data(X, envir = environment())
>
> test <- test.l(X = X, l = seq(from = 0, to = 0.6, by = 0.1))
>
>
>
>
>
> dev.off()
null device
1
>