Description Details Author(s) References Examples

The Rmixmod package provides model-based clustering by fitting a mixture model (e.g. Gaussian components for quantitative continuous data) to the data and identifying each cluster with one of its components. The number of components can be determined from the data, typically using the BIC criterion. In practice, however, individual clusters can be poorly fitted by Gaussian distributions, and in that case model-based clustering tends to represent one non-Gaussian cluster by a mixture of two or more Gaussian components. This package, RmixmodCombi, following the article in the references, combines the components of the EM/BIC solution (provided by the Rmixmod package) hierarchically according to an entropy criterion. This yields a clustering for each number of clusters less than or equal to K. These clusterings can be compared on substantive grounds, and we also provide a way of selecting the number of clusters via a piecewise linear regression fit to the (possibly rescaled) entropy plot.

Package: | RmixmodCombi |

Type: | Package |

Version: | 1.0 |

Date: | 2013-09-20 |

Depends: | Rmixmod, Rcpp |

See the Rmixmod package documentation for more details about fitting mixture models. See the cited article for more details about the combining methodology.

J.-P. Baudry and G. Celeux Maintainer: J.-P. Baudry <[email protected]>

J.-P. Baudry, A. E. Raftery, G. Celeux, K. Lo and R. Gottardo (2010). Combining mixture components for clustering. *Journal of Computational and Graphical Statistics, 19(2):332-353.*

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 | ```
##### Example of quantitative data #####
set.seed(1)
data(Baudry_etal_2010_JCGS_examples)
res <- mixmodCombi(ex4.1, nbCluster = 1:8)
res # is of class MixmodCombi
res@mixmodOutput # is the initial EM/BIC solution (provided by mixmodCluster or by the user as a
# [\code{\linkS4class{MixmodCluster}}] object) from which the hierarchy is computed
res@hierarchy[[3]] # is the 3-cluster solution obtained by hierarchically combining the initial
# EM/BIC solution
## Not run:
plot(res) # This is a simulated example where the clusters should obviously not be identified to
# Gaussian components...
hist(res, nbCluster = 4)
## End(Not run)
##### Example of qualtitative data #####
set.seed(1)
data(car)
res <- mixmodCombi(car[1:300,], nbCluster = 1:10) # Only the 300 first observations for a
# quick example
res # is of class MixmodCombi
res@mixmodOutput # is the initial EM/BIC solution (provided by mixmodCluster or by the user as a
# [\code{\linkS4class{MixmodCluster}}] object) from which the hierarchy is computed
res@hierarchy[[res@ICLNbCluster]] # is the solution obtained by hierarchically combining the initial
# EM/BIC solution for the number of clusters selected with ICL
## Not run:
plot(res)
barplot(res)
## End(Not run)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.