doi:10.1016/j.ejor.2003.12.027
Copyright © 2004 Published by Elsevier B.V.
Stochastics and Statistics
Using the sum-of-uniforms method to generate correlated random variates with certain marginal distribution
Jung-Tai Chen
, 
Department of Asia-Pacific Industrial and Business Management, National University of Kaohsiung, 700 Kaohsiung University Road, Nantz District, Kaohsiung 811, Taiwan, ROC
Received 3 May 2003;
accepted 11 December 2003.
Available online 24 June 2004.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
This paper extends the sum-of-uniforms method to generate correlated random variables with certain marginal distributions. We first use the transformation method to derive the joint probability density function of the correlated uniform random variables. We also demonstrate that the sum-of-uniforms method can be extended to generate correlated random variables with certain marginal distributions including uniform, exponential, Erlang, Bernoulli, binomial, geometric, and negative binomial. Finally, this paper presents the exact correlation coefficients of such correlated random variables.
Author Keywords: Simulation; Multivariate generation; Correlated random variates
Fig. 1. The domain of (Un−1,Wn) for c
1 with positive correlation.
Fig. 2. The domain of (Un−1,Un): (a) c
1 with positive correlation; (b) c<1 with positive correlation; (c) c
1 with negative correlation; (d) c<1 with negative correlation.
Fig. 3. The joint pdf of (Un−1,Un): (a) c=1.5 with positive correlation; (b) c=0.5 with positive correlation; (c) c=1.5 with negative correlation; (d) c=0.5 with negative correlation.
Fig. 4. The scatter diagram of (Un−1,Un), based on 10,000 pairs of observations: (a) c=1.5 with positive correlation; (b) c=0.5 with positive correlation; (c) c=1.5 with negative correlation; (d) c=0.5 with negative correlation.
Fig. 5. The joint pdf of (Yn−1,Yn) with μ=1: (a) c=1.5 with positive correlation; (b) c=0.5 with positive correlation; (c) c=1.5 with negative correlation; (d) c=0.5 with negative correlation.
Fig. 6. The scatter diagram of (Yn−1,Yn) with μ=1, based on 10,000 pairs of observations: (a) c=1.5 with positive correlation; (b) c=0.5 with positive correlation; (c) c=1.5 with negative correlation; (d) c=0.5 with negative correlation.
Table 1. The exact value of ρ(c) for correlated uniform random variables

Table 2. Numerical values of ρ(c) for correlated exponential random variables

Table 3. Numerical values of ρ(c,p) for positively correlated Bernoulli random variables

Table 4. Numerical values of ρ(c,p) for negatively correlated Bernoulli random variables

LB denotes the lower bound.
Table 5. Numerical values of ρ(c,p) for positively correlated geometric random variables

Table 6. Numerical values of ρ(c,p) for negatively correlated geometric random variables

LB denotes the lower bound.