Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-26T14:53:51.809Z Has data issue: false hasContentIssue false

On spherical Monte Carlo simulations for multivariate normal probabilities

Published online by Cambridge University Press:  21 March 2016

Huei-Wen Teng*
Affiliation:
National Central University
Ming-Hsuan Kang*
Affiliation:
National Chiao Tung University
Cheng-Der Fuh*
Affiliation:
National Central University
*
Postal address: Graduate Institute of Statistics, National Central University, 300 Zhongda Road, Zhongli District, Taoyuan City, 32001, Taiwan, R.O.C.
∗∗∗ Postal address: Department of Applied Mathematics, National Chiao Tung University, 1001 University Road, Hsinchu, 30010, Taiwan, R.O.C.
Postal address: Graduate Institute of Statistics, National Central University, 300 Zhongda Road, Zhongli District, Taoyuan City, 32001, Taiwan, R.O.C.
Rights & Permissions [Opens in a new window]

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the ‘Save PDF’ action button.

The calculation of multivariate normal probabilities is of great importance in many statistical and economic applications. In this paper we propose a spherical Monte Carlo method with both theoretical analysis and numerical simulation. We start by writing the multivariate normal probability via an inner radial integral and an outer spherical integral using the spherical transformation. For the outer spherical integral, we apply an integration rule by randomly rotating a predetermined set of well-located points. To find the desired set, we derive an upper bound for the variance of the Monte Carlo estimator and propose a set which is related to the kissing number problem in sphere packings. For the inner radial integral, we employ the idea of antithetic variates and identify certain conditions so that variance reduction is guaranteed. Extensive Monte Carlo simulations on some probabilities confirm these claims.

Type
General Applied Probability
Copyright
Copyright © Applied Probability Trust 2015 

References

Anderson, T. W., Olkin, I. and Underhill, L. G. (1987). Generation of random orthogonal matrices. SIAM J. Sci. Statist. Comput. 8, 625629.CrossRefGoogle Scholar
Asmussen, S. and Glynn, P. W. (2007). Stochastic Simulation: Algorithms and Analysis. Springer, New York.CrossRefGoogle Scholar
Conway, J. H. and Sloane, N. J. A. (1999). Sphere Packings, Lattices and Groups, 3rd edn. Springer, New York.CrossRefGoogle Scholar
Craig, P. (2008). A new reconstruction of multivariate normal orthant probabilities. J. R. Statist. Soc. B 70, 227243.CrossRefGoogle Scholar
Davis, P. J. and Rabinowitz, P. (1984). Methods of Numerical Integration, 2nd edn. Academic Press, Orlando, FL.Google Scholar
Deák, I. (1980). Three digit accurate multiple normal probabilities. Numer. Math. 35, 369380.CrossRefGoogle Scholar
Deák, I. (2000). Subroutines for computing normal probabilities of sets—computer experiences. Ann. Operat. Res. 100, 103122.CrossRefGoogle Scholar
Diaconis, P. and Shahshahani, M. (1987). The subgroup algorithm for generating uniform random variables. Prob. Eng. Inf. Sci. 1, 1532.CrossRefGoogle Scholar
Fang, K.-T. and Wang, Y. (1994). Number-Theoretic Methods in Statistics (Monogr. Statist. Appl. Prob. 51), Chapman & Hall, London.Google Scholar
Genz, A. (1992). Numerical computation of multivariate normal probabilities. J. Comput. Graphical Statist. 1, 141-150.Google Scholar
Genz, A. (1993). Comparison of methods for the computation of multivariate normal probabilities. Comput. Sci. Statist. 25, 400405.Google Scholar
Genz, A. and Bretz, F. (2002). Methods for the computation of multivariate t probabilities. J. Comput. Graphical Statist. 11, 950971.CrossRefGoogle Scholar
Genz, A. and Bretz, F. (2009). Computation of Multivariate Normal and t Probabilities (Lecture Notes Statist. 195), Springer, Dordrecht.CrossRefGoogle Scholar
Glasserman, P., Heidelberger, P. and Shahabuddin, P. (2000). Variance reduction techniques for estimating value-at-risk. Manag. Sci. 46, 13491364.CrossRefGoogle Scholar
Gupton, G. M., Finger, C. C. and Bhatia, M. (1997). CreditMetrics – Technical Document. Morgan, New York.Google Scholar
Hajivassiliou, V., McFadden, D. and Ruud, P. (1996). Simulation of multivariate normal rectangle probabilities and their derivatives: theoretical and computational results. J. Econometrics 72, 85134.Google Scholar
Heiberger, R. M. (1978). Algorithm AS 127: generation of random orthogonal matrices. J. R. Statist. Soc. C 27, 199206.Google Scholar
Hsu, J. C. (1996). Multiple Comparisons. Chapman & Hall, London.Google Scholar
Miwa, T., Hayter, A. J. and Kuriki, S. (2003). The evaluation of general non-centred orthant probabilities. J. R. Statist. Soc. B 65, 223234.CrossRefGoogle Scholar
Monahan, J. and Genz, A. (1997). Spherical-radial integration rules for Bayesian computation. J. Amer. Statist. Assoc. 92, 664674.CrossRefGoogle Scholar
Nebe, G. and Sloane, N. J. A. A catalogue of lattices. Available at http://www.math.rwth-aachen.de/∼Gabriele.Nebe/LATTICES/.Google Scholar
Ross, S. M. (2013). Simulation, 5th edn. Elsevier, Amsterdam.Google Scholar
Sándor, Z. and András, P. (2004). Alternative sampling methods for estimating multivariate normal probabilities.break J. Econometrics 120, 207234.CrossRefGoogle Scholar
Sloane, N. J. A. Spherical Codes. Available at http://neilsloane.com/packings/index.html#I.Google Scholar
Somerville, P. N. (2001). Numerical computation of multivariate normal and multivariate t probabilities over ellipsoidal regions. J. Statist. Software 6, 10pp.Google Scholar
Stewart, G. W. (1980). The efficient generation of random orthogonal matrices with an application to condition estimators. SIAM J. Numerical Anal. 17, 403409.CrossRefGoogle Scholar
Vijverberg, W. P. M. (1997). Monte Carlo evaluation of multivariate normal probabilities. J. Econometrics 76, 281-307.Google Scholar