# Inverse gamma distribution definition

Wikipedia says the pdf for the gamma function is:

$X \sim \operatorname{Gamma}(\alpha,\beta) \implies \Pr(X=x) \propto x^{\alpha-1}e^{-\beta x}$

If $Y = 1/X$, then

$\Pr(Y=y) = \Pr(X=1/y) \propto (1\mathbin/y)^{\alpha-1}e^{-\beta/y} = y^{-(\alpha-1)}e^{-\beta\mathbin/y} = y^{-\alpha+1}e^{-\beta\mathbin/y}$

However, wikipedia says of the inverse gamma distribution that

$\Pr(Y=y) \propto y^{-\alpha\mathbf{-1}}e^{-\beta/y}$

Eg, -1, not +1.

Is this just a parameterization difference, or a mistake in wikipedia, or a mistake in my derivation?

You have made a mistake assuming that you can derive the pdf of $1/X$ the way you would do it for the discrete case. But remember this is the continuous case so solving for the variable and substituting back is not enough. You still have to multiplty by the Jacobian. And also some people might take issue with the equality sign within the probability as in the continuous setting the probability of any value is zero. But you can say something like
$$P\left( y<Y <y +\Delta y\right) = f_Y(y)$$
for sufficiently small $\Delta y$.
• Ah, of course. Should have been obvious had I started with $\Pr(Y \le y) = \Pr(X \ge 1\mathbin/y)$. – Alec Dec 28 '15 at 23:36