Implements a beta-binomial pair of conjugate distributions.
In this model, samples are distributed according to \[ \pi \sim \mathrm {Beta}(\alpha, \beta) \\ x \sim \mathrm {Bernoulli} (\pi) \] The nextDouble() method returns a sample of \(x\) and the nextMean() returns a sample of \(\pi\).