The plots in the first post were at n_experiments=100, while the second are at n_experiments=1000? Between them, the variance in the posterior of both P and N has shrunk considerably, and the mean of the posterior distribution is zeroing in on the true parameter value. I don’t think you’re doing anything wrong.