A Comparison of R esampling Methods for Clustering Ensembles

1. A Comparison of Resampling Methods for Clustering Ensembles Behrouz Minaei, Alexander Topchyand William Punch Department of Computer Science and Engineering MLMTA 2004, Las Vegas, June 22th 2004

2. Outline • Clustering Ensemble • How to generate different partitions? • How to combine multiple partitions? • Resampling Methods • Bootstrap vs. Subsampling • Experimental study • Methods • Results • Conclusion

3. Ensemble Benefits • Combinations of classifiers proved to be very effective in supervised learning framework, e.g. bagging and boosting algorithms • Distributed data mining requires efficient algorithms capable to integrate the solutions obtained from multiple sources of data and features • Ensembles of clusterings can provide novel, robust, and stable solutions

4. 2 clusters 2 clusters 3 clusters 6 3 clusters 4 2 0 -2 -4 -6 -5 0 5 Is Meaningful Clustering Combination Possible? “Combination” of 4 different partitions can lead to true clusters!

5. Pattern Matrix, Distance matrix

6. Representation of Multiple Partitions • Combination of partitions can be viewed as another clustering problem, where each Pi represents a new feature with categorical values • Cluster membership of a pattern in different partitions is regarded as a new feature vector • Combining the partitions is equivalent to clustering these tuples 7 objects clustered by 4 algorithms

7. Re-labeling and Voting

8. Co-association As Consensus Function • Similarity between objects can be estimated by the number of clusters shared by two objects in all the partitions of an ensemble • This similarity definition expresses the strength of co-association of n objects by an n x n matrix • xi: the i-th pattern; pk(xi): cluster label of xi in the k-th partition; I(): Indicator function; N = no. of different partitions • This consensus function eliminates the need for solving the label correspondence problem

9. Clustering Combination Approaches Generative mechanism Consensus function Single link Different algorithms Comp. link Co-association-based Avg. link Different initialization for one algorithm Voting approach Others … Deterministic Information Theory approach Different subsets of objects Resampling CSPA Different subsets of features Hypergraph methods HGPA Project to 1D MCLA Projection to subspaces Mixture Model (EM) Rand. cuts/plane Taxonomy of Clustering Combination Approaches

10. Resampling Methods • Bootstrapping (Sampling with replacement) • Create an artificial list by randomly drawing N elements from that list. Some elements will be picked more than once. • Statistically on average 37% of elements are repeated • Subsampling (Sampling without replacement) • Control over the size of subsample

11. Experiment: Data sets

12. Subsampling results on Star/Galaxy

13. Subsampling results on Halfrings

14. Subsampling on Halfrings

15. Error Rate for Individual Clustering

17. Taxonomy of Clustering Methods Clustering • Clustering methods call for particular clustering algorithms • Hundreds of clustering algorithms proposed Model-based Density-based Heuristic Pattern matrix Proximity matrix Mixture model (Gaussian mixture, Latent class) Spatial Clustering Linkage methods (single-link, complete-link, CHAMELEON) Graph-theoretic (MST, spectral clustering, Min-cut) Prototype-based (k-means, k-medoid) Mode seeking (mean-shift) Kernel-based (DENCLUE) Grid-based (Wave-Cluster, STING)

18. Clustering Combination • Observations • Each clustering algorithm addresses differently the issues of number of clusters, and the structure imposed on the data, and produces different data partitions • A single clustering algorithm can produce distinct results on the same data set depending on parameter values • Instead of a single clustering algorithm, use multiple algorithms • Success of the ensemble of classifiers in supervised learning (bagging and boosting) is the motivation

19. 1.2 1.2 1 1 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 0 0.5 1 -0.2 0 0.5 1 Are There Any Clusters in Data? • Most clustering algorithms find clusters, even if the data is uniform! 100 uniform data points k-means, k=3

20. Clustering in new space • Essentially consensus must be found by clustering objects in new space – space of cluster labels • In fact, features of new space are extracted from “old” features by different clustering algorithms • New space can possibly be extended by adding original “old” features • Problem becomes very rich with infinite ways to solve

21. How to Generate the Ensemble? • Apply different clustering algorithms • Single-link, EM, complete-link, spectral clustering, k-means • Use random initializations of the same algorithm • Output of k-means and EM depends on initial cluster centers • Use different parameter settings of the same algorithm • number of clusters, “width” parameter in spectral clustering, thresholds in linkage methods • Use different feature subsets or data projections • Re-sample data with or without replacement

22. Consensus Functions • Co-association matrix (Fred & Jain 2002) • Similarity between 2 patterns estimated by counting the number of shared clusters • Single-link with max life-time for finding the consensus partition • Hyper-graph methods (Strehl & Ghosh 2002) • Clusters in different partitions represented by hyper-edges • Consensus partition found by a k-way min-cut of the hyper-graph • Re-labeling and voting (Fridlyand & Dudoit 2001) • If the label correspondence problem is solved for the given partitions, a simple voting procedure can be used to assign objects in clusters • Mutual Information (Topchy, Jain & Punch 2003) • Maximize the mutual information between the individual partitions and the target consensus partition • Finite mixture model (Topchy, Jain & Punch 2003) • Maximum likelihood solution to latent class analysis problem in the space of cluster labels via EM

23. Consensus functions: Co-association approach (Fred 2001, Fred & Jain 2002) • Similarity between objects can be estimated by the number of clusters shared by two objects in all the partitions of an ensemble. • This similarity definition expresses the strength of co-association of objects by a matrix containing the values: • Thus, one can use numerous similarity-based clustering algorithms by applying them to the matrix of co-association values.

24. Consensus functions: Hypergraph approach (Strehl & Ghosh 2002) • All the clusters in the ensemble partitions can be represented as hyperedges on a graph with N vertices. • Each hyperedge describes a set of objects belonging to the same clusters. • A consensus function is formulated as a solution to k-way min-cut hypergraph partitioning problem. Each connected component after the cut corresponds to a cluster in the consensus partition. • Hypergraph partitioning problem is NP-hard, but very efficient heuristics are developed for its solution with complexity proportional to the number of hyperedges O(|E|).

25. Consensus functions: Mutual information(Topchy et al. 2003a) • Mutual information (MI) between the empirical probability distribution of labels in the consensus partition and the labels in the ensemble must be maximized. Under the assumption of independence of partitions, mutual information can be written as sum of pair-wise MIs between target and given partitions. • An elegant solution can be obtained from a generalized definition of MI. Quadratic MI information can be effectively maximized by the k-means algorithm in the space of specially transformed cluster labels of given ensemble. Computational complexity of the algorithm is O(kNH)

26. Consensus functions: Finite Mixture model(Topchy et al. 2003b) • The main assumption is that the labels yi are modeled as random variables drawn from a probability distribution described as a mixture of multivariate multinomial component densities:

27. Consensus functions: Finite Mixture model The objective of consensus clustering is formulated as a maximum likelihood estimation problem. To find the best fitting mixture density for a given data Y we must maximize the likelihood function with respect to the unknown parameters : EM algorithm is used to solve this maximum likelihood problem

28. Consensus functions: Re-labeling and Voting(Fridlyand & Dudoit 2001) • If a label correspondence problem is solved for all given partitions, then a simple voting procedure can be used to assign objects in clusters. However, label correspondence is exactly what makes unsupervised combination difficult. • A heuristic approximation to consistent labeling is possible. All the partitions in the ensemble can be re-labeled according to their best agreement with some chosen reference partition. • The reference partition can be taken as one from the ensemble, or from a new clustering of the dataset.

29. Subsampling results on Galaxy/Star