Chi-squared distribution: Difference between revisions

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search
No edit summary
 
imported>Nrco0e
rmv repeated link of gamma distribution
 
Line 10: Line 10:
   | parameters = <math>k \in \mathbb{N}^{*}~~</math> (known as "degrees of freedom")
   | parameters = <math>k \in \mathbb{N}^{*}~~</math> (known as "degrees of freedom")
   | support    = <math>x \in (0, +\infty)\;</math>
   | support    = <math>x \in (0, +\infty)\;</math>
   | pdf        = <math>\frac{1}{2^{k/2}\Gamma(k/2)}\; x^{k/2-1} e^{-x/2}\; </math>
   | pdf        = <math>\frac{1}{2^{k/2}\Gamma(k/2)}\; x^{(k/2)-1} e^{-x/2}\; </math>
   | cdf        = <math>\frac{1}{\Gamma(k/2 )} \; \gamma\left(\frac{k}{2},\,\frac{x}{2}\right)\;</math>
   | cdf        = <math>\frac{1}{\Gamma(k/2 )} \; \gamma{\left(\frac{k}{2},\,\frac{x}{2}\right)}\;</math>
   | mean      = <math>k</math>
   | mean      = <math>k</math>
   | median    = <math>\approx k\bigg(1-\frac{2}{9k}\bigg)^3\;</math>
   | median    = <math>\approx k\bigg(1-\frac{2}{9k}\bigg)^3\;</math>
   | mode      = <math>\max(k-2,0)\;</math>
   | mode      = <math>\max(k-2,0)\;</math>
   | variance  = <math>2k\;</math>
   | variance  = <math>2k\;</math>
   | skewness  = <math>\sqrt{8/k}\,</math>
   | skewness  = <math display="inline">\sqrt{8/k}\,</math>
   | kurtosis  = <math>\frac{12}{k}</math>
   | kurtosis  = <math>\frac{12}{k}</math>
   | entropy    = <math>\begin{align}\frac{k}{2}&+\log\left(2\Gamma\Bigl(\frac{k}{2}\Bigr)\right) \\ &\!+\left(1-\frac{k}{2}\right)\psi\left(\frac{k}{2}\right)\end{align}</math>
   | entropy    = <math>\begin{align} \frac{k}{2}& + \log\left(2\Gamma{\left(\frac{k}{2}\right)}\right) \\ &\!+\left(1-\frac{k}{2}\right) \psi{\left(\frac{k}{2}\right)} \end{align}</math>
   | mgf        = <math>(1-2t)^{-k/2} \text{ for } t < \frac{1}{2}\;</math>
   | mgf        = <math>(1-2t)^{-k/2} </math>{{quad}} for <math> t < \tfrac{1}{2}\;</math>
   | char      = <math>(1-2it)^{-k/2}</math><ref>{{cite web | url=http://www.planetmathematics.com/CentralChiDistr.pdf | title=Characteristic function of the central chi-square distribution | author=M.A. Sanders | access-date=2009-03-06 | archive-url=https://web.archive.org/web/20110715091705/http://www.planetmathematics.com/CentralChiDistr.pdf | archive-date=2011-07-15 | url-status=dead }}</ref>
   | char      = <math>(1-2it)^{-k/2}</math><ref>{{cite web | url=http://www.planetmathematics.com/CentralChiDistr.pdf | title = Characteristic function of the central chi-square distribution | author=M.A. Sanders | access-date=2009-03-06 | archive-url=https://web.archive.org/web/20110715091705/http://www.planetmathematics.com/CentralChiDistr.pdf | archive-date=2011-07-15 }}</ref>
|pgf=<math>(1-2\ln t)^{-k/2} \text{ for } 0<t<\sqrt{e}\;</math>}}
  | pgf       = <math>(1-2\ln t)^{-k/2} </math>{{quad}} for <math> 0 < t < \sqrt{e}\;</math>
}}


