Gaussian mixture models for compositional data {Compositional}

Gaussian mixture models for compositional data

Description

Gaussian mixture models for compositional data.

Usage

mix.compnorm(x, g, model, type = "alr", veo = FALSE)

Arguments

x

A matrix with the compositional data.

g

How many clusters to create.

model

The type of model to be used.

1. "EII": All groups have the same diagonal covariance matrix, with the same variance for all variables.

2. "VII": Different diagonal covariance matrices, with the same variance for all variables within each group.

3. "EEI": All groups have the same diagonal covariance matrix.

4. "VEI": Different diagonal covariance matrices. If we make all covariance matrices have determinant 1, (divide the matrix with the $p$-th root of its determinant) then all covariance matrices will be the same.

5. "EVI": Different diagonal covariance matrices with the same determinant.

6. "VVI": Different diagonal covariance matrices, with nothing in common.

7. "EEE": All covariance matrices are the same.

8. "EEV": Different covariance matrices, but with the same determinant and in addition, if we make them have determinant 1, they will have the same trace.

9. "VEV": Different covariance matrices but if we make the matrices have determinant 1, then they will have the same trace.

10. "VVV": Different covariance matrices with nothing in common.

11. "EVE": Different covariance matrices, but with the same determinant. In addition, calculate the eigenvectors for each covariance matrix and you will see the extra similarities.

12. "VVE": Different covariance matrices, but they have something in common with their directions. Calculate the eigenvectors of each covariance matrix and you will see the similarities.

13. "VEE": Different covariance matrices, but if we make the matrices have determinant 1, then they will have the same trace. In addition, calculate the eigenvectors for each covariance matrix and you will see the extra similarities.

14. "EVV": Different covariance matrices, but with the same determinant.

type

The type of trasformation to be used, either the additive log-ratio ("alr"), the isometric log-ratio ("ilr") or the pivot coordinate ("pivot") transformation.

veo

Stands for "Variables exceed observations". If TRUE then if the number variablesin the model exceeds the number of observations, but the model is still fitted.

Details

A log-ratio transformation is applied and then a Gaussian mixture model is constructed.

Value

A list including:

mu

A matrix where each row corresponds to the mean vector of each cluster.

su

An array containing the covariance matrix of each cluster.

prob

The estimated mixing probabilities.

est

The estimated cluster membership values.

Author(s)

Michail Tsagris.

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

References

Ryan P. Browne, Aisha ElSherbiny and Paul D. McNicholas (2015). R package mixture: Mixture Models for Clustering and Classification.

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

See Also

bic.mixcompnorm, rmixcomp, mix.compnorm.contour, alfa.mix.norm, alfa.knn, alfa.rda, comp.nb

Examples

x <- as.matrix(iris[, 1:4])
x <- x/ rowSums(x)
mod1 <- mix.compnorm(x, 3, model = "EII" )
mod2 <- mix.compnorm(x, 4, model = "VII")