Hotellings T-kwadraat

Hotellings T², kort voor de Hotellings ${\displaystyle T^{2}}$-toets, is een statistische toets, genoemd naar Harold Hotelling, die men kan zien als een multivariate generalisatie van de t-toets. Met Hotellings ${\displaystyle T^{2}}$ wordt eigenlijk de toetsingsgrootheid bedoeld, gedefinieerd door:

${\displaystyle T^{2}=n({\overline {X}}-\mu _{0})'S^{-1}({\overline {X}}-\mu _{0}),}$

gebaseerd op een (aselecte) steekproef ${\displaystyle X_{1},\ldots ,X_{n}}$ uit een ${\displaystyle p}$-dimensionale verdeling.

In de formule stelt op de gebruikelijke manier ${\displaystyle {\overline {X}}}$ het steekproefgemiddelde voor en

${\displaystyle S={\frac {1}{n-1}}\sum _{i=1}^{n}(X_{i}-{\overline {X}})(X_{i}-{\overline {X}})'}$

de ${\displaystyle p\times p}$-steekproefcovariantiematrix.

Als de ${\displaystyle p}$-dimensionale verdeling een normale verdeling is met vector van verwachtingswaarden ${\displaystyle \mu _{0}}$, geldt dat:

${\displaystyle {\frac {n-p}{p(n-1)}}T^{2}}$

een ${\displaystyle F(p,n-p)p}$-verdeling heeft, dus een F-verdeling met ${\displaystyle p}$ vrijheidsgraden in de teller en ${\displaystyle n-p}$ in de noemer.

Hotellings ${\displaystyle T^{2}}$ wordt vaak gebruikt om uitbijters in multivariate data te detecteren. De Hotelling ${\displaystyle T^{2}}$-toets is een multivariate t-toets die gebruikmaakt van de covariantiematrix ${\displaystyle S}$ van de data:

${\displaystyle S={\frac {1}{n-1}}\sum _{i=1}^{k}\sum _{j=1}^{n_{i}}(X_{ij}-{\overline {X}})(X_{ij}-{\overline {X}})'}$

Daarin is

• ${\displaystyle n}$ het totaal aantal stalen
• ${\displaystyle n_{i}}$ het aantal stalen van de genomen subgroep i
• ${\displaystyle k}$ het aantal gemeten variabelen per staal
• ${\displaystyle X_{ij}}$ de vector met het ${\displaystyle j}$-de staal uit de ${\displaystyle i}$-de groep
• ${\displaystyle {\overline {X}}}$ de vector met de gemiddelde meetwaarden over alle stalen en subgroepen

De waarden van Hotellings ${\displaystyle T^{2}}$ worden berekend afhankelijk van het feit of er subgroepen zijn en of het populatiegemiddelde bekend is. Men moet subgroepen definiëren, wanneer men merkt dat de meetgegevens in groepen voorkomen. Bijvoorbeeld, wanneer een industriële productie niet volcontinu, maar per ketel, per ton, per vrachtwagen binnengekomen grondstof, ... (men zegt dan ook per batch) gebeurt, dan vormt elke geproduceerde batch een afzonderlijke eenheid of subgroep in de meetgegevens. Van elke zo geproduceerde eenheid of subgroep neemt men dan 1 of meerdere stalen waarvan meetgegevens verzameld worden.

De waarden van Hotellings ${\displaystyle T^{2}}$ worden berekend volgens de vierde kolom van onderstaande tabel. Deze waarden volgen een F-verdeling, waarvan het aantal vrijheidsgraden wordt weergegeven tussen haakjes in de formules uit de vijfde kolom. De grenswaarde waarmee vergeleken wordt, is het ${\displaystyle p}$-de percentiel (de waarde van de verdeling waarvoor p% van de stalen een hogere waarde heeft) van deze verdeling. In de Hotellings ${\displaystyle T^{2}}$-verdeling wordt verondersteld dat de populaties komen uit multivariate normaal verdeelde populaties. Zijn de varianties van de subgroepen niet gelijk, dan gebruikt men een gewogen covariantiematrix ${\displaystyle S_{W}}$:

${\displaystyle S_{W}={\frac {1}{n-k}}\sum _{i=1}^{k}(n_{i}-1)S_{i}}$,

met ${\displaystyle S_{i}}$ de covariantiematrix voor de groep ${\displaystyle i}$.

Grootte subgroep (${\displaystyle m}$)	Populatie-gemiddelde ${\displaystyle \mu _{T}}$	Covariantie	${\displaystyle T^{2}}$-statistiek	kritieke waarde voor ${\displaystyle T^{2}}$
1	Bekend	${\displaystyle S}$	${\displaystyle (X-\mu _{T})'S^{-1}(X-\mu _{T})}$	${\displaystyle {\frac {p(n-1)}{n-p}}F_{(p,n-p)}}$
		${\displaystyle S_{W}}$	${\displaystyle (X-\mu _{T})'S_{W}^{-1}(X-\mu _{T})}$	${\displaystyle {\frac {p(n-k)}{n-k-p+1}}F_{(p,n-k-p+1)}}$
	Onbekend	${\displaystyle S}$	${\displaystyle (X-{\overline {X}})'S^{-1}(X-{\overline {X}})}$	${\displaystyle {\frac {p(n+1)(n-1)}{n(n-p)}}F_{(p,n-p)}}$
		${\displaystyle S_{W}}$	${\displaystyle (X-{\overline {X}})'S_{W}^{-1}(X-{\overline {X}})}$	${\displaystyle {\frac {p(n-k)(n+1)}{n(n-k-p+1)}}F_{(p,n-k-p+1)}}$
> 1	Bekend	${\displaystyle S}$	${\displaystyle ({\overline {X}}_{s}-\mu _{T})'S^{-1}({\overline {X}}_{s}-\mu _{T})}$	${\displaystyle {\frac {p(n-1)}{m(n-p)}}F_{(p,n-p)}}$
		${\displaystyle S_{W}}$	${\displaystyle ({\overline {X}}_{s}-\mu _{T})'S_{W}^{-1}({\overline {X}}_{s}-\mu _{T})}$	${\displaystyle {\frac {p(n-k)}{m(n-k-p+1)}}F_{(p,n-k-p+1)}}$
	Onbekend	${\displaystyle S}$	${\displaystyle ({\overline {X}}_{s}-{\overline {X}})'S^{-1}({\overline {X}}_{s}-{\overline {X}})}$	${\displaystyle {\frac {p(m+n)(n-1)}{mn(n-p)}}F_{(p,n-p)}}$
		${\displaystyle S_{W}}$	${\displaystyle ({\overline {X}}_{s}-{\overline {X}})'S_{W}^{-1}({\overline {X}}_{s}-{\overline {X}})}$	${\displaystyle {\frac {p(m+n)(n-k)}{mn(n-k-p+1)}}F_{(p,n-k-p+1)}}$

Hierbij is:

${\displaystyle \mu _{T}}$ het gemiddelde van de totale populatie

${\displaystyle {\overline {X}}_{s}}$: het gemiddelde voor de groep stalen die men met de populatie vergelijkt

Verder hebben alle andere parameters in deze tabel dezelfde betekenis als vermeld bij de andere vergelijkingen.

Literatuur

• Mason, R.L.; Chou, Y.-M.; Young, J.C. Journal of Quality Technology, 2001, 33, 466-478

Zie ook

• Uitbijter