-1
$\begingroup$

Let $X$ be a Gaussian vector in dimension $n$, with $0$ mean and covariance identity. Let $A$ be a square matrix of size $n$, and $Y = A X$. Let $N$ be the square of $Y$ euclidean norm: $N = \sum Y_i^2$. One computes easily the mean of $N$: $E[N] = \text{Tr}(A A')$. But what about its variance?

$\endgroup$
7
  • 2
    $\begingroup$ This is simplified a lot by using the singular value decomposition to reduce to the case in which $A$ is diagonal. If my arithmetic is right the variance is $2 \sum \sigma_i^4$, where $\sigma_i$ are the singular values of $A$. $\endgroup$ Sep 16, 2013 at 12:36
  • 1
    $\begingroup$ @Ricardo Andrade: Is there a reason you felt it necessary to capitalize Gaussian but not euclidean? $\endgroup$ Sep 16, 2013 at 12:37
  • $\begingroup$ Dear @Mark: That is a very good point. I am afraid the capitalization was a very silly lapse on my part. I was editing the tags and noticed the occurrence of 'gaussian' in the title. Unfortunately, I completely missed the occurrence of 'euclidean'. Now the title looks unbalanced. I apologize for my poor edit. /// Dear msfr: I am sorry for my poor edit. Please undo or correct my partial capitalization. $\endgroup$ Sep 16, 2013 at 17:22
  • $\begingroup$ Dear Mark and Ricardo. Actually, I wanted to avoid svd ! $\endgroup$
    – msfr
    Sep 16, 2013 at 19:22
  • 1
    $\begingroup$ @msfr: If you just want, say, an upper or lower bound then one can find more easily computable expressions which may or may not be sharp enough for your purposes, but the variance is what it is. $\endgroup$ Sep 16, 2013 at 20:46

1 Answer 1

0
$\begingroup$

Apparently, one possible answer is $$ E[N^2] = 3\ \text{Tr} (A A' A A') + 2 \sum_{\text{all } A_i 2\times 2 \text{ submatrices}} \text{Det} (A_i)^2. $$ Does this sound familiar to anyone ?

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.