Deep Unsupervised Clustering with Gaussian Mixture Variational Autoencoders
We study a variant of the variational autoencoder model (VAE) with a Gaussian mixture as a prior distribution, with the goal of performing unsupervised clustering through deep generative models. We observe that the known problem of over-regularisation that has been shown to arise in regular VAEs also manifests itself in our model and leads to cluster degeneracy. We show that a heuristic called minimum information constraint that has been shown to mitigate this effect in VAEs can also be applied to improve unsupervised clustering performance with our model. Furthermore we analyse the effect of this heuristic and provide an intuition of the various processes with the help of visualizations. Finally, we demonstrate the performance of our model on synthetic data, MNIST and SVHN, showing that the obtained clusters are distinct, interpretable and result in achieving competitive performance on unsupervised clustering to the state-of-the-art results.READ FULL TEXT VIEW PDF
Clustering high-dimensional data, such as images or biological measureme...
Recent advances in deep learning have shown their ability to learn stron...
Outcome regressed on class labels identified by unsupervised clustering ...
Resting-state functional connectivity states are often identified as clu...
Detecting stellar clusters have always been an important research proble...
Extracting insight from the enormous quantity of data generated from
Quantitative medical image computing (radiomics) has been widely applied...
Deep Unsupervised Clustering with Gaussian Mixture Variational Autoencoders
Denoising Variational Autoencoder
Unsupervised clustering remains a fundamental challenge in machine learning research. While long-established methods such as
-means and Gaussian mixture models (GMMs)(Bishop, 2006) still lie at the core of numerous applications (Aggarwal & Reddy, 2013), their similarity measures are limited to local relations in the data space and are thus unable to capture hidden, hierarchical dependencies in latent spaces. Alternatively, deep generative models can encode rich latent structures. While they are not often applied directly to unsupervised clustering problems, they can be used for dimensionality reduction, with classical clustering techniques applied to the resulting low-dimensional space (Xie et al., 2015). This is an unsatisfactory approach as the assumptions underlying the dimensionality reduction techniques are generally independent of the assumptions of the clustering techniques.
Deep generative models try to estimate the density of observed data under some assumptions about its latent structure, i.e., its hidden causes. They allow us to reason about data in more complex ways than in models trained purely through supervised learning. However, inference in models with complicated latent structures can be difficult. Recent breakthroughs in approximate inference have provided tools for constructing tractable inference algorithms. As a result of combining differentiable models with variational inference, it is possible to scale up inference to datasets of sizes that would not have been possible with earlier inference methods(Rezende et al., 2014). One popular algorithm under this framework is the variational autoencoder (VAE) (Kingma & Welling, 2013; Rezende et al., 2014).
In this paper, we propose an algorithm to perform unsupervised clustering within the VAE framework. To do so, we postulate that generative models can be tuned for unsupervised clustering by making the assumption that the observed data is generated from a multimodal prior distribution, and, correspondingly, construct an inference model that can be directly optimised using the reparameterization trick. We also show that the problem of over-regularisation in VAEs can severely effect the performance of clustering, and that it can be mitigated with the minimum information constraint introduced by Kingma et al. (2016).
Unsupervised clustering can be considered a subset of the problem of disentangling latent variables, which aims to find structure in the latent space in an unsupervised manner. Recent efforts have moved towards training models with disentangled latent variables corresponding to different factors of variation in the data. Inspired by the learning pressure in the ventral visual stream, Higgins et al. (2016) were able to extract disentangled features from images by adding a regularisation coefficient to the lower bound of the VAE. As with VAEs, there is also effort going into obtaining disentangled features from generative adversarial networks (GANs) (Goodfellow et al., 2014). This has been recently achieved with InfoGANs (Chen et al., 2016a)
, where structured latent variables are included as part of the noise vector, and the mutual information between these latent variables and the generator distribution is then maximised as a mini-max game between the two networks. Similarly, Tagger(Greff et al., 2016), which combines iterative amortized grouping and ladder networks, aims to perceptually group objects in images by iteratively denoising its inputs and assigning parts of the reconstruction to different groups. Johnson et al. (2016) introduced a way to combine amortized inference with stochastic variational inference in an algorithm called structured VAEs. Structured VAEs are capable of training deep models with GMM as prior distribution. Shu et al. (2016) introduced a VAE with a multimodal prior where they optimize the variational approximation to the standard variational objective showing its performance in video prediction task.
The work that is most closely related to ours is the stacked generative semi-supervised model (M1+M2) by Kingma et al. (2014)
. One of the main differences is the fact that their prior distribution is a neural network transformation of both continuous and discrete variables, with Gaussian and categorical priors respectively. The prior for our model, on the other hand, is a neural network transformation of Gaussian variables, which parametrise the means and variances of a mixture of Gaussians, with categorical variables for the mixture components. Crucially,Kingma et al. (2014) apply their model to semi-supervised classification tasks, whereas we focus on unsupervised clustering. Therefore, our inference algorithm is more specific to the latter.
VAEs are the result of combining variational Bayesian methods with the flexibility and scalability provided by neural networks (Kingma & Welling, 2013; Rezende et al., 2014). Using variational inference it is possible to turn intractable inference problems into optimisation problems (Wainwright & Jordan, 2008), and thus expand the set of available tools for inference to include optimisation techniques as well. Despite this, a key limitation of classical variational inference is the need for the likelihood and the prior to be conjugate in order for most problems to be tractably optimised, which in turn can limit the applicability of such algorithms. Variational autoencoders introduce the use of neural networks to output the conditional posterior (Kingma & Welling, 2013)
and thus allow the variational inference objective to be tractably optimised via stochastic gradient descent and standard backpropagation. This technique, known as the reparametrisation trick, was proposed to enable backpropagation through continuous stochastic variables. While under normal circumstances backpropagation through stochastic variables would not be possible without Monte Carlo methods, this is bypassed by constructing the latent variables through the combination of a deterministic function and a separate source of noise. We refer the reader toKingma & Welling (2013) for more details.
In regular VAEs, the prior over the latent variables is commonly an isotropic Gaussian. This choice of prior causes each dimension of the multivariate Gaussian to be pushed towards learning a separate continuous factor of variation from the data, which can result in learned representations that are structured and disentangled. While this allows for more interpretable latent variables (Higgins et al., 2016), the Gaussian prior is limited because the learnt representation can only be unimodal and does not allow for more complex representations. As a result, numerous extensions to the VAE have been developed, where more complicated latent representations can be learned by specifying increasingly complex priors (Chung et al., 2015; Gregor et al., 2015; Eslami et al., 2016).
In this paper we choose a mixture of Gaussians as our prior, as it is an intuitive extension of the unimodal Gaussian prior. If we assume that the observed data is generated from a mixture of Gaussians, inferring the class of a data point is equivalent to inferring which mode of the latent distribution the data point was generated from. While this gives us the possibility to segregate our latent space into distinct classes, inference in this model is non-trivial. It is well known that the reparametrisation trick which is generally used for VAEs cannot be directly applied to discrete variables. Several possibilities for estimating the gradient of discrete variables have been proposed (Glynn, 1990; Titsias & Lázaro-Gredilla, 2015). Graves (2016) also suggested an algorithm for backpropagation through GMMs. Instead, we show that by adjusting the architecture of the standard VAE, our estimator of the variational lower bound of our Gaussian mixture variational autoencoder (GMVAE) can be optimised with standard backpropagation through the reparametrisation trick, thus keeping the inference model simple.
Consider the generative model , where an observed sample is generated from a set of latent variables , and under the following process:
where is a predefined number of components in the mixture, and , and are given by neural networks with parameters and , respectively. That is, the observed sample is generated from a neural network observation model parametrised by and the continuous latent variable . Furthermore, the distribution of is a Gaussian mixture with means and variances specified by another neural network model parametrised by and with input .
More specifically, the neural network parameterised by outputs a set of means and variances , given as input. A one-hot vector
is sampled from the mixing probability, which chooses one component from the Gaussian mixture. We set the parameter to make uniformly distributed. The generative and variational views of this model are depicted in Fig. 1.
The generative model is trained with the variational inference objective, i.e. the log-evidence lower bound (ELBO), which can be written as
We assume the mean-field variational family as a proxy to the posterior which factorises as , where indexes over data points. To simplify further notation, we will drop and consider one data point at a time. We parametrise each variational factor with the recognition networks and that output the parameters of the variational distributions and specify their form to be Gaussian posteriors. We derived the -posterior, , as:
The lower bound can then be written as,
We refer to the terms in the lower bound as the reconstruction term, conditional prior term, -prior term and -prior term respectively.
The reconstruction term can be estimated by drawing Monte Carlo samples from , where the gradient can be backpropagated with the standard reparameterisation trick (Kingma & Welling, 2013). The -prior term can be calculated analytically.
Importantly, by constructing the model this way, the conditional prior term can be estimated using Eqn. 5 without the need to sample from the discrete distribution .
Since can be computed for all with one forward pass, the expectation over it can be calculated in a straightforward manner and backpropagated as usual. The expectation over can be estimated with Monte Carlo samples and the gradients can be backpropagated via the reparameterisation trick. This method of calculating the expectation is similar to the marginalisation approach of Kingma et al. (2014), with a subtle difference. Kingma et al. (2014) need multiple forward passes to obtain each component of the -posterior. Our method requires wider output layers of the neural network parameterised by , but only need one forward pass. Both methods scale up linearly with the number of clusters.
The most unusual term in our ELBO is the -prior term. The -posterior calculates the clustering assignment probability directly from the value of and , by asking how far is from each of the cluster positions generated by . Therefore, the -prior term can reduce the KL divergence between the -posterior and the uniform prior by concurrently manipulating the position of the clusters and the encoded point . Intuitively, it would try to merge the clusters by maximising the overlap between them, and moving the means closer together. This term, similar to other KL-regularisation terms, is in tension with the reconstruction term, and is expected to be over-powered as the amount of training data increases.
The possible overpowering effect of the regularisation term on VAE training has been described numerous times in the VAE literature (Bowman et al., 2015; Sønderby et al., 2016; Kingma et al., 2016; Chen et al., 2016b). As a result of the strong influence of the prior, the obtained latent representations are often overly simplified and poorly represent the underlying structure of the data. So far there have been two main approaches to overcome this effect: one solution is to anneal the KL term during training by allowing the reconstruction term to train the autoencoder network before slowly incorporating the regularization from the KL term (Sønderby et al., 2016). The other main approach involves modifying the objective function by setting a cut-off value that removes the effect of the KL term when it is below a certain threshold (Kingma et al., 2016). As we show in the experimental section below, this problem of over-regularisation is also prevalent in the assignment of the GMVAE clusters and manifests itself in large degenerate clusters. While we show that the second approach suggested by Kingma et al. (2016) does indeed alleviate this merging phenomenon, finding solutions to the over-regularization problem remains a challenging open problem.
The main objective of our experiments is not only to evaluate the accuracy of our proposed model, but also to understand the optimisation dynamics involved in the construction of meaningful, differentiated latent representations of the data. This section is divided in three parts:
We first study the inference process in a low-dimensional synthetic dataset, and focus in particular on how the over-regularisation problem affects the clustering performance of the GMVAE and how to alleviate the problem;
We then evaluate our model on an MNIST unsupervised clustering task; and
We finally show generated images from our model, conditioned on different values of the latent variables, which illustrate that the GMVAE can learn disentangled, interpretable latent representations.
Throughout this section we make use of the following datasets:
Synthetic data: We create a synthetic dataset mimicking the presentation of Johnson et al. (2016), which is a 2D dataset with 10,000 data points created from the arcs of 5 circles.
MNIST: The standard handwritten digits dataset, composed of 28x28 grayscale images and consisting of 60,000 training samples and 10,000 testing samples (LeCun et al., 1998).
SVHN: A collection of 32x32 images of house numbers (Netzer et al., 2011). We use the cropped version of the standard and the extra training sets, adding up to a total of approximately 600,000 images.
We quantify clustering performance by plotting the magnitude of the -prior term described in Eqn. 6 during training. This quantity can be thought of as a measure of how much different clusters overlap. Since our goal is to achieve meaningful clustering in the latent space, we would expect this quantity to go down as the model learns the separate clusters.
Empirically, however, we have found this not to be the case. The latent representations that our model converges to merges all classes into the same large cluster instead of representing information about the different clusters, as can be seen in Figs. 1(d) and 2(a). As a result, each data point is equally likely to belong to any of clusters, rendering our latent representations completely uninformative with respect to the class structure.
We argue that this phenomenon can be interpreted as the result of over-regularisation by the -prior term. Given that this quantity is driven up by the optimisation of KL term in the lower bound, it reaches its maximum possible value of zero, as opposed to decreasing with training to ensure encoding of information about the classes. We suspect that the prior has too strong of an influence in the initial training phase and drives the model parameters into a poor local optimum that is hard to be driven out off by the reconstruction term later on.
This observation is conceptually very similar to the over-regularisation problem encountered in regular VAEs and we thus hypothesize that applying similar heuristics should help alleviate the problem. We show in Fig. 1(f) that by using the previously mentioned modification to the lower-bound proposed by Kingma et al. (2016), we can avoid the over-regularisation caused by the -prior. This is achieved by maintaining the cost from the -prior at a constant value until it exceeds that threshold. Formally, the modified -prior term is written as:
This modification suppresses the initial effect of the -prior to merge all clusters thus allowing them to spread out until the cost from the -prior cost is high enough. At that point its effect is significantly reduced and is mostly limited to merging individual clusters that are overlapping sufficiently. This can be seen clearly in Figs. 1(e) and 1(f). The former shows the clusters before the -prior cost is taken into consideration, and as such the clusters have been able to spread out. Once the -prior is activated, clusters that are very close together will be merged as seen in Fig. 1(f).
Finally, in order to illustrate the benefits of using neural networks for the transformation of the distributions, we compare the density observed by our model (Fig. 1(c)) with a regular GMM (Fig. 1(c)) in data space. As illustrated by the figures, the GMVAE allows for a much richer, and thus more accurate representations than regular GMMs, and is therefore more successful at modelling non-Gaussian data.
: (a) Data is distributed with 5 modes on the 2 dimensional data space. (b) GMVAE learns the density model that can model data using a mixture of non-Gaussian distributions in the data space. (c) GMM cannot represent the data as well because of the restrictive Gaussian assumption. (d) GMVAE, however, suffers from over-regularisation and can result in poor minima when looking at the latent space. (e) Using the modification to the ELBO(Kingma et al., 2016) allows the clusters to spread out. (f) As the model converges the -prior term is activated and regularises the clusters in the final stage by merging excessive clusters.
We now assess the model’s ability to represent discrete information present in the data on an image clustering task. We train a GMVAE on the MNIST training dataset and evaluate its clustering performance on the test dataset. To compare the cluster assignments given by the GMVAE with the true image labels we follow the evaluation protocol of Makhzani et al. (2015), which we summarise here for clarity. In this method, we find the element of the test set with the highest probability of belonging to cluster and assign that label to all other test samples belonging to . This is then repeated for all clusters , and the assigned labels are compared with the true labels to obtain an unsupervised classification error rate.
While we observe the cluster degeneracy problem when training the GMVAE on the synthetic dataset, the problem does not arise with the MNIST dataset. We thus optimise the GMVAE using the ELBO directly, without the need for any modifications. A summary of the results obtained on the MNIST benchmark with the GMVAE as well as other recent methods is shown in Table 1. We achieve classification scores that are competitive with the state-of-the-art techniques111It is worth noting that shortly after our initial submission, Rui Shu published a blog post (http://ruishu.io/2016/12/25/gmvae/) with an analysis on Gaussian mixture VAEs. In addition to providing insightful comparisons to the aforementioned M2 algorithm, he implements a version that achieves competitive clustering scores using a comparably simple network architecture. Crucially, he shows that model M2 does not use discrete latent variables when trained without labels. The reason this problem is not as severe in the GMVAE might possibly be the more restrictive assumptions in the generative process, which helps the optimisation, as argued in his blog., except for adversarial autoencoders (AAE). We suspect the reason for this is, again, related to the KL terms in the VAE’s objective. As indicated by Hoffman et al., the key difference in the adversarial autoencoders objective is the replacement of the KL term in the ELBO by an adversarial loss that allows the latent space to be manipulated more carefully (Hoffman & Johnson, 2016). Details of the network architecture used in these experiments can be found in Appendix A.
Empirically, we observe that increasing the number of Monte Carlo samples and the number of clusters makes the GMVAE more robust to initialisation and more stable as shown in Fig. 4. If fewer samples or clusters are used then the GMVAE can occasionally converge faster to poor local minima, missing some of the modes of the data distribution.
|Method||K||Best Run||Average Run|
|CatGAN (Springenberg, 2015)||20||90.30||-|
|AAE (Makhzani et al., 2015)||16||-||90.45 2.05|
|AAE (Makhzani et al., 2015)||30||-||95.90 1.13|
|DEC (Xie et al., 2015)||10||84.30||-|
|GMVAE (M = 1)||10||87.31||77.78 5.75|
|GMVAE (M = 10)||10||88.54||82.31 3.75|
|GMVAE (M = 1)||16||89.01||85.09 1.99|
|GMVAE (M = 10)||16||96.92||87.82 5.33|
|GMVAE (M = 1)||30||95.84||92.77 1.60|
|GMVAE (M = 10)||30||93.22||89.27 2.50|
So far we have argued that the GMVAE picks up natural clusters in the dataset, and that these clusters share some structure with the actual classes of the images. Now we train the GMVAE with on MNIST to show that the learnt components in the distribution of the latent space actually represent meaningful properties of the data. First, we note that there are two sources of stochasticity in play when sampling from the GMVAE, namely
Sampling from its prior, which will generate the means and variances of through a neural network ; and
Sampling from the Gaussian mixture determined by and , which will generate the image through a neural network .
In Fig. 4(a) we explore the latter option by setting and sampling multiple times from the resulting Gaussian mixture. Each row in Fig. 4(a) corresponds to samples from a different component of the Gaussian mixture, and it can be clearly seen that samples from the same component consistently result in images from the same class of digit. This confirms that the learned latent representation contains well differentiated clusters, and exactly one per digit. Additionally, in Fig. 4(b) we explore the sensitivity of the generated image to the Gaussian mixture components by smoothly varying and sampling from the same component. We see that while reliably controls the class of the generated image, sets the “style” of the digit.
Finally, in Fig. 6 we show images sampled from a GMVAE trained on SVHN, showing that the GMVAE clusters visually similar images together.
We have introduced a class of variational autoencoders in which one level of the latent encoding space has the form of a Gaussian mixture model, and specified a generative process that allows us to formulate a variational Bayes optimisation objective. We then discuss the problem of over-regularisation in VAEs. In the context of our model, we show that this problem manifests itself in the form of cluster degeneracy. Crucially, we show that this specific manifestation of the problem can be solved with standard heuristics.
We evaluate our model on unsupervised clustering tasks using popular datasets and achieving competitive results compared to the current state of the art. Finally, we show via sampling from the generative model that the learned clusters in the latent representation correspond to meaningful features of the visible data. Images generated from the same cluster in latent space share relevant high-level features (e.g. correspond to the same MNIST digit) while being trained in an entirely unsupervised manner.
It is worth noting that GMVAEs can be stacked by allowing the prior on
to be a Gaussian mixture distribution as well. A deep GMVAE could scale much better with number of clusters given that it would be combinatorial with regards to both number of layers and number of clusters per layer. As such, while future research on deep GMVAEs for hierarchical clustering is a possibility, it is crucial to also address the enduring optimisation challenges associated with VAEs in order to do so.
We would like to acknowledge the NVIDIA Corporation for the donation of a GeForce GTX Titan Z used in our experiments. We would like to thank Jason Rolfe, Rui Shu and the reviewers for useful comments. Importantly, we would also like to acknowledge that the variational family which we used throughout this version of the paper was suggested by an anonymous reviewer.
Draw: A recurrent neural network for image generation.In Proceedings of The 32nd International Conference on Machine Learning, pp. 1462–1471, 2015.
Early Visual Concept Learning with Unsupervised Deep Learning.ArXiv e-prints, June 2016.
Workshop in Advances in Approximate Bayesian Inference, NIPS, 2016.
European Conference on Computer Vision (ECCV) Workshop on Action and Anticipation for Visual Learning, 2016.
For optimisation, we use Adam (Kingma & Ba, 2014) with a learning rate of
and standard hyperparameter values, and . The model architectures used in our experiments are shown in Tables A.1, A.2 and A.3.
fc 120 ReLU 120 ReLU
|= 2, = 2 (Exp),|
|= 2, = 2 (Exp)|
|MNIST||28x28||conv 16x6x6* (1-0) 32x6x6* (1-0)||= 150, = 150 (Exp),|
|64x4x4* (2-1) 500*||= 200, = 200 (Exp)|
|SVHN||32x32||conv 64x4x4* (2-1) 128x4x4* (2-1)||= 150, = 150 (Exp),|
|246x4x4* (2-1) 500*||= 200, = 200 (Exp)|
. We exponentiate the variance components to keep their value positive. An asterisk (*) indicates the use of batch normalization and a ReLU nonlinearity. For convolutional layers, the numbers in parentheses indicate stride-padding.
|Synthetic||2||fc 120 Tanh|
|MNIST||150||fc 500 Tanh|
|SVHN||150||fc 500 Tanh|
|Synthetic||2||fc 120 ReLU 120 ReLU|
|MNIST||200||500* full-conv 64x4x4* (2-1) 32x6x6* (1-0)||28x28 (Sigmoid)|
|SVHN||200||500* full-conv 246x4x4* (2-1) 128x4x4* (2-1)||32x32 (Sigmoid)|