Equivalent ways of parametrizing Gamma distribution From scipy documentation at https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.gamma.html, Gamma distribution is written as
$$f(x, \alpha) = \frac{x^{\alpha - 1} e^{-x}}{\Gamma(\alpha)}$$
The doc also says this is equivalent to the more common way of parametrizing Gamma distribution
$$f(x, \alpha, \beta) = \frac{\beta^\alpha x^{\alpha - 1} e^{-\beta x}}{\Gamma(\alpha)}$$
but with a scale of $\frac{1}{\beta}$.
Can anyone provide more details to show how the two are equivalent?
I'm not sure if this is the right way to do it, but if I substitute $x = \beta y$ into the first equation, it seems there would be a factor $\beta$ missing compared to the 2nd equation.

A change of variables in the density requires more than substitution.   In particular you need to multiply by the absolute value of the derivative of the inverse function.  This would be more obvious if you considered the cumulative distribution function and then differentiated to get the density: the extra multiplicative factor would come through the chain rule.
If $g(x)$ is strictly increasing, you could consider $$f_Y(y)=\tfrac{d}{dy} F_Y(y) = \tfrac{d}{dy} F_X\big(g^{-1}(y)\big) = f_X\big(g^{-1}(y)\big) \tfrac{d}{dy} \big(g^{-1}(y)\big)$$ and something similar if $g(x)$ was strictly decreasing.  Combining these two results, a typical statement is that if you have $f_X(x)$ and want to consider the density of $Y=g(X)$, then $$f_Y(y) =f_X\big(g^{-1}(y)\big)  \left| \tfrac{d}{dy} \big(g^{-1}(y)\big) \right|$$ though it gets more complicated when $g(x)$ is not a bijection
In your example $g(x)=\frac x\beta$ so $g^{-1}(y)=\beta y$ and thus you need to multiply by $\left| \tfrac{d}{dy} \big(g^{-1}(y)\big) \right| = \beta$
