Assume a group of data points, such as one plotted here (this graph isn’t specific to my problem, but just used as a suitable example):
Inspecting the scatter graph visually, it’s fairly obvious the data points form two ‘groups’, with some random points that do not obviously belong to either.
I’m looking for an algorithm, that would allow me to:
- start with a data set of two or more dimensions.
- detect such groups from the dataset without prior knowledge on how many (or if any) might be there
- once the groups have been detected, ‘ask’ the model of groups, if a new sample point seems to fit to any of the groups
There are many choices, but if you are interested in the probability that a new data point belongs to a particular mixture, I would use a probabilistic approach such as Gaussian mixture modeling either estimated by maximum likelihood or Bayes.
Maximum likelihood estimation of mixtures models is implemented in Matlab.
Your requirement that the number of components is unknown makes your model more complex. The dominant probabilistic approach is to place a Dirichlet Process prior on the mixture distribution and estimate by some Bayesian method. For instance, see this paper on infinite Gaussian mixture models. The DP mixture model will give you inference over the number of components and the components each elements belong to, which is exactly what you want. Alternatively you could perform model selection on the number of components, but this is generally less elegant.
There are many implementation of DP mixture models models, but they may not be as convenient. For instance, here’s a Matlab implementation.
Your graph suggests you are an R user. In that case, if you are looking for prepacked solutions, the answer to your question lies on this Task View for cluster analysis.