El cambio de variable $\;u:=\dfrac {\beta}{\theta}\;$ debería ayudar mucho aquí dando $\;\theta:=\dfrac {\beta}u,\;d\theta=-\dfrac{\beta}{u^2}du\;$ y :
\begin {align} I( \alpha , \beta )&:= \dfrac { \beta ^ \alpha }{ \Gamma ( \alpha )} \int_0 ^{ \infty } { \theta ^{-( \alpha +1)}}\;e^{- \large { \frac { \beta }{ \theta }}}d \theta\\ &= \dfrac { \beta ^ \alpha }{ \Gamma ( \alpha )} \int_ { \infty }^0 \left ( \frac { \beta }u \right )^{-( \alpha +1)};e^{-u} \frac {- \beta {{u^2}\a, du \\ &= \dfrac { \beta ^ \alpha }{ \Gamma ( \alpha )} \int_0 ^{ \infty } \beta\left ( \frac {u}{ \beta } \right )^{ \alpha +1};e^{-u} \frac 1{u^2}\N-, du \\ &= \dfrac {1}{ \Gamma ( \alpha )} \int_0 ^{ \infty } u^{ \alpha -1};e^{-u}\N-,du \\ &= \dfrac { \Gamma ( \alpha )}{ \Gamma ( \alpha )}, \quad (*) \\ &=1 \end {align} $(*)$ de la definición de la $\Gamma$ función para $\,\alpha>0,\;\beta>0$ .