robCompositions (version 1.9.1)

isomLR: Isometric log-ratio transformation

Description

An isometric log-ratio transformation with a special choice of the balances according to Hron et al. (2010).

Usage

isomLR(x, fast=FALSE)

Arguments

x

object of class data.frame or matrix with positive entries

fast

if TRUE, approx. 10 times faster but with numerical problems for high-dimensional data.

Value

• The isomLR transformed data.

Details

The isomLR transformation moves D-part compositional data from the simplex into a (D-1)-dimensional real space isometrically. From this choice of the balances, all the relative information of the part $x_1$ from the remaining parts is separated. It is useful for estimating missing values in $x_1$ by regression of the remaining variables.

References

Egozcue J.J., V. Pawlowsky-Glahn, G. Mateu-Figueras and C. Barcel'o-Vidal (2003) Isometric logratio transformations for compositional data analysis. Mathematical Geology, 35(3) 279-300. \

Hron, K. and Templ, M. and Filzmoser, P. (2010) Imputation of missing values for compositional data using classical and robust methods Computational Statistics and Data Analysis, vol 54 (12), pages 3095-3107.

See Also

isomLRinv

Examples

require(MASS)
Sigma <- matrix(c(5.05,4.95,4.95,5.05), ncol=2, byrow=TRUE)
z <- isomLRinv(mvrnorm(100, mu=c(0,2), Sigma=Sigma))

data(expenditures)
isomLR(expenditures)

x <- exp(mvrnorm(2000, mu=rep(1,10), diag(10)))
system.time(isomLR(x))
system.time(isomLR(x, fast=TRUE))