def log_marginal_likelihood(self, params): return np.sum(log_beta(params.a, params.b) - log_beta(self.priors.a, self.priors.b))