EM_algorithm