log_likelihood function

Calculate log-likelihood