Pearson chi-squared statistic for contingency tables ($X^2={\sum }_{ij}(O_{ij}-E_{ij}{)}^2/E_{ij}$)

For observed counts $O_{ij}$ and positive expected counts $E_{ij}$ in a contingency table, the statistic is $X^2={\sum }_{ij}(O_{ij}-E_{ij}{)}^2/E_{ij}$. Under common category probabilities or statistical independence, estimate $E_{ij}$ by the product of the row and column totals divided by the grand total. With large expected counts, this gives the test statistic for the Pearson chi-squared test of homogeneity or the Pearson chi-squared test of independence, with asymptotic chi-squared distribution on $(r-1)(c-1)$ degrees of freedom.

Consider independent binomial counts $X_i\sim \mathrm{Bin}(n_i,p)$ under the common-probability null hypothesis, with fixed positive row sizes $n_i$. Let $N={\sum }_i{n}_i$, $K={\sum }_i{X}_i$ and assume $0<K<N$. Conditional on $K$, the probability mass function is the displayed multivariate hypergeometric distribution: in the joint likelihood, the factor $p^K(1-p{)}^{N-K}$ is common to all admissible tables and cancels. The Vandermonde identity gives the denominator, so the conditional distribution is free of the unknown nuisance parameter $p$.

For a table $x$, its Pearson chi-squared statistic isAn exact conditional p-value sums the displayed probabilities over all tables with the same margins and $Q(x)\ge Q(x_{\mathrm{obs}})$. The ordering is by the Pearson chi-squared statistic, rather than by the probability of an individual table; different exact-test orderings can give different answers. Enumeration is practical when $K$ and the number of rows are small. Monte Carlo tests can approximate the same tail by sampling tables with fixed margins.^[ref]

An exact arithmetic comparison is easy too. Let $M$ be the least common multiple of the row sizes. Compare the integersince every table has the same positive denominator $MK(N-K)$. This includes boundary ties without a floating-point tolerance. Conditional p-values are discrete and can be conservative; “exact” concerns their null calibration, not equality of rejection probability with every chosen significance level.

For a two-by-two contingency table, subtracting one half from each absolute observed-minus-expected count difference, with truncation at zero, gives the displayed corrected Pearson chi-squared statistic. Expected counts use the independence model $E_{ij}=O_{i+}{O}_{+j}/N$. The subtraction compensates approximately for using a continuous chi-squared distribution for discrete counts. It can make the test conservative. The corrected statistic should be distinguished from the ordinary uncorrected Pearson chi-squared statistic.