Beta binomial prior used here (no choose constant)

dbb_prior(x, N, u, v)