Skip to main content
2 of 3
added 255 characters in body; deleted 3 characters in body
Peter
  • 21
  • 4

How to approx. decompose a sym. p.d. matrix M into X'X?

M: pxp symmetric p.d. matrix with unit diagonals n: number much smaller than p

Want a nonrandom nxp matrix X such that X'X is close to M element-wise. If n gets larger, hopefully difference gets smaller.

I have a no so good method. Get the eigendecomposition M = D'VD, and rearrange the diagonals of V in decreasing order and so the rows of D. Take only the first n rows of D as X. This method works only if the eigenvalues decrease very fast.

Any input is appreciated. Thanks a lot!

Peter
  • 21
  • 4