Compositional (version 2.0)

bMixture model selection via BIC: Mixture model selection via BIC

Description

Mixture model selection via BIC.

Usage

bic.mixcompnorm(x, A, type = "alr")

Arguments

x

A matrix with the compositional data.

A

The maximum number of components, clusters, to be considered.

type

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

Value

a plot with the BIC of the best model for each number of components versus the number of components. A list including:

Details

The alr or the ilr-transformation is applied to the compositional data first and then mixtures of multivariate Gaussian distributions are fitted. BIC is used to decide on the optimal model and number of components.

References

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

Ryan P. Browne and Paul D. McNicholas (2014). Estimating Common Principal Components in High Dimensions. Advances in Data Analysis and Classification, 8(2), 217-226.

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

See Also

mix.compnorm, mixnorm.contour, rmixcomp

Examples

library(MASS)
x <- iris[, 1:4]
bic.mixcompnorm(x, 6, type = "alr")
bic.mixcompnorm(x, 6, type = "ilr")