An Algorithm for Fast Generation of Bivariate Poisson Random Vectors

Files
Main article (1.06 MB)
Downloads: 1395
TR Number
Date
2010
Journal Title
Journal ISSN
Volume Title
Publisher
INFORMS
Abstract

We present the "trivariate reduction extension" (TREx)-an exact algorithm for the fast generation of bivariate Poisson random vectors. Like the normal-to-anything (NORTA) procedure, TREx has two phases: a preprocessing phase when the required algorithm parameters are identified, and a generation phase when the parameters identified during the preprocessing phase are used to generate the desired Poisson vector. We prove that the proposed algorithm covers the entire range of theoretically feasible correlations, and we provide efficient-computation directives and rigorous bounds for truncation error control. We demonstrate through extensive numerical tests that TREx, being a specialized algorithm for Poisson vectors, has a preprocessing phase that is uniformly a hundred to a thousand times faster than a fast implementation of NORTA. The generation phases of TREx and NORTA are comparable in speed, with that of TREx being marginally faster. All code is publicly available.

Description
Keywords
Statistics, Simulation, Random variable generation, Multivariate, Distribution, Correlation, Marginals
Citation
Shin, Kaeyoung; Pasupathy, Raghu. An Algorithm for Fast Generation of Bivariate Poisson Random Vectors. INFORMS Journal on Computing 2010 22:1, 81-92. doi: 10.1287/ijoc.1090.0332