Para outras páxinas con títulos homónimos véxase: Distribución.

**Normal multivariante**| Función de multivariante |
| Función de distribución |
| **Parámetros** | {\displaystyle \mu =[\mu \_{1},\dots ,\mu \_{n}]^{T}} (vector real){\displaystyle \Sigma } matriz de covarianza (matriz real definida positiva de dimensión {\displaystyle n\times n} ) |
| **Soporte** | {\displaystyle x\in \mathbb {R} ^{n}\!} |
| **Función de densidade** | {\displaystyle f\_{X}(x\_{1},\dots ,x\_{n})={\frac {1}{(2\pi )^{n/2}\left|\Sigma \right|^{1/2}}}}{\displaystyle \exp \left(-{\frac {1}{2}}(x-\mu )^{\top }\Sigma ^{-1}(x-\mu )\right)} |
| **Función de distribución** | Sen expresión analítica |
| **Media** | {\displaystyle \mu \,\!} |
| **Mediana** | {\displaystyle \mu \,\!} |
| **Moda** | {\displaystyle \mu \,\!} |
| **Varianza** | {\displaystyle \sigma ^{2}\,\!} |
| **Asimetría** | {\displaystyle 0\,} |
| **Curtose** | {\displaystyle 0\,} |
| **Entropía** | {\displaystyle \ln \left({\sqrt {(2\,\pi \,e)^{n}\left|\Sigma \right|}}\right)\!} |
| **F. xeradora de momentos** | {\displaystyle M\_{X}(t)=\exp \left(\mu ^{\top }t+{\frac {1}{2}}t^{\top }\Sigma t\right)} |
| **Func. caract.** | {\displaystyle \phi \_{X}(t;\mu ,\Sigma )=\exp \left(i\mu ^{\top }t-{\frac {1}{2}}t^{\top }\Sigma t\right)} |

En probabilidade e estatística, unha **distribución normal multivariante**, tamén chamada **distribución gaussiana multivariante**, é unha xeneralización da distribución normal unidimensional a dimensións superiores.

Caso xeral
----------

Un vector aleatorio {\displaystyle \ X=[X\_{1},\dots ,X\_{n}]^{T}} segue unha distribución normal multivariante se satisfai as seguintes condicións equivalentes:

* Toda combinación linear {\displaystyle \ Y=a\_{1}X\_{1}+\cdots +a\_{n}X\_{n}} está normalmente distribuída.
* Hai un vector aleatorio {\displaystyle \ Z=[Z\_{1},\dots ,Z\_{m}]^{T}}, con compoñentes que son variables aleatorias independentes distribuídas segundo a normal estándar, un vector {\displaystyle \ \mu =[\mu \_{1},\dots ,\mu \_{n}]^{T}} e unha matriz {\displaystyle n\times m} {\displaystyle \ A} tal que {\displaystyle \ X=AZ+\mu }.
* Hai un vector \mu  e unha matriz semidefinida positiva simétrica {\displaystyle \ \Sigma } tal que a función característica de *X* es

> 
> 
>   
> 
> 
> 
> {\displaystyle \phi \_{X}\left(u;\mu ,\Sigma \right)=\exp \left(i\mu ^{\top }u-{\frac {1}{2}}u^{\top }\Sigma u\right).}
> 

Se {\displaystyle \ \Sigma } é unha matriz non singular, entón a distribución pode describirse pola seguinte función de densidade:

> 
> 
>   
> 
> {\displaystyle f\_{X}(x\_{1},\dots ,x\_{n})={\frac {1}{(2\pi )^{n/2}|\Sigma |^{1/2}}}\exp \left(-{\frac {1}{2}}(\mathbf {x} -\mathbf {\mu } )^{\top }\Sigma ^{-1}(\mathbf {x} -\mathbf {\mu } )\right)}
> 
> 
> 

onde {\displaystyle \ \left|\Sigma \right|} é o determinante de {\displaystyle \ \Sigma }. A ecuación redúcese á distribución normal se {\displaystyle \ \Sigma } é un escalar (é dicir, unha matriz 1x1).

O vector μ nestas circunstancias é a esperanza de *X* e a matriz {\displaystyle \ \Sigma =AA^{T}} é a matriz de covarianza das compoñentes *X**i*.

