Index: The Book of Statistical ProofsProbability Distributions ▷ Univariate continuous distributions ▷ Gamma distribution ▷ Quantile function

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:

Metadata: ID: P194 | shortcut: gam-qf | author: JoramSoch | date: 2020-11-19, 07:31.