simulate.bbdml function

Simulate from beta-binomial model