É importante comprender que a matriz de covarianza pode ser singular (aínda que non estea así descrita pola fórmula de arriba, para a que {\displaystyle \ \Sigma ^{-1}} está definida).

Este caso aparece con frecuencia en estatística; por exemplo, na distribución do vector de residuos en problemas ordinarios de regresión linear. Ademais, os *X**i* son en xeral *non* independentes; poden verse como o resultado de aplicar a transformación linear *A* a unha colección de variables normais *Z*.

Esta distribución dun vector aleatorio *X* que segue unha distribución normal multivariante pode ser descrita coa seguinte notación:

> 
> 
>   
> 
> {\displaystyle X\ \sim {\mathcal {N}}(\mu ,\Sigma ),}
> 
> 
> 

ou facer explícito que *X* é *n*-dimensional,

> 
> 
>   
> 
> {\displaystyle X\ \sim {\mathcal {N}}\_{N}(\mu ,\Sigma ).}
> 
> 
> 

### Función de distribución

A función de distribución {\displaystyle F(x)} defínese como a probabilidade de que todos os valores dun vector aleatorio X sexan menores ou iguais que os valores correspondentes dun vector x. Aínda que F non teña unha fórmula, hai unha serie de algoritmos que permiten estimala numericamente.

### Contraexemplo

O feito de que dúas variables aleatorias *X* e *Y* sigan unha distribución normal, cada unha, non implica que o par (*X*, *Y*) siga unha distribución normal conxunta. Un exemplo simple dáse con *X* Normal(0,1), *Y* = *X* se |*X*| > 1 e *Y* = −*X* se |*X*| < 1. isto tamén é certo para máis de dúas variables aleatorias.

### Normalmente distribuídas e independencia

Se *X* e *Y* están normalmente distribuídas e son independentes, a súa distribución conxunta tamén está normalmente distribuída, é dicir, o par (*X*, *Y*) debe ter unha distribución normal bivariante. En calquera caso, un par de variables aleatorias normalmente distribuídas non teñen por que ser independentes ao seren consideradas de forma conxunta.

Caso bivariante
---------------

No caso particular de dúas dimensións, a función de densidade (con media (0, 0)) é

{\displaystyle f(x,y)={\frac {1}{2\pi \sigma \_{x}\sigma \_{y}{\sqrt {1-\rho ^{2}}}}}\exp \left(-{\frac {1}{2(1-\rho ^{2})}}\left({\frac {x^{2}}{\sigma \_{x}^{2}}}+{\frac {y^{2}}{\sigma \_{y}^{2}}}-{\frac {2\rho xy}{(\sigma \_{x}\sigma \_{y})}}\right)\right)}
onde \rho  é o coeficiente de correlación entre X e {\displaystyle Y}. Neste caso,

{\displaystyle \Sigma ={\begin{bmatrix}\sigma \_{x}^{2}&\rho \sigma \_{x}\sigma \_{y}\\\rho \sigma \_{x}\sigma \_{y}&\sigma \_{y}^{2}\end{bmatrix}}.}

Transformación afín
-------------------

Se {\displaystyle Y=c+BX\,} é unha transformación afín de {\displaystyle X\ \sim {\mathcal {N}}(\mu ,\Sigma ),} onde c\, é un vector de constantes {\displaystyle M\times 1} e {\displaystyle B\,} unha matriz {\displaystyle M\times N}, entón {\displaystyle Y\,} ten unha distribución normal multivariante con esperanza {\displaystyle c+B\mu \,} e varianza {\displaystyle B\Sigma B^{T}\,} é dicir, {\displaystyle Y\sim {\mathcal {N}}\left(c+B\mu ,B\Sigma B^{T}\right)}. En particular, calquera subconxunto dos {\displaystyle X\_{i}\,} ten unha distribución marxinal que é tamén unha normal multivariante.

Para ver isto, pode considerarse o seguinte exemplo: para extraer o subconxunto {\displaystyle (X\_{1},X\_{2},X\_{4})^{T}\,}, emprégase

