Noncentral t-distribution

Short description: Probability distribution

Noncentral Student's t

Probability density function File:325px
Parameters	ν > 0 degrees of freedom $\mu \in \mathrm{\Re }$ noncentrality parameter
Support	$x\in (-\mathrm{\infty };+\mathrm{\infty })$
PDF	see text
CDF	see text
Mean	see text
Mode	see text
Variance	see text
Skewness	see text
Kurtosis	see text

The noncentral t-distribution generalizes Student's t-distribution using a noncentrality parameter. Whereas the central probability distribution describes how a test statistic t is distributed when the difference tested is null, the noncentral distribution describes how t is distributed when the null is false. This leads to its use in statistics, especially calculating statistical power. The noncentral t-distribution is also known as the singly noncentral t-distribution, and in addition to its primary use in statistical inference, is also used in robust modeling for data.

Definitions

If Z is a standard normal random variable, and V is a chi-squared distributed random variable with ν degrees of freedom that is independent of Z, then

$T=\frac{Z+\mu }{\sqrt{V/\nu }}$

is a noncentral t-distributed random variable with ν degrees of freedom and noncentrality parameter μ ≠ 0. Note that the noncentrality parameter may be negative.

Cumulative distribution function

The cumulative distribution function of noncentral t-distribution with ν degrees of freedom and noncentrality parameter μ can be defined using the integral

$F_{\nu ,\mu }(x)=\frac{1}{2^{\nu /2-1}\mathrm{\Gamma }(\nu /2)}\underset{0}{\overset{\mathrm{\infty }}{\int }}\mathrm{\Phi }(\frac{ux}{\sqrt{\nu }}-\mu )u^{\nu -1}{e}^{-\frac{1}{2}{u}^2}du,$

where $\mathrm{\Phi }(x)$ is the cumulative distribution function of the standard normal distribution. A more practical approach is to express it as

