add EM estimation