Last data update: 2014.03.03

R: Function to test maximum valid l value.
test.lR Documentation

Function to test maximum valid l value.

Description

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 
>