Consider the following probabilistic scheme:
$\alpha \sim \mbox{Gamma}(a,b)$
$\rho \sim \mbox{Beta}(\alpha, \beta)$
$z_i \sim \mbox{Bernoully}(\rho) \quad$ for $i=1..N$
The variables $z_i, a, b, \beta$ are all known. What is the posterior distribution over $\alpha$?
That is, how is $\alpha | \mathbf{z}, a, b, \beta$ distributed?

