Eigendecomposition
Posted: Wed Nov 26, 2014 4:15 pm
I am trying to estimate an eigenvalue decomposition for Matrix A using the makepcomp command. However the eigenvectors (say vector M) and the eigenvectors (say matrix U) estimated do not give back the matrix A when I am trying to estimate UxMx@inverse(U). The commands I use are the following:
A.makepcomp(cov=cov, cpnorm, eigval=M, eigvec=U)
Matrix X=Ux@makediagonal(M)x@inverse(U)
A.makepcomp(cov=cov, cpnorm, eigval=M, eigvec=U)
Matrix X=Ux@makediagonal(M)x@inverse(U)