12

It is well-known that a covariance matrix must be semi-positive definite, however, is the converse true?

That is, does every semi-positive definite matrix correspond to a covariance matrix?

Glen_b
  • 257,508
  • 32
  • 553
  • 939
Jingjings
  • 1,173
  • 1
  • 9
  • 13

1 Answers1

12

Going by the definitions of PD and PSD here, yes, I think so, since we can do this by construction. I'll assume for a slightly simpler argument that you mean for matrices with real elements, but with appropriate changes it would extend to complex matrices.

Let $A$ be some real PSD matrix; from the definition I linked to, it will be symmetric. Any real symmetric positive definite matrix $A$ can be written as $A = LL^T$. This can be done by $L=Q\sqrt{D}Q^T$ if $A=QDQ^T$ with orthogonal $Q$ and diagonal $D$ and $\sqrt{D}$ as matrix of component wise square roots of $D$. Thus, it needn't be full rank.

Let $Z$ be some vector random variable, of the appropriate dimension, with covariance matrix $I$ (which is easy to create).

Then $LZ$ has covariance matrix $A$.

[At least that's in theory. In practice there'd be various numerical issues to deal with if you wanted good results, and - because of the usual problems with floating point calculation - you'd only approximately get what you need; that is, the population variance of a computed $LZ$ usually wouldn't be exactly $A$. But this sort of thing is always an issue when we come to actually calculate things]

amoeba
  • 93,463
  • 28
  • 275
  • 317
Glen_b
  • 257,508
  • 32
  • 553
  • 939
  • While it is true that a decomposition $A=LL'$ is possible without full rank, the Cholesky algorithm only works with regular $A$. So without full rank, it cannot be a Cholesky decomposition. Computationally, one could do this decomposition in the singular case by diagonalization. (Although this is far more expensive) – Horst Grünbusch Aug 13 '14 at 12:14
  • @Horst: Why would $L=Q\sqrt{D}Q^T$ be lower triangular? – amoeba Aug 13 '14 at 12:41
  • 1
    @amoeba While one could organize it so it is, it doesn't have to be lower triangular for the argument to work - it's a feature of the Cholesky but it's not required for the result to work. – Glen_b Aug 13 '14 at 13:02
  • @Glen_b: Precisely. So the answer, as it reads now (after the edit), contains a mistake. – amoeba Aug 13 '14 at 13:05
  • 1
    @Glen Is being symmetric a necessary condition for being PSD or is that definition one of many? – 114 Aug 13 '14 at 13:42
  • @amoeba: One lesson in proper answering and editing to me: Omit as much as possible. Then it will be a great answer. – Horst Grünbusch Aug 13 '14 at 16:58
  • Thanks for the edits @amoeba and Horst. I came back this morning (my time) to fix it up and ... it's all done. – Glen_b Aug 13 '14 at 21:18
  • 1
    @114 for the relation between symmetric and PSD see http://math.stackexchange.com/questions/516533/symmetric-vs-positive-semidefinite – Frank Aug 13 '14 at 22:13