Dagum distribution

Definition

The cumulative distribution function of the Dagum distribution (Type I) is given by :F(x;a,b,p)= \left( 1+\left(\frac{x}{b}\right)^{-a} \right)^{-p} \text{ for } x 0 \text{ where } a, b, p 0 .

The corresponding probability density function is given by :f(x;a,b,p)= \frac{ap}{x} \left( \frac{(\tfrac{x}{b})^{a p}}{\left((\tfrac{x}{b})^a + 1 \right)^{p+1}} \right) .

The quantile function is given by

:Q(u;a,b,p)= b(u^{-1/p}-1)^{-1/a}=\frac{bu^{\frac{1}{pa}}}{\left(1-u^{\frac{1}{p}}\right)^\frac{1}{a}}

The Dagum distribution can be derived as a special case of the generalized Beta II (GB2) distribution (a generalization of the Beta prime distribution):

: X \sim D(a,b,p) \iff X \sim GB2(a,b,p, 1)

There is also an intimate relationship between the Dagum and Singh–Maddala / Burr distribution. : X \sim D(a,b,p) \iff \frac{1}{X} \sim SM(a,\tfrac{1}{b},p)

The cumulative distribution function of the Dagum (Type II) distribution adds a point mass at the origin and then follows a Dagum (Type I) distribution over the rest of the support (i.e. over the positive halfline) :F(x;a,b,p,\delta)= \delta + (1-\delta) \left( 1+\left(\frac{x}{b}\right)^{-a} \right)^{-p}.

The quantile function of the Dagum (Type II) distribution is :Q(u;a,b,p,\delta)= b \left[\left(\frac{u-\delta}{1-\delta}\right)^{-\frac{1}{p}}-1\right]^{-\frac{1}{a}}.

Use in economics

The Dagum distribution is often used to model income and wealth distribution. The relation between the Dagum Type I and the Gini coefficient is summarized in the formula below:

:G=\frac{\Gamma (p) \Gamma (2p+1/a)}{\Gamma (2p) \Gamma (p+1/a)}-1 ,

where \Gamma(\cdot) is the gamma function. Note that this value is independent from the scale-parameter, b. The Gini coefficient for a continuous probability distribution takes the form:

:G = {1\over{\mu}}\int_{0}^{\infty}F(1-F)dx

where F is the CDF of the distribution and \mu is the expected value. For the Dagum distribution, assuming a 1, the formula for the Gini coefficient becomes:

:G = \frac{\Gamma(p)}{b \Gamma(1-1/a)\Gamma(p+1/a)} \int_{0}^{\infty} \left[ 1 + \left( \frac{x}{b} \right)^{-a} \right]^{-p} \left{ 1 - \left[ 1 + \left( \frac{x}{b} \right)^{-a} \right]^{-p} \right} dx

Defining the substitution u = (x/b)^{-a} leads to the simpler equation:

:G = \frac{\Gamma(p)}{a \Gamma(1-1/a)\Gamma(p+1/a)} \int_{0}^{\infty} u^{-1/a-1}\left( 1 + u \right)^{-p} \left[ 1 - \left( 1 + u \right)^{-p} \right] du

And making the substitution z = u/(1+u) further changes the Gini coefficient formula to:

:G = \frac{\Gamma(p)}{a \Gamma(1-1/a)\Gamma(p+1/a)} \int_{0}^{1} z^{-1/a-1}\left( 1 - z \right)^{p + 1/a - 1} \left[ 1 - \left( 1 - z \right)^{p} \right] dz

Each of the integral components is equivalent to the standard beta function, which may be written in terms of the gamma function:

:\text{B}(x,y) = {\Gamma(x)\Gamma(y)\over{\Gamma(x+y)}}

Then using this definition, and the fact that a\Gamma(1-1/a) = -\Gamma(-1/a) , the Gini coefficient may be written as:

:\begin{aligned} G &= \frac{\Gamma(p)}{ \Gamma(-1/a)\Gamma(p+1/a)} \left[ \frac{\Gamma(-1/a)\Gamma(2p+1/a)}{\Gamma(2p)} - \frac{\Gamma(-1/a)\Gamma(p+1/a)}{\Gamma(p)} \right] \ &= \frac{\Gamma(p)\Gamma(2p+1/a)}{\Gamma(2p)\Gamma(p+1/a)} - 1 \end{aligned}

with the last term being the Gini coefficient for the Dagum distribution. Although the Dagum distribution is not the only three-parameter distribution used to model income distribution, one study found it to usually be a better fit than other three-parameter models.

The Dagum distribution has been extended to model net wealth distribution, accounting for the observed frequencies of negative and null net wealth. This generalized model, known as the Dagum General Model of Net Wealth Distribution, is a mixture model consisting of an atomic distribution at zero (representing economic units with no wealth) with two continuous distributions for negative and positive net wealth.

References

References

1. (2018). "Using the GB2 Income Distribution". Econometrics.
2. Dagum, Camilo. (1975). "A model of income distribution and the conditions of existence of moments of finite order". [[Bulletin of the International Statistical Institute]].
3. Dagum, Camilo. (1977). "A new model of personal income distribution: Specification and estimation". Économie Appliquée.
4. Kleiber, Christian. (2008). "Modeling Income Distributions and Lorenz Curves". Springer.
5. (2003). "Statistical Size Distributions in Economics and Actuarial Sciences". Wiley.
6. Kleiber, Christian. (2007). "A Guide to the Dagum Distributions". Working Paper.
7. (2002). "A Comparison of Parametric Models of Income Distribution Across Countries and Over Time". Luxembourg Income Study Working Paper No. 305.
8. Dagum, Camilo. (2006). "Wealth distribution models: analysis and applications". Statistica.