PyMC doesn’t know how to compute convolutions of random variables (yet!) – this is why you are getting an error when you try to add n1 + n2
. Remember that addition of random variables is convolution.
In your specific case, you end up with a poisson-binomal distribution. This discussion is relevant.