Okay, I understand. I also noticed that you are setting alpha in the NegativeBinomial likelihood to a static value of 1. Is there a particular reason why you are doing this? I would try setting a prior on the alpha something like pm.Exponential(lam = 1/10)
1 Like