## 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 |

## Abstract

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