Perturbation operation {Compositional}

Perturbation operation

Description

Perturbation operation.

Usage

perturbation(x, y, oper = "+")

Arguments

x

A matrix with the compositional data.

y

Either a matrix with compositional data or a vector with compositional data. In either case, the data may not be compositional data, as long as they non negative.

oper

For the summation this must be "*" and for the negation it must be "/". According to Aitchison (1986), multiplication is equal to summation in the log-space, and division is equal to negation.

Details

This is the perturbation operation defined by Aitchison (1986).

Value

A matrix with the perturbed compositional data.

Author(s)

Michail Tsagris.

R implementation and documentation: Michail Tsagris [e-mail].

References

Aitchison J. (1986). The statistical analysis of compositional data. Chapman & Hall.

See Also

power

Examples

x <- as.matrix(iris[1:15, 1:4])
y <- as.matrix(iris[21:35, 1:4])
perturbation(x, y)
perturbation(x, y[1, ])