Let $f(x;\theta)$ be the poisson frequency function with mean $\lambda$. and $p(\lambda)$ the Gamma distribution with mean $\mu$, and variance $\mu^2/\alpha$.
I have to show that $g(x)=\frac{\Gamma(x+\alpha)}{x!\Gamma(\alpha)}\left(\frac{\alpha}{\alpha+\mu}\right)^{\alpha}\left(\frac{\mu}{\alpha+\mu}\right)^{x}$.
However, from my calculations:
$\displaystyle g(x)=\int f(x;\theta)p(\lambda) d\lambda=\int\frac{e^{-\lambda}\lambda^x}{x!}\cdot\frac{\mu^{\alpha}\lambda^{\alpha-1}}{\Gamma(\alpha)}e^{-\mu\lambda}d\lambda = \frac{\Gamma(\alpha+x)}{x!\Gamma(\alpha)}\left(1/(\mu+1\right))^{x+\alpha}\mu^{\alpha}$, where the last equality is obtained by doing a coordinate change to $z=\lambda(\mu+1)$
So where did I go wrong?
Any help will be appreciated.