CS229 Lecture notes

Andrew Ng

Mixtures of Gaussians and the EM algorithm

In this set of notes, we discuss the EM (Expectation-Maximization) for den-
sity estimation.
Suppose that we are given a training set {x(1) , . . . , x(m) } as usual. Since
we are in the unsupervised learning setting, these points do not come with
any labels.
We wish to model the data by specifying a joint distribution P p(x(i) , z (i) ) =
p(x |z )p(z ). Here, z Multinomial() (where j 0, kj=1 j = 1,
(i) (i) (i) (i)

and the parameter j gives p(z (i) = j),), and x(i) |z (i) = j N (j , j ). We
let k denote the number of values that the z (i) s can take on. Thus, our
model posits that each x(i) was generated by randomly choosing z (i) from
{1, . . . , k}, and then x(i) was drawn from one of k Gaussians depending on
z (i) . This is called the mixture of Gaussians model. Also, note that the
z (i) s are latent random variables, meaning that theyre hidden/unobserved.
This is what will make our estimation problem difficult.
The parameters of our model are thus , and . To estimate them, we
can write down the likelihood of our data:
X
m
(, , ) = log p(x(i) ; , , )
i=1
X
m X
k
= log p(x(i) |z (i) ; , )p(z (i) ; ).
i=1 z (i) =1

However, if we set to zero the derivatives of this formula with respect to

the parameters and try to solve, well find that it is not possible to find the
maximum likelihood estimates of the parameters in closed form. (Try this
yourself at home.)
The random variables z (i) indicate which of the k Gaussians each x(i)
had come from. Note that if we knew what the z (i) s were, the maximum

1
 2

likelihood problem would have been easy. Specifically, we could then write
down the likelihood as
X
m
(, , ) = log p(x(i) |z (i) ; , ) + log p(z (i) ; ).
i=1

Maximizing this with respect to , and gives the parameters:

1 X
m
j = 1{z (i) = j},
m i=1
Pm (i)
i=1 1{z = j}x(i)
j = P m (i) = j}
,
i=1 1{z
Pm (i)
i=1 1{z P = j}(x(i) j )(x(i) j )T
j = m (i) = j}
.
i=1 1{z

Indeed, we see that if the z (i) s were known, then maximum likelihood
estimation becomes nearly identical to what we had when estimating the
parameters of the Gaussian discriminant analysis model, except that here
the z (i) s playing the role of the class labels.1
However, in our density estimation problem, the z (i) s are not known.
What can we do?
The EM algorithm is an iterative algorithm that has two main steps.
Applied to our problem, in the E-step, it tries to guess the values of the
z (i) s. In the M-step, it updates the parameters of our model based on our
guesses. Since in the M-step we are pretending that the guesses in the first
part were correct, the maximization becomes easy. Heres the algorithm:

Repeat until convergence: {

(E-step) For each i, j, set

(i)
wj := p(z (i) = j|x(i) ; , , )
1
There are other minor differences in the formulas here from what wed obtained in
PS1 with Gaussian discriminant analysis, first because weve generalized the z (i) s to be
multinomial rather than Bernoulli, and second because here we are using a different j
for each Gaussian.
 3

(M-step) Update the parameters:

1 X (i)
m
j := w ,
m i=1 j
Pm (i) (i)
i=1 wj x
j := Pm (i) ,
i=1 wj
Pm (i) (i)
i=1 wj (x j )(x(i) j )T
j := Pm (i)
i=1 wj

In the E-step, we calculate the posterior probability of our parameters

the z (i) s, given the x(i) and using the current setting of our parameters. I.e.,
using Bayes rule, we obtain:

p(x(i) |z (i) = j; , )p(z (i) = j; )

p(z (i) = j|x(i) ; , , ) = Pk
(i) (i) = l; , )p(z (i) = l; )
l=1 p(x |z

Here, p(x(i) |z (i) = j; , ) is given by evaluating the density of a Gaussian

with mean j and covariance j at x(i) ; p(z (i) = j; ) is given by j , and so
(i)
on. The values wj calculated in the E-step represent our soft guesses2 for
the values of z (i) .
Also, you should contrast the updates in the M-step with the formulas we
had when the z (i) s were known exactly. They are identical, except that in-
stead of the indicator functions 1{z (i) = j} indicating from which Gaussian
(i)
each datapoint had come, we now instead have the wj s.
The EM-algorithm is also reminiscent of the K-means clustering algo-
rithm, except that instead of the hard cluster assignments c(i), we instead
(i)
have the soft assignments wj . Similar to K-means, it is also susceptible
to local optima, so reinitializing at several different initial parameters may
be a good idea.
Its clear that the EM algorithm has a very natural interpretation of
repeatedly trying to guess the unknown z (i) s; but how did it come about,
and can we make any guarantees about it, such as regarding its convergence?
In the next set of notes, we will describe a more general view of EM, one
2
The term soft refers to our guesses being probabilities and taking values in [0, 1]; in
contrast, a hard guess is one that represents a single best guess (such as taking values
in {0, 1} or {1, . . . , k}).
 4

that will allow us to easily apply it to other estimation problems in which

there are also latent variables, and which will allow us to give a convergence
guarantee.