{\displaystyle B={\begin{bmatrix}1&0&0&0&0&\ldots &0\\0&1&0&0&0&\ldots &0\\0&0&0&1&0&\ldots &0\end{bmatrix}}}
o que extrae directamente os elementos desexados.

Outro corolario sería que a distribución de {\displaystyle Z=b\cdot X}, onde b é un vector da mesma lonxitude que X e o punto indica un produto vectorial, é unha distribución gaussiana unidimensional con {\displaystyle Z\sim {\mathcal {N}}\left(b\cdot \mu ,b^{T}\Sigma b\right)}. Este resultado obtense empregando

{\displaystyle B={\begin{bmatrix}b\_{1}&b\_{2}&\ldots &b\_{n}\\0&0&\ldots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\ldots &0\end{bmatrix}}}
e considerando só a primeira compoñente do produto (a primeira fila de B é o vector b). A definición positiva de {\displaystyle \Sigma } implica que a varianza do produto vectorial debería ser positiva.

Interpretación xeométrica
-------------------------

As curvas de equidensidade dunha distribución normal multivariante son elipsoides (é dicir, transformacións lineares de hiperesferas) centrados na media. as direccións dos eixes principais dos elipsoides veñen dadas polos vectores propios da matriz de covarianza {\displaystyle \Sigma }. As lonxitudes relativas dos cadrados dos eixes principais veñen dadas polos correspondentes vectores propios.

Se {\displaystyle \Sigma =U\Lambda U^{T}=U\Lambda ^{1/2}(U\Lambda ^{1/2})^{T}} é unha descomposición espectral onde as columnas de U son vectores propios unitarios e \Lambda  é unha matriz diagonal de valores propios, entón temos

{\displaystyle X\ \sim N(\mu ,\Sigma )\iff X\ \sim \mu +U\Lambda ^{1/2}N(0,I)\iff X\ \sim \mu +UN(0,\Lambda ).}
Ademais, *U* pode escollerse de tal modo que sexa unha matriz de rotación, tal que invertendo un eixe non teña ningún efecto en {\displaystyle N(0,\Lambda )}, pero invertendo unha columna, cambie o signo do determinante de *U'*. A distribución {\displaystyle N(\mu ,\Sigma )} é en efecto {\displaystyle N(0,I)} escalada por {\displaystyle \Lambda ^{1/2}}, rotada por *U* e trasladada por \mu .

Reciprocamente, calquera escolla de \mu , matriz de rango completo *U*, e valores diagonais positivos {\displaystyle \Lambda \_{i}} cede o paso a unha distribución normal non singular multivariante. Se calquera {\displaystyle \Lambda \_{i}} é cero e *U* é cadrada, a matriz de covarianza {\displaystyle U\Lambda U^{T}} é singular. Xeometricamente isto significa que cada curva elipsoide é infinitamente delgada e ten volume cero nun espazo *n*-dimensional, así como, polo menos, un dos principais eixes ten lonxitude cero.

Correlacións e independencia
----------------------------

En xeral, as variables aleatorias poden ser incorreladas, pero altamente dependentes, mais se un vector aleatorio ten unha distribución normal multivariante, entón calquera dúas ou máis das súas compoñentes que sexan incorreladas, son independentes.

Porén, non é certo que dúas variables aleatorias que están (separadamente, marxinalmente) normalmente distribuídas e incorreladas sexan independentes. Dúas variables aleatorias que están normalmente distribuídas poden que non o estean conxuntamente.

Momentos máis altos
-------------------

O momento estándar de *k*-ésima orde de *X* defínese como

{\displaystyle \mu \_{1,\dots ,N}(X)\ {\stackrel {\mathrm {def} }{=}}\ \mu \_{r\_{1},\dots ,r\_{N}}(X)\ {\stackrel {\mathrm {def} }{=}}\ E\left[\prod \limits \_{j=1}^{N}X\_{j}^{r\_{j}}\right]}
onde {\displaystyle r\_{1}+r\_{2}+\cdots +r\_{N}=k.}

Os momentos centrais de orde *k* veñen dados como segue:

(a) Se *k* é impar, {\displaystyle \mu \_{1,\dots ,N}(X-\mu )=0}.

