Mixture model selection via BIC

Description

Mixture model selection via BIC.

Usage

1
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").

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.

Value

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

mod

A message informing the user about the best model.

BIC

The BIC values for every possible model and number of components.

Author(s)

Michail Tsagris

R implementation and documentation: Michail Tsagris <mtsagris@yahoo.gr> and Giorgos Athineou <athineou@csd.uoc.gr>

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

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

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.