Generalized multivariate log-gamma distribution explained

In probability theory and statistics, the generalized multivariate log-gamma (G-MVLG) distribution is a multivariate distribution introduced by Demirhan and Hamurkaroglu^[1] in 2011. The G-MVLG is a flexible distribution. Skewness and kurtosis are well controlled by the parameters of the distribution. This enables one to control dispersion of the distribution. Because of this property, the distribution is effectively used as a joint prior distribution in Bayesian analysis, especially when the likelihood is not from the location-scale family of distributions such as normal distribution.

Joint probability density function

If

\boldsymbol{Y} \sim G - MVLG (\delta,\nu,\boldsymbol{ λ },\boldsymbol{\mu}) , the joint probability density function (pdf) of \boldsymbol{Y} = (Y _{1} , ... ,Y _{k} ) is given as the following:f(y _{1,...,y} _{k)=} \delta ^{\nu} \sum ^{inf ty} _{n =0} (( 1- \delta) ^{n} \prod ^{k} _{i =1} \mu _{i} λ ^{- \nu - n} _{i})/([\Gamma(\nu + n)] ^{k-1} \Gamma(\nu)n!) \exp \ {(\nu + n)\sum ^{k} _{i =1} \mu _{i} y _{i} - \sum ^{k} _{i =1} (1)/(λ _{i}) \exp\{\mu _{i} y _{i\}\},}

where

\boldsymbol{y}\in R ^{k} , \nu> 0 , λ _{j} > 0 , \mu _{j} > 0 for j =1 , ... ,k, \delta = \det(\boldsymbol{\Omega}) ^{(1)/(k -1)} , and\boldsymbol{\Omega} = \left(\begin{array}{cccc} 1 & \sqrt{ abs (\rho _{12} )} & … & \sqrt{ abs (\rho _{1k} )} \\ \sqrt{ abs (\rho _{12} )} & 1 & … & \sqrt{ abs (\rho _{2k} )} \\ \vdots & \vdots & \ddots & \vdots \\ \sqrt{ abs (\rho _{1k} )} & \sqrt{ abs (\rho _{2k} )} & … & 1 \end{array} \right), \rho _{ij} is the correlation between Y _{i} and Y _{j} , \det( ⋅ ) and abs ( ⋅ ) denote determinant and absolute value of inner expression, respectively, and \boldsymbol{g} = (\delta,\nu,\boldsymbol{ λ } ^{T,\boldsymbol{\mu}} ^{T)} includes parameters of the distribution.

Properties

Joint moment generating function

The joint moment generating function of G-MVLG distribution is as the following:

M _{\boldsymbol{Y} }(\boldsymbol) =\delta^\nu \bigg(\prod_^k\lambda_i^\bigg)\sum_^\infty \frac(1-\delta)^n \prod_^k \frac.

Marginal central moments

r ^{th} marginal central moment of Y _{i} is as the following:{\mu _{i}'} ^{r = \left[ (( λ ^{t _{i/\mu} _{i}} _{i/\delta)})/(\Gamma(\nu))} \sum ^{r} _{k =0} \binom{r}{k}\left[ (ln ( λ _{i/\delta)})/(\mu _{i}) \right] ^{r-k} (\partial ^{k\Gamma(\nu + t} _{i/\mu} _{i)})/(\partial t ^{k} _{i}) \right] _{t _{i=0}} .

Marginal expected value and variance

Marginal expected value

Y _{i} is as the following:\operatorname{E}(Y _{i} ) = (1)/(\mu _{i}) [ ln ( λ _{i/\delta)+\digamma(\nu)],} \operatorname{var}(Z ^{[ 1 ]} _{i) = \digamma} (\nu)/(\mu ^{2} _{i)}

where

\digamma(\nu) and \digamma ^{[1]} (\nu) are values of digamma and trigamma functions at \nu , respectively.

Related distributions

Demirhan and Hamurkaroglu establish a relation between the G-MVLG distribution and the Gumbel distribution (type I extreme value distribution) and gives a multivariate form of the Gumbel distribution, namely the generalized multivariate Gumbel (G-MVGB) distribution. The joint probability density function of

\boldsymbol{T}\sim G - MVGB (\delta,\nu,\boldsymbol{ λ },\boldsymbol{\mu}) is the following:f(t _{1,...,t} _{k;} \delta,\nu,\boldsymbol{ λ },\boldsymbol{\mu})) = \delta ^{\nu} \sum ^{inf ty} _{n =0} (( 1- \delta) ^{k} _{i =1} \mu _{i} λ ^{- \nu - n} _{i})/([\Gamma(\nu + n)] ^{k-1} \Gamma(\nu)n!) \exp \ { - (\nu + n)\sum ^{k} _{i =1} \mu _{i} t _{i} - \sum ^{k} _{i =1} (1)/(λ _{i}) \exp\{ - \mu _{i} t _{i\}\},} t _{i\in} R .

The Gumbel distribution has a broad range of applications in the field of risk analysis. Therefore, the G-MVGB distribution should be beneficial when it is applied to these types of problems..

Notes and References

1. Demirhan. Haydar. Hamurkaroglu, Canan. On a multivariate log-gamma distribution and the use of the distribution in the Bayesian analysis. Journal of Statistical Planning and Inference. 2011. 141. 1141–1152. 10.1016/j.jspi.2010.09.015. 3.