Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down on April 23, and its replacement, is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: RE: matrix decomposition

From   "Feiveson, Alan H. (JSC-SK311)" <>
To   "" <>
Subject   st: RE: matrix decomposition
Date   Tue, 10 Aug 2010 08:18:26 -0500

David if V = A*A' where A is the Cholesky decomposition, you can get any other decomposition of V by multiplying A by an arbitrary orthogonal matrix, say P. Thus

(AP)*(AP)' = APP'A' = AIA' = AA' = V

since by definition, PP' = I

Of course, given A, you would have to figure out how construct P to do what you want and still be orthogonal.

Al Feiveson

-----Original Message-----
From: [] On Behalf Of
Sent: Tuesday, August 10, 2010 5:46 AM
Subject: st: matrix decomposition

Dear all,

I want to get Stata to decompose a variance covariance (pos definite)
matrix, say V, into A*A', where A is not triangular as in the cholesky
decomposition (by the command cholesky(V)), but has other linear
restrictions e.g. zeros in different places. Anyone have any tips/examples
of how to program/command this?

Many thanks for any help,
*   For searches and help try:

*   For searches and help try:

© Copyright 1996–2015 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index