[ADMB Users] Beta-binomial random number generation

dave fournier davef at otter-rsch.com
Mon Nov 4 10:00:27 PST 2013

If you need to generate beta distributed random numbers you could try

double inv_cumd_beta_stable(double a, double b, double y, double eps =

and adjust the 1.e-7 as needed.  It is in there to keep the function 
from blowing up
in mixed models.

More information about the Users mailing list