nllik function

Calculate the predictive negative log-likelihood