12 votos

¿Qué es $\mathbb E \lVert X \rVert$ para una normal multivariante $X \sim \mathcal N(\mu, \Sigma)$?

$\DeclareMathOperator\E{\mathbb E} \DeclareMathOperator\Var{\mathrm{Var}} \newcommand\R{\mathbb R} \DeclareMathOperator\N{\mathcal N} \DeclareMathOperator\tr{\mathrm{tr}}$Suponga $X \sim \N(\mu, \Sigma)$ donde $\mathrm{supp}(X) = \R^n$. Podemos suponer $\Sigma$ a ser nonsingular. ¿Qué es $\E \lVert X \rVert$?

Si $\Sigma = \sigma I$, $\lVert X \rVert / \sigma$ sigue un noncentral chi distribución, que ha conocido la media (en términos de un polinomio generalizado de Laguerre). Que se da en este caso $$ \E \lVert X \rVert = \sigma \sqrt{\frac{\pi}{2}} \; L_{{1/2}}^{{(n/2-1)}}\left({\frac {-\lVert \mu \rVert^{2}}{2 \sigma}}\right)\, ;$$ si asumimos, además, que $\mu = 0$, entonces se convierte en una distribución chi y $$ \E \lVert X \rVert = \sigma \sqrt {2} \, \frac{\Gamma((n+1)/2)}{\Gamma(n/2)} .$$

En el caso general, tenemos un fácil límite superior a través de la desigualdad de Jensen. Dejando $C^T C = \Sigma$$Z \sim \N(0, I)$: \begin{align} \left( \E \lVert X \rVert \right)^2 &< \E \lVert X \rVert^2 \\&= \E \lVert \mu + C Z \rVert^2 \\&= \lVert \mu \rVert^2 + 2 \mu^T C \E Z + \E \tr( Z^T C^T C Z ) \\&= \lVert \mu \rVert^2 + \tr( C^T C \,\E Z Z^T ) \\&= \lVert \mu \rVert^2 + \tr( \Sigma ) ,\end{align} donde la desigualdad es estricta, ya hemos asumido $\lVert X \rVert$ es no degenerada.

También tenemos un límite inferior de la misma manera: desde $\lVert \cdot \rVert$ es convexa, $\E \lVert X \rVert \ge \lVert \E X \rVert = \lVert \mu \rVert$. (Yo previamente había algo más complicado aquí, basado en la delimitación de la varianza de la función de $z \mapsto \lVert \mu + C z \rVert$, pero resultó en una peor obligado.)

Podemos apretar estos límites o, preferiblemente, encontrar una expresión exacta para $\E \lVert X \rVert$ en el caso general? Si no, ¿qué acerca de la diagonal, pero no esférica caso, o $\mu = 0$ o de otras interesantes subcases?


Ha habido cierta discusión en este sitio de la distribución de $Y = \lVert X \rVert^2$, sobre todo aquí y aquí. Si tenemos la distribución completa de la $Y = \lVert X \rVert^2$, podemos ser capaces de encontrar $\E \sqrt{Y}$. Que la distribución es bruto, a pesar de que, analíticas y de respuestas parecen difíciles.

9voto

kjetil b halvorsen Puntos 7012

No voy a escribir una respuesta completa, que va a tomar mucho tiempo ahora. La respuesta completa a la pregunta se puede encontrar en Mathai Y Provost: "Cuadráticas formas en las Variables Aleatorias", véase, por ejemplo, mi respuesta: http://math.stackexchange.com/questions/442472/sum-of-squares-of-dependent-gaussian-random-variables/442916#442916

La respuesta es largo, complicado e implica funciones especiales, tales como "la Lauricella de la función". No sé si hay un paquete de R para que (podía encontrar uno). La respuesta se puede encontrar en ese libro en la sección 3.2.b.5: Inversa y la fracción de momentos en el Nonsingular Caso Normal, entre las páginas 56--61.

La distribución de la forma cuadrática en las variables normales puede ser calculada con el paquete R CompQuadForm, tal vez usted podría tratar de integración numérica con la ayuda de ese paquete.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X