run_estimations_with_gibbs function

Run all EM-based estimations with Gibbs sampling (CRAN-safe)