naginterfaces.library.contab.chisq

naginterfaces.library.contab.chisq(nobs)

chisq computes ${\chi }^2$ statistics for a two-way contingency table. For a $2\times 2$ table with a small number of observations exact probabilities are computed.

For full information please refer to the NAG Library document for g11aa

https://support.nag.com/numeric/nl/nagdoc_30.3/flhtml/g11/g11aaf.html

Parameters

nobsint, array-like, shape $(\text{nrow},\text{ncol})$

The contingency table $nobs[\text{i}-1,\text{j}-1]$ must contain $n_{\text{i}\text{j}}$, for $\text{j}=1,2,\dots ,c$, for $\text{i}=1,2,\dots ,r$.

Returns

exptfloat, ndarray, shape $(\text{nrow},\text{ncol})$

The table of expected values. $expt[\text{i}-1,\text{j}-1]$ contains $f_{\text{i}\text{j}}$, for $\text{j}=1,2,\dots ,c$, for $\text{i}=1,2,\dots ,r$.

chistfloat, ndarray, shape $(\text{nrow},\text{ncol})$

The table of ${\chi }^2$ contributions. $chist[\text{i}-1,\text{j}-1]$ contains $\frac{{(n_{\text{i}\text{j}}-f_{\text{i}\text{j}})}^2}{f_{\text{i}\text{j}}}$, for $\text{j}=1,2,\dots ,c$, for $\text{i}=1,2,\dots ,r$.

probfloat

If $c=2$, $r=2$ and $n\le 40$ then $prob$ contains the two tail significance level for Fisher’s exact test, otherwise $prob$ contains the significance level from the Pearson ${\chi }^2$ statistic.

chifloat

The Pearson ${\chi }^2$ statistic.

gfloat

The likelihood ratio test statistic.

dffloat

The degrees of freedom for the statistics.

Raises

NagValueError

(errno $1$)

On entry, $\text{ncol}=⟨value⟩$.

Constraint: $\text{ncol}\ge 2$.

(errno $1$)

On entry, $\text{nrow}=⟨value⟩$.

Constraint: $\text{nrow}\ge 2$.

(errno $2$)

On entry, all elements of $nobs=0$.

(errno $2$)

On entry, $i=⟨value⟩$, $j=⟨value⟩$ and $nobs[i-1,j-1]=⟨value⟩$.

Constraint: $nobs[i-1,j-1]\ge 0$.

(errno $3$)

On entry, a $2\times 2$ table has a row or column with both values zero.

Warns

NagAlgorithmicWarning

(errno $4$)

At least one cell has an expected frequency, $f_{ij}\le 0.5$.

Notes

In the NAG Library the traditional C interface for this routine uses a different algorithmic base. Please contact NAG if you have any questions about compatibility.

For a set of $n$ observations classified by two variables, with $r$ and $c$ levels respectively, a two-way table of frequencies with $r$ rows and $c$ columns can be computed.

$$\begin{array}{c}\begin{array}{ccccc}{n}_{11}& n_{12}& \dots & n_{1c}& n_{1.}\\ n_{21}& n_{22}& \dots & n_{2c}& n_{2.}\\ ⋮& ⋮& ⋮& ⋮& ⋮\\ n_{r1}& n_{r2}& \dots & n_{rc}& n_{r.}\\ n_{.1}& n_{.2}& \dots & n_{.c}& n\end{array}\end{array}$$

To measure the association between the two classification variables two statistics that can be used are, the Pearson ${\chi }^2$ statistic, ${\sum }_{i=1}^r{\sum }_{j=1}^c\frac{{(n_{ij}-f_{ij})}^2}{f_{ij}}$, and the likelihood ratio test statistic, $2{\sum }_{i=1}^r{\sum }_{j=1}^c{n}_{ij}\times \log \left(n_{ij}/f_{ij}\right)$, where $f_{ij}$ are the fitted values from the model that assumes the effects due to the classification variables are additive, i.e., there is no association. These values are the expected cell frequencies and are given by

$$f_{ij}=n_{i.}{n}_{.j}/n\text{.}$$

Under the hypothesis of no association between the two classification variables, both these statistics have, approximately, a ${\chi }^2$-distribution with $(c-1)(r-1)$ degrees of freedom. This distribution is arrived at under the assumption that the expected cell frequencies, $f_{ij}$, are not too small. For a discussion of this point see Everitt (1977). He concludes by saying, ‘… in the majority of cases the chi-square criterion may be used for tables with expectations in excess of $0.5$ in the smallest cell’.

In the case of the $2\times 2$ table, i.e., $c=2$ and $r=2$, the ${\chi }^2$ approximation can be improved by using Yates’ continuity correction factor. This decreases the absolute value of $(n_{ij}-f_{ij})$ by $\frac{1}{2}$. For $2\times 2$ tables with a small value of $n$ the exact probabilities from Fisher’s test are computed. These are based on the hypergeometric distribution and are computed using stat.prob_hypergeom. A two tail probability is computed as $min(1,2p_u,2p_l)$, where $p_u$ and $p_l$ are the upper and lower one-tail probabilities from the hypergeometric distribution.

References

Everitt, B S, 1977, The Analysis of Contingency Tables, Chapman and Hall

Kendall, M G and Stuart, A, 1973, The Advanced Theory of Statistics (Volume 2), (3rd Edition), Griffin