def calc_logpdf_marginal(N, x_sum, alpha, beta): return betaln(x_sum + alpha, N - x_sum + beta) - betaln(alpha, beta)