Approximating the Laplace transform of the sum of dependent lognormals
Let (X1,…,Xn) be multivariate normal, with mean vector μ and covariance matrix Σ, and Sn=eX1+⋯+eXn. The Laplace transform L(θ)=Ee−θSn∝∫exp{−hθ(x)}dx is represented as L~(θ)I(θ), where L~(θ) is given in closed-form and I(θ) is the error factor (≈1). We obtain L~(θ) by replacing hθ(x) with a second order Taylor expansion around its minimiser x∗. An algorithm for calculating the asymptotic expansion of x∗ is presented, and it is shown that I(θ)→1 as θ→∞. A variety of numerical methods for evaluating I(θ) are discussed, including Monte Carlo with importance sampling and quasi-Monte Carlo. Numerical examples (including Laplace transform inversion for the density of Sn) are also given.