Loi T² d'Hotelling

En statistiques, plus particulièrement dans les tests d'hypothèses, la loi TModèle:2 de Hotelling, proposée par Harold Hotelling, est une loi de probabilité multivariée qui est étroitement liée à la loi de Fisher et qui se distingue surtout par le fait qu'elle apparaît comme la loi d'un ensemble de statistiques d'échantillon qui sont des généralisations naturelles des statistiques sous-jacentes à la loi t de Student. La statistique tModèle:2 de Hotelling est une généralisation de la statistique t de Student qui est utilisée dans les tests d'hypothèses multivariées.

Motivation

La loi apparaît dans les statistiques multivariées lors de la réalisation de tests sur les différences entre les moyennes (multivariées) de différentes populations, où les tests pour les problèmes univariés utiliseraient un test t . La loi porte le nom de Harold Hotelling, qui l'a développée comme une généralisation de la loi t de Student.

Définition

Si le vecteur $d$ est une loi normale multivariée avec une moyenne nulle et une matrice de covariance unitaire $𝒩({\mathrm{𝟎}}_p,{𝐈}_{p,p})$ et $M$ est une matrice $p\times p$ aléatoire suivant un loi de Wishart $W({𝐈}_{p,p},m)$ avec une matrice d'échelle unitaire et à m degrés de liberté, et d et M sont indépendants l'un de l'autre, alors la forme quadratique $X$ suit une loi de Hotelling (de paramètres $p$ et $m$)

$\frac{m-p+1}{pm}X\sim F_{p,m-p+1}$

On peut montrer que si une variable aléatoire X suit la loi TModèle:2 de Hotelling, $X\sim T_{p,m}^2$, alors:

$t^2\sim T_{p,n-1}^2=\frac{p(n-1)}{n-p}{F}_{p,n-p},$

où $F_{p,m-p+1}$ est la loi de Fisher de paramètres p et m − p + 1.

Statistique tModèle:2 d'Hotelling

Soit $\hat{\mathrm{\Sigma }}$ la covariance de l'échantillon :

$\hat{\mathrm{\Sigma }}=\frac{1}{n-1}\sum _{i=1}^n({𝐱}_i-\overline{𝐱})({𝐱}_i-\overline{𝐱}{)}^T$

où l'exposant T désigne la transposition. On peut montrer que $\hat{\mathrm{\Sigma }}$ est une matrice semi-définie positive et $(n-1)\hat{\mathrm{\Sigma }}$ suit une loi de Wishart p-variée à n − 1 degrés de liberté. L'échantillon de matrice de covariance de la moyenne s'obtient par ${\hat{\mathrm{\Sigma }}}_{\overline{𝐱}}=\hat{\mathrm{\Sigma }}/n$ .

La statistique tModèle:2 de Hotelling est alors définie comme:

$t^2=n(\overline{𝐱}-\mathit{\mu }{)}^{\mathrm{T}}{\hat{\mathrm{\Sigma }}}_{\overline{𝐱}}^{-1}(\overline{𝐱}-\mathit{\mu }),$

qui est proportionnelle à la distance de Mahalanobis entre la moyenne de l'échantillon et $\mathit{\mu }$ . Pour cette raison, on devrait s’attendre à ce que la statistique prenne des valeurs faibles si $\overline{𝐱}\approx \mathit{\mu }$, et des valeurs élevées si elles sont différentes.

Par la définition de la loi,

$X=m{d}^T{M}^{-1}d\sim T^2(p,m).$

où $F_{p,n-p}$ est la loi de Fisher avec les paramètres p et n − p .

Afin de calculer une valeur p (sans rapport avec la variable p ici), il faut remarquer que la loi de $t^2$ implique de manière équivalente que

$\frac{n-p}{p(n-1)}{t}^2\sim F_{p,n-p}.$

Ensuite, on utilise la quantité sur le côté gauche pour évaluer la valeur p correspondant à l'échantillon, qui provient de la loi de Fisher. Un ellipsoïde de confiance peut également être déterminé en utilisant une logique similaire.

Motivation

Soit ${𝒩}_p(\mathit{\mu },\mathrm{\Sigma })$ une loi normale p-variée de vecteur moyenne $\mathit{\mu }$ et matrice de covariance connue $\mathrm{\Sigma }$. Soient

${𝐱}_1,\dots ,{𝐱}_n\sim {𝒩}_p(\mathit{\mu },\mathrm{\Sigma })$

n variables aléatoires indépendantes identiquement distribuées (iid), qui peuvent être représentées comme $p\times 1$ vecteurs colonnes de nombres réels. On pose

$\overline{𝐱}=\frac{1}{n_x}\sum _{i=1}^{n_x}{𝐱}_i\overline{𝐲}=\frac{1}{n_y}\sum _{i=1}^{n_y}{𝐲}_i$

comme la moyenne de l'échantillon avec de covariance ${\mathrm{\Sigma }}_{\overline{𝐱}}=\mathrm{\Sigma }/n$ . On peut montrer que

$(\overline{𝐱}-\mathit{\mu }{)}^{\mathrm{T}}{\mathrm{\Sigma }}_{\overline{𝐱}}^{-1}(\overline{𝐱}-\mathit{\mu })\sim {\chi }_p^2,$

