I understood Mitra to want a matrix with
the eigenvalues along the diagonal. That's
what was asked for. It's true, as I understand
it, that diagonalization as such is a different
thing.
Of course, if the eigenvalues have imaginary
parts it is rather difficult to do anything
with that matrix in Stata, but I didn't think it
necessary to spell that out. Either way, -geneigen-
would seem to be the first step in Stata, even
if one can't go beyond it.
Nick
[email protected]
> -----Original Message-----
> Clyde Schechter
>
> Yesterday, Mitra Pinaki asked for help in diagonalizing a
> non-symmetric
> matrix. Nick Cox pointed him to Kit Baum's -geneigen-
> which will, indeed,
> calculate the eigenvalues of a non-symmetric matrix. But that won't
> diagonalize the matrix. In fact, it's fairly simple linear
> algebra to show
> that any matrix of real numbers which _can_ be diagonalized
> is necessarily
> symmetric! So it's unclear what M.P. is trying to do.
>
> The most that can be done with a non-symmetric matrix is to
> identify its
> eigenvalues and eigenvectors (which, due to non-symmetry, will have
> cardinality strictly less than the rank of the matrix) and
> then create a
> diagonal matrix out of those. But that diagonal matrix
> will be of smaller
> rank than the original, and is not a conjugate of the
> original matrix. Any
> conjugate of the original which does diagonalize as much as
> possible will
> still have a piece which has non-zero elements off the diagonal.
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/