lk_obs_score_between function

Compute observed log-likelihood and score