où ${\chi }_p^2$ est la loi du chi carré avec p degrés de liberté.

Statistique sur deux échantillons

Si on définit ${𝐱}_1,\dots ,{𝐱}_{n_x}\sim N_p(\mathit{\mu },\mathrm{\Sigma })$ et ${𝐲}_1,\dots ,{𝐲}_{n_y}\sim N_p(\mathit{\mu },\mathrm{\Sigma })$, avec les échantillons tirés indépendamment de deux lois normales multivariées indépendantes avec la même moyenne et la même covariance, et on pose

${\hat{\mathrm{\Sigma }}}_{𝐱}=\frac{1}{n_x-1}\sum _{i=1}^{n_x}({𝐱}_i-\overline{𝐱})({𝐱}_i-\overline{𝐱}{)}^T$

${\hat{\mathrm{\Sigma }}}_{𝐲}=\frac{1}{n_y-1}\sum _{i=1}^{n_y}({𝐲}_i-\overline{𝐲})({𝐲}_i-\overline{𝐲}{)}^T$

comme les moyennes des échantillons, et

$\hat{\mathrm{\Sigma }}=\frac{(n_x-1){\hat{\mathrm{\Sigma }}}_{𝐱}+(n_y-1){\hat{\mathrm{\Sigma }}}_{𝐲}}{n_x+n_y-2}$

comme les matrices de covariance d'échantillon respectives. Alors

$t^2=\frac{n_x{n}_y}{n_x+n_y}(\overline{𝐱}-\overline{𝐲}{)}^T{\hat{\mathrm{\Sigma }}}^{-1}(\overline{𝐱}-\overline{𝐲})\sim T^2(p,n_x+n_y-2)$

est l'estimateur non biaisé de la matrice de covariance groupée (une extension de la variance composite).

Enfin, la statistique tModèle:2 à deux échantillons de Hotelling est

$\frac{n_x+n_y-p-1}{(n_x+n_y-2)p}{t}^2\sim F(p,n_x+n_y-1-p).$

Notions associées

Cela peut être lié à la loi de Fisher par

$\frac{n_x+n_y-p-1}{(n_x+n_y-2)p}{t}^2\sim F(p,n_x+n_y-1-p;\delta ),$

La loi non nulle de cette statistique est la loi de Fisher non centrée (le rapport d'une variable aléatoire suivant la loi du χ² non centrée et d'une variable aléatoire centrale indépendante suivant une loi du χ² )

$\delta =\frac{n_x{n}_y}{n_x+n_y}{𝒅}^T{\mathrm{\Sigma }}^{-1}𝒅,$

avec

$d_1={\bar{x}}_1-{\bar{y}}_1,d_2={\bar{x}}_2-{\bar{y}}_2$

où $𝒅=\overline{𝐱}\mathbf{-}\overline{𝐲}$ est le vecteur de différence entre les moyennes de la population.

Dans le cas à deux variables, la formule se simplifie bien, permettant d'apprécier comment la corrélation $\rho$ entre les variables affecte $t^2$ . Si l'on définit

$d_1={\bar{x}}_1-{\bar{y}}_1,d_2={\bar{x}}_2-{\bar{y}}_2$

et

$s_1=\sqrt{{\mathrm{\Sigma }}_{11}}{s}_2=\sqrt{{\mathrm{\Sigma }}_{22}}\rho ={\mathrm{\Sigma }}_{12}/(s_1{s}_2)={\mathrm{\Sigma }}_{21}/(s_1{s}_2)$

alors

$t^2=\frac{n_x{n}_y}{(n_x+n_y)(1-{\rho }^2)}[{\left(\frac{d_1}{s_1}\right)}^2+{\left(\frac{d_2}{s_2}\right)}^2-2\rho \left(\frac{d_1}{s_1}\right)\left(\frac{d_2}{s_2}\right)]$

Ainsi, si les différences entre les deux lignes du vecteur $𝐝=\overline{𝐱}-\overline{𝐲}$ sont du même signe, en général, $t^2$ devient plus petit à mesure $\rho$ devient plus positif. Si les différences sont de signes opposés, $t^2$ devient plus grand à mesure $\rho$ devient plus positif.

Un cas particulier univarié peut être trouvé dans le test t de Welch.

Des tests plus robustes et plus puissants que le test à deux échantillons de Hotelling ont été proposés dans la littérature, voir par exemple les tests basés sur la distance entre points qui peuvent également être appliqués lorsque le nombre de variables est comparable, voire supérieur, au nombre de sujetsModèle:,.

Voir également

• Test de Student pour les statistiques univariées
• Loi de Student pour les probabilités univariées
• Loi de Student multivariée
• Loi de Fisher (généralement tabulée ou disponible dans les bibliothèques numériques, et donc utilisée pour tester la statistique TModèle:2 en utilisant la relation donnée ci-dessus)
• Loi du lambda de Wilks (dans les statistiques multivariées, la Modèle:Lien de Wilks est au T ² de Hotelling ce que la loi de Snedecor est à la loi de Student dans les statistiques univariées)

Liens externes

• Modèle:EncycloMath
• Modèle:Mathworld