Excellent analysis, Thomas. An alternative to looking at EISPACK
documentation is to do the following computations:
(1) Values %*% diag(vectors) %*% solve(Values)
(2) solve(Values) %*% diag(vectors) %*% Values
One of these two should return the original matrix; the other will
likely be very
Mikael -
The matrix PA9900 is not a symmetric matrix. Eigen() will
automatically detect this. help("eigen") says explicitly:
"For `eigen( , symmetric = FALSE)' the choice of length of
the eigenvectors is not defined by LINPACK. In all other
cases the vectors are normalized to unit length."
I
Dear R-listers
Is there anyone who knows why I get different eigenvectors when I run
MatLab and R? I run both programs in Windows Me. Can I make R to produce
the same vectors as MatLab?
#R Matrix
PA9900<-c(11/24 ,10/53 ,0/1 ,0/1 ,29/43 ,1/24 ,27/53 ,0/1 ,0/1 ,13/43
,14/24 ,178/53 ,146/244 ,17/23
