On the gamma difference distribution

Peter J. Forrester

2023-09-25

Abstract:The gamma difference distribution is defined as the difference of two gamma distributions, with in general different shape and rate parameters. Starting with knowledge of the corresponding characteristic function, a second order linear differential equation characterisation of the probability density function is given. This is used to derive a Stein-type differential identity relating to the expectation with respect to the gamma difference distribution of a general twice differentiable function $g(x)$. Choosing $g(x) = x^k$ gives a second order recurrence for the positive integer moments, which are also shown to permit evaluations in terms of ${}_2 F_1$ hypergeometric polynomials. A hypergeometric function evaluation is given for the absolute continuous moments. Specialising the gamma difference distribution gives the variance gamma distribution. Results of the type obtained herein have previously been obtained for this distribution, allowing for comparisons to be made.

Statistics Theory,Probability

What problem does this paper attempt to address?

This paper aims to solve several key problems related to the Gamma Difference Distribution. Specifically, the author attempts to:

1. **Derive the differential equation characteristics of the probability density function of the Gamma Difference Distribution**:
- Through the known characteristic function, the author gives a second - order linear differential equation description of the probability density function \(p_{X_1 - X_2}(x)\). This differential equation is:
\[
x p''_{X_1 - X_2}(x)+(x(\beta_1-\beta_2)+(2 - \alpha_1-\alpha_2))p'_{X_1 - X_2}(x)+((\beta_1-\beta_2+\alpha_1\beta_2-\alpha_2\beta_1)-x\beta_1\beta_2)p_{X_1 - X_2}(x) = 0
\]
This equation does not hold when \(\alpha_1+\alpha_2 < 1\), because at this time \(p_{X_1 - X_2}(x)\) diverges at \(x = 0\).

2. **Establish Stein - type differential identities**:
- Based on the above differential equation, the author derives a Stein - type differential identity for the expected value of a general twice - differentiable function \(g(x)\):
\[
E\left(Xg''(X)+(X(\beta_1-\beta_2)+(\alpha_1+\alpha_2))g'(X)+(\alpha_1\beta_2-\alpha_2\beta_1 - X\beta_1\beta_2)g(X)\right)=0
\]

3. **Derive the recurrence relation of positive integer moments**:
- Starting from the Stein identity, the author obtains a third - order recurrence relation of positive integer moments \(m_k(X)\):
\[
\beta_1\beta_2m_{k + 1}(X)=\left(k(\beta_2-\beta_1)+(\alpha_1\beta_2-\alpha_2\beta_1)\right)m_k(X)+k(k - 1+\alpha_1+\alpha_2)m_{k-1}(X),\quad k\geq0
\]
where the initial condition is \(m_0(X)=1\).

4. **Give the closed - form expressions of positive integer moments and absolutely continuous moments**:
- The positive integer moment \(m_k(X)\) can be represented by the hypergeometric polynomial \({}_2F_1\):
\[
m_k(X)=\frac{1}{\beta_1^k}\left(\prod_{i = 1}^k(\alpha_1 + i-1)\right){}_2F_1\left(-k,\alpha_2;-( \alpha_1 + k - 1);-\frac{\beta_1}{\beta_2}\right)
\]
- The absolutely continuous moment \(E(|X|^{b - 1})\) can also be represented by the hypergeometric function:
\[
E(|X|^{b-1})=\frac{\beta_1^{\alpha_1}\beta_2^{\alpha_2}}{(\beta_1+\beta_2)^{\alpha_1+\alpha_2 + b - 1}}\frac{\Gamma(