lk_obs_score function

Compute observed log-likelihood and score