I am trying to do this using the Poisson-Gamma mixture parameterization.
I think the derivation is related to the derivation of the expected value. But I have only seen how to derive the mean using conditioning, not from summation. Could someone help?