Add Beta-Bernoulli model to C++