Exact beta-binomial density using sums.

dBetaBinom_cpp_old calculates the beta-binomial density via a number of sums, which is slow for high-value data but fast for low-value data. This function is called by dBetaBinom if necessary, and should not be called outside of this.

dBetaBinom_cpp_old(ms, ns, pi, theta, LOG)

Arguments

ms

Numeric vector. Vector of number of successes

ns

Numeric vector. Vector of number of trials

pi

Number. Probability of success; 0 >= pi >= 1

theta

Number. Dispersion parameter; 0 >= theta > +Inf

LOG

Logical. if TRUE, return log-densities

Value

A numeric vector of the same length as ms and ns, containing (log-)beta-binomial densities