likelihood_saem function

likelihood_saem