Changed run gibbs sampler to run em algorithm