\operatorname{P}(X=k) = \int\limits_{0}^{\infty} \frac{\lambda^{k}}{k!}e^{-\lambda} \,\,\pi(\lambda)\,\mathrm d\lambda