Behavior of the NORTA method for correlated random vector generation as the dimension increases

The NORTA method is a fast general-purpose method for generating samples of a random vector with given marginal distributions and given correlation matrix. It is known that there exist marginal distributions and correlation matrices that the NORTA method cannot match, even though a random vector with the prescribed qualities exists. We investigate this problem as the dimension of the random vector increases. Simulation results show that the problem rapidly becomes acute, in the sense that NORTA fails to work with an increasingly large proportion of correlation matrices. Simulation results also show that if one is willing to settle for a correlation matrix that is "close" to the desired one, then NORTA performs well with increasing dimension. As part of our analysis, we develop a method for sampling correlation matrices uniformly (in a certain precise sense) from the set of all such matrices. This procedure can be used more generally for sampling uniformly from the space of all symmetric positive definite matrices with diagonal elements fixed at positive values.

[1]  R. Iman,et al.  A distribution-free approach to inducing rank correlation among input variables , 1982 .

[2]  G HendersonShane,et al.  Behavior of the NORTA method for correlated random vector generation as the dimension increases , 2003 .

[3]  R. Clemen,et al.  Correlations and Copulas for Decision and Risk Analysis , 1999 .

[4]  R. M. Cooke,et al.  CONDITIONAL , PARTIAL AND RANK CORRELATION FOR THE ELLIPTICAL COPULA ; DEPENDENCE MODELLING IN UNCERTAINTY ANALYSIS , 2001 .

[5]  Louis Guttman,et al.  Enlargement Methods for Computing the Inverse Matrix , 1946 .

[6]  L. Devroye Non-Uniform Random Variate Generation , 1986 .

[7]  Shane G. Henderson,et al.  Chessboard Distributions and Random Vectors with Specified Marginals and Covariance Matrix , 2002, Oper. Res..

[8]  Diane Valérie Ouellette Schur complements and statistics , 1981 .

[9]  Mark E. Johnson Multivariate Statistical Simulation: Johnson/Multivariate , 1987 .

[10]  C. Richard Johnson,et al.  Matrix Completion Problems: A Survey , 1990 .

[11]  H. Wolkowicz,et al.  Matrix Completion Problems , 2000 .

[12]  W. Whitt Bivariate Distributions with Given Marginals , 1976 .

[13]  Joseph L. Hammond,et al.  Generation of Pseudorandom Numbers with Specified Univariate Distributions and Correlation Coefficients , 1975, IEEE Transactions on Systems, Man, and Cybernetics.

[14]  Philip M. Lurie,et al.  An Approximate Method for Sampling Correlated Random Variables From Partially-Specified Distributions , 1998 .

[15]  Bruce W. Schmeiser,et al.  Bivariate Gamma Random Vectors , 1982, Oper. Res..

[16]  W. Kruskal Ordinal Measures of Association , 1958 .

[17]  Roger M. Cooke,et al.  Generating "dependent" quasi-random numbers , 2000, 2000 Winter Simulation Conference Proceedings (Cat. No.00CH37165).

[18]  Averill M. Law,et al.  Simulation Modeling and Analysis , 1982 .

[19]  I. Olkin,et al.  Generating Correlation Matrices , 1984 .

[20]  A. Edelman Eigenvalues and condition numbers of random matrices , 1988 .

[21]  Shane G. Henderson,et al.  Chessboard distributions , 2001, Proceeding of the 2001 Winter Simulation Conference (Cat. No.01CH37304).

[22]  R. Saigal,et al.  Handbook of semidefinite programming : theory, algorithms, and applications , 2000 .

[23]  H. Gifford Fong,et al.  Impact of Correlated Default Risk on Credit Portfolios , 2001 .

[24]  Shane G. Henderson,et al.  Mathematics for simulation , 2000, Proceeding of the 2001 Winter Simulation Conference (Cat. No.01CH37304).

[25]  H. Joe Multivariate models and dependence concepts , 1998 .

[26]  Satishs Iyengar,et al.  Multivariate Models and Dependence Concepts , 1998 .

[27]  Charles H. Reilly,et al.  The Effects of Coefficient Correlation Structure in Two-Dimensional Knapsack Problems on Solution Procedure Performance , 2000 .

[28]  Henry Wolkowicz,et al.  Handbook of Semidefinite Programming , 2000 .

[29]  Maurice G. Kendall,et al.  A Course in the Geometry of n Dimensions , 1962 .

[30]  Shane G. Henderson,et al.  Properties of the NORTA method in higher dimensions , 2002, Proceedings of the Winter Simulation Conference.

[31]  Shane G. Henderson,et al.  Simulation mathematics and random number generation: mathematics for simulation , 2001, WSC '01.