Noncentral chi-square distribution

Noncentral chi-square

Probability density function 325px
Cumulative distribution function 325px
Parameters	<math>k > 0\,</math> degrees of freedom <math>\lambda > 0\,</math> non-centrality parameter
Support	<math>x \in [0; +\infty)\,</math>
Probability density function (pdf)	<math>\frac{1}{2}e^{-(x+\lambda)/2}\left (\frac{x}{\lambda} \right)^{k/4-1/2} I_{k/2-1}(\sqrt{\lambda x})</math>
Cumulative distribution function (cdf)	:<math>\sum_{j=0}^\infty e^{-\lambda/2} \frac{(\lambda/2)^j}{j!} \frac{\gamma(j+k/2,x/2)}{\Gamma(j+k/2)}\,</math>
Mean	<math>k+\lambda\,</math>
Median
Mode
Variance	<math>2(k+2\lambda)\,</math>
Skewness	<math>\frac{2^{3/2}(k+3\lambda)}{(k+2\lambda)^{3/2
Excess kurtosis	{{{kurtosis}}}
Entropy	{{{entropy}}}
Moment-generating function (mgf)	{{{mgf}}}
Characteristic function	{{{char}}}

</math>|

 kurtosis   =<math>\frac{12(k+4\lambda)}{(k+2\lambda)^2}</math>|
 entropy    =|
 mgf        =<math>\frac{\exp\left(\frac{ \lambda t}{1-2t }\right)}{(1-2 t)^{k/2}}</math>|
 char       =<math>\frac{\exp\left(\frac{i\lambda t}{1-2it}\right)}{(1-2it)^{k/2}}</math>

}}

In probability theory and statistics, the noncentral chi-square or noncentral <math>\chi^2</math> distribution is a generalization of the chi-square distribution. If <math>X_i</math> are k independent, normally distributed random variables with means <math>\mu_i</math> and variances <math>\sigma_i^2</math>, then the random variable

<math>\sum_1^k \left(\frac{X_i}{\sigma_i}\right)^2</math>

is distributed according to the noncentral chi-square distribution. The noncentral chi-square distribution has two parameters: <math>k</math> which specifies the number of degrees of freedom (i.e. the number of <math>X_i</math>), and <math>\lambda</math> which is related to the mean of the random variables <math>X_i</math> by:

<math>\lambda=\sum_1^k \left(\frac{\mu_i}{\sigma_i}\right)^2.</math>

Note that some references define <math>\lambda</math> as one half of the above sum.

Properties

The noncentral chi-square distribution is equivalent to a (central) chi-square distribution with <math>k + 2P</math> degrees of freedom, where <math>P</math> is a Poisson random variable with parameter <math>\lambda/2</math>. Thus, the probability density function is given by

<math>

f_X(x; k,\lambda) = \sum_{i=0}^\infty \frac{e^{-\lambda/2} (\lambda/2)^i}{i!} f_{Y_{k+2i}}(x), </math> where <math>Y_q</math> is distributed as chi-square with <math>q</math> degrees of freedom.

Alternatively, the pdf can be written as

<math>f_X(x;k,\lambda)=\frac{1}{2} e^{-(x+\lambda)/2} \left (\frac{x}{\lambda}\right)^{k/4-1/2} I_{k/2-1}(\sqrt{\lambda x})</math>

where <math>I_\nu(z)</math> is a modified Bessel function of the first kind given by

<math> I_a(y) := (y/2)^a \sum_{j=0}^\infty \frac{ (y^2/4)^j}{j! \Gamma(a+j+1)}</math>

The moment generating function is given by

<math>M(t;k,\lambda)=\frac{\exp\left(\frac{ \lambda t}{1-2t }\right)}{(1-2 t)^{k/2}}.</math>

The first few raw moments are:

<math>\mu^'_1=k+\lambda</math>

<math>\mu^'_2=(k+\lambda)^2 + 2(k + 2\lambda) </math>

<math>\mu^'_3=(k+\lambda)^3 + 6(k+\lambda)(k+2\lambda)+8(k+3\lambda)</math>

