You are here: MIMS > EPrints
MIMS EPrints

2006.146: Numerically stable generation of correlation matrices and their factors

2006.146: Philip I. Davies and Nicholas J. Higham (2000) Numerically stable generation of correlation matrices and their factors. BIT Numerical Mathematics, 40 (4). pp. 640-651. ISSN 1572-9125

Full text available as:

PDF - Archive staff only - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
217 Kb

DOI: 10.1023/A:1022384216930


Correlation matrices—symmetric positive semidefinite matrices with unit diagonal—are important in statistics and in numerical linear algebra. For simulation and testing it is desirable to be able to generate random correlation matrices with specified eigenvalues (which must be nonnegative and sum to the dimension of the matrix). A popular algorithm of Bendel and Mickey takes a matrix having the specified eigenvalues and uses a finite sequence of Givens rotations to introduce 1s on the diagonal. We give improved formulae for computing the rotations and prove that the resulting algorithm is numerically stable. We show by example that the formulae originally proposed, which are used in certain existing Fortran implementations, can lead to serious instability. We also show how to modify the algorithm to generate a rectangular matrix with columns of unit 2-norm. Such a matrix represents a correlation matrix in factored form, which can be preferable to representing the matrix itself, for example when the correlation matrix is nearly singular to working precision.

Item Type:Article
Uncontrolled Keywords:Random correlation matrix - Bendel-Mickey algorithm - eigenvalues - singular value decomposition - test matrices - forward error bounds - relative error bounds - IMSL - NAG Library - Jacobi method
Subjects:MSC 2000 > 15 Linear and multilinear algebra; matrix theory
MSC 2000 > 65 Numerical analysis
MIMS number:2006.146
Deposited By:Miss Louise Stait
Deposited On:27 June 2006

Download Statistics: last 4 weeks
Repository Staff Only: edit this item