Multivariate t distribution - Alchetron, the free social encyclopedia


Notation t ν ( μ , Σ ) {displaystyle t_{ u }({oldsymbol {mu }},{oldsymbol {Sigma }})} Parameters μ = [ μ 1 , … , μ p ] T {displaystyle {oldsymbol {mu }}=[mu _{1},dots ,mu _{p}]^{T}} location (real p × 1 {displaystyle p imes 1} vector) Σ {displaystyle {oldsymbol {Sigma }}} covariance matrix (positive-definite real p × p {displaystyle p imes p} matrix) ν {displaystyle u } is the degrees of freedom Support x ∈ R p {displaystyle mathbf {x} in mathbb {R} ^{p}!} PDF Γ [ ( ν + p ) / 2 ] Γ ( ν / 2 ) ν p / 2 π p / 2 \| Σ \| 1 / 2 [ 1 + 1 ν ( x − μ ) T Σ − 1 ( x − μ ) ] − ( ν + p ) / 2 {displaystyle {rac {Gamma left[( u +p)/2 ight]}{Gamma ( u /2) u ^{p/2}pi ^{p/2}left\|{oldsymbol {Sigma }} ight\|^{1/2}}}left[1+{rac {1}{ u }}({mathbf {x} }-{oldsymbol {mu }})^{ m {T}}{oldsymbol {Sigma }}^{-1}({mathbf {x} }-{oldsymbol {mu }}) ight]^{-( u +p)/2}} CDF No analytic expression, but see text for approximations Mean μ {displaystyle {oldsymbol {mu }}} if ν > 1 {displaystyle u >1} ; else undefined

In statistics, the multivariate t-distribution (or multivariate Student distribution) is a multivariate probability distribution. It is a generalization to random vectors of the Student's t-distribution, which is a distribution applicable to univariate random variables. While the case of a random matrix could be treated within this structure, the matrix t-distribution is distinct and makes particular use of the matrix structure.

Definition

One common method of construction of a multivariate t-distribution, for the case of p dimensions, is based on the observation that if y and u are independent and distributed as N ( 0 , Σ ) and χ ν 2 (i.e. multivariate normal and chi-squared distributions) respectively, the covariance Σ is a p × p matrix, and y / u / ν = x − μ , then x has the density

Γ [ ( ν + p ) / 2 ] Γ ( ν / 2 ) ν p / 2 π p / 2 | Σ | 1 / 2 [ 1 + 1 ν ( x − μ ) T Σ − 1 ( x − μ ) ] − ( ν + p ) / 2

and is said to be distributed as a multivariate t-distribution with parameters Σ , μ , ν .

In the special case ν = 1 , the distribution is a multivariate Cauchy distribution.

Derivation

There are in fact many candidates for the multivariate generalization of Student's t-distribution. An extensive survey of the field has been given by Kotz and Nadarajah (2004). The essential issue is to define a probability density function of several variables that is the appropriate generalization of the formula for the univariate case. In one dimension ( p = 1 ), with t = x − μ and Σ = 1 , we have the probability density function

f ( t ) = Γ [ ( ν + 1 ) / 2 ] ν π Γ [ ν / 2 ] ( 1 + t 2 / ν ) − ( ν + 1 ) / 2

and one approach is to write down a corresponding function of several variables. This is the basic idea of elliptical distribution theory, where one writes down a corresponding function of p variables t i that replaces t 2 by a quadratic function of all the t i . It is clear that this only makes sense when all the marginal distributions have the same degrees of freedom ν . With A = Σ − 1 , one has a simple choice of multivariate density function

f ( t ) = Γ ( ( ν + p ) / 2 ) | A | 1 / 2 ν p π p Γ ( ν / 2 ) ( 1 + ∑ i , j = 1 p , p A i j t i t j / ν ) − ( ν + p ) / 2

which is the standard but not the only choice.

An important special case is the standard bivariate t-distribution, p = 2:

f ( t 1 , t 2 ) = | A | 1 / 2 2 π ( 1 + ∑ i , j = 1 2 , 2 A i j t i t j / ν ) − ( ν + 2 ) / 2

Note that Γ ( ν + 2 2 ) π ν Γ ( ν 2 ) = 1 2 π .

Now, if A is the identity matrix, the density is

f ( t 1 , t 2 ) = 1 2 π ( 1 + ( t 1 2 + t 2 2 ) / ν ) − ( ν + 2 ) / 2 .

The difficulty with the standard representation is revealed by this formula, which does not factorize into the product of the marginal one-dimensional distributions. When Σ is diagonal the standard representation can be shown to have zero correlation but the marginal distributions do not agree with statistical independence.

Cumulative distribution function

The definition of the cumulative distribution function (cdf) in one dimension can be extended to multiple dimensions by defining the following probability (here x is a real vector):

F ( x ) = P ( X ≤ x ) , where X ∼ t ν ( μ , Σ ) .

There is no simple formula for F ( x ) , but it can be approximated numerically via Monte Carlo integration.

Further theory

Many such distributions may be constructed by considering the quotients of normal random variables with the square root of a sample from a chi-squared distribution. These are surveyed in the references and links below.

Copulas based on the multivariate t

The use of such distributions is enjoying renewed interest due to applications in mathematical finance, especially through the use of the Student's t copula.

In univariate statistics, the Student's t-test makes use of Student's t-distribution. Hotelling's T-squared distribution is a distribution that arises in multivariate statistics. The matrix t-distribution is a distribution for random variables arranged in a matrix structure.

References

Multivariate t-distribution Wikipedia

(Text) CC BY-SA

Contents