How to infer parameters for a continuous mixture model: Poisson-Beta model

Thank you so much for your time and help! The code works perfectly!