This optimization problem can be solved using a Generalized Expectation-Maximization (GEM) algorithm as described in (Mei et al., 2008a).
in the E-step of the algorithm, the distribution of hidden variables is computed using Eq.