(b) Se *k* é par, con {\displaystyle k=2\lambda }, entón

{\displaystyle \mu \_{1,\dots ,2\lambda }(X-\mu )=\sum \left(\sigma \_{ij}\sigma \_{k\ell }\cdots \sigma \_{XZ}\right)}
onda a suma se toma sobre todas as disposicións de conxuntos {\displaystyle \left\{1,\dots ,2\lambda \right\}} en \lambda  pares (non ordenados). É dicir, se se ten un *k*-ésimo ({\displaystyle {}=2\lambda =6}) momento central, estaranse a sumar os produtos de {\displaystyle \lambda =3} covarianzas (a notación -\mu  suprimiuse para facilitar a lectura):

{\displaystyle {\begin{aligned}&{}E[X\_{1}X\_{2}X\_{3}X\_{4}X\_{5}X\_{6}]\\&{}=E[X\_{1}X\_{2}]E[X\_{3}X\_{4}]E[X\_{5}X\_{6}]+E[X\_{1}X\_{2}]E[X\_{3}X\_{5}]E[X\_{4}X\_{6}]+E[X\_{1}X\_{2}]E[X\_{3}X\_{6}]E[X\_{4}X\_{5}]\\&{}+E[X\_{1}X\_{3}]E[X\_{2}X\_{4}]E[X\_{5}X\_{6}]+E[X\_{1}X\_{3}]E[X\_{2}X\_{5}]E[X\_{4}X\_{6}]+E[X\_{1}X\_{3}]E[X\_{2}X\_{6}]E[X\_{4}X\_{5}]\\&+E[X\_{1}X\_{4}]E[X\_{2}X\_{3}]E[X\_{5}X\_{6}]+E[X\_{1}X\_{4}]E[X\_{2}X\_{5}]E[X\_{3}X\_{6}]+E[X\_{1}X\_{4}]E[X\_{2}X\_{6}]E[X\_{3}X\_{5}]\\&+E[X\_{1}X\_{5}]E[X\_{2}X\_{3}]E[X\_{4}X\_{6}]+E[X\_{1}X\_{5}]E[X\_{2}X\_{4}]E[X\_{3}X\_{6}]+E[X\_{1}X\_{5}]E[X\_{2}X\_{6}]E[X\_{3}X\_{4}]\\&+E[X\_{1}X\_{6}]E[X\_{2}X\_{3}]E[X\_{4}X\_{5}]+E[X\_{1}X\_{6}]E[X\_{2}X\_{4}]E[X\_{3}X\_{5}]+E[X\_{1}X\_{6}]E[X\_{2}X\_{5}]E[X\_{3}X\_{4}].\end{aligned}}}
Isto dá lugar a {\displaystyle (2\lambda -1)!/(2^{\lambda -1}(\lambda -1)!)} termos na suma (15 no caso superior), nos que cada un é o produto de \lambda  (3 neste caso) covarianzas. Para momentos de orde cuarta (catro variables) hai tres termos. Para momentos de orde sexta hai 3 × 5 = 15 termos, e para momentos de orde oitava hai 3 × 5 × 7 = 105 termos.

As covarianzas son entón determinadas mediante a substitución dos termos da lista {\displaystyle \left[1,\dots ,2\lambda \right]} polos termos correspondentes da lista que consiste en {\displaystyle r\_{1}} uns, logo {\displaystyle r\_{2}} douses etc... Para ilustrar isto, pode examinarse o seguinte caso de momento central de orde cuarta:

