Imagine I have $x\stackrel{d}{=}\vert\xi A_i U\vert$,
where $\xi$ is 1-dimensional, non negative random variable with finite 2nd moment and independent of $U$, $A_i$ is a fixed $n$-dimensional vector and $U$ is uniformly distributed on the $n$-dimensional unit sphere $\mathbb{S}_{n-1}$ (see for a definition here: https://en.wikipedia.org/wiki/Unit_sphere#Unit_spheres_and_balls_in_Euclidean_space - the first equation)
Now I want to calculate the median of it; Or more generally I don't need to calculate it explicitly, but I want to show the following:
Let $x\stackrel{d}{=}\vert\xi A U\vert$, where $\xi$ is defined the same way, U as well but $A$ denotes a symmetric positive definite $n \times n$ matrix;
Then $x$ is also $n$-dimensional; What I want to prove now is:
$$\frac{\text{med}(\vert\xi A_i U\vert)}{dA_iA_i^\top }=\frac{\text{med}(\vert\xi A_j U\vert)}{dA_jA_j^\top }$$ for all $i,j \in \{1, ...,n\}$; d denotes here: $\sqrt{\mathbb E(\xi^2)}$;
But I really have no clue on how to do it - calculating a median is tricky because even if two random variables are independent, the median of the product can't be separated generally;
If you have any ideas - even if they don't lead to the solution immediately, even if they later turn out to be wrong - let me know; I really appreciate any help.
Thank you
For Marcus I try to extend my question so you can see where it is coming from:
I have an elliptical distributed random variable- so I have:
$x\stackrel{d}{=}\mu+\xi A U$,
where $\xi$ is 1-dimensional, non negative random variable with finite 2nd moment and independent of $U$, $A$ is a fixed $n \times n$, symmetric positive definite matrix and $U$ is uniformly distributed on the $n$-dimensional unit sphere $\mathbb{S}_{n-1}$ and is a column vector;
Then $x$ is an $n$-dimensional vector and we have:
$x_i\stackrel{d}{=}\mu+\xi A_i U$,
where $A_i$ denotes the ith row; Now I want the median absolute deviation and the standard error; more precisely I want to prove that:
$$\frac{\text{med}(\vert x_i-\text{med}(x_i)\vert)}{sd(x_i)}=\frac{\text{med}(\vert x_j-\text{med}(x_j)\vert)}{sd(x_j)}$$
Here sd means the square root of the variance and the upper term is just the definition of the median absolute deviation; According to google the median of $x$ is $\mu$, that's why the first term reduces to the median of $\vert\xi A_i U\vert$ - for the standard deviation we use the definition of the variance:
Var$(x_i)=\mathbb E((x_i-\mathbb E(x_i))(x_i-\mathbb E(x_i))^\top)=E((\xi A_i U)(\xi A_i U)^\top)=\mathbb E(\xi^2)E(A_i UU^\top A_i^\top)=\mathbb E(\xi)^2A_iE(UU^\top )A_i^\top=\mathbb E(\xi)^2A_iA_i^\top$
where we used the fact that 1. $\mathbb E(x_i)=\mu_i$ ($\mu$ is median and expectation) 2.independence of $\xi$ and $U$ and that $U*U^\top=1$ by definition of being on the unit sphere;
Then I really need to correct myself; I want to show that:
$$\frac{\text{med}(\vert\xi A_i U\vert)}{d\sqrt{A_iA_i^\top} }=\frac{\text{med}(\vert\xi A_j U\vert)}{d\sqrt{A_jA_j^\top} }$$ for all $i,j \in \{1, ...,n\}$; d denotes here: $\sqrt{\mathbb E(\xi^2)}$;