Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

mix.compnorm

Gaussian mixture models for compositional data

Description

Gaussian mixture models for compositional data.

Usage

mix.compnorm(x, g, model = NULL, 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.If NULL all 14 of the ones listed below are fitted. Otherwise, choose one from below.

  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

Either the additive ("alr") or the isometric ("ilr") log-ratio transformation is to be used..

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 eahc 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 mtsagris@uoc.gr and Giorgos Athineou <gioathineou@gmail.com>.

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

Examples

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

## End(Not run)
Compositional

Compositional Data Analysis

v4.6
GPL (>= 2)
Authors
Michail Tsagris [aut, cre], Giorgos Athineou [aut], Abdulaziz Alenazi [ctb]
Initial release
2021-04-27

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.