<math>\mu^'_4=(k+\lambda)^4+12(k+\lambda)^2(k+2\lambda)+4(11k^2+44k\lambda+36\lambda^2)+48(k+4\lambda)</math>

The first few central moments are:

<math>\mu_2=2(k+2\lambda)\,</math>

<math>\mu_3=8(k+3\lambda)\,</math>

<math>\mu_4=12(k+2\lambda)^2+48(k+4\lambda)\,</math>

The nth cumulant is

<math>K_n=2^{n-1}(n-1)!(k+n\lambda).\,</math>

Hence

<math>\mu^'_n = 2^{n-1}(n-1)!(k+n\lambda)+\sum_{j=1}^{n-1} \frac{(n-1)!2^{j-1}}{(n-j)!}(k+j\lambda )\mu^'_{n-j} </math>

Again using the relation between the central and noncentral chi-square distributions, the cumulative distribution function (cdf) can be written as

<math>P(x; k, \lambda ) = \sum_{j=0}^\infty e^{-\lambda/2} \frac{(\lambda/2)^j}{j!} Q(x; k+2j)</math>

where <math>Q(x; k)</math> is the cumulative density of the central chi-squared distribution which is given by

<math>Q(x;k)=\frac{\gamma(k/2,x/2)}{\Gamma(k/2)}\,</math>

where <math>\gamma(k,z)</math> is the lower incomplete Gamma function.

Derivation of the pdf

The derivation of the probability density function is most easily done by performing the following steps:

1. Start with the joint PDF of two independent non-zero mean Gaussian distributions, <math>X</math> and <math>Y</math>.
2. Convert the joint density <math>f(X, Y)</math> to polar: <math>f(R, A)</math> where <math>R^2 = (X^2+Y^2)</math>, <math>tan(A) = Y/X</math>.
3. Integrate over the angular variable <math>A</math>.
4. Convert from R to r where <math>r^2 = R</math>. This will yield a series expansion in r one factor of which matches the modified Bessel function <math>I_0</math>.
5. Take the Laplace (Fourier) transform term-by-term and the special case K = 2 and the MGF will result.
6. For the general case, take the K = 2 MGF and raise it to the <math>K/2</math> power.
7. The final trick to hide the K-dependence in the numerator of the MGF is to note that <math>\lambda</math> is a function of K; that is,

<math>\lambda_2=\sum_1^2 \left(\frac{\mu_i}{\sigma_i}\right)^2</math>

<math>\lambda_K=\sum_1^k \left(\frac{\mu_i}{\sigma_i}\right)^2 = \lambda</math>

and therefore, <math>\lambda</math> is not explicitly a function of K in the above table.

Related distributions

• If <math>Z</math> is chi-square distributed <math>Z \sim \chi_k^2</math> then <math>Z</math> is also non-central chi-square distributed: <math>Z \sim {\chi'}^2_k(0)</math>

• If <math>J \sim Poisson(\lambda/2)</math>, then <math>{\chi'}_k^2(\lambda) \sim \chi_{k+2J}^2</math>

Various chi and chi-square distributions

Name	Statistic
chi-square distribution	<math>\sum_1^k \left(\frac{X_i-\mu_i}{\sigma_i}\right)^2</math>
noncentral chi-square distribution	<math>\sum_1^k \left(\frac{X_i}{\sigma_i}\right)^2</math>
chi distribution	<math>\sqrt{\sum_1^k \left(\frac{X_i-\mu_i}{\sigma_i}\right)^2}</math>
noncentral chi distribution	<math>\sqrt{\sum_1^k \left(\frac{X_i}{\sigma_i}\right)^2}</math>

References

• Abramowitz, M. and Stegun, I.A. (1972), Handbook of Mathematical Functions, Dover. Section 26.4.25.
• Johnson, N. L. and Kotz, S., (1970), Continuous Univariate Distributions, vol. 2, Houghton-Mifflin.

Related Links

This distribution (and many others) is available in the free interactive statistical tables program, STATTAB. The cumulative distribution function, its inverse, and parameters of the distribution can be calculated with these packages. A free Fortran library for these distributions is in CDFLIB. The URL for download is [1]