In [[probability theory]] and [[statistics]], the '''<math>\chi^2</math>-distribution''' with <math>k</math> [[Degrees of freedom (statistics)|degrees of freedom]] is the distribution of a sum of the squares of <math>k</math> [[Independence (probability theory)|independent]] [[standard normal]] random variables.<ref>{{Cite web |last=Weisstein |first=Eric W. |title=Chi-Squared Distribution |url=https://mathworld.wolfram.com/Chi-SquaredDistribution.html |access-date=2024-10-11 |website=mathworld.wolfram.com |language=en}}</ref>
In [[probability theory]] and [[statistics]], the '''<math>\chi^2</math>-distribution''' with <math>k</math> [[Degrees of freedom (statistics)|degrees of freedom]] is the distribution of a sum of the squares of <math>k</math> [[Independence (probability theory)|independent]] [[standard normal]] random variables.<ref>{{Cite web |last=Weisstein |first=Eric W. |title=Chi-Squared Distribution |url=https://mathworld.wolfram.com/Chi-SquaredDistribution.html |access-date=2024-10-11 |website=mathworld.wolfram.com |language=en}}</ref>


The chi-squared distribution <math> \chi^2_k </math> is a special case of the [[gamma distribution]] and the univariate [[Wishart distribution]].  Specifically if
The chi-squared distribution <math> \chi^2_k </math> is a special case of the [[gamma distribution]] and the univariate [[Wishart distribution]].  Specifically if <math> X \sim \chi^2_k </math> then <math display="inline"> X \sim \text{Gamma}(\alpha=\frac{k}{2}, \theta=2) </math> (where <math>\alpha</math> is the shape parameter and <math>\theta</math> the scale parameter of the gamma distribution) and <math> X \sim \text{W}_1(1,k) </math>.
<math> X \sim \chi^2_k </math> then <math> X \sim \text{Gamma}(\alpha=\frac{k}{2}, \theta=2) </math> (where <math>\alpha</math> is the shape parameter and <math>\theta</math> the scale parameter of the gamma distribution) and <math> X \sim \text{W}_1(1,k) </math>.


The '''scaled chi-squared distribution''' <math>s^2 \chi^2_k </math> is a reparametrization of the [[gamma distribution]] and the univariate [[Wishart distribution]].  Specifically if
The '''scaled chi-squared distribution''' <math>s^2 \chi^2_k </math> is a reparametrization of the gamma distribution and the univariate [[Wishart distribution]].  Specifically if <math> X  \sim s^2 \chi^2_k </math> then <math display="inline"> X \sim \text{Gamma}(\alpha=\frac{k}{2}, \theta=2 s^2) </math> and <math> X \sim \text{W}_1(s^2,k) </math>.
<math> X  \sim s^2 \chi^2_k </math> then <math> X \sim \text{Gamma}(\alpha=\frac{k}{2}, \theta=2 s^2) </math> and <math> X \sim \text{W}_1(s^2,k) </math>.


The chi-squared distribution is one of the most widely used [[probability distribution]]s in [[inferential statistics]], notably in [[hypothesis testing]] and in construction of [[confidence interval]]s.<ref name="United States Department of Commerce, National Bureau of Standards; Dover Publications-1983">{{Abramowitz Stegun ref|26|940}}</ref><ref>NIST (2006). [http://www.itl.nist.gov/div898/handbook/eda/section3/eda3666.htm Engineering Statistics Handbook – Chi-Squared Distribution]</ref><ref name="Johnson-1994">{{cite book
The chi-squared distribution is one of the most widely used [[probability distribution]]s in [[inferential statistics]], notably in [[hypothesis testing]] and in construction of [[confidence interval]]s.<ref name="United States Department of Commerce, National Bureau of Standards; Dover Publications-1983">{{Abramowitz Stegun ref|26|940}}</ref><ref>NIST (2006). [http://www.itl.nist.gov/div898/handbook/eda/section3/eda3666.htm Engineering Statistics Handbook – Chi-Squared Distribution]</ref><ref name="Johnson-1994">{{cite book | last1 = Johnson | first1 = N. L. | first2 = S. | last2 = Kotz | first3 = N. | last3 = Balakrishnan | title = Continuous Univariate Distributions | edition = Second | volume = 1 | chapter = Chi-Square Distributions including Chi and Rayleigh | pages = 415–493 | publisher = John Wiley and Sons | year = 1994 | isbn = 978-0-471-58495-7 }}</ref><ref>{{cite book
  | last1 = Johnson
  | first1 = N. L.
  | first2 = S. |last2=Kotz |first3=N. |last3=Balakrishnan
  | title = Continuous Univariate Distributions |edition=Second |volume=1 |chapter=Chi-Square Distributions including Chi and Rayleigh |pages=415–493
  | publisher = John Wiley and Sons
  | year = 1994
  | isbn = 978-0-471-58495-7
}}</ref><ref>{{cite book
   | last1 = Mood
   | last1 = Mood
   | first1 = Alexander
   | first1 = Alexander
Line 53: Line 44:
== Definitions ==
== Definitions ==
If {{math|''Z''<sub>1</sub>, ..., ''Z''<sub>''k''</sub>}} are [[independence (probability theory)|independent]], [[standard normal]] random variables, then the sum of their squares,
If {{math|''Z''<sub>1</sub>, ..., ''Z''<sub>''k''</sub>}} are [[independence (probability theory)|independent]], [[standard normal]] random variables, then the sum of their squares,
: <math>X\ = \sum_{i=1}^k Z_i^2,</math>
<math display="block">X\ = \sum_{i=1}^k Z_i^2,</math>
is distributed according to the chi-squared distribution with {{mvar|k}} degrees of freedom. This is usually denoted as
is distributed according to the chi-squared distribution with {{mvar|k}} degrees of freedom. This is usually denoted as
: <math> X\ \sim\ \chi^2(k)\ \ \text{or}\ \ X\ \sim\ \chi^2_k.</math>
<math display="block"> X\ \sim\ \chi^2(k)\ \ \text{or}\ \ X\ \sim\ \chi^2_k.</math>


The chi-squared distribution has one parameter: a positive integer {{mvar|k}} that specifies the number of [[degrees of freedom (statistics)|degrees of freedom]] (the number of random variables being summed, ''Z''<sub>''i''</sub> s).
The chi-squared distribution has one parameter: a positive integer {{mvar|k}} that specifies the number of [[degrees of freedom (statistics)|degrees of freedom]] (the number of random variables ''Z''<sub>''i''</sub> being summed).


=== Introduction ===
=== Introduction ===
Line 78: Line 69:
An additional reason that the chi-squared distribution is widely used is that it turns up as the large sample distribution of generalized [[Likelihood-ratio test|likelihood ratio tests]] (LRT).<ref name="Westfall-2013">{{cite book|last1=Westfall|first1=Peter H.|title=Understanding Advanced Statistical Methods|date=2013|publisher=CRC Press|location=Boca Raton, FL|isbn=978-1-4665-1210-8}}</ref> LRTs have several desirable properties; in particular, simple LRTs commonly provide the highest power to reject the null hypothesis ([[Neyman–Pearson lemma]]) and this leads also to optimality properties of generalised LRTs. However, the normal and chi-squared approximations are only valid asymptotically. For this reason, it is preferable to use the ''t'' distribution rather than the normal approximation or the chi-squared approximation for a small sample size. Similarly, in analyses of contingency tables, the chi-squared approximation will be poor for a small sample size, and it is preferable to use [[Fisher's exact test]]. Ramsey shows that the exact [[binomial test]] is always more powerful than the normal approximation.<ref name="Ramsey-1988">{{cite journal|last1=Ramsey|first1=PH|title=Evaluating the Normal Approximation to the Binomial Test|journal=Journal of Educational Statistics|date=1988|volume=13|issue=2|pages=173–82|doi=10.2307/1164752|jstor=1164752}}</ref>
An additional reason that the chi-squared distribution is widely used is that it turns up as the large sample distribution of generalized [[Likelihood-ratio test|likelihood ratio tests]] (LRT).<ref name="Westfall-2013">{{cite book|last1=Westfall|first1=Peter H.|title=Understanding Advanced Statistical Methods|date=2013|publisher=CRC Press|location=Boca Raton, FL|isbn=978-1-4665-1210-8}}</ref> LRTs have several desirable properties; in particular, simple LRTs commonly provide the highest power to reject the null hypothesis ([[Neyman–Pearson lemma]]) and this leads also to optimality properties of generalised LRTs. However, the normal and chi-squared approximations are only valid asymptotically. For this reason, it is preferable to use the ''t'' distribution rather than the normal approximation or the chi-squared approximation for a small sample size. Similarly, in analyses of contingency tables, the chi-squared approximation will be poor for a small sample size, and it is preferable to use [[Fisher's exact test]]. Ramsey shows that the exact [[binomial test]] is always more powerful than the normal approximation.<ref name="Ramsey-1988">{{cite journal|last1=Ramsey|first1=PH|title=Evaluating the Normal Approximation to the Binomial Test|journal=Journal of Educational Statistics|date=1988|volume=13|issue=2|pages=173–82|doi=10.2307/1164752|jstor=1164752}}</ref>


Lancaster shows the connections among the binomial, normal, and chi-squared distributions, as follows.<ref name="Lancaster-1969">{{Citation
Lancaster shows the connections among the binomial, normal, and chi-squared distributions, as follows.<ref name="Lancaster-1969">{{ citation | last = Lancaster | first = H.O. | title = The Chi-squared Distribution | year = 1969 | publisher = Wiley }}</ref> De Moivre and Laplace established that a binomial distribution could be approximated by a normal distribution. Specifically they showed the asymptotic normality of the random variable
|last=Lancaster
|first=H.O.
|title=The Chi-squared Distribution
|year=1969
|publisher=Wiley
}}</ref> De Moivre and Laplace established that a binomial distribution could be approximated by a normal distribution. Specifically they showed the asymptotic normality of the random variable


:<math> \chi = {m - Np \over \sqrt{Npq}} </math>
<math display="block"> \chi = \frac{m - Np}{\sqrt{Npq}} </math>


where <math>m</math> is the observed number of successes in <math>N</math> trials, where the probability of success is <math>p</math>, and <math>q = 1 - p</math>.
where <math>m</math> is the observed number of successes in <math>N</math> trials, where the probability of success is <math>p</math>, and <math>q = 1 - p</math>.
Line 92: Line 77:
Squaring both sides of the equation gives
Squaring both sides of the equation gives


: <math style="block"> \chi^2 = {(m - Np)^2\over Npq} </math>
<math display="block" style="block"> \chi^2 = \frac{\left(m - Np\right)^2}{Npq} </math>


Using <math>N = Np + N(1 - p)</math>, <math>N = m + (N - m)</math>, and <math>q = 1 - p</math>, this equation can be rewritten as
Using <math>N = Np + N(1 - p)</math>, <math>N = m + (N - m)</math>, and <math>q = 1 - p</math>, this equation can be rewritten as


: <math style="block"> \chi^2 = {(m - Np)^2\over Np} + {(N - m - Nq)^2\over Nq} </math>
<math display="block" style="block"> \chi^2 = \frac{\left(m - Np\right)^2}{Np} + \frac{\left(N - m - Nq\right)^2}{Nq} </math>


The expression on the right is of the form that [[Karl Pearson]] would generalize to the form
The expression on the right is of the form that [[Karl Pearson]] would generalize to the form


: <math style="block"> \chi^2 = \sum_{i=1}^n \frac{(O_i - E_i)^2}{E_i} </math>
<math display="block" style="block"> \chi^2 = \sum_{i=1}^n \frac{\left(O_i - E_i\right)^2}{E_i} </math>


where
where


<math style="block"> \chi^2</math> = Pearson's cumulative test statistic, which asymptotically approaches a <math>\chi^2</math> distribution;
* <math style="block"> \chi^2</math> = Pearson's cumulative test statistic, which asymptotically approaches a <math>\chi^2</math> distribution;
<math style="block">O_i</math> = the number of observations of type <math>i</math>;
* <math style="block">O_i</math> = the number of observations of type <math>i</math>;
<math style="block">E_i = N p_i</math> = the expected (theoretical) frequency of type <math>i</math>, asserted by the null hypothesis that the fraction of type <math>i</math> in the population is <math> p_i</math>; and
* <math style="block">E_i = N p_i</math> = the expected (theoretical) frequency of type <math>i</math>, asserted by the null hypothesis that the fraction of type <math>i</math> in the population is <math> p_i</math>; and
<math style="block">n</math> = the number of cells in the table.{{cn|date=November 2023}}
* <math style="block">n</math> = the number of cells in the table.{{cn|date=November 2023}}


In the case of a binomial outcome (flipping a coin), the binomial distribution may be approximated by a normal distribution (for sufficiently large <math>n</math>). Because the square of a standard normal distribution is the chi-squared distribution with one degree of freedom, the probability of a result such as 1 heads in 10 trials can be approximated either by using the normal distribution directly, or the chi-squared distribution for the normalised, squared difference between observed and expected value. However, many problems involve more than the two possible outcomes of a binomial, and instead require 3 or more categories, which leads to the multinomial distribution. Just as de Moivre and Laplace sought for and found the normal approximation to the binomial, Pearson sought for and found a degenerate multivariate normal approximation to the multinomial distribution (the numbers in each category add up to the total sample size, which is considered fixed). Pearson showed that the chi-squared distribution arose from such a multivariate normal approximation to the multinomial distribution, taking careful account of the statistical dependence (negative correlations) between numbers of observations in different categories.<ref name="Lancaster-1969" />
In the case of a binomial outcome (flipping a coin), the binomial distribution may be approximated by a normal distribution (for sufficiently large <math>n</math>). Because the square of a standard normal distribution is the chi-squared distribution with one degree of freedom, the probability of a result such as 1 heads in 10 trials can be approximated either by using the normal distribution directly, or the chi-squared distribution for the normalised, squared difference between observed and expected value. However, many problems involve more than the two possible outcomes of a binomial, and instead require 3 or more categories, which leads to the multinomial distribution. Just as de Moivre and Laplace sought for and found the normal approximation to the binomial, Pearson sought for and found a degenerate multivariate normal approximation to the multinomial distribution (the numbers in each category add up to the total sample size, which is considered fixed). Pearson showed that the chi-squared distribution arose from such a multivariate normal approximation to the multinomial distribution, taking careful account of the statistical dependence (negative correlations) between numbers of observations in different categories.<ref name="Lancaster-1969" />
Line 113: Line 98:
=== Probability density function ===
=== Probability density function ===
The [[probability density function]] (pdf) of the chi-squared distribution is
The [[probability density function]] (pdf) of the chi-squared distribution is
:<math>
<math display="block">
f(x;\,k) =
f(x;\,k) = \begin{cases}
\begin{cases}
   \dfrac{x^{k/2 -1} e^{-x/2}}{2^{k/2} \, \Gamma{\left(\frac k 2 \right)}}, & x > 0; \\
   \dfrac{x^{k/2 -1} e^{-x/2}}{2^{k/2} \Gamma\left(\frac k 2 \right)}, & x > 0; \\ 0, & \text{otherwise}.
  0, & \text{otherwise}.
\end{cases}
\end{cases}
</math>
</math>
Line 128: Line 113:


Its [[cumulative distribution function]] is:
Its [[cumulative distribution function]] is:
: <math>
<math display="block">
     F(x;\,k) = \frac{\gamma(\frac{k}{2},\,\frac{x}{2})}{\Gamma(\frac{k}{2})} = P\left(\frac{k}{2},\,\frac{x}{2}\right),
     F(x;\,k) = \frac{\gamma{\left(\frac{k}{2},\,\frac{x}{2}\right)}}{\Gamma{\left(\frac{k}{2}\right)}} = P{\left(\frac{k}{2},\,\frac{x}{2}\right)},
   </math>
   </math>
where <math>\gamma(s,t)</math> is the [[lower incomplete gamma function]] and <math display="inline">P(s,t)</math> is the [[Incomplete gamma function#Regularized gamma functions and Poisson random variables|regularized gamma function]].
where <math>\gamma(s,t)</math> is the [[lower incomplete gamma function]] and <math display="inline">P(s,t)</math> is the [[Incomplete gamma function#Regularized gamma functions and Poisson random variables|regularized gamma function]].


In a special case of <math>k = 2</math> this function has the simple form:
In a special case of <math>k = 2</math> this function has the simple form:
: <math>
<math display="block">
     F(x;\,2) = 1 - e^{-x/2}
     F(x;\,2) = 1 - e^{-x/2}
   </math>
   </math>
which can be easily derived by integrating <math>f(x;\,2)=\frac{1}{2}e^{-x/2}</math> directly. The integer recurrence of the gamma function makes it easy to compute <math>F(x;\,k)</math> for other small, even <math>k</math>.
which can be easily derived by integrating <math display="inline">f(x;\,2) = \frac{1}{2}e^{-x/2}</math> directly. The integer recurrence of the gamma function makes it easy to compute <math>F(x;\,k)</math> for other small, even <math>k</math>.


Tables of the chi-squared cumulative distribution function are widely available and the function is included in many [[spreadsheet]]s and all [[statistical packages]].
Tables of the chi-squared cumulative distribution function are widely available and the function is included in many [[spreadsheet]]s and all [[statistical packages]].
Line 145: Line 130:


The tail bound for the cases when <math>z > 1</math>, similarly, is
The tail bound for the cases when <math>z > 1</math>, similarly, is
: <math>
<math display="block">
     1-F(z k;\,k) \leq (z e^{1-z})^{k/2}.
     1-F(z k;\,k) \leq (z e^{1-z})^{k/2}.
   </math>
   </math>
Line 158: Line 143:


'''Theorem.''' If <math>Z_1,...,Z_n</math> are [[independence (probability theory)|independent]] identically distributed (i.i.d.), [[standard normal]] random variables, then
'''Theorem.''' If <math>Z_1,...,Z_n</math> are [[independence (probability theory)|independent]] identically distributed (i.i.d.), [[standard normal]] random variables, then
<math>\sum_{t=1}^n(Z_t - \bar Z)^2 \sim \chi^2_{n-1}</math>
<math display="inline">\sum_{t=1}^n \left(Z_t - \bar Z\right)^2 \sim \chi^2_{n-1}</math>
where <math>\bar Z = \frac{1}{n} \sum_{t=1}^n Z_t.</math>
where <math display="inline">\bar Z = \frac{1}{n} \sum_{t=1}^n Z_t.</math>


{{hidden begin|style=width:100%|ta1=center|border=1px #aaa solid|title=[Proof]}}
{{hidden begin|style=width:100%|ta1=center|border=1px #aaa solid|title=[Proof]}}
Line 185: Line 170:
=== Sample mean ===
=== Sample mean ===
The sample mean of <math>n</math> [[i.i.d.]] chi-squared variables of degree <math>k</math> is distributed according to a gamma distribution with shape <math>\alpha</math> and scale <math>\theta</math> parameters:
The sample mean of <math>n</math> [[i.i.d.]] chi-squared variables of degree <math>k</math> is distributed according to a gamma distribution with shape <math>\alpha</math> and scale <math>\theta</math> parameters:
:<math> \overline X = \frac{1}{n} \sum_{i=1}^n X_i \sim \operatorname{Gamma}\left(\alpha=n\, k /2, \theta= 2/n \right) \qquad \text{where } X_i \sim \chi^2(k)</math>
<math display="block"> \overline X = \frac{1}{n} \sum_{i=1}^n X_i \sim \operatorname{Gamma}\left(\alpha{=}\tfrac{n k}{2}, \,\theta{=}\tfrac{2}{n} \right) \qquad \text{where } X_i \sim \chi^2(k)</math>


[[#Asymptotic properties|Asymptotically]], given that for a shape parameter <math> \alpha </math> going to infinity, a Gamma distribution converges towards a normal distribution with expectation <math> \mu = \alpha\cdot \theta </math> and variance <math> \sigma^2 = \alpha\, \theta^2 </math>, the sample mean converges towards:
[[#Asymptotic properties|Asymptotically]], given that for a shape parameter <math> \alpha </math> going to infinity, a Gamma distribution converges towards a normal distribution with expectation <math> \mu = \alpha \theta </math> and variance {{nowrap|<math> \sigma^2 = \alpha \theta^2 </math>,}} the sample mean converges towards:


<math style="block"> \overline X \xrightarrow{n \to \infty} N(\mu = k, \sigma^2 = 2\, k /n ) </math>
<math display="block">\overline X \xrightarrow{n \to \infty} N{\left(\mu{=}k, \, \sigma^2{=}\tfrac{2k}{n} \right)} </math>


Note that we would have obtained the same result invoking instead the [[central limit theorem]], noting that for each chi-squared variable of degree <math>k</math> the expectation is <math> k </math> , and its variance <math> 2\,k </math> (and hence the variance of the sample mean <math> \overline{X}</math> being <math> \sigma^2 = \frac{2k}{n} </math>).
Note that we would have obtained the same result invoking instead the [[central limit theorem]], noting that for each chi-squared variable of degree <math>k</math> the expectation is {{nowrap|<math> k </math>,}} and its variance <math> 2k </math> (and hence the variance of the sample mean <math> \overline{X}</math> being <math display="inline"> \sigma^2 = \tfrac{2k}{n} </math>).


=== Entropy ===
=== Entropy ===
The [[differential entropy]] is given by
The [[differential entropy]] is given by
: <math>
<math display="block">
     h = \int_{0}^\infty f(x;\,k)\ln f(x;\,k) \, dx
     \begin{align}
      = \frac k 2 + \ln \left[2\,\Gamma \left(\frac k 2 \right)\right] + \left(1-\frac k 2 \right)\, \psi\!\left(\frac k 2 \right),
h &= \int_0^\infty f(x;\,k) \ln f(x;\,k) \, dx \\
  &= \frac k 2 + \ln \left[2\,\Gamma{\left(\frac k 2 \right)}\right] + \left(1-\frac k 2 \right) \psi\!\left(\frac k 2 \right),
\end{align}
   </math>
   </math>
where <math>\psi(x)</math> is the [[Digamma function]].
where <math>\psi(x)</math> is the [[Digamma function]].
Line 205: Line 192:
=== Noncentral moments ===
=== Noncentral moments ===
The noncentral moments (raw moments) of a chi-squared distribution with <math>k</math> degrees of freedom are given by<ref>[http://mathworld.wolfram.com/Chi-SquaredDistribution.html Chi-squared distribution], from [[MathWorld]], retrieved Feb. 11, 2009</ref><ref>M. K. Simon, ''Probability Distributions Involving Gaussian Random Variables'', New York: Springer, 2002, eq. (2.35), {{ISBN|978-0-387-34657-1}}</ref>
The noncentral moments (raw moments) of a chi-squared distribution with <math>k</math> degrees of freedom are given by<ref>[http://mathworld.wolfram.com/Chi-SquaredDistribution.html Chi-squared distribution], from [[MathWorld]], retrieved Feb. 11, 2009</ref><ref>M. K. Simon, ''Probability Distributions Involving Gaussian Random Variables'', New York: Springer, 2002, eq. (2.35), {{ISBN|978-0-387-34657-1}}</ref>
: <math>
<math display="block">
     \operatorname{E}(X^m) = k (k+2) (k+4) \cdots (k+2m-2) = 2^m \frac{\Gamma\left(m+\frac{k}{2}\right)}{\Gamma\left(\frac{k}{2}\right)}.
\begin{align}
  </math>
     \operatorname{E}(X^m) &= k (k+2) (k+4) \cdots (k+2m-2) \\[1ex]
&= 2^m \frac{\Gamma{\left(m+\frac{k}{2}\right)}}{\Gamma{\left(\frac{k}{2}\right)}}.
\end{align} </math>


=== Cumulants ===
=== Cumulants ===
The [[cumulant]]s are readily obtained by a [[power series]] expansion of the logarithm of the characteristic function:
The [[cumulant]]s are readily obtained by a [[power series]] expansion of the logarithm of the characteristic function:
: <math>\kappa_n = 2^{n-1}(n-1)!\,k</math>
<math display="block">\kappa_n = 2^{n-1}(n-1)!\,k</math>
with [[cumulant generating function]] <math>\ln E[e^{tX}] = - \frac k2 \ln(1-2t) </math>.
with [[cumulant generating function]] <math display="inline">\ln \operatorname{E}[e^{tX}] = - \frac{k}{2} \ln(1-2t) </math>.


=== Concentration ===
=== Concentration ===


The chi-squared distribution exhibits strong concentration around its mean. The standard Laurent-Massart<ref>{{Cite journal |last1=Laurent |first1=B. |last2=Massart |first2=P. |date=2000-10-01 |title=Adaptive estimation of a quadratic functional by model selection |journal=The Annals of Statistics |volume=28 |issue=5 |doi=10.1214/aos/1015957395 |s2cid=116945590 |issn=0090-5364|doi-access=free }}</ref> bounds are:
The chi-squared distribution exhibits strong concentration around its mean. The standard Laurent-Massart<ref>{{Cite journal |last1=Laurent |first1=B. |last2=Massart |first2=P. |date=2000-10-01 |title=Adaptive estimation of a quadratic functional by model selection |journal=The Annals of Statistics |volume=28 |issue=5 |doi=10.1214/aos/1015957395 |s2cid=116945590 |issn=0090-5364|doi-access=free }}</ref> bounds are:
: <math>\operatorname{P}(X - k \ge 2 \sqrt{k x} + 2x) \le \exp(-x)</math>
<math display="block">\Pr(X - k \ge 2 \sqrt{k x} + 2x) \le e^{-x}</math>
: <math>\operatorname{P}(k - X \ge 2 \sqrt{k x}) \le \exp(-x)</math>
<math display="block">\Pr(k - X \ge 2 \sqrt{k x}) \le e^{-x}</math>
One consequence is that, if <math>Z \sim N(0, 1)^k</math> is a gaussian random vector in <math>\R^k</math>, then as the dimension <math>k</math> grows, the squared length of the vector is concentrated tightly around <math>k</math> with a width <math>k^{1/2 + \alpha}</math>:<math display="block">Pr(\|Z\|^2 \in [k - 2k^{1/2+\alpha}, k + 2k^{1/2+\alpha} + 2k^{\alpha}]) \geq 1-e^{-k^\alpha}</math>where the exponent <math>\alpha</math> can be chosen as any value in <math>\R</math>.
One consequence is that, if <math>Z \sim N(0, 1)^k</math> is a Gaussian random vector in <math>\R^k</math>, then as the dimension <math>k</math> grows, the squared length of the vector is concentrated tightly around <math>k</math> with a width <math>k^{1/2 + \alpha}</math>:<math display="block">\Pr\left(\left\|Z\right\|^2 \in \left[k - 2k^{1/2+\alpha}, \; k + 2k^{1/2+\alpha} + 2k^{\alpha}\right]\right) \geq 1-e^{-k^\alpha}</math>where the exponent <math>\alpha</math> can be chosen as any value in <math>\R</math>.


Since the cumulant generating function for <math>\chi^2(k)</math> is <math>K(t) = -\frac k2 \ln(1-2t) </math>, and its [[Convex conjugate|convex dual]] is <math>K^*(q) = \frac 12 (q-k + k\ln\frac kq) </math>, the standard [[Chernoff bound]] yields<math display="block">\begin{aligned}
Since the cumulant generating function for <math>\chi^2(k)</math> is <math display="inline">K(t) = -\frac k2 \ln(1-2t) </math>, and its [[Convex conjugate|convex dual]] is <math display="inline">K^*(q) = \frac{1}{2} \left(q - k + k\ln\frac{k}{q}\right) </math>, the standard [[Chernoff bound]] yields<math display="block">\begin{aligned}
\ln Pr(X \geq (1 + \epsilon) k) &\leq -\frac k2 ( \epsilon - \ln(1+\epsilon)) \\
\ln \Pr(X \geq (1 + \varepsilon) k) &\leq -\frac{k}{2} \left( \varepsilon - \ln(1+\varepsilon)\right) \\
\ln Pr(X \leq (1 - \epsilon) k) &\leq -\frac k2 ( -\epsilon - \ln(1-\epsilon))
\ln \Pr(X \leq (1 - \varepsilon) k) &\leq -\frac{k}{2} \left(-\varepsilon - \ln(1-\varepsilon)\right)
\end{aligned}</math>where <math>0< \epsilon < 1</math>. By the union bound,<math display="block">Pr(X \in (1\pm \epsilon ) k ) \geq 1 - 2e^{-\frac k2 (\frac 12 \epsilon^2 - \frac 13 \epsilon^3)} </math>This result is used in proving the [[Johnson–Lindenstrauss lemma]].<ref>[https://ocw.mit.edu/courses/18-s096-topics-in-mathematics-of-data-science-fall-2015/f9261308512f6b90e284599f94055bb4_MIT18_S096F15_Ses15_16.pdf MIT 18.S096 (Fall 2015): Topics in Mathematics of Data Science, Lecture 5, Johnson-Lindenstrauss Lemma and Gordons Theorem]</ref>
\end{aligned}</math>where <math>0 < \varepsilon < 1</math>. By the union bound,<math display="block">Pr(X \in (1\pm \varepsilon ) k ) \geq 1 - 2 e^{-\frac k2 (\frac{1}{2} \varepsilon^2 - \frac 13 \varepsilon^3)} </math>This result is used in proving the [[Johnson–Lindenstrauss lemma]].<ref>[https://ocw.mit.edu/courses/18-s096-topics-in-mathematics-of-data-science-fall-2015/f9261308512f6b90e284599f94055bb4_MIT18_S096F15_Ses15_16.pdf MIT 18.S096 (Fall 2015): Topics in Mathematics of Data Science, Lecture 5, Johnson-Lindenstrauss Lemma and Gordons Theorem]</ref>


=== Asymptotic properties ===
=== Asymptotic properties ===
[[File:Chi-square median approx.png|thumb|upright=1.818|Approximate formula for median (from the Wilson–Hilferty transformation) compared with numerical quantile (top); and difference ({{font color|blue|blue}}) and relative difference ({{font color|red|red}}) between numerical quantile and approximate formula (bottom). For the chi-squared distribution, only the positive integer numbers of degrees of freedom (circles) are meaningful.]]
[[File:Chi-square median approx.png|thumb|upright=1.818|Approximate formula for median (from the Wilson–Hilferty transformation) compared with numerical quantile (top); and difference ({{font color|blue|blue}}) and relative difference ({{font color|red|red}}) between numerical quantile and approximate formula (bottom). For the chi-squared distribution, only the positive integer numbers of degrees of freedom (circles) are meaningful.]]


By the [[central limit theorem]], because the chi-squared distribution is the sum of <math>k</math> independent random variables with finite mean and variance, it converges to a normal distribution for large <math>k</math>. For many practical purposes, for <math>k>50</math> the distribution is sufficiently close to a [[normal distribution]], so the difference is ignorable.<ref>{{cite book|title=Statistics for experimenters|author=Box, Hunter and Hunter|publisher=Wiley|year=1978|isbn=978-0-471-09315-2|page=[https://archive.org/details/statisticsforexp00geor/page/118 118]|url-access=registration|url=https://archive.org/details/statisticsforexp00geor/page/118}}</ref> Specifically, if <math>X \sim \chi^2(k)</math>, then as <math>k</math> tends to infinity, the distribution of <math>(X-k)/\sqrt{2k}</math> [[convergence of random variables#Convergence in distribution|tends]] to a standard normal distribution. However, convergence is slow as the [[skewness]] is <math>\sqrt{8/k}</math> and the [[excess kurtosis]] is <math>12/k</math>.
By the [[central limit theorem]], because the chi-squared distribution is the sum of <math>k</math> independent random variables with finite mean and variance, it converges to a normal distribution for large <math>k</math>. For many practical purposes, for <math>k>50</math> the distribution is sufficiently close to a [[normal distribution]], so the difference is ignorable.<ref>{{cite book|title=Statistics for experimenters|author=Box, Hunter and Hunter|publisher=Wiley|year=1978|isbn=978-0-471-09315-2|page=[https://archive.org/details/statisticsforexp00geor/page/118 118]|url-access=registration|url=https://archive.org/details/statisticsforexp00geor/page/118}}</ref> Specifically, if <math>X \sim \chi^2(k)</math>, then as <math>k</math> tends to infinity, the distribution of <math>(X-k)/\sqrt{2k}</math> [[convergence of random variables#Convergence in distribution|tends]] to a standard normal distribution. However, convergence is slow as the [[skewness]] is <math display="inline">\sqrt{8/k}</math> and the [[excess kurtosis]] is <math>12/k</math>.


The sampling distribution of <math>\ln(\chi^2)</math> converges to normality much faster than the sampling distribution of <math>\chi^2</math>,<ref>{{cite journal |first1=M. S. |last1=Bartlett |first2=D. G. |last2=Kendall |title=The Statistical Analysis of Variance-Heterogeneity and the Logarithmic Transformation |journal=Supplement to the Journal of the Royal Statistical Society |volume=8 |issue=1 |year=1946 |pages=128–138 |jstor=2983618 |doi=10.2307/2983618 }}</ref> as the [[logarithmic transformation|logarithmic transform]] removes much of the asymmetry.<ref name="Pillai-2016">{{Cite journal|last=Pillai|first=Natesh S.|year=2016|title=An unexpected encounter with Cauchy and Lévy|journal=[[Annals of Statistics]]|volume=44|issue=5|pages=2089–2097|doi=10.1214/15-aos1407|arxiv=1505.01957|s2cid=31582370}}</ref>
The sampling distribution of <math>\ln(\chi^2)</math> converges to normality much faster than the sampling distribution of <math>\chi^2</math>,<ref>{{cite journal |first1=M. S. |last1=Bartlett |first2=D. G. |last2=Kendall |title=The Statistical Analysis of Variance-Heterogeneity and the Logarithmic Transformation |journal=Supplement to the Journal of the Royal Statistical Society |volume=8 |issue=1 |year=1946 |pages=128–138 |jstor=2983618 |doi=10.2307/2983618 }}</ref> as the [[logarithmic transformation|logarithmic transform]] removes much of the asymmetry.<ref name="Pillai-2016">{{Cite journal|last=Pillai|first=Natesh S.|year=2016|title=An unexpected encounter with Cauchy and Lévy|journal=[[Annals of Statistics]]|volume=44|issue=5|pages=2089–2097|doi=10.1214/15-aos1407|arxiv=1505.01957|s2cid=31582370}}</ref>


Other functions of the chi-squared distribution converge more rapidly to a normal distribution. Some examples are:
Other functions of the chi-squared distribution converge more rapidly to a normal distribution. Some examples are:
* If <math>X \sim \chi^2(k)</math> then <math>\sqrt{2X}</math> is approximately normally distributed with mean <math>\sqrt{2k-1}</math> and unit variance (1922, by [[R. A. Fisher]], see (18.23), p.&nbsp;426 of Johnson.<ref name="Johnson-1994" />
* If <math>X \sim \chi^2(k)</math> then <math>\sqrt{2X}</math> is approximately normally distributed with mean <math>\sqrt{2k-1}</math> and unit variance (1922, by [[R. A. Fisher]], see (18.23), p.&nbsp;426 of Johnson).<ref name="Johnson-1994" />
* If <math>X \sim \chi^2(k)</math> then <math>\sqrt[3]{X/k}</math> is approximately normally distributed with mean <math> 1-\frac{2}{9k}</math> and variance <math>\frac{2}{9k} .</math><ref>{{cite journal |last1=Wilson |first1=E. B. |last2=Hilferty |first2=M. M. |year=1931 |title=The distribution of chi-squared |journal=[[Proc. Natl. Acad. Sci. USA]] |volume=17 |issue=12 |pages=684–688 |bibcode=1931PNAS...17..684W |doi=10.1073/pnas.17.12.684 |pmid=16577411 |pmc=1076144 |doi-access=free }}</ref> This is known as the '''Wilson–Hilferty transformation''', see (18.24), p.&nbsp;426 of Johnson.<ref name="Johnson-1994" />
* If <math>X \sim \chi^2(k)</math> then <math display="inline">\sqrt[3]{X/k}</math> is approximately normally distributed with mean <math> 1-\frac{2}{9k}</math> and variance <math>\frac{2}{9k} .</math><ref>{{cite journal |last1=Wilson |first1=E. B. |last2=Hilferty |first2=M. M. |year=1931 |title=The distribution of chi-squared |journal=[[Proc. Natl. Acad. Sci. USA]] |volume=17 |issue=12 |pages=684–688 |bibcode=1931PNAS...17..684W |doi=10.1073/pnas.17.12.684 |pmid=16577411 |pmc=1076144 |doi-access=free }}</ref> This is known as the '''Wilson–Hilferty transformation''', see (18.24), p.&nbsp;426 of Johnson.<ref name="Johnson-1994" />
** This [[Data transformation (statistics)#Transforming to normality|normalizing transformation]] leads directly to the commonly used median approximation <math>k\bigg(1-\frac{2}{9k}\bigg)^3\;</math> by back-transforming from the mean, which is also the median, of the normal distribution.
** This [[Data transformation (statistics)#Transforming to normality|normalizing transformation]] leads directly to the commonly used median approximation <math>k\bigg(1-\frac{2}{9k}\bigg)^3\;</math> by back-transforming from the mean, which is also the median, of the normal distribution.


Line 243: Line 232:
* <math> \chi_k^2 \sim {\chi'}^2_k(0)</math> ([[noncentral chi-squared distribution]] with non-centrality parameter <math> \lambda = 0 </math>)
* <math> \chi_k^2 \sim {\chi'}^2_k(0)</math> ([[noncentral chi-squared distribution]] with non-centrality parameter <math> \lambda = 0 </math>)
* If <math>Y \sim \mathrm{F}(\nu_1, \nu_2)</math> then <math>X = \lim_{\nu_2 \to \infty} \nu_1 Y</math> has the chi-squared distribution <math>\chi^2_{\nu_{1}}</math>
* If <math>Y \sim \mathrm{F}(\nu_1, \nu_2)</math> then <math>X = \lim_{\nu_2 \to \infty} \nu_1 Y</math> has the chi-squared distribution <math>\chi^2_{\nu_{1}}</math>
:*As a special case, if <math>Y \sim \mathrm{F}(1, \nu_2)\,</math> then <math>X = \lim_{\nu_2 \to \infty} Y\,</math> has the chi-squared distribution <math>\chi^2_{1}</math>
**As a special case, if <math>Y \sim \mathrm{F}(1, \nu_2)\,</math> then <math>X = \lim_{\nu_2 \to \infty} Y\,</math> has the chi-squared distribution <math>\chi^2_{1}</math>
* <math> \|\boldsymbol{N}_{i=1,\ldots,k} (0,1) \|^2 \sim \chi^2_k </math> (The squared [[Norm (mathematics)|norm]] of ''k'' standard normally distributed variables is a chi-squared distribution with ''k'' [[degrees of freedom (statistics)|degrees of freedom]])
* <math> \left\|\boldsymbol{N}_{i=1,\ldots,k} (0,1) \right\|^2 \sim \chi^2_k </math> (The squared [[Norm (mathematics)|norm]] of ''k'' standard normally distributed variables is a chi-squared distribution with ''k'' [[degrees of freedom (statistics)|degrees of freedom]])
* If <math>X \sim \chi^2_\nu\,</math> and <math>c>0 \,</math>, then <math>cX \sim \Gamma(k=\nu/2, \theta=2c)\,</math>. ([[gamma distribution]])
* If <math>X \sim \chi^2_\nu\,</math> and <math>c>0 \,</math>, then <math>cX \sim \Gamma(k=\nu/2, \theta=2c)\,</math>. ([[gamma distribution]])
* If <math>X \sim \chi^2_k</math> then <math>\sqrt{X} \sim \chi_k</math> ([[chi distribution]])
* If <math>X \sim \chi^2_k</math> then <math>\sqrt{X} \sim \chi_k</math> ([[chi distribution]])
Line 257: Line 246:
* If <math> X \sim \operatorname{U}(0,1)\, </math> ([[Uniform distribution (continuous)|uniform distribution]]) then <math> -2\log(X) \sim \chi^2_2\,</math>
* If <math> X \sim \operatorname{U}(0,1)\, </math> ([[Uniform distribution (continuous)|uniform distribution]]) then <math> -2\log(X) \sim \chi^2_2\,</math>
* If <math>X_i \sim \operatorname{Laplace}(\mu,\beta)\,</math> then <math>\sum_{i=1}^n \frac{2 |X_i-\mu|}{\beta} \sim \chi^2_{2n}\,</math>
* If <math>X_i \sim \operatorname{Laplace}(\mu,\beta)\,</math> then <math>\sum_{i=1}^n \frac{2 |X_i-\mu|}{\beta} \sim \chi^2_{2n}\,</math>
* If <math>X_i</math> follows the [[generalized normal distribution]] (version 1) with parameters <math>\mu,\alpha,\beta</math> then <math>\sum_{i=1}^n \frac{2 |X_i-\mu|^\beta}{\alpha} \sim \chi^2_{2n/\beta}\,</math> <ref>{{cite journal |last= Bäckström |first= T. |author2=Fischer, J. |date=January 2018|title= Fast Randomization for Distributed Low-Bitrate Coding of Speech and Audio|journal= IEEE/ACM Transactions on Audio, Speech, and Language Processing |volume= 26|issue= 1|pages= 19–30|doi= 10.1109/TASLP.2017.2757601|s2cid= 19777585 |url= https://research.aalto.fi/files/27158975/ELEC_backstrom_et_al_Fast_randomization.pdf }}</ref>
* If <math>X_i</math> follows the [[generalized normal distribution]] (version 1) with parameters <math>\mu,\alpha,\beta</math> then <math>\sum_{i=1}^n \frac{2 |X_i-\mu|^\beta}{\alpha} \sim \chi^2_{2n/\beta}\,</math><ref>{{cite journal |last= Bäckström |first= T. |author2=Fischer, J. |date=January 2018|title= Fast Randomization for Distributed Low-Bitrate Coding of Speech and Audio|journal= IEEE/ACM Transactions on Audio, Speech, and Language Processing |volume= 26|issue= 1|pages= 19–30|doi= 10.1109/TASLP.2017.2757601|s2cid= 19777585 |url= https://research.aalto.fi/files/27158975/ELEC_backstrom_et_al_Fast_randomization.pdf }}</ref>
* The chi-squared distribution is a transformation of [[Pareto distribution]]
* The chi-squared distribution is a transformation of [[Pareto distribution]]
* [[Student's t-distribution]] is a transformation of chi-squared distribution
* [[Student's t-distribution]] is a transformation of chi-squared distribution
Line 266: Line 255:
A chi-squared variable with <math>k</math> degrees of freedom is defined as the sum of the squares of <math>k</math> independent [[standard normal]] random variables.
A chi-squared variable with <math>k</math> degrees of freedom is defined as the sum of the squares of <math>k</math> independent [[standard normal]] random variables.


If <math>Y</math> is a <math>k</math>-dimensional Gaussian random vector with mean vector <math>\mu</math> and rank <math>k</math> covariance matrix <math>C</math>, then <math>X = (Y-\mu )^{T}C^{-1}(Y-\mu)</math> is chi-squared distributed with <math>k</math> degrees of freedom.
If <math>Y</math> is a <math>k</math>-dimensional Gaussian random vector with mean vector <math>\mu</math> and rank <math>k</math> covariance matrix <math>C</math>, then <math>X = (Y-\mu )^\mathsf{T}C^{-1}(Y-\mu)</math> is chi-squared distributed with <math>k</math> degrees of freedom.


The sum of squares of [[statistically independent]] unit-variance Gaussian variables which do ''not'' have mean zero yields a generalization of the chi-squared distribution called the [[noncentral chi-squared distribution]].
The sum of squares of [[statistically independent]] unit-variance Gaussian variables which do ''not'' have mean zero yields a generalization of the chi-squared distribution called the [[noncentral chi-squared distribution]].


If <math>Y</math> is a vector of <math>k</math> [[i.i.d.]] standard normal random variables and <math>A</math> is a <math>k\times k</math> [[symmetric matrix|symmetric]], [[idempotent matrix]] with [[rank (linear algebra)|rank]] <math>k-n</math>, then the [[quadratic form]] <math>Y^TAY</math> is chi-square distributed with <math>k-n</math> degrees of freedom.
If <math>Y</math> is a vector of <math>k</math> [[i.i.d.]] standard normal random variables and <math>A</math> is a <math>k\times k</math> [[symmetric matrix|symmetric]], [[idempotent matrix]] with [[rank (linear algebra)|rank]] <math>k-n</math>, then the [[quadratic form]] <math>Y^\mathsf{T}\!AY</math> is chi-square distributed with <math>k-n</math> degrees of freedom.


If <math>\Sigma</math> is a <math>p\times p</math> positive-semidefinite covariance matrix with strictly positive diagonal entries, then for <math>X\sim N(0,\Sigma)</math> and <math>w</math> a random <math>p</math>-vector independent of <math>X</math> such that <math>w_1+\cdots+w_p=1</math> and <math>w_i\geq 0, i=1,\ldots,p,</math> then
If <math>\Sigma</math> is a <math>p\times p</math> positive-semidefinite covariance matrix with strictly positive diagonal entries, then for <math>X\sim N(0,\Sigma)</math> and <math>w</math> a random <math>p</math>-vector independent of <math>X</math> such that <math>w_1+\cdots+w_p=1</math> and <math>w_i\geq 0, i=1,\ldots,p,</math> then<ref name="Pillai-2016" />


: <math>\frac{1}{\left(\frac{w_1}{X_1},\ldots,\frac{w_p}{X_p}\right)\Sigma\left(\frac{w_1}{X_1},\ldots,\frac{w_p}{X_p}\right)^\top} \sim \chi_1^2.</math><ref name="Pillai-2016" />
<math display="block">\frac{1}{\tilde{w}^\mathsf{T} \Sigma \tilde{w}} \sim \chi_1^2.</math>where <math>\tilde{w} = (w_1/X_1, \dots, w_p/X_p)</math>.


The chi-squared distribution is also naturally related to other distributions arising from the Gaussian. In particular,
The chi-squared distribution is also naturally related to other distributions arising from the Gaussian. In particular,
Line 284: Line 273:


=== Linear combination ===
=== Linear combination ===
If <math>X_1,\ldots,X_n</math> are chi square random variables and <math>a_1,\ldots,a_n\in\mathbb{R}_{>0}</math>, then the distribution of <math>X=\sum_{i=1}^n a_i X_i</math> is a special case of a [[Generalized chi-squared distribution|Generalized Chi-squared Distribution]].
If <math>X_1,\ldots,X_n</math> are chi square random variables and <math>a_1,\ldots,a_n\in\mathbb{R}_{>0}</math>, then the distribution of <math display="inline">X = \sum_{i=1}^n a_i X_i</math> is a special case of the [[generalized chi-squared distribution]].
A closed expression for this distribution is not known. It may be, however, approximated efficiently using the [[Characteristic function (probability theory)#Properties|property of characteristic functions]] of chi-square random variables.<ref>{{cite journal
A closed expression for this distribution is not known. It may be, however, approximated efficiently using the [[Characteristic function (probability theory)#Properties|property of characteristic functions]] of chi-square random variables.<ref>{{cite journal
|first=J.
|first=J.
Line 293: Line 282:
|issue=50
|issue=50
|year=2013
|year=2013
|pages=505202
|article-number=505202
|doi=10.1088/1751-8113/46/50/505202 |bibcode=2013JPhA...46X5202B
|doi=10.1088/1751-8113/46/50/505202 |bibcode=2013JPhA...46X5202B
|arxiv=1208.2691
|arxiv=1208.2691
Line 310: Line 299:


=== Gamma, exponential, and related distributions ===
=== Gamma, exponential, and related distributions ===
The chi-squared distribution <math>X \sim \chi_k^2</math> is a special case of the [[gamma distribution]], in that <math>X \sim \Gamma \left(\frac{k}2,\frac{1}2\right)</math> using the rate parameterization of the gamma distribution (or
The chi-squared distribution <math>X \sim \chi_k^2</math> is a special case of the [[gamma distribution]], in that <math display="inline">X \sim \Gamma{\left(\tfrac{k}{2},\tfrac{1}{2}\right)}</math> using the rate parameterization of the gamma distribution (or <math display="inline">X \sim \Gamma {\left(\tfrac{k}{2},2 \right)}</math> using the scale parameterization of the gamma distribution)
<math>X \sim \Gamma \left(\frac{k}2,2 \right)</math> using the scale parameterization of the gamma distribution)
where {{mvar|k}} is an integer.
where {{mvar|k}} is an integer.


Because the [[exponential distribution]] is also a special case of the gamma distribution, we also have that if <math>X \sim \chi_2^2</math>, then <math>X\sim \operatorname{exp}\left(\frac 1 2\right)</math> is an [[exponential distribution]].
Because the [[exponential distribution]] is also a special case of the gamma distribution, we also have that if <math>X \sim \chi_2^2</math>, then <math display="inline">X\sim \operatorname{exp}\left(\tfrac 1 2\right)</math> is an [[exponential distribution]].


The [[Erlang distribution]] is also a special case of the gamma distribution and thus we also have that if <math>X \sim\chi_k^2</math> with even <math>k</math>, then <math>X</math> is Erlang distributed with shape parameter <math>k/2</math> and scale parameter <math>1/2</math>.
The [[Erlang distribution]] is also a special case of the gamma distribution and thus we also have that if <math>X \sim\chi_k^2</math> with even <math>k</math>, then <math>X</math> is Erlang distributed with shape parameter <math>k/2</math> and scale parameter <math>1/2</math>.
Line 322: Line 310:


Following are some of the most common situations in which the chi-squared distribution arises from a Gaussian-distributed sample.
Following are some of the most common situations in which the chi-squared distribution arises from a Gaussian-distributed sample.
* if <math>X_1, ..., X_n</math> are [[i.i.d.]] <math>N(\mu, \sigma^2)</math> [[random variable]]s, then <math>\sum_{i=1}^n(X_i - \overline{X})^2 \sim \sigma^2 \chi^2_{n-1}</math> where <math>\overline{X} = \frac{1}{n} \sum_{i=1}^n X_i</math>.
* if <math>X_1, ..., X_n</math> are [[i.i.d.]] <math>N(\mu, \sigma^2)</math> [[random variable]]s, then <math display="inline">\sum_{i=1}^n \left(X_i - \bar{X}\right)^2 \sim \sigma^2 \chi^2_{n-1}</math> where <math display="inline">\bar{X} = \frac{1}{n} \sum_{i=1}^n X_i</math>.
* The box below shows some [[statistics]] based on <math>X_i \sim N(\mu_i, \sigma^2_i), i= 1, \ldots, k</math> independent random variables that have probability distributions related to the chi-squared distribution:
* The box below shows some [[statistics]] based on <math>X_i \sim N(\mu_i, \sigma^2_i), i= 1, \ldots, k</math> independent random variables that have probability distributions related to the chi-squared distribution:
{| class="wikitable"  style="margin:1em auto;" align="center"
{| class="wikitable"  style="margin:1em auto;" align="center"
Line 336: Line 324:
| [[noncentral chi distribution]] || <math>\sqrt{\sum_{i=1}^k \left(\frac{X_i}{\sigma_i}\right)^2}</math>
| [[noncentral chi distribution]] || <math>\sqrt{\sum_{i=1}^k \left(\frac{X_i}{\sigma_i}\right)^2}</math>
|}
|}
The chi-squared distribution is also often encountered in [[magnetic resonance imaging]].<ref>den Dekker A. J., Sijbers J., (2014) "Data distributions in magnetic resonance images: a review", ''Physica Medica'', [https://dx.doi.org/10.1016/j.ejmp.2014.05.002]</ref>
The chi-squared distribution is also often encountered in [[magnetic resonance imaging]].<ref>{{cite journal | last = den Dekker | first = A. J. | last2 = Sijbers | first2 = J. | date = 2014 | title = Data distributions in magnetic resonance images: A review |journal = Physica Medica | doi = 10.1016/j.ejmp.2014.05.002 | url = https://linkinghub.elsevier.com/retrieve/pii/S1120179714000829 | language = en |volume = 30 | issue = 7 | pages = 725–741| url-access = subscription }}</ref>


== Computational methods ==
== Computational methods ==
Line 344: Line 332:
The table below gives a number of ''p''-values matching to <math> \chi^2 </math> for the first 10 degrees of freedom.
The table below gives a number of ''p''-values matching to <math> \chi^2 </math> for the first 10 degrees of freedom.
{| class="wikitable"
{| class="wikitable"
! Degrees of freedom (df)
! Degrees of{{pb}}freedom (df)
!colspan=11| <math> \chi^2 </math> value<ref>[http://www2.lv.psu.edu/jxm57/irp/chisquar.html Chi-Squared Test] {{Webarchive|url=https://web.archive.org/web/20131118011437/http://www2.lv.psu.edu/jxm57/irp/chisquar.html |date=2013-11-18 }} Table B.2. Dr. Jacqueline S. McLaughlin at The Pennsylvania State University. In turn citing: R. A. Fisher and F. Yates, Statistical Tables for Biological Agricultural and Medical Research, 6th ed., Table IV. Two values have been corrected, 7.82 with 7.81 and 4.60 with 4.61</ref>
!colspan=11| <math> \chi^2 </math> value<ref>[http://www2.lv.psu.edu/jxm57/irp/chisquar.html Chi-Squared Test] {{Webarchive|url=https://web.archive.org/web/20131118011437/http://www2.lv.psu.edu/jxm57/irp/chisquar.html |date=2013-11-18 }} Table B.2. Dr. Jacqueline S. McLaughlin at The Pennsylvania State University. In turn citing: R. A. Fisher and F. Yates, Statistical Tables for Biological Agricultural and Medical Research, 6th ed., Table IV. Two values have been corrected, 7.82 with 7.81 and 4.60 with 4.61</ref>
|-
|-
| style="text-align:center;" | 1
| style="text-align:center;" | 1
| 0.004
| 0.004 || 0.02 || 0.06 || 0.15 || 0.46 || 1.07 || 1.64 |2.71 |3.84 |6.63 || 10.83
| 0.02
| 0.06
| 0.15
| 0.46
| 1.07
| 1.64
| 2.71
| 3.84
| 6.63
| 10.83
|-
|-
| style="text-align:center;" | 2
| style="text-align:center;" | 2
| 0.10
| 0.10 || 0.21 || 0.45 || 0.71 || 1.39 |2.41 || 3.22 |4.61 |5.99 |9.21 || 13.82
| 0.21
| 0.45
| 0.71
| 1.39
| 2.41
| 3.22
| 4.61
| 5.99
| 9.21
| 13.82
|-
|-
| style="text-align:center;" | 3
| style="text-align:center;" | 3
| 0.35
| 0.35 || 0.58 || 1.01 || 1.42 || 2.37 |3.66 || 4.64 |6.25 |7.81 || 11.34 || 16.27
| 0.58
| 1.01
| 1.42
| 2.37
| 3.66
| 4.64
| 6.25
| 7.81
| 11.34
| 16.27
|-
|-
| style="text-align:center;" | 4
| style="text-align:center;" | 4
| 0.71
| 0.71 || 1.06 || 1.65 || 2.20 || 3.36 |4.88 || 5.99 |7.78 |9.49 || 13.28 || 18.47
| 1.06
| 1.65
| 2.20
| 3.36
| 4.88
| 5.99
| 7.78
| 9.49
| 13.28
| 18.47
|-
|-
| style="text-align:center;" | 5
| style="text-align:center;" | 5
| 1.14
| 1.14 || 1.61 || 2.34 || 3.00 || 4.35 |6.06 || 7.29 |9.24 || 11.07 || 15.09 || 20.52
| 1.61
| 2.34
| 3.00
| 4.35
| 6.06
| 7.29
| 9.24
| 11.07
| 15.09
| 20.52
|-
|-
| style="text-align:center;" | 6
| style="text-align:center;" | 6
| 1.63
| 1.63 || 2.20 || 3.07 || 3.83 || 5.35 |7.23 || 8.56 || 10.64 || 12.59 || 16.81 || 22.46
| 2.20
| 3.07
| 3.83
| 5.35
| 7.23
| 8.56
| 10.64
| 12.59
| 16.81
| 22.46
|-
|-
| style="text-align:center;" | 7
| style="text-align:center;" | 7
| 2.17
| 2.17 || 2.83 || 3.82 || 4.67 || 6.35 |8.38 || 9.80 || 12.02 || 14.07 || 18.48 || 24.32
| 2.83
| 3.82
| 4.67
| 6.35
| 8.38
| 9.80
| 12.02
| 14.07
| 18.48
| 24.32
|-
|-
| style="text-align:center;" | 8
| style="text-align:center;" | 8
| 2.73
| 2.73 || 3.49 || 4.59 || 5.53 || 7.34 |9.52 || 11.03 || 13.36 || 15.51 || 20.09 || 26.12
| 3.49
| 4.59
| 5.53
| 7.34
| 9.52
| 11.03
| 13.36
| 15.51
| 20.09
| 26.12
|-
|-
| style="text-align:center;" | 9
| style="text-align:center;" | 9
| 3.32
| 3.32 || 4.17 || 5.38 || 6.39 || 8.34 || 10.66 || 12.24 || 14.68 || 16.92 || 21.67 || 27.88
| 4.17
| 5.38
| 6.39
| 8.34
| 10.66
| 12.24
| 14.68
| 16.92
| 21.67
| 27.88
|-
|-
| style="text-align:center;" | 10
| style="text-align:center;" | 10
| 3.94
| 3.94 || 4.87 || 6.18 || 7.27 || 9.34 || 11.78 || 13.44 || 15.99 || 18.31 || 23.21 || 29.59
| 4.87
| 6.18
| 7.27
| 9.34
| 11.78
| 13.44
| 15.99
| 18.31
| 23.21
| 29.59
13.25
|-
|-
! scope="row" style="text-align:right;" | ''p''-value (probability)
! scope="row" style="text-align: center;" | ''p''-value{{pb}}(probability)
| style="background: #ffa2aa" | 0.95
| style="background: #ffa2aa" | 0.95
| style="background: #efaaaa" | 0.90
| style="background: #efaaaa" | 0.90
Line 498: Line 385:


The distribution was independently rediscovered by the English mathematician [[Karl Pearson]] in the context of [[goodness of fit]], for which he developed his [[Pearson's chi-squared test]], published in 1900, with computed table of values published in {{Harv|Elderton|1902}}, collected in {{Harv|Pearson|1914|pp=xxxi–xxxiii, 26–28|loc=Table XII}}.
The distribution was independently rediscovered by the English mathematician [[Karl Pearson]] in the context of [[goodness of fit]], for which he developed his [[Pearson's chi-squared test]], published in 1900, with computed table of values published in {{Harv|Elderton|1902}}, collected in {{Harv|Pearson|1914|pp=xxxi–xxxiii, 26–28|loc=Table XII}}.
The name "chi-square" ultimately derives from Pearson's shorthand for the exponent in a [[multivariate normal distribution]] with the Greek letter [[Chi (letter)|Chi]], writing {{mvar|−½χ<sup>2</sup>}} for what would appear in modern notation as {{math|−½'''x'''<sup>T</sup>Σ<sup>−1</sup>'''x'''}} (Σ being the [[covariance matrix]]).<ref>R. L. Plackett, ''Karl Pearson and the Chi-Squared Test'', International Statistical Review, 1983, [https://www.jstor.org/stable/1402731?seq=3 61f.]
The name "chi-square" ultimately derives from Pearson's shorthand for the exponent in a [[multivariate normal distribution]] with the Greek letter [[Chi (letter)|Chi]], writing {{math|−{{1/2}}''χ''<sup>2</sup>}} for what would appear in modern notation as {{math|−{{1/2}}'''x'''<sup>T</sup>Σ<sup>−1</sup>'''x'''}} ({{math|Σ}} being the [[covariance matrix]]).<ref>R. L. Plackett, ''Karl Pearson and the Chi-Squared Test'', International Statistical Review, 1983, [https://www.jstor.org/stable/1402731?seq=3 61f.]
See also Jeff Miller, [http://jeff560.tripod.com/c.html Earliest Known Uses of Some of the Words of Mathematics].</ref> The idea of a family of "chi-squared distributions", however, is not due to Pearson but arose as a further development due to Fisher in the 1920s.{{sfn|Hald|1998|pp=633–692|loc=27. Sampling Distributions under Normality}}
See also Jeff Miller, [http://jeff560.tripod.com/c.html Earliest Known Uses of Some of the Words of Mathematics].</ref> The idea of a family of "chi-squared distributions", however, is not due to Pearson but arose as a further development due to Fisher in the 1920s.{{sfn|Hald|1998|pp=633–692|loc=27. Sampling Distributions under Normality}}


Line 512: Line 399:
* [[Reduced chi-squared statistic]]
* [[Reduced chi-squared statistic]]
* [[Wilks's lambda distribution]]
* [[Wilks's lambda distribution]]
* [[Modified half-normal distribution]]<ref name="Sun-2021">{{cite journal |last1=Sun |first1=Jingchao |last2=Kong |first2=Maiying |last3=Pal |first3=Subhadip |title=The Modified-Half-Normal distribution: Properties and an efficient sampling scheme |journal=Communications in Statistics - Theory and Methods |date=22 June 2021 |volume=52 |issue=5 |pages=1591–1613 |doi=10.1080/03610926.2021.1934700 |s2cid=237919587 |url=https://figshare.com/articles/journal_contribution/The_Modified-Half-Normal_distribution_Properties_and_an_efficient_sampling_scheme/14825266/1/files/28535884.pdf |issn=0361-0926}}</ref> with the pdf on <math>(0, \infty)</math> is given as <math> f(x)= \frac{2\beta^{\alpha/2} x^{\alpha-1} \exp(-\beta x^2+ \gamma x )}{\Psi{\left(\frac{\alpha}{2}, \frac{ \gamma}{\sqrt{\beta}}\right)}}</math>, where <math>\Psi(\alpha,z)={}_1\Psi_1\left(\begin{matrix}\left(\alpha,\frac{1}{2}\right)\\(1,0)\end{matrix};z \right)</math> denotes the [[Fox–Wright Psi function]].
* [[Modified half-normal distribution]]<ref name="Sun-2021">{{cite journal |last1=Sun |first1=Jingchao |last2=Kong |first2=Maiying |last3=Pal |first3=Subhadip |title=The Modified-Half-Normal distribution: Properties and an efficient sampling scheme | journal = Communications in Statistics - Theory and Methods |date=22 June 2021 |volume=52 |issue=5 |pages=1591–1613 | doi = 10.1080/03610926.2021.1934700 |s2cid=237919587 |url=https://figshare.com/articles/journal_contribution/The_Modified-Half-Normal_distribution_Properties_and_an_efficient_sampling_scheme/14825266/1/files/28535884.pdf |issn=0361-0926}}</ref> with the pdf on <math>(0, \infty)</math> is given as {{nowrap|<math> f(x) = \frac{2\beta^{\alpha/2} x^{\alpha-1}}{\Psi{\left(\frac{\alpha}{2}, \frac{ \gamma}{\sqrt{\beta}}\right)}} e^{-\beta x^2+ \gamma x}</math>,}} where <math>\Psi(\alpha,z) = {}_1\Psi_1{\left(\begin{matrix}\left(\alpha,\frac{1}{2}\right)\\(1,0)\end{matrix};z \right)}</math> denotes the [[Fox–Wright Psi function]].
{{Colend}}
{{Colend}}


Line 519: Line 406:


==Sources==
==Sources==
{{refbegin}}
* {{cite book |title=A history of mathematical statistics from 1750 to 1930 |last=Hald |first=Anders |author-link=Anders Hald |year=1998 |publisher=Wiley |location=New York |isbn=978-0-471-17912-2 }}
* {{cite book |title=A history of mathematical statistics from 1750 to 1930 |last=Hald |first=Anders |author-link=Anders Hald |year=1998 |publisher=Wiley |location=New York |isbn=978-0-471-17912-2 }}
* {{Cite journal |last=Elderton |first=William Palin |author-link=William Palin Elderton |title=Tables for Testing the Goodness of Fit of Theory to Observation |doi=10.1093/biomet/1.2.155 |journal=Biometrika |volume=1 |issue=2 |pages=155–163 |year=1902 |url=https://zenodo.org/record/1431595}}
* {{Cite journal |last=Elderton |first=William Palin |author-link=William Palin Elderton |title=Tables for Testing the Goodness of Fit of Theory to Observation |doi=10.1093/biomet/1.2.155 |journal=Biometrika |volume=1 |issue=2 |pages=155–163 |year=1902 |url=https://zenodo.org/record/1431595}}
* {{cite journal |last=Pearson |first=Karl |title=On the probability that two independent distributions of frequency are really samples of the same population, with special reference to recent work on the identity of Trypanosome strains |date=1914 |journal=Biometrika |volume=10 |pages=85–154 |doi=10.1093/biomet/10.1.85}}
* {{cite journal |last=Pearson |first=Karl |title=On the probability that two independent distributions of frequency are really samples of the same population, with special reference to recent work on the identity of Trypanosome strains |date=1914 |journal=Biometrika |volume=10 |pages=85–154 |doi=10.1093/biomet/10.1.85}}
{{refend}}


== Further reading ==
== Further reading ==
Line 541: Line 430:
[[Category:Normal distribution]]
[[Category:Normal distribution]]
[[Category:Infinitely divisible probability distributions]]
[[Category:Infinitely divisible probability distributions]]
[[Category:Exponential family distributions]]

Latest revision as of 19:22, 18 November 2025

Template:Short description Script error: No such module "about".

Template:Probability distribution

In probability theory and statistics, the χ2-distribution with k degrees of freedom is the distribution of a sum of the squares of k independent standard normal random variables.[1]

The chi-squared distribution χk2 is a special case of the gamma distribution and the univariate Wishart distribution. Specifically if Xχk2 then XGamma(α=k2,θ=2) (where α is the shape parameter and θ the scale parameter of the gamma distribution) and XW1(1,k).

The scaled chi-squared distribution s2χk2 is a reparametrization of the gamma distribution and the univariate Wishart distribution. Specifically if Xs2χk2 then XGamma(α=k2,θ=2s2) and XW1(s2,k).

The chi-squared distribution is one of the most widely used probability distributions in inferential statistics, notably in hypothesis testing and in construction of confidence intervals.[2][3][4][5] This distribution is sometimes called the central chi-squared distribution, a special case of the more general noncentral chi-squared distribution.[6]

The chi-squared distribution is used in the common chi-squared tests for goodness of fit of an observed distribution to a theoretical one, the independence of two criteria of classification of qualitative data, and in finding the confidence interval for estimating the population standard deviation of a normal distribution from a sample standard deviation. Many other statistical tests also use this distribution, such as Friedman's analysis of variance by ranks.

Definitions

If Template:Math are independent, standard normal random variables, then the sum of their squares, X =i=1kZi2, is distributed according to the chi-squared distribution with Template:Mvar degrees of freedom. This is usually denoted as X  χ2(k)  or  X  χk2.

The chi-squared distribution has one parameter: a positive integer Template:Mvar that specifies the number of degrees of freedom (the number of random variables Zi being summed).

Introduction

The chi-squared distribution is used primarily in hypothesis testing, and to a lesser extent for confidence intervals for population variance when the underlying distribution is normal. Unlike more widely known distributions such as the normal distribution and the exponential distribution, the chi-squared distribution is not as often applied in the direct modeling of natural phenomena. It arises in the following hypothesis tests, among others:

It is also a component of the definition of the t-distribution and the F-distribution used in t-tests, analysis of variance, and regression analysis.

The primary reason for which the chi-squared distribution is extensively used in hypothesis testing is its relationship to the normal distribution. Many hypothesis tests use a test statistic, such as the t-statistic in a t-test. For these hypothesis tests, as the sample size, Template:Mvar, increases, the sampling distribution of the test statistic approaches the normal distribution (central limit theorem). Because the test statistic (such as Template:Mvar) is asymptotically normally distributed, provided the sample size is sufficiently large, the distribution used for hypothesis testing may be approximated by a normal distribution. Testing hypotheses using a normal distribution is well understood and relatively easy. The simplest chi-squared distribution is the square of a standard normal distribution. So wherever a normal distribution could be used for a hypothesis test, a chi-squared distribution could be used.

Suppose that Z is a random variable sampled from the standard normal distribution, where the mean is 0 and the variance is 1: ZN(0,1). Now, consider the random variable X=Z2. The distribution of the random variable X is an example of a chi-squared distribution:  X  χ12. The subscript 1 indicates that this particular chi-squared distribution is constructed from only 1 standard normal distribution. A chi-squared distribution constructed by squaring a single standard normal distribution is said to have 1 degree of freedom. Thus, as the sample size for a hypothesis test increases, the distribution of the test statistic approaches a normal distribution. Just as extreme values of the normal distribution have low probability (and give small p-values), extreme values of the chi-squared distribution have low probability.

An additional reason that the chi-squared distribution is widely used is that it turns up as the large sample distribution of generalized likelihood ratio tests (LRT).[7] LRTs have several desirable properties; in particular, simple LRTs commonly provide the highest power to reject the null hypothesis (Neyman–Pearson lemma) and this leads also to optimality properties of generalised LRTs. However, the normal and chi-squared approximations are only valid asymptotically. For this reason, it is preferable to use the t distribution rather than the normal approximation or the chi-squared approximation for a small sample size. Similarly, in analyses of contingency tables, the chi-squared approximation will be poor for a small sample size, and it is preferable to use Fisher's exact test. Ramsey shows that the exact binomial test is always more powerful than the normal approximation.[8]

Lancaster shows the connections among the binomial, normal, and chi-squared distributions, as follows.[9] De Moivre and Laplace established that a binomial distribution could be approximated by a normal distribution. Specifically they showed the asymptotic normality of the random variable

χ=mNpNpq

where m is the observed number of successes in N trials, where the probability of success is p, and q=1p.

Squaring both sides of the equation gives

χ2=(mNp)2Npq

Using N=Np+N(1p), N=m+(Nm), and q=1p, this equation can be rewritten as

χ2=(mNp)2Np+(NmNq)2Nq

The expression on the right is of the form that Karl Pearson would generalize to the form

χ2=i=1n(OiEi)2Ei

where

  • χ2 = Pearson's cumulative test statistic, which asymptotically approaches a χ2 distribution;
  • Oi = the number of observations of type i;
  • Ei=Npi = the expected (theoretical) frequency of type i, asserted by the null hypothesis that the fraction of type i in the population is pi; and
  • n = the number of cells in the table.Script error: No such module "Unsubst".

In the case of a binomial outcome (flipping a coin), the binomial distribution may be approximated by a normal distribution (for sufficiently large n). Because the square of a standard normal distribution is the chi-squared distribution with one degree of freedom, the probability of a result such as 1 heads in 10 trials can be approximated either by using the normal distribution directly, or the chi-squared distribution for the normalised, squared difference between observed and expected value. However, many problems involve more than the two possible outcomes of a binomial, and instead require 3 or more categories, which leads to the multinomial distribution. Just as de Moivre and Laplace sought for and found the normal approximation to the binomial, Pearson sought for and found a degenerate multivariate normal approximation to the multinomial distribution (the numbers in each category add up to the total sample size, which is considered fixed). Pearson showed that the chi-squared distribution arose from such a multivariate normal approximation to the multinomial distribution, taking careful account of the statistical dependence (negative correlations) between numbers of observations in different categories.[9]

Probability density function

The probability density function (pdf) of the chi-squared distribution is f(x;k)={xk/21ex/22k/2Γ(k2),x>0;0,otherwise. where Γ(k/2) denotes the gamma function, which has closed-form values for integer k.

For derivations of the pdf in the cases of one, two and k degrees of freedom, see Proofs related to chi-squared distribution.

Cumulative distribution function

File:Chernoff-bound.svg
Chernoff bound for the CDF and tail (1-CDF) of a chi-squared random variable with ten degrees of freedom (k=10)

Its cumulative distribution function is: F(x;k)=γ(k2,x2)Γ(k2)=P(k2,x2), where γ(s,t) is the lower incomplete gamma function and P(s,t) is the regularized gamma function.

In a special case of k=2 this function has the simple form: F(x;2)=1ex/2 which can be easily derived by integrating f(x;2)=12ex/2 directly. The integer recurrence of the gamma function makes it easy to compute F(x;k) for other small, even k.

Tables of the chi-squared cumulative distribution function are widely available and the function is included in many spreadsheets and all statistical packages.

Letting zx/k, Chernoff bounds on the lower and upper tails of the CDF may be obtained.[10] For the cases when 0<z<1 (which include all of the cases when this CDF is less than half): F(zk;k)(ze1z)k/2.

The tail bound for the cases when z>1, similarly, is 1F(zk;k)(ze1z)k/2.

For another approximation for the CDF modeled after the cube of a Gaussian, see under Noncentral chi-squared distribution.

Properties

Cochran's theorem

Script error: No such module "Labelled list hatnote". The following is a special case of Cochran's theorem.

Theorem. If Z1,...,Zn are independent identically distributed (i.i.d.), standard normal random variables, then t=1n(ZtZ¯)2χn12 where Z¯=1nt=1nZt.

<templatestyles src="Template:Hidden begin/styles.css"/>

[Proof]

Proof. Let Z𝒩(0¯,11) be a vector of n independent normally distributed random variables, and Z¯ their average. Then t=1n(ZtZ¯)2=t=1nZt2nZ¯2=Z[111n1¯1¯]Z=:ZMZ where 11 is the identity matrix and 1¯ the all ones vector. M has one eigenvector b1:=1n1¯ with eigenvalue 0, and n1 eigenvectors b2,...,bn (all orthogonal to b1) with eigenvalue 1, which can be chosen so that Q:=(b1,...,bn) is an orthogonal matrix. Since also X:=QZ𝒩(0¯,Q11Q)=𝒩(0¯,11), we have t=1n(ZtZ¯)2=ZMZ=XQMQX=X22+...+Xn2χn12, which proves the claim.

Additivity

It follows from the definition of the chi-squared distribution that the sum of independent chi-squared variables is also chi-squared distributed. Specifically, if Xi,i=1,n are independent chi-squared variables with ki, i=1,n degrees of freedom, respectively, then Y=X1++Xn is chi-squared distributed with k1++kn degrees of freedom.

Sample mean

The sample mean of n i.i.d. chi-squared variables of degree k is distributed according to a gamma distribution with shape α and scale θ parameters: X=1ni=1nXiGamma(α=nk2,θ=2n)where Xiχ2(k)

Asymptotically, given that for a shape parameter α going to infinity, a Gamma distribution converges towards a normal distribution with expectation μ=αθ and variance σ2=αθ2, the sample mean converges towards:

XnN(μ=k,σ2=2kn)

Note that we would have obtained the same result invoking instead the central limit theorem, noting that for each chi-squared variable of degree k the expectation is k, and its variance 2k (and hence the variance of the sample mean X being σ2=2kn).

Entropy

The differential entropy is given by h=0f(x;k)lnf(x;k)dx=k2+ln[2Γ(k2)]+(1k2)ψ(k2), where ψ(x) is the Digamma function.

The chi-squared distribution is the maximum entropy probability distribution for a random variate X for which E(X)=k and E(ln(X))=ψ(k/2)+ln(2) are fixed. Since the chi-squared is in the family of gamma distributions, this can be derived by substituting appropriate values in the Expectation of the log moment of gamma. For derivation from more basic principles, see the derivation in moment-generating function of the sufficient statistic.

Noncentral moments

The noncentral moments (raw moments) of a chi-squared distribution with k degrees of freedom are given by[11][12] E(Xm)=k(k+2)(k+4)(k+2m2)=2mΓ(m+k2)Γ(k2).

Cumulants

The cumulants are readily obtained by a power series expansion of the logarithm of the characteristic function: κn=2n1(n1)!k with cumulant generating function lnE[etX]=k2ln(12t).

Concentration

The chi-squared distribution exhibits strong concentration around its mean. The standard Laurent-Massart[13] bounds are: Pr(Xk2kx+2x)ex Pr(kX2kx)ex One consequence is that, if ZN(0,1)k is a Gaussian random vector in k, then as the dimension k grows, the squared length of the vector is concentrated tightly around k with a width k1/2+α:Pr(Z2[k2k1/2+α,k+2k1/2+α+2kα])1ekαwhere the exponent α can be chosen as any value in .

Since the cumulant generating function for χ2(k) is K(t)=k2ln(12t), and its convex dual is K*(q)=12(qk+klnkq), the standard Chernoff bound yieldslnPr(X(1+ε)k)k2(εln(1+ε))lnPr(X(1ε)k)k2(εln(1ε))where 0<ε<1. By the union bound,Pr(X(1±ε)k)12ek2(12ε213ε3)This result is used in proving the Johnson–Lindenstrauss lemma.[14]

Asymptotic properties

File:Chi-square median approx.png
Approximate formula for median (from the Wilson–Hilferty transformation) compared with numerical quantile (top); and difference (blue) and relative difference (red) between numerical quantile and approximate formula (bottom). For the chi-squared distribution, only the positive integer numbers of degrees of freedom (circles) are meaningful.

By the central limit theorem, because the chi-squared distribution is the sum of k independent random variables with finite mean and variance, it converges to a normal distribution for large k. For many practical purposes, for k>50 the distribution is sufficiently close to a normal distribution, so the difference is ignorable.[15] Specifically, if Xχ2(k), then as k tends to infinity, the distribution of (Xk)/2k tends to a standard normal distribution. However, convergence is slow as the skewness is 8/k and the excess kurtosis is 12/k.

The sampling distribution of ln(χ2) converges to normality much faster than the sampling distribution of χ2,[16] as the logarithmic transform removes much of the asymmetry.[17]

Other functions of the chi-squared distribution converge more rapidly to a normal distribution. Some examples are:

  • If Xχ2(k) then 2X is approximately normally distributed with mean 2k1 and unit variance (1922, by R. A. Fisher, see (18.23), p. 426 of Johnson).[4]
  • If Xχ2(k) then X/k3 is approximately normally distributed with mean 129k and variance 29k.[18] This is known as the Wilson–Hilferty transformation, see (18.24), p. 426 of Johnson.[4]
    • This normalizing transformation leads directly to the commonly used median approximation k(129k)3 by back-transforming from the mean, which is also the median, of the normal distribution.

Related distributions

Template:More citations needed section

A chi-squared variable with k degrees of freedom is defined as the sum of the squares of k independent standard normal random variables.

If Y is a k-dimensional Gaussian random vector with mean vector μ and rank k covariance matrix C, then X=(Yμ)TC1(Yμ) is chi-squared distributed with k degrees of freedom.

The sum of squares of statistically independent unit-variance Gaussian variables which do not have mean zero yields a generalization of the chi-squared distribution called the noncentral chi-squared distribution.

If Y is a vector of k i.i.d. standard normal random variables and A is a k×k symmetric, idempotent matrix with rank kn, then the quadratic form YTAY is chi-square distributed with kn degrees of freedom.

If Σ is a p×p positive-semidefinite covariance matrix with strictly positive diagonal entries, then for XN(0,Σ) and w a random p-vector independent of X such that w1++wp=1 and wi0,i=1,,p, then[17]

1w~TΣw~χ12.where w~=(w1/X1,,wp/Xp).

The chi-squared distribution is also naturally related to other distributions arising from the Gaussian. In particular,

  • Y is F-distributed, YF(k1,k2) if Y=X1/k1X2/k2, where X1χk12 and X2χk22 are statistically independent.
  • If X1χk12 and X2χk22 are statistically independent, then X1+X2χk1+k22. If X1 and X2 are not independent, then X1+X2 is not chi-square distributed.

Generalizations

The chi-squared distribution is obtained as the sum of the squares of Template:Mvar independent, zero-mean, unit-variance Gaussian random variables. Generalizations of this distribution can be obtained by summing the squares of other types of Gaussian random variables. Several such distributions are described below.

Linear combination

If X1,,Xn are chi square random variables and a1,,an>0, then the distribution of X=i=1naiXi is a special case of the generalized chi-squared distribution. A closed expression for this distribution is not known. It may be, however, approximated efficiently using the property of characteristic functions of chi-square random variables.[20]

Chi-squared distributions

Noncentral chi-squared distribution

Script error: No such module "Labelled list hatnote". The noncentral chi-squared distribution is obtained from the sum of the squares of independent Gaussian random variables having unit variance and nonzero means.

Generalized chi-squared distribution

Script error: No such module "Labelled list hatnote". The generalized chi-squared distribution is obtained from the quadratic form Template:Math where Template:Mvar is a zero-mean Gaussian vector having an arbitrary covariance matrix, and Template:Mvar is an arbitrary matrix.

Gamma, exponential, and related distributions

The chi-squared distribution Xχk2 is a special case of the gamma distribution, in that XΓ(k2,12) using the rate parameterization of the gamma distribution (or XΓ(k2,2) using the scale parameterization of the gamma distribution) where Template:Mvar is an integer.

Because the exponential distribution is also a special case of the gamma distribution, we also have that if Xχ22, then Xexp(12) is an exponential distribution.

The Erlang distribution is also a special case of the gamma distribution and thus we also have that if Xχk2 with even k, then X is Erlang distributed with shape parameter k/2 and scale parameter 1/2.

Occurrence and applicationsScript error: No such module "anchor".

The chi-squared distribution has numerous applications in inferential statistics, for instance in chi-squared tests and in estimating variances. It enters the problem of estimating the mean of a normally distributed population and the problem of estimating the slope of a regression line via its role in Student's t-distribution. It enters all analysis of variance problems via its role in the F-distribution, which is the distribution of the ratio of two independent chi-squared random variables, each divided by their respective degrees of freedom.

Following are some of the most common situations in which the chi-squared distribution arises from a Gaussian-distributed sample.

  • if X1,...,Xn are i.i.d. N(μ,σ2) random variables, then i=1n(XiX¯)2σ2χn12 where X¯=1ni=1nXi.
  • The box below shows some statistics based on XiN(μi,σi2),i=1,,k independent random variables that have probability distributions related to the chi-squared distribution:
Name Statistic
chi-squared distribution i=1k(Xiμiσi)2
noncentral chi-squared distribution i=1k(Xiσi)2
chi distribution i=1k(Xiμiσi)2
noncentral chi distribution i=1k(Xiσi)2

The chi-squared distribution is also often encountered in magnetic resonance imaging.[21]

Computational methods

Table of Template:Math values vs Template:Math-values

The p-value is the probability of observing a test statistic at least as extreme in a chi-squared distribution. Accordingly, since the cumulative distribution function (CDF) for the appropriate degrees of freedom (df) gives the probability of having obtained a value less extreme than this point, subtracting the CDF value from 1 gives the p-value. A low p-value, below the chosen significance level, indicates statistical significance, i.e., sufficient evidence to reject the null hypothesis. A significance level of 0.05 is often used as the cutoff between significant and non-significant results.

The table below gives a number of p-values matching to χ2 for the first 10 degrees of freedom.

Degrees ofTemplate:Pbfreedom (df) χ2 value[22]
1 0.004 0.02 0.06 0.15 0.46 1.07 1.64 2.71 3.84 6.63 10.83
2 0.10 0.21 0.45 0.71 1.39 2.41 3.22 4.61 5.99 9.21 13.82
3 0.35 0.58 1.01 1.42 2.37 3.66 4.64 6.25 7.81 11.34 16.27
4 0.71 1.06 1.65 2.20 3.36 4.88 5.99 7.78 9.49 13.28 18.47
5 1.14 1.61 2.34 3.00 4.35 6.06 7.29 9.24 11.07 15.09 20.52
6 1.63 2.20 3.07 3.83 5.35 7.23 8.56 10.64 12.59 16.81 22.46
7 2.17 2.83 3.82 4.67 6.35 8.38 9.80 12.02 14.07 18.48 24.32
8 2.73 3.49 4.59 5.53 7.34 9.52 11.03 13.36 15.51 20.09 26.12
9 3.32 4.17 5.38 6.39 8.34 10.66 12.24 14.68 16.92 21.67 27.88
10 3.94 4.87 6.18 7.27 9.34 11.78 13.44 15.99 18.31 23.21 29.59
p-valueTemplate:Pb(probability) 0.95 0.90 0.80 0.70 0.50 0.30 0.20 0.10 0.05 0.01 0.001

These values can be calculated evaluating the quantile function (also known as "inverse CDF" or "ICDF") of the chi-squared distribution;[23] e. g., the Template:Math ICDF for Template:Math and Template:Math yields Template:Math as in the table above, noticing that Template:Math is the p-value from the table.

History

This distribution was first described by the German geodesist and statistician Friedrich Robert Helmert in papers of 1875–6,Template:Sfn[24] where he computed the sampling distribution of the sample variance of a normal population. Thus in German this was traditionally known as the Helmert'sche ("Helmertian") or "Helmert distribution".

The distribution was independently rediscovered by the English mathematician Karl Pearson in the context of goodness of fit, for which he developed his Pearson's chi-squared test, published in 1900, with computed table of values published in Script error: No such module "Footnotes"., collected in Script error: No such module "Footnotes".. The name "chi-square" ultimately derives from Pearson's shorthand for the exponent in a multivariate normal distribution with the Greek letter Chi, writing Template:Math for what would appear in modern notation as Template:Math (Template:Math being the covariance matrix).[25] The idea of a family of "chi-squared distributions", however, is not due to Pearson but arose as a further development due to Fisher in the 1920s.Template:Sfn

See also

Script error: No such module "Portal". Template:Colbegin

Template:Colend

References

Template:Reflist

Sources

Template:Refbegin

  • Script error: No such module "citation/CS1".
  • Script error: No such module "Citation/CS1".
  • Script error: No such module "Citation/CS1".

Template:Refend

Further reading

Template:Refbegin

Template:Refend

External links

Template:ProbDistributions Template:Authority control

  1. Script error: No such module "citation/CS1".
  2. Template:Abramowitz Stegun ref
  3. NIST (2006). Engineering Statistics Handbook – Chi-Squared Distribution
  4. a b c Script error: No such module "citation/CS1".
  5. Script error: No such module "citation/CS1".
  6. Script error: No such module "citation/CS1".
  7. Script error: No such module "citation/CS1".
  8. Script error: No such module "Citation/CS1".
  9. a b Script error: No such module "citation/CS1".
  10. Script error: No such module "Citation/CS1".
  11. Chi-squared distribution, from MathWorld, retrieved Feb. 11, 2009
  12. M. K. Simon, Probability Distributions Involving Gaussian Random Variables, New York: Springer, 2002, eq. (2.35), Template:ISBN
  13. Script error: No such module "Citation/CS1".
  14. MIT 18.S096 (Fall 2015): Topics in Mathematics of Data Science, Lecture 5, Johnson-Lindenstrauss Lemma and Gordons Theorem
  15. Script error: No such module "citation/CS1".
  16. Script error: No such module "Citation/CS1".
  17. a b Script error: No such module "Citation/CS1".
  18. Script error: No such module "Citation/CS1".
  19. Script error: No such module "Citation/CS1".
  20. Script error: No such module "Citation/CS1".
  21. Script error: No such module "Citation/CS1".
  22. Chi-Squared Test Template:Webarchive Table B.2. Dr. Jacqueline S. McLaughlin at The Pennsylvania State University. In turn citing: R. A. Fisher and F. Yates, Statistical Tables for Biological Agricultural and Medical Research, 6th ed., Table IV. Two values have been corrected, 7.82 with 7.81 and 4.60 with 4.61
  23. Script error: No such module "citation/CS1".
  24. F. R. Helmert, "Ueber die Wahrscheinlichkeit der Potenzsummen der Beobachtungsfehler und über einige damit im Zusammenhange stehende Fragen", Zeitschrift für Mathematik und Physik 21, 1876, pp. 192–219
  25. R. L. Plackett, Karl Pearson and the Chi-Squared Test, International Statistical Review, 1983, 61f. See also Jeff Miller, Earliest Known Uses of Some of the Words of Mathematics.
  26. Script error: No such module "Citation/CS1".