Compositional (version 5.0)

Perturbation operation: 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.

Value

A matrix with the perturbed compositional data.

Details

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

References

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

See Also

power

Examples

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