Hotelling's T-squared distribution
In statistics Hotelling's T-squared distribution is a univariate distribution proportional to the F-distribution and arises importantly as the distribution of a set of statistics which are natural generalizations of the statistics underlying Student's t-distribution. In particular, the distribution arises in multivariate statistics in undertaking tests of the differences between the (multivariate) means of different populations, where tests for univariate problems would make use of a t-test.
The distribution is named for Harold Hotelling, who developed it[1] as a generalization of Student's t-distribution.
Contents
The distribution
If the vector pd1 is Gaussian multivariate-distributed with zero mean and unit covariance matrix N(p01,pIp) and mMp is a p x p matrix with a Wishart distribution with unit scale matrix and m degrees of freedom W(pIp,m) then m(1d' pM−1pd1) has a Hotelling T2 distribution with dimensionality parameter p and m degrees of freedom.[2]
If the notation is used to denote a random variable having Hotelling's T-squared distribution with parameters p and m then, if a random variable X has Hotelling's T-squared distribution,
then[1]
where is the F-distribution with parameters p and m−p+1.
Hotelling's T-squared statistic
Hotelling's T-squared statistic is a generalization of Student's t statistic that is used in multivariate hypothesis testing, and is defined as follows.[1]
Let denote a p-variate normal distribution with location and covariance . Let
be n independent random variables, which may be represented as column vectors of real numbers. Define
to be the sample mean. It can be shown that
where is the chi-squared distribution with p degrees of freedom. To show this use the fact that and then derive the characteristic function of the random variable . This is done below,
However, is often unknown and we wish to do hypothesis testing on the location .
Sum of p squared t's
Define
to be the sample covariance. Here we denote transpose by an apostrophe. It can be shown that is a positive (semi) definite matrix and follows a p-variate Wishart distribution with n−1 degrees of freedom.[3] Hotelling's T-squared statistic is then defined[4] to be
and, also from above,
i.e.
where is the F-distribution with parameters p and n−p. In order to calculate a p value, multiply the t2 statistic by the above constant and use the F-distribution.
Hotelling's two-sample T-squared statistic
If and , with the samples independently drawn from two independent multivariate normal distributions with the same mean and covariance, and we define
as the sample means, and
as the unbiased pooled covariance matrix estimate, then Hotelling's two-sample T-squared statistic is
and it can be related to the F-distribution by[3]
The non-null distribution of this statistic is the noncentral F-distribution (the ratio of a non-central Chi-squared random variable and an independent central Chi-squared random variable)
with
where is the difference vector between the population means.
More robust and powerful tests than Hotelling's two-sample test have been proposed in the literature, see for example the interpoint distance based tests which can be applied also when the number of variables is comparable with, or even larger than, the number of subjects.[5][6]
In the two variable case, the formula simplifies nicely allowing appreciation of how the correlation, , between the variables affects . If we define
and
then
Thus, if the differences in the two rows of the vector are of the same sign, in general, becomes smaller as becomes more positive. If the differences are of opposite sign becomes larger as becomes more positive.
See also
- Student's t-test in univariate statistics
- Student's t-distribution in univariate probability theory
- Multivariate Student distribution.
- F-distribution (commonly tabulated or available in software libraries, and hence used for testing the T-squared statistic using the relationship given above)
- Wilks' lambda distribution (in multivariate statistics Wilks's Λ is to Hotelling's T2 as Snedecor's F is to Student's t in univariate statistics).
References
- ↑ 1.0 1.1 1.2 Lua error in package.lua at line 80: module 'strict' not found.
- ↑ Eric W. Weisstein, CRC Concise Encyclopedia of Mathematics, Second Edition, Chapman & Hall/CRC, 2003, p. 1408
- ↑ 3.0 3.1 Lua error in package.lua at line 80: module 'strict' not found.
- ↑ http://www.itl.nist.gov/div898/handbook/pmc/section5/pmc543.htm
- ↑ Lua error in package.lua at line 80: module 'strict' not found.
- ↑ Lua error in package.lua at line 80: module 'strict' not found.
External links
- Lua error in package.lua at line 80: module 'strict' not found.