mclust (version 1.1-7)

mhtree.VI: Classification tree for hierarchical clustering for Gaussian models with diagonal variance.

Description

Computes a classification tree for agglomerative hierarchical clustering using a Gaussian model in which clusters are spherical and may vary in volume volume.

Usage

mhtree.VI(data, partition, min.clusters = 1, alpha = 1)

Arguments

data

matrix of observations.

partition

initial classification of the data. The default puts every observation in a singleton cluster.

min.clusters

minimum number of clusters desired. The default is to carry out agglomerative hierarchical clustering until termination, that is, until all observations belong to a single group. The default value is 1.

alpha

The value alpha times the trace of the sample crossproduct matrix of all the observations divided by the product of the data dimensions, is used for the purpose of initalization. The default value is 1.

Value

• an object of class "mhtree", which consists of a classification tree with the following attributes:
• call
  a copy of the call to mhtree.VI.
• change
  value of the optimal change in likelihood at each stage.
• dimensions
  the data dimensions.
• initial.partition
  the partition at which agglomerative hierarchical clustering is initiated.

NOTES

The value alpha is needed because the criterion is not defined for singleton clusters or clusters consisting only of points that coincide.

References

J. D. Banfield and A. E. Raftery, Model-based Gaussian and non-Gaussian Clustering, Biometrics,49:803-821 (1993).

C. Fraley, Algorithms for Model-based Gaussian Hierarchical Clustering,Technical Report No. 311, Department of Statistics, University of Washington (October 1996), to appear in SIAM Journal on Scientific Computing.

See Also

mhtree, mhclass, awe, partuniq

Examples

data(iris)
mhtree.VI(iris[,1:4])