Exact beta-binomial density using a multiprecision library.

dBetaBinom_MP calculates the beta-binomial density while avoiding numerical mistakes (catastrophic cancellations) due to extreme parameter values. This function is called by dBetaBinom if necessary, and should not be called outside of this.

dBetaBinom_MP(ms, ns, piX, thetaX, LOG, NecPres)

Arguments

ms

Numeric vector. Vector of number of successes

ns

Numeric vector. Vector of number of trials

piX

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

thetaX

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

LOG

Logical. if TRUE, return log-densities

NecPres

Number. Necessary Precision, i.e. number of bits, for an accurate density calculatation, as determined by the function dBetaBinom

Value

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