The Frechet mean for compositional data {Compositional}

The Frechet mean for compositional data

Description

Mean vector or matrix with mean vectors of compositional data using the \alpha-transformation.

Usage

frechet(x, a)

Arguments

x

A matrix with the compositional data.

a

The value of the power transformation, it has to be between -1 and 1. If zero values are present it has to be greater than 0. If \alpha=0 the isometric log-ratio transformation is applied and the closed geometric mean is calculated. You can also provide a sequence of values of alpha and in this case a matrix of Frechet means will be returned.

Details

The power transformation is applied to the compositional data and the mean vector is calculated. Then the inverse of it is calculated and the inverse of the power transformation applied to the last vector is the Frechet mean.

Value

If \alpha is a single value, the function will return a vector with the Frechet mean for the given value of \alpha. Otherwise the function will return a matrix with the Frechet means for each value of \alpha.

Author(s)

Michail Tsagris.

R implementation and documentation: Michail Tsagris [e-mail] and Giorgos Athineou <[e-mail]>.

References

Tsagris M.T., Preston S. and Wood A.T.A. (2011). A data-based power transformation for compositional data. In Proceedings of the 4th Compositional Data Analysis Workshop, Girona, Spain. https://arxiv.org/pdf/1106.1451.pdf

See Also

alfa, alfainv, profile

Examples

library(MASS)
x <- as.matrix(fgl[, 2:9])
x <- x / rowSums(x)
frechet(x, 0.2)
frechet(x, 1)