The model is modularized in such a way that running EM on all variables simultaneously would be very difficult.
we adopt a variational approach which optimizes various subgroups of the variables in a round-robin fashion, holding approximations to the others fixed.