I have data that has been used to fit a Bayesian non-parametric mixture model. One of the things I can do with that is sample conditional posterior cluster assignments. Given a sufficiently large sample, I’ve used that to calculate a similarity matrix–that is, a (n x n) matrix of pairwise probabilities of falling into the same cluster.