Posterior distributions of modal parameters are obtained using the VB approach. The VB approach tends to slightly underestimate uncertainties in modal parameters. The uncertainties in modal parameters can be used to remove spurious modes. The VB approach is recommended comparing to the EM algorithm and Gibbs sampler.