def negative_log_predictive_density(dist, y):
    """
    Negative log predictive density of model (normalized by number of observations).
    """
    return -dist.log_prob(y.ravel()).item() / y.shape[0]