{\displaystyle E\left[X\_{i}^{4}\right]=3\sigma \_{ii}^{2}}
{\displaystyle E\left[X\_{i}^{3}X\_{j}\right]=3\sigma \_{ii}\sigma \_{ij}}
{\displaystyle E\left[X\_{i}^{2}X\_{j}^{2}\right]=\sigma \_{ii}\sigma \_{jj}+2\left(\sigma \_{ij}\right)^{2}}
{\displaystyle E\left[X\_{i}^{2}X\_{j}X\_{k}\right]=\sigma \_{ii}\sigma \_{jk}+2\sigma \_{ij}\sigma \_{ik}}
{\displaystyle E\left[X\_{i}X\_{j}X\_{k}X\_{n}\right]=\sigma \_{ij}\sigma \_{kn}+\sigma \_{ik}\sigma \_{jn}+\sigma \_{in}\sigma \_{jk}.}
onde {\displaystyle \sigma \_{ij}} é a covarianza de {\displaystyle X\_{i}} e {\displaystyle X\_{j}}. A idea do método descrito é que primeiro se atopa o caso xeral para o momento k-ésimo, onde se teñen k diferentes variables X - {\displaystyle E\left[X\_{i}X\_{j}X\_{k}X\_{n}\right]} e entón pódense simplificar apropiadamente. Se se ten {\displaystyle E\left[X\_{i}^{2}X\_{k}X\_{n}\right]} entón, simplemente se considera {\displaystyle X\_{i}=X\_{j}} e dedúcese que {\displaystyle \sigma \_{ii}=\sigma \_{i}^{2}}.

Distribucións condicionais
--------------------------

Se \mu  y {\displaystyle \Sigma } se dividen como segue:

{\displaystyle \mu ={\begin{bmatrix}\mu \_{1}\\\mu \_{2}\end{bmatrix}}\quad } con tamaños {\displaystyle {\begin{bmatrix}q\times 1\\(N-q)\times 1\end{bmatrix}}}
{\displaystyle \Sigma ={\begin{bmatrix}\Sigma \_{11}&\Sigma \_{12}\\\Sigma \_{21}&\Sigma \_{22}\end{bmatrix}}\quad } con tamaños {\displaystyle {\begin{bmatrix}q\times q&q\times (N-q)\\(N-q)\times q&(N-q)\times (N-q)\end{bmatrix}}}
entón a distribución de {\displaystyle x\_{1}} condicionada a {\displaystyle x\_{2}=a} é unha normal multivariante {\displaystyle (X\_{1}|X\_{2}=a)\sim N({\bar {\mu }},{\overline {\Sigma }})} onde

{\displaystyle {\bar {\mu }}=\mu \_{1}+\Sigma \_{12}\Sigma \_{22}^{-1}\left(a-\mu \_{2}\right)}
e matriz de covarianza

{\displaystyle {\overline {\Sigma }}=\Sigma \_{11}-\Sigma \_{12}\Sigma \_{22}^{-1}\Sigma \_{21}.}
Esta matriz é o complemento de Schur de {\displaystyle {\mathbf {\Sigma } \_{22}}} en {\displaystyle {\mathbf {\Sigma } }}. Isto significa que para calcular a matriz condicional de covarianza se inverte a matriz global de covarianza, desprézanse as filas e as columnas correspondentes ás variables baixo as que está condicionada e entón invértese de novo para conseguir a matriz condicional de covarianza.

Sábese que {\displaystyle x\_{2}=a} altera a varianza, aínda que a nova varianza non dependa do valor específico de a; quizais resulta máis sorprendente que a media se cambia por {\displaystyle \Sigma \_{12}\Sigma \_{22}^{-1}\left(a-\mu \_{2}\right)}; pode compararse isto coa situación na que non se coñece o valor de a, caso no que {\displaystyle x\_{1}} tería como distribución

{\displaystyle N\_{q}\left(\mu \_{1},\Sigma \_{11}\right)}.

A matriz {\displaystyle \Sigma \_{12}\Sigma \_{22}^{-1}} coñécese como a matriz de coeficientes de regresión.

### Esperanza condicional bivariante

No caso

{\displaystyle {\begin{pmatrix}X\_{1}\\X\_{2}\end{pmatrix}}\sim {\mathcal {N}}\left({\begin{pmatrix}0\\0\end{pmatrix}},{\begin{pmatrix}1&\rho \\\rho &1\end{pmatrix}}\right)}

entón

{\displaystyle E(X\_{1}|X\_{2}>z)=\rho {\phi (z) \over \Phi (-z)}}

onde esta última razón adoita chamarse razón inversa de Mills.

Matriz de información de Fisher
-------------------------------

A matriz de información de Fisher (MIF) para unha distribución normal toma unha formulación especial. O elemento {\displaystyle (m,n)} da MIF para {\displaystyle X\sim N(\mu (\theta ),\Sigma (\theta ))} é

