Estimating the number of clusters

For a collection of points in R^n, is there a statistic that I can compute which will estimate the number of clusters with some level of confidence?

-

This is an age old question. Which actually does (I think even cannot have) definite answer. Because, first you need to define what you mean by a cluster and so on. Famous saying in this regard is that "cluster is in the eye of a beholder", it is easy to construct examples where somebody could see one cluster but somebody else more than one.

This being said, MDL (minimum description length) principle would lead you to device a (IMHO) a most principled way clustering cost function. Which by optimizing you could the find the cluster assignments and number of clusters simultaneously. For multinomial data you can see following: P.Kontkanen, P.Myllymäki, W.Buntine, J.Rissanen, H.Tirri, An MDL Framework for Data Clustering. In Advances in Minimum Description Length: Theory and Applications, edited by P. Grünwald, I.J. Myung and M. Pitt. The MIT Press, 2005.

The intuitively appealing idea behind MDL clustering is that by clustering you create a model of the data. So assumption is that very good model is such which lets you compress the data well.

Anyway MDL might not be easy to apply, if you are looking for a practical way to detect the number of clusters. BIC (Bayesian information criteria) and F-ratio have proven to work ok in practice.

BIC formulation can be found in: http://staff.utia.cas.cz/nagy/skola/Projekty/Classification/Xmeans.pdf

-

I don't know if this is the most effective answer, but it has to be one of the coolest:

Barcodes: the persistent topology of data

It lets you estimate not only the number of connected components but also the higher Betti numbers. (For instance, if you take a random sample of enough points from an annulus, you will be able to measure both H0 = Z and H1 = Z.)

-
 I don't think persistent homology actually would help with this. The question doesn't seem to be about data analysis, but clustering of data. – Josh Roberts Oct 21 2009 at 11:37 Isn't persistent 0th homology pretty relevant to the OP's problem? – Darsh Ranjan Oct 21 2009 at 16:06 Yes, that's true. – Josh Roberts Oct 21 2009 at 17:48 Actually, this "Barcodes" paper points to projection of high-dimensional datasets into lower dimensions and then recognizing shapes. The projection vector used to decrease the dimensionality of the data points has a strong effect on the possible separability of the clusters. I actually don't know if there's a good way to pre-estimate the number of expected clusters without considering a priori factors about what the data set it. But the paper does not necessarily address what the best projective vector might be. – ABh Aug 28 2010 at 21:43

Carlsson has developed methods from applied topology to do clustering work. Robert Ghrist called a talk about this "Clusterfunctor" since it involves a functor from metric spaces to "persistent sets". It's talked about in Topology and Data, a survey article about using topology to do data analysis.

-

I think that for computing a level of confidence you need a model of clusters, for example a gaussian mixture model.

If the level of confidence is not an issue and you need simply the "right" number of clusters I would repeat what Ville said plus the following article:

Estimating the number of clusters in a dataset via the Gap statistic (2000) Tibshirani, Walther, Hastie http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.34.6664

-

There is a whole research field on this, and the answers so far have only highlighted a few of the specific algorithms developed. Check out this Wikipedia article for an overview and some references.

-