Welcome to MathOverflow

Question

1

Recently i got stuck into a problem and couldn't find its satisfactory answer anywhere.

My question is simple. Suppose i have a fat random matrix (i,e $R$ has dimensions $k\times d$ where $k<d$ ) whose elements are chosen from a i.i.d standard normal distribution N(0,1).
If i find its pseudo inverse, given by: $R^+ = (R' R)^{-1} R'$.

Will the this pseudo inverse matrix will still remain random ?
If yes, will it contain elements distributed with normal distribution?
If yes, what would be the mean and variance of this this normal distribution?

I am asking these questions because i have experimented with a lot of random matrices (with elements distributed with N(0,1). When in plot a histogram of pseudo inverse elements, it comes a normal distribution with mean = 0 and variance = 1/(variance of $R$ $\times d^2$) ; where d are the columns in R.)

I have tried to find PDF using Jacobian transform but i could not figure out how will it shape up the variance.

I would be thankful if you could guide me or clarify my problem. Thanks,

edited Dec 13 at 14:41

Mark Meckes
5,972●8●31

asked Dec 13 at 13:55

Salman
11●1

I cleaned up your math by putting it into LaTeX. I think though that your pseudoinverse formula is wrong. – Mark Meckes Dec 13 at 14:46

2

1. Yes, of course. 2. No. Think about the case k=d: the inverse of a single normal random variable is not normal. This paper looks like a good place to find the kind of information you want: ugr.es/~ramongs/articulos%20en%20pdf/cimat1.pdf – Mark Meckes Dec 13 at 14:49

Thanks Mark. But can you tell a condition under which inverse of a single normal random variable becomes normal. When i plot histogram of the inverse of \mathbf{R} with k and d very large, i get a nearly normal distribution. – Salman Dec 15 at 8:01

Distribution of Inverse of a Random Matrix

Remember to vote up questions/answers you find interesting or helpful (requires 15 reputation points)

You can accept an answer to one of your own questions by clicking the check mark next to it. This awards 15 reputation points to the person who answered and 2 reputation points to you.