alfa.mix.norm: Gaussian mixture models for compositional data using the...

In Compositional: Compositional Data Analysis

Gaussian mixture models for compositional data using the \alpha-transformation

Description

Gaussian mixture models for compositional data using the \alpha-transformation.

Usage

alfa.mix.norm(x, g, a, model, veo = FALSE)

Arguments

x	A matrix with the compositional data.
g	How many clusters to create.
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.
model	The type of model to be used. "EII": All groups have the same diagonal covariance matrix, with the same variance for all variables. "VII": Different diagonal covariance matrices, with the same variance for all variables within each group. "EEI": All groups have the same diagonal covariance matrix. "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. "EVI": Different diagonal covariance matrices with the same determinant. "VVI": Different diagonal covariance matrices, with nothing in common. "EEE": All covariance matrices are the same. "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. "VEV": Different covariance matrices but if we make the matrices have determinant 1, then they will have the same trace. "VVV": Different covariance matrices with nothing in common. "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. "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. "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. "EVV": Different covariance matrices, but with the same determinant.
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.

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

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

Examples

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

Compositional documentation built on Oct. 9, 2024, 5:10 p.m.