{\displaystyle {\mathcal {I}}\_{m,n}={\frac {\partial \mu }{\partial \theta \_{m}}}\Sigma ^{-1}{\frac {\partial \mu ^{\top }}{\partial \theta \_{n}}}+{\frac {1}{2}}\mathrm {tr} \left(\Sigma ^{-1}{\frac {\partial \Sigma }{\partial \theta \_{m}}}\Sigma ^{-1}{\frac {\partial \Sigma }{\partial \theta \_{n}}}\right)}
onde

* {\displaystyle {\frac {\partial \mu }{\partial \theta \_{m}}}={\begin{bmatrix}{\frac {\partial \mu \_{1}}{\partial \theta \_{m}}}&{\frac {\partial \mu \_{2}}{\partial \theta \_{m}}}&\cdots &{\frac {\partial \mu \_{N}}{\partial \theta \_{m}}}&\end{bmatrix}}}
* {\displaystyle {\frac {\partial \mu ^{\top }}{\partial \theta \_{m}}}=\left({\frac {\partial \mu }{\partial \theta \_{m}}}\right)^{\top }={\begin{bmatrix}{\frac {\partial \mu \_{1}}{\partial \theta \_{m}}}\\\\{\frac {\partial \mu \_{2}}{\partial \theta \_{m}}}\\\\\vdots \\\\{\frac {\partial \mu \_{N}}{\partial \theta \_{m}}}\\\\\end{bmatrix}}}
* {\displaystyle {\frac {\partial \Sigma }{\partial \theta \_{m}}}={\begin{bmatrix}{\frac {\partial \Sigma \_{1,1}}{\partial \theta \_{m}}}&{\frac {\partial \Sigma \_{1,2}}{\partial \theta \_{m}}}&\cdots &{\frac {\partial \Sigma \_{1,N}}{\partial \theta \_{m}}}\\\\{\frac {\partial \Sigma \_{2,1}}{\partial \theta \_{m}}}&{\frac {\partial \Sigma \_{2,2}}{\partial \theta \_{m}}}&\cdots &{\frac {\partial \Sigma \_{2,N}}{\partial \theta \_{m}}}\\\\\vdots &\vdots &\ddots &\vdots \\\\{\frac {\partial \Sigma \_{N,1}}{\partial \theta \_{m}}}&{\frac {\partial \Sigma \_{N,2}}{\partial \theta \_{m}}}&\cdots &{\frac {\partial \Sigma \_{N,N}}{\partial \theta \_{m}}}\end{bmatrix}}}
* {\displaystyle \mathrm {tr} } é a función traza dunha matriz.

Diverxencia de Kullback-Leibler
-------------------------------

A diverxencia de Kullback-Leibler de {\displaystyle N0\_{N}(\mu \_{0},\Sigma \_{0})} a {\displaystyle N1\_{N}(\mu \_{1},\Sigma \_{1})} é:

{\displaystyle D\_{\text{KL}}(N0\|N1)={1 \over 2}\left(\log \_{e}\left({\det \Sigma \_{1} \over \det \Sigma \_{0}}\right)+\mathrm {tr} \left(\Sigma \_{1}^{-1}\Sigma \_{0}\right)+\left(\mu \_{1}-\mu \_{0}\right)^{\top }\Sigma \_{1}^{-1}(\mu \_{1}-\mu \_{0})-N\right).}
O logaritmo debe tomarse con base *e* nos dous termos (logaritmos neperianos); seguindo o logaritmo están os logaritmos neperianos das expresións que son ambos os factores da función de densidade ou se non, xorden naturalmente. A diverxencia de arriba mídese en nats. Dividindo a expresión de arriba por log*e* 2 dáse paso a a diverxencia en bits.

Estimación de parámetros
------------------------

A función de densidade de probabilidade dunha normal multivariante *N*-dimensional é

{\displaystyle f(x)=(2\pi )^{-N/2}\det(\Sigma )^{-1/2}\exp \left(-{1 \over 2}(x-\mu )^{T}\Sigma ^{-1}(x-\mu )\right)}
e o estimador de máxima verosimilitude da matriz de covarianza para unha mostra de *n* observacións é

