Simulations of multivariate gamma distributions and multifactor gamma distributions
Philippe Bernardoff, Bénédicte Puig
TL;DR
This work extends simulations from bivariate to multivariate and multifactor gamma distributions by developing a general mgd framework defined via $\mathbb{E}\{\exp(-\boldsymbol{\theta}^{\top}\mathbf{X})\}=[P(\boldsymbol{\theta})]^{-\lambda}$ with an affine $P_n$, and by constructing $\mathbf{X}=\mathbf{Y}+\mathbf{Z}$ to enable practical simulation. It derives a general joint distribution form $\boldsymbol{\gamma}_{(P_n,\lambda)}$, marginal gamma laws, and multifactor extensions with conditional Laplace transforms, including explicit infinite-divisibility criteria and special functions (e.g., Lauricella/Horn) for densities in low dimensions. The paper then provides concrete conditional-Laplace-transform expressions and detailed simulation algorithms for key cases ($n=2,3,4$) in the $k=1$ setting, and discusses a Markovian simplification that yields an efficient sequential simulation scheme for general $n$. It culminates with extensive simulations in dimensions 2–5 illustrating bgd, tgd, quadrivariate gamma, and Mmgd constructions, validating the methods and clarifying parameter choices and correlation structures. The results significantly extend the toolbox for simulating high-dimensional infinitely divisible gamma-type models, with potential applications in multivariate risk, classification within exponential families, and stochastic modeling of dependent gamma processes.
Abstract
This article provides a general expression for infinitely divisible multivariate gamma distributions defined by their Laplace transforms, as well as the conditional Laplace transform of infinitely divisible multivariate gamma distributions.We give algorithms for simulating infinitely divisible gamma distributions and infinitely divisible multifactor gamma distributions in dimension 2,3,4 and for all dimensions greater than 2 in the Markovian case. We give examples of simulations in dimension 2,3,4 and in dimension 5 in the Markovian case.
