Could you be a bit more specific with what you’re trying to achieve? Glancing at the code snippet, it looks like you’re implementing a Beta-Binomial model to model click-through rates (CTR). The Beta-Binomial model is actually a conjugate model, so you wouldn’t even need PyMC3 for that: there’s a closed-form expression for the posterior!
Perhaps I’m misunderstanding?