log_mlik_all function

Fit models and compute the log marginal likelihood for all clusters