Before revision:  In a simulation study, we benchmark the performance of the proposed VB method against MSL and MCMC in terms of parameter recovery, predictive accuracy and computational efficiency.
After revision:  In a simulation study, we benchmark the performance of the proposed VB method against maximum simulated likelihood (MSL) and Markov chain Monte Carlo (MCMC) methods in terms of parameter recovery, predictive accuracy and computational efficiency.