Added EM estimation