def log_predictive_likelihood(self, data_point, params): ll = log_beta(params.a + data_point, params.b + (1 - data_point)) ll -= log_beta(params.a, params.b) return np.sum(ll)