$F_{\nu ,\mu }(x)=\{\begin{array}{ll}{\tilde{F}}_{\nu ,\mu }(x),& \text{if }x\ge 0;\\ 1-{\tilde{F}}_{\nu ,-\mu }(x),& \text{if }x<0,\end{array}$

where

${\tilde{F}}_{\nu ,\mu }(x)=\mathrm{\Phi }(-\mu )+\frac{1}{2}\sum _{j=0}^{\mathrm{\infty }}[p_j{I}_y(j+\frac{1}{2},\frac{\nu }{2})+q_j{I}_y(j+1,\frac{\nu }{2})],$

$I_y(a,b)$ is the regularized incomplete beta function,

$y=\frac{x^2}{x^2+\nu },$

$p_j=\frac{1}{j!}\exp (-\frac{{\mu }^2}{2}){\left(\frac{{\mu }^2}{2}\right)}^j,$

$q_j=\frac{\mu }{\sqrt{2}\mathrm{\Gamma }(j+3/2)}\exp (-\frac{{\mu }^2}{2}){\left(\frac{{\mu }^2}{2}\right)}^j,$

Although this formula is more complicated, it is computationally more efficient, as each series term can be calculated from the previous term using simple arithmetic operations. In statistical programming language R, the cumulative distribution function is implemented as pt. In the Python library SciPy, it is implemented as scipy.stats.nct.cdf.

Probability density function

The probability density function (pdf) for the noncentral t-distribution with ν > 0 degrees of freedom and noncentrality parameter μ can be expressed in several forms. A simple integral form is given by:

$f(x)=\frac{2^{\frac{1-\nu }{2}}}{\sqrt{\nu \pi }\mathrm{\Gamma }(\nu /2)}\underset{0}{\overset{\mathrm{\infty }}{\int }}{u}^{\nu }\exp (-\frac{1}{2}{(\frac{ux}{\sqrt{\nu }}-\mu )}^2-\frac{1}{2}{u}^2)du$

The confluent hypergeometric function form of the density function is

$f(x)=\underset{\text{StudentT}(x;\mu =0)}{\underbrace{\frac{\mathrm{\Gamma }(\frac{\nu +1}{2})}{\sqrt{\nu \pi }\mathrm{\Gamma }(\frac{\nu }{2})}{(1+\frac{x^2}{\nu })}^{-\frac{\nu +1}{2}}}}\exp (-\frac{{\mu }^2}{2})(A_{\nu }(x;\mu )+B_{\nu }(x;\mu )),$

where

$\begin{array}{rl}{A}_{\nu }(x;\mu )& ={{}_{1}F}_1(\frac{\nu +1}{2};\frac{1}{2};\frac{{\mu }^2{x}^2}{2(x^2+\nu )}),\\ B_{\nu }(x;\mu )& =\frac{\sqrt{2}\mu x}{\sqrt{x^2+\nu }}\frac{\mathrm{\Gamma }(\frac{\nu }{2}+1)}{\mathrm{\Gamma }(\frac{\nu +1}{2})}{{}_{1}F}_1(\frac{\nu }{2}+1;\frac{3}{2};\frac{{\mu }^2{x}^2}{2(x^2+\nu )}),\end{array}$

and where ₁F₁ is a confluent hypergeometric function.

An alternative integral form is

$f(x)=\frac{{\nu }^{\frac{\nu }{2}}\exp (-\frac{\nu {\mu }^2}{2(x^2+\nu )})}{\sqrt{\pi }\mathrm{\Gamma }(\frac{\nu }{2})2^{\frac{\nu -1}{2}}(x^2+\nu {)}^{\frac{\nu +1}{2}}}\underset{0}{\overset{\mathrm{\infty }}{\int }}{y}^{\nu }\exp (-\frac{1}{2}{(y-\frac{\mu x}{\sqrt{x^2+\nu }})}^2)dy.$

A third form of the density is obtained using its cumulative distribution functions, as follows.

$f(x)=\{\begin{array}{ll}\frac{\nu }{x}\{F_{\nu +2,\mu }\left(x\sqrt{1+\frac{2}{\nu }}\right)-F_{\nu ,\mu }(x)\},& \text{if }x\ne 0;\\ \frac{\mathrm{\Gamma }(\frac{\nu +1}{2})}{\sqrt{\pi \nu }\mathrm{\Gamma }(\frac{\nu }{2})}\exp (-\frac{{\mu }^2}{2}),& \text{if }x=0.\end{array}$

This is the approach implemented by the dt function in R.

Derivation

The integral form of the CDF can be derived directly from the definition of a noncentral t variable:

$F_{\nu ,\mu }(x)=\mathrm{P}\{\frac{Z+\mu }{\sqrt{V/\nu }}<x\}=\mathrm{P}\{\frac{Z+\mu }{\chi /\sqrt{\nu }}<x\},$

where $\chi$ represents a chi random variable with $\nu$ degrees of freedom. Since $Z$ and $\chi$ are independent random variables, the above probability can be quantified as the double integral of the product of their probability density functions over the region where the inequality holds:

$\mathrm{P}\{\frac{Z+\mu }{\chi /\sqrt{\nu }}<x\}=\underset{R}{\iint }{f}_Z(z)f_{\chi }(u)dzdu=\underset{R}{\iint }\left(\frac{1}{\sqrt{2\pi }}{e}^{-\frac{1}{2}{z}^2}\right)\left(\frac{1}{2^{\nu /2-1}\mathrm{\Gamma }(\nu /2)}{u}^{\nu -1}{e}^{-\frac{1}{2}{u}^2}\right)dzdu,$

$\mathrm{w}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{e}R=\{(u,z)|u>0,\frac{z+\mu }{u/\sqrt{\nu }}<x\}=\{(u,z)|u>0,z<\frac{x}{\sqrt{\nu }}u-\mu \}$

Here, region $R$ represents the area in the $(u,z)$-plane lying to the right of the line $u=0$, and below the line $z=\frac{x}{\sqrt{\nu }}u-\mu$. Thus, the integral can be written with specified limits:

$F_{\nu ,\mu }(x)=\underset{0}{\overset{\mathrm{\infty }}{\int }}\underset{-\mathrm{\infty }}{\overset{\frac{x}{\sqrt{\nu }}u-\mu }{\int }}\left(\frac{1}{\sqrt{2\pi }}{e}^{-\frac{1}{2}{z}^2}\right)dz\left(\frac{1}{2^{\nu /2-1}\mathrm{\Gamma }(\nu /2)}{u}^{\nu -1}{e}^{-\frac{1}{2}{u}^2}\right)du$

Finally, integrating with respect to $z$ and separating the constant term yields:

$F_{\nu ,\mu }(x)=\frac{1}{2^{\nu /2-1}\mathrm{\Gamma }(\nu /2)}\underset{0}{\overset{\mathrm{\infty }}{\int }}\mathrm{\Phi }(\frac{ux}{\sqrt{\nu }}-\mu )u^{\nu -1}{e}^{-\frac{1}{2}{u}^2}du$.

The formula for the PDF can be found by differentiating the CDF with respect to $x$:

$f(x)=\frac{d}{dx}{F}_{\nu ,\mu }(x)=\frac{1}{2^{\nu /2-1}\mathrm{\Gamma }(\nu /2)}\underset{0}{\overset{\mathrm{\infty }}{\int }}\frac{d}{dx}\left[\mathrm{\Phi }(\frac{ux}{\sqrt{\nu }}-\mu )\right]u^{\nu -1}{e}^{-\frac{1}{2}{u}^2}du$

This simplifies to the first integral form of the PDF given above:

$f(x)=\frac{2^{\frac{1-\nu }{2}}}{\sqrt{\nu \pi }\mathrm{\Gamma }(\nu /2)}\underset{0}{\overset{\mathrm{\infty }}{\int }}{u}^{\nu }\exp (-\frac{1}{2}{(\frac{ux}{\sqrt{\nu }}-\mu )}^2-\frac{1}{2}{u}^2)du$

Properties

Moments of the noncentral t-distribution

In general, the kth raw moment of the noncentral t-distribution is

$\text{E}\left[T^k\right]=\{\begin{array}{ll}{\left(\frac{\nu }{2}\right)}^{\frac{k}{2}}\frac{\mathrm{\Gamma }\left(\frac{\nu -k}{2}\right)}{\mathrm{\Gamma }\left(\frac{\nu }{2}\right)}\text{exp}(-\frac{{\mu }^2}{2})\frac{d^k}{d{\mu }^k}\text{exp}\left(\frac{{\mu }^2}{2}\right),& \text{if }\nu >k;\\ \text{Does not exist},& \text{if }\nu \le k.\\ \end{array}$

In particular, the mean and variance of the noncentral t-distribution are

$\begin{array}{rl}\text{E}\left[T\right]& =\{\begin{array}{ll}\mu \sqrt{\frac{\nu }{2}}\frac{\mathrm{\Gamma }((\nu -1)/2)}{\mathrm{\Gamma }(\nu /2)},& \text{if }\nu >1;\\ \text{Does not exist},& \text{if }\nu \le 1,\\ \end{array}\\ \text{Var}\left[T\right]& =\{\begin{array}{ll}\frac{\nu (1+{\mu }^2)}{\nu -2}-\frac{{\mu }^2\nu }{2}{\left(\frac{\mathrm{\Gamma }((\nu -1)/2)}{\mathrm{\Gamma }(\nu /2)}\right)}^2,& \text{if }\nu >2;\\ \text{Does not exist},& \text{if }\nu \le 2.\\ \end{array}\end{array}$

An excellent approximation to $\sqrt{\frac{\nu }{2}}\frac{\mathrm{\Gamma }((\nu -1)/2)}{\mathrm{\Gamma }(\nu /2)}$ is ${(1-\frac{3}{4\nu -1})}^{-1}$, which can be used in both formulas.

Asymmetry

The non-central t-distribution is asymmetric unless μ is zero, i.e., a central t-distribution. In addition, the asymmetry becomes smaller the larger degree of freedom. The right tail will be heavier than the left when μ > 0, and vice versa. However, the usual skewness is not generally a good measure of asymmetry for this distribution, because if the degrees of freedom is not larger than 3, the third moment does not exist at all. Even if the degrees of freedom is greater than 3, the sample estimate of the skewness is still very unstable unless the sample size is very large.

Mode

The noncentral t-distribution is always unimodal and bell shaped, but the mode is not analytically available, although for μ ≠ 0 we have

$\sqrt{\frac{\nu }{\nu +(5/2)}}<\frac{\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{e}}{\mu }<\sqrt{\frac{\nu }{\nu +1}}$

In particular, the mode always has the same sign as the noncentrality parameter μ. Moreover, the negative of the mode is exactly the mode for a noncentral t-distribution with the same number of degrees of freedom ν but noncentrality parameter −μ.

The mode is strictly increasing with μ (it always moves in the same direction as μ is adjusted in). In the limit, when μ → 0, the mode is approximated by

$\sqrt{\frac{\nu }{2}}\frac{\mathrm{\Gamma }\left(\frac{\nu +2}{2}\right)}{\mathrm{\Gamma }\left(\frac{\nu +3}{2}\right)}\mu ;$

and when μ → ∞, the mode is approximated by

$\sqrt{\frac{\nu }{\nu +1}}\mu .$

Related distributions

• Central t-distribution: the central t-distribution can be converted into a location/scale family. This family of distributions is used in data modeling to capture various tail behaviors. The location/scale generalization of the central t-distribution is a different distribution from the noncentral t-distribution discussed in this article. In particular, this approximation does not respect the asymmetry of the noncentral t-distribution. However, the central t-distribution can be used as an approximation to the noncentral t-distribution.
• If T is noncentral t-distributed with ν degrees of freedom and noncentrality parameter μ and F = T², then F has a noncentral F-distribution with 1 numerator degree of freedom, ν denominator degrees of freedom, and noncentrality parameter μ².
• If T is noncentral t-distributed with ν degrees of freedom and noncentrality parameter μ and $Z=\underset{\nu \to \mathrm{\infty }}{\lim }T$, then Z has a normal distribution with mean μ and unit variance.
• When the denominator noncentrality parameter of a doubly noncentral t-distribution is zero, then it becomes a noncentral t-distribution.

Special cases

• When μ = 0, the noncentral t-distribution becomes the central (Student's) t-distribution with the same degrees of freedom.

Occurrence and applications

Use in power analysis

Suppose we have an independent and identically distributed sample X₁, ..., X_n each of which is normally distributed with mean θ and variance σ², and we are interested in testing the null hypothesis θ = 0 vs. the alternative hypothesis θ ≠ 0. We can perform a one sample t-test using the test statistic

$T=\frac{\overline{X}}{\hat{\sigma }/\sqrt{n}}=\frac{\frac{\overline{X}-\theta }{(\sigma /\sqrt{n})}+\frac{\theta }{(\sigma /\sqrt{n})}}{\sqrt{\left(\frac{{\hat{\sigma }}^2}{{\sigma }^2/(n-1)}\right)/(n-1)}}$

where $\overline{X}$ is the sample mean and ${\hat{\sigma }}^2$ is the unbiased sample variance. Since the right hand side of the second equality exactly matches the characterization of a noncentral t-distribution as described above, T has a noncentral t-distribution with n−1 degrees of freedom and noncentrality parameter $\sqrt{n}\theta /\sigma$.

If the test procedure rejects the null hypothesis whenever $|T|>t_{1-\alpha /2}$, where $t_{1-\alpha /2}$ is the upper α/2 quantile of the (central) Student's t-distribution for a pre-specified α ∈ (0, 1), then the power of this test is given by

$1-F_{n-1,\sqrt{n}\theta /\sigma }(t_{1-\alpha /2})+F_{n-1,\sqrt{n}\theta /\sigma }(-t_{1-\alpha /2}).$

Similar applications of the noncentral t-distribution can be found in the power analysis of the general normal-theory linear models, which includes the above one sample t-test as a special case.

Use in tolerance intervals

One-sided normal tolerance intervals have an exact solution in terms of the sample mean and sample variance based on the noncentral t-distribution. This enables the calculation of a statistical tolerance limit below which, with some confidence level, a specified proportion of a sampled population falls.

Let $X$ be a random normal variable with unknown population mean and variance, and let a sample of size $n$ yield sample mean $\bar{x}$ and sample standard deviation $s$. The probability that the proportion of $X$ values lying below limit $L$ is at least $p$ is given by the noncentral t-CDF:

$\mathrm{P}\{\mathrm{P}\{X<L\}>p\}=F_{\nu ,\mu }\left(\frac{L-\bar{x}}{s/\sqrt{n}}\right),\mathrm{w}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{e}\nu =n-1,\mu ={\mathrm{\Phi }}^{-1}(p)\sqrt{n}$

Tolerance bounds are typically expressed in the form $\bar{x}+ks$, where k is referred to as its k-factor. This factor is analogous to a z-score, but uses sample statistics rather than population parameters. In this framing, the probability that the proportion of $X$ lying below the limit $\bar{x}+ks$ is at least $p$ is given by

$\mathrm{P}\{\mathrm{P}\{X<\bar{x}+ks\}>p\}=F_{\nu ,\mu }(\sqrt{n}k),\mathrm{w}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{e}\nu =n-1,\mu ={\mathrm{\Phi }}^{-1}(p)\sqrt{n}$

Thus, calculating the k-factor corresponding to a given proportion and confidence level requires using the inverse cumulative distribution function of the noncentral t-distribution, which is implemented as qt in R, and as scipy.stats.nct.ppf in SciPy.

See also

• Noncentral F-distribution

External links

• Eric W. Weisstein. "Noncentral Student's t-Distribution." From MathWorld—A Wolfram Web Resource
• High accuracy calculation for life or science.: Noncentral t-distribution From Casio company.