{\displaystyle {\widehat {\Sigma }}={1 \over n}\sum \_{i=1}^{n}(X\_{i}-{\overline {X}})(X\_{i}-{\overline {X}})^{T}}
o que é, simplemente, a matriz da mostra de covarianza. Este é un estimador con nesgo que ten como esperanza

{\displaystyle E[{\widehat {\Sigma }}]={n-1 \over n}\Sigma .}
Unha covarianza da mostra sen nesgo é

{\displaystyle {\widehat {\Sigma }}={1 \over n-1}\sum \_{i=1}^{n}(X\_{i}-{\overline {X}})(X\_{i}-{\overline {X}})^{T}.}

Entropía
--------

A entropía diferencial da distribución normal multivariante é

{\displaystyle {\begin{aligned}h\left(f\right)&=-\int \_{-\infty }^{\infty }\int \_{-\infty }^{\infty }\cdots \int \_{-\infty }^{\infty }f(x)\ln f(x)\,dx\\&={\frac {1}{2}}\left(N+N\ln \left(2\pi \right)+\ln \left|\Sigma \right|\right)\\&={\frac {1}{2}}\ln\{(2\pi e)^{N}\left|\Sigma \right|\}\end{aligned}}}
onde {\displaystyle \left|\Sigma \right|} é o determinante da matriz de covarianza {\displaystyle \Sigma }.

Tests de normalidade multivariante
----------------------------------

Os tests de normalidade multivariante comproban a similitude dun conxunto dado de datos coa distribución normal multivariante. A hipótese nula é que o conxunto de datos é similar á distribución normal e polo tanto un p-valor suficientemente pequeno indica datos non normais. Os tests de normalidade multivariante inclúen o test de Cox-Small e a adaptación de Smith y Jain
 do test de Friedman-Rafsky.

Simulando valores da distribución
---------------------------------

Un método amplamente empregado para simular un vector aleatorio X da distribución normal multivariada {\displaystyle N}-dimensional con vector de medias \mu  e matriz de covarianza {\displaystyle \Sigma } (requirida para que sexa simétrica e definida positiva) funciona como segue:

1. Calcúlase a descomposición de Cholesky de {\displaystyle \Sigma }, é dicir, atópase a única matriz triangular inferior A tal que {\displaystyle A\,A^{T}=\Sigma }. Podería empregarse calquera outra matriz A que satisfaga esta condición, ou sexa, que é un a raíz cadrada de {\displaystyle \Sigma }, mais habitualmente atopar esa matriz, distinta da da descomposición de Cholesky, sería bastante máis custoso en termos de computación.
2. Sexa {\displaystyle Z=(z\_{1},\dots ,z\_{N})^{T}} un vector con compoñentes {\displaystyle N} normais e independentes que varían e que pode xerarse, por exemplo, empregando o método de Box-Muller.
3. Sexa {\displaystyle X=\,} {\displaystyle \mu +AZ\,.}

Notas
-----

1. ↑  Pode verse MVNDST en , que inclúe código FORTRAN ou , con código MATLAB.
2. ↑  Nikolaus Hansen. "The CMA Evolution Strategy: A Tutorial" (PDF). Arquivado dende o orixinal (PDF) o 31 de marzo de 2010. Consultado o 19 de maio de 2016.
3. ↑  Gokhale, DV; Ahmed, NA; Res, BC; Piscataway, NJ (1989). "Entropy Expressions and Their Estimators for Multivariate Distributions". *Information Theory, IEEE Transactions on* **35** (3): 688–692. doi:10.1109/18.30996.
4. ↑  Cox, D. R.; Small, N. J. H. (1978). "Testing multivariate normality". *Biometrika* **65** (2): 263–272. doi:10.1093/biomet/65.2.263.
5. ↑  Smith, Stephen P.; Jain, Anil K. (1988). "A test to determine the multivariate normality of a dataset". *IEEE Transactions on Pattern Analysis and Machine Intelligence* **10** (5): 757–761. doi:10.1109/34.6789.

|  |  |
| --- | --- |
| Control de autoridades | 
* Wd: Q1149000
* GND: 4227589-1
* MAG: 177384507

 |