A simple generalization of the Box–Muller method for obtaining a pair of correlated standard normal variables

Abstract

In the simple and widely used method of Box–Muller [G. Box and M. Muller, A note on the generation of random normal deviates, Ann. Math. Statist. 29 (1958), pp. 610–611], from a pair of uniform and independent random variables in (0,1), a pair of standard and independent normal variables is obtained. In this article, we present a very simple and elegant generalization of this method to obtain a pair of correlated standard normal variables with a given coefficient of correlation. This generalized method, which is computationally very easy, is interpreted in geometric terms, considering a translation of the uniform interval (0,1) and a rotation of a defined angle, both related to the coefficient of correlation. Some numerical results are simulated and statistically analysed, proving that the generalization is extremely simple and powerful.