Mixture of two beta geometric with beta shape > 1

Hi @luca_giacomel, sorry about the delay response :sweat_smile: I tried your code some time ago and didnt manage to make it work, and it looks like a bug. I think in this particular example, it might be straightforward to write down the mixture logp directly. You can have a look of an example here:
http://docs.pymc.io/notebooks/gaussian-mixture-model-advi.html