[R] Spectral Decomposition

2007-06-29 Thread Doran, Harold
All of my resources for numerical analysis show that the spectral decomposition is A = CBC' Where C are the eigenvectors and B is a diagonal matrix of eigen values. Now, using the eigen function in R # Original matrix aa - matrix(c(1,-1,-1,1), ncol=2) ss - eigen(aa) # This results yields back

Re: [R] Spectral Decomposition

2007-06-29 Thread Richard M. Heiberger
For a general square matrix A, the eigenvalue decomposition is A = CBC^{-1} For the special case of symmetric A, C^{-1} = C' __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Spectral Decomposition [OOPS]

2007-06-29 Thread Ted Harding
On 29-Jun-07 13:23:05, Ted Harding wrote: [Sorry -- a silly typo in my previous]: If A is not symmetric, you have left eigenvectors: x'*A = lambda*x' and right eigenvectors: A*x = lambda*x and the left eigenvectors are not the same as the right eigenvectors, though you have the

Re: [R] Spectral Decomposition

2007-06-29 Thread Ted Harding
On 29-Jun-07 12:29:31, Doran, Harold wrote: All of my resources for numerical analysis show that the spectral decomposition is A = CBC' Where C are the eigenvectors and B is a diagonal matrix of eigen values. Now, using the eigen function in R # Original matrix aa -