Proof: Quantile function of the gamma distribution
Index:
The Book of Statistical Proofs ▷
Probability Distributions ▷
Univariate continuous distributions ▷
Gamma distribution ▷
Quantile function
Metadata: ID: P194 | shortcut: gam-qf | author: JoramSoch | date: 2020-11-19, 07:31.
Theorem: Let $X$ be a random variable following a gamma distribution:
\[\label{eq:gam} X \sim \mathrm{Gam}(a,b) \; .\]Then, the quantile function of $X$ is
\[\label{eq:gam-qf} Q_X(p) = \left\{ \begin{array}{rl} -\infty \; , & \text{if} \; p = 0 \\ \gamma^{-1}(a, \Gamma(a) \cdot p)/b \; , & \text{if} \; p > 0 \end{array} \right.\]where $\gamma^{-1}(s, y)$ is the inverse of the lower incomplete gamma function $\gamma(s, x)$
Proof: The cumulative distribution function of the gamma distribution is:
\[\label{eq:gam-cdf} F_X(x) = \left\{ \begin{array}{rl} 0 \; , & \text{if} \; x < 0 \\ \frac{\gamma(a,bx)}{\Gamma(a)} \; , & \text{if} \; x \geq 0 \; . \end{array} \right.\]The quantile function $Q_X(p)$ is defined as the smallest $x$, such that $F_X(x) = p$:
\[\label{eq:qf} Q_X(p) = \min \left\lbrace x \in \mathbb{R} \, \vert \, F_X(x) = p \right\rbrace \; .\]Thus, we have $Q_X(p) = -\infty$, if $p = 0$. When $p > 0$, it holds that
\[\label{eq:gam-qf-s1} Q_X(p) = F_X^{-1}(x) \; .\]This can be derived by rearranging equation \eqref{eq:gam-cdf}:
\[\label{eq:gam-qf-s2} \begin{split} p &= \frac{\gamma(a,bx)}{\Gamma(a)} \\ \Gamma(a) \cdot p &= \gamma(a,bx) \\ \gamma^{-1}(a, \Gamma(a) \cdot p) &= bx \\ x &= \frac{\gamma^{-1}(a, \Gamma(a) \cdot p)}{b} \; . \end{split}\]∎
Sources: - Wikipedia (2020): "Incomplete gamma function"; in: Wikipedia, the free encyclopedia, retrieved on 2020-11-19; URL: https://en.wikipedia.org/wiki/Incomplete_gamma_function#Definition.
Metadata: ID: P194 | shortcut: gam-qf | author: JoramSoch | date: 2020-11-19, 07:31.