Chi-Square Tests

In this section, we will study a number of important hypothesis tests that fall under the general term chi-square tests. These are named, as you might guess, because in each case the test statistics has (in the limit) a chi-square distribution. Although there are several different tests in this general category, they all share some common themes:

Topics

The One-Sample Bernoulli Model

Suppose that

X X 1 X 2 X n

is a random sample from the Bernoulli distribution with unknown success parameter

p 01

. Thus, these are independent random variables taking the values 1 and 0 with probabilities

p

and

1 p

respectively. We want to test

p p 0

versus

p p 0

, where

p 0 01

is specified. Of course, we have already studied such tests in the Bernoulli model. But keep in mind that our methods in this section will generalize to a variety of new models that we have not yet studied.

Let

O 1 j 1 n X j

and

O 0 n O 1 j 1 n 1 X j

. These statistics give the number of times (frequency) that outcomes 1 and 0 occur, respectively. Moreover, we know that each has a binomial distribution;

O 1

has parameters

n

and

p

, while

O 0

has parameters

n

and

1 p

. In particular,

O 1 n p

O 0 n 1 p

, and

O 1 O 0 n p 1 p

. Moreover, recall that

O 1

is sufficient for

p

. Thus, any good test statistic should be a function of

O 1

. Next, recall that when

n

is large, the distribution of

O 1

is approximately normal, by the central limit theorem. Let

Z O 1 n p 0 n p 0 1 p 0

Note that

Z

is the standard score of

O 1

under

. Hence if

n

is large,

Z

has approximately the standard normal distribution under

, and therefore

V Z 2

has approximately the chi-square distribution with 1 degree of freedom under

. As usual, let

k

denote the quantile function of the chi-square distribution with

k

degrees of freedom.

Show that an approximate test of versus at the $α$ level of significance is to reject if and only if $V 1 1 α$ .

Show that the test in Exercise 1 is equivalent to the unbiased test with test statistic $Z$ (the approximate normal test) derived in the section on Tests in the Bernoulli model.

For purposes of generalization, the critical result in the next exercise is a special representation of

V

. Let

e 0 n 1 p 0

and let

e 1 n p 0

. Note that these are the expected frequencies of the outcomes 0 and 1, respectively, under

Show that $V O 0 e 0 2 e 0 O 1 e 1 2 e 1$

This representation shows that our test statistic

V

measures the discrepancy between the expected frequencies, under

, and the observed frequencies. Of course, large values of

V

are evidence in favor of

. Finally, note that although there are two terms in the expansion of

V

in Exercise 3, there is only one degree of freedom since

O 0 O 1 n

. The observed and expected frequencies could be stored in a

12

table.

The Multi-Sample Bernoulli Model

Suppose now that we have samples from several (possibly) different, independent Bernoulli trials processes. Specifically, suppose that

X i X i 1 X i 2 X i n i

is a random sample of size

n i

from the Bernoulli distribution with unknown success parameter

p i 01

for each

i 12 m

. Moreover, the samples

X 1 X 2 X m

are independent. We want to test hypotheses about the unknown parameter vector

p p 1 p 2 p m

. There are two common cases that we consider below, but first let's set up the essential notation that we will need for both cases. For

i 12 m

, and

j 01

let

O i j

denote the number of times that outcome

j

occurs in sample

X i

. The observed frequency

O i j

has a binomial distribution;

O i 1

has parameters

n i

and

p i

while

O i 0

has parameters

n i

and

1 p i

The Completely Specified Case

Consider a given parameter vector

p 0 p 01 p 02 p 0 m 01 m

. We want to test the null hypothesis

p p 0

, versus

p p 0

. Since the null hypothesis specifies the value of

p i

for each

i

, this is called the completely specified case. Now let

e i 0 n i 1 p i 0

and let

e i 1 n i p i 0

. These are the expected frequencies of the outcomes 0 and 1, respectively, from sample

X i

under

Use Exercise 3 and independence to show that if $n i$ is large for each $i$ , then $V i 1 m j 01 O i j e i j 2 e i j$ has approximately the chi-square distribution with $m$ degrees of freedom.

As a rule of thumb, large means that we need

e i j 5

for each

i 12 m

and

j 01

. But of course, the larger these expected frequencies the better.

Under the large sample assumption, show that an approximate test of versus at the $α$ level of significance is to reject if and only if $V m 1 α$ .

Once again, note that the test statistic

V

measures the discrepancy between the expected and observed frequencies, over all outcomes and all samples. There are

2 m

terms in the expansion of

V

in Exercise 4, but only

m

degrees of freedom, since

O i 0 O i 1 n i

for each

i 12 m

. The observed and expected frequencies could be stored in a

m 2

table.

The Equal Probability Case

Suppose now that we want to test the null hypothesis

p 1 p 2 p m

that all of the success probabilities are the same, versus the complementary alternative hypothesis

that the probabilities are not all the same. Note, in contrast to the previous model, that the null hypothesis does not specify the value of the common success probability

p

. But note also that under the null hypothesis, the

m

samples can be combined to form one large sample of Bernoulli trials with success probability

p

. Thus, a natural approach is to estimate

p

and then define the test statistic that measures the discrepancy between the expected and observed frequencies, just as before. The challenge will be to find the distribution of the test statistic.

Let

n i 1 m n i

denote the total sample size when the samples are combined. Then the overall sample mean, which in this context is the overall sample proportion of successes, is

P 1 n i 1 m j 1 n i X i j 1 n i 1 m O i 1

The sample proportion

P

is the best estimate of

p

, in just about any sense of the word. Next, let

E i 0 n i 1 P

and let

E i 1 n i P

. These are the estimated expected frequencies of 0 and 1, respectively, from sample

i

under

. Of course these estimated frequencies are now statistics (and hence random) rather than parameters. Just as before, we define our test statistic

V i 1 m j 01 O i j E i j 2 E i j

It turns out that under

, the distribution of

V

converges to the chi-square distribution with

m 1

degrees of freedom as

n

Show that an approximate test of versus at the $α$ level of significance is to reject if and only if $V m 1 1 α$ .

Intuitively, we lost a degree of freedom over the completely specified case because we had to estimate the unknown common success probability

p

. Again, the observed and expected frequencies could be stored in a

m 2

table.

The One-Sample Multinomial Model

Our next model generalizes the one-sample Bernoulli model in a different direction. Suppose that

X X 1 X 2 X n

is a sequence of multinomial trials. Thus, these are independent, identically distributed random variables, each taking values in a set

S

with

k

elements. If we want, we can assume that

S 01 k 1

; the one-sample Bernoulli model then corresponds to

k 2

. Let

f

denote the common probability density function of the sample variables on

S

, so that

f j X i j

for

i 12 n

and

j S

. The values of

f

are assumed unknown, but of course we must have

j S f j 1

, so there are really only

k 1

unknown parameters. For a given probability density function

f 0

S

we want to test

f f 0

versus

f f 0

By this time, our general approach should be clear. We let

O j

denote the number of times that outcome

j

occurs in sample

X

. Note that

O j

has the binomial distribution with parameters

n

and

f j

. Thus,

e j n f 0 j

is the expected number of times that outcome

j

occurs, under

. Out test statistic, of course, is

V j S O j e j 2 e j

It turns out that under

, the distribution of

V

converges to the chi-square distribution with

k 1

degrees of freedom as

n

. Note that there are

k

terms in the expansion of

V

, but only

k 1

degrees of freedom since

i 1 n O i n

Show that an approximate test of versus at the $α$ level of significance is to reject if and only if $V k 1 1 α$ .

Again, as a rule of thumb, we need

e j 5

for each

j S

, but the larger the expected frequencies the better.

The Multi-Sample Multinomial Model

As you might guess, our final generalization is to the multi-sample multinomial model. Specifically, suppose that

X i X i 1 X i 2 X i n i

is a random sample of size

n i

from a distribution on a set

S

with

k

elements, for each

i 12 m

. Moreover, we assume that the samples

X 1 X 2 X m

are independent. Again there is no loss in generality if we take

S 01 k 1

. Then

k 2

reduces to the multi-sample Bernoulli model, and

m 1

corresponds to the one-sample multinomial model.

Let

f i

denote the common probability density function of the variables in sample

X i

, so that

f i j X i l j

for

i 12 m

l 12 n i

, and

j S

. These are generally unknown, so that our vector of parameters is the vector of probability density functions:

f f 1 f 2 f m

. Of course,

j S f i j 1

for

i 12 m

, so there are actually

m k 1

unknown parameters. We are interested in testing hypotheses about

f

. As in the multi-sample Bernoulli model, there are two common cases that we consider below, but first let's set up the essential notation that we will need for both cases. For

i 12 m

and

j S

, let

O i j

denote the number of times that outcome

j

occurs in sample

X i

. The observed frequency

O i j

has a binomial distribution with parameters

n i

and

f i j

The Completely Specified Case

Consider a given vector of probability density functions on

S

, denoted

f 0 f 01 f 02 f 0 m

. We want to test the null hypothesis

f f 0

, versus

f f 0

. Since the null hypothesis specifies the value of

f i j

for each

i

and

j

, this is called the completely specified case. Let

e i j n i f 0 i j

. This is the expected frequency of outcome

j

in sample

X i

under

Use the result from the one-sample multinomial case and independence to show that if $n i$ is large for each $i$ , then $V i 1 m j S O i j e i j 2 e i j$ has approximately the chi-square distribution with $m k 1$ degrees of freedom.

As usual, our rule of thumb is that we need

e i j 5

for each

i 12 m

and

j S

. But of course, the larger these expected frequencies the better.

Under the large sample assumption, show that an approximate test of versus at the $α$ level of significance is to reject if and only if $V m k 1 1 α$ .

As always, the test statistic

V

measures the discrepancy between the expected and observed frequencies, over all outcomes and all samples. There are

m k

terms in the expansion of

V

in Exercise 8, but we lose

m

degrees of freedom, since

j S O i j n i

for each

i 12 m

The Equal PDF Case

Suppose now that we want to test the null hypothesis

f 1 f 2 f m

that all of the probability density functions are the same, versus the complementary alternative hypothesis

that the probability density functions are not all the same. Note, in contrast to the previous model, that the null hypothesis does not specify the value of the common success probability density function

f

. But note also that under the null hypothesis, the

m

samples can be combined to form one large sample of multinomial trials with probability density function

f

. Thus, a natural approach is to estimate the values of

f

and then define the test statistic that measures the discrepancy between the expected and observed frequencies, just as before.

Let

n i 1 m n i

denote the total sample size when the samples are combined. Under

, our best estimate of

f j

P j 1 n i 1 m O i j

Hence our estimate of the expected frequency of outcome

j

in sample

X i

under

E i j n i P j

. Again, this estimated frequency is now a statistic (and hence random) rather than a parameter. Just as before, we define our test statistic

V i 1 m j S O i j E i j 2 E i j

As you no doubt expect by now, it turns out that under

, the distribution of

V

converges to a chi-square distribution as

n

. But let's see if we can determine the degrees of freedom heuristically.

Argue that the limiting distribution of $V$ has $k 1 m 1$ degrees of freedom.

There are $k m$ terms in the expansion of $V$ .
We lose $m$ degrees of freedom since $j S O i j n i$ for each $i 12 m$ .
We must estimate all but one of the probabilities $f j$ for $j S$ , thus losing $k 1$ degrees of freedom.

Show that an approximate test of versus at the $α$ level of significance is to reject if and only if $V k 1 m 1 1 α$ .

A Goodness of Fit Test

A goodness of fit test is an hypothesis test that an unknown sampling distribution is a particular, specified distribution or belongs to a parametric family of distributions. Such tests are clearly fundamental and important. The one-sample multinomial model leads to a quite general goodness of fit test.

To set the stage, suppose that we have an observable random variable

X

for an experiment, taking values in a general set

S

. Random variable

X

might have a continuous or discrete distribution, and might be single-variable or multi-variable. We want to test the null hypothesis that

X

has a given, completely specified distribution, or that the distribution of

X

belongs to a particular parametric family.

Our first step, in either case, is to sample from the distribution of

X

to obtain a sequence of independent, identically distributed variables

X X 1 X 2 X n

. Next, we select

k

and partition

S

into

k

(disjoint) subsets. We will denote the partition by

A j j J

where

J k

. Next, we define the sequence of random variables

Y Y 1 Y 2 Y n

Y i j

if and only if

X i A j

for

i 12 n

and

j J

Argue that $Y$ is a multinomial trials sequence with parameters $n$ and $f$ , where $f j X A j$ for $j J$ .

The Completely Specified Case

Let

H

denote the statement that

X

has a given, completely specified distribution. Let

f 0

denote the probability density function on

J

defined by

f 0 j X A j H

for

j J

. To test hypothesis

H

, we can formally test

f f 0

versus

f f 0

, which of course, is precisely the problem we solved in the one-sample multinomial model.

Generally, we would partition the space

S

into as many subsets as possible, subject to the restriction that the expected frequencies all be at least 5.

The Partially Specified Case

Often we don't really want to test whether

X

has a completely specified distribution (such as the normal distribution with mean 5 and variance 9), but rather whether the distribution of

X

belongs to a specified parametric family (such as the normal). A natural course of action in this case would be to estimate the unknown parameters and then proceed just as above. As we have seen before, the expected frequencies would be statistics

E j

because they would be based on the estimated parameters. As a rule of thumb, we lose a degree of freedom in the chi-square statistic

V

for each parameter that we estimate, although the precise mathematics can be complicated.

A Test of Independence

Suppose that we have observable random variables

X

and

Y

for an experiment, where

X

takes values in a set

S

with

k

elements, and

Y

takes values in a set

T

with

m

elements. Let

f

denote the joint probability density function of

X Y

, so that

f i j X i Y j

for

i S

and

j T

. Recall that the marginal probability density functions of

X

and

Y

are the functions

g

and

h

respectively, where

g i j T f i j, i S

h j i S f i j, j T

Usually, of course,

f

g

, and

h

are unknown. In this section, we are interested in testing whether

X

and

Y

are independent, a basic and important test. Formally then we want to test the null hypothesis

i j S T f i j g i h j

versus the complementary alternative

Our first step, of course, is to draw a random sample

X Y X 1 Y 1 X 2 Y 2 X n Y n

from the distribution of

X Y

. Since the state spaces are finite, this sample forms a sequence of multinomial trials. Thus, with our usual notation, let

O i j

denote the number of times that

i j

occurs in the sample, for each

i j S T

. This statistic has the binomial distribution with trial parameter

n

and success parameter

f i j

. Under

, the success parameter is

g i h j

. However, since we don't know the success parameters, we must estimate them in order to compute the expected frequencies. Our best estimate of

f i j

is the sample proportion

1 n O i j

. Thus, our best estimates of

g i

and

h j

are

1 n N i

and

1 n M j

, respectively, where

N i

is the number of times that

i

occurs in sample

X

and

M j

is the number of times that

j

occurs in sample

Y

N i j T O i j

M j i S O i j

Thus, our estimate of the expected frequency of

i j

under

E i j n 1 n N i 1 n M j 1 n N i M j

Of course, we define our test statistic by

V i S j T O i j E i j 2 E i j

As you now expect, the distribution of

V

converges to a chi-square distribution as

n

. But let's see if we can determine the appropriate degrees of freedom on heuristic grounds.

Argue that the limiting distribution of $V$ has $k 1 m 1$ degrees of freedom.

There are $k m$ terms in the expansion of $V$ .
We lose one degree of freedom since $i S j T O i j n$ ,
We must estimate all but one of the probabilities $g i$ for $i S$ , thus losing $k 1$ degrees of freedom.
We must estimate all but one of the probabilities $h j$ for $j T$ , thus losing $m 1$ degrees of freedom.

Show that an approximate test of versus at the $α$ level of significance is to reject if and only if $V k 1 m 1 1 α$ .

The observed frequencies are often recorded in a

k m

table, known as a contingency table, so that

O i j

is the number in row

i

and column

j

. In this setting, note that

N i

is the sum of the frequencies in the

i

row and

M j

is the sum of the frequencies in the

j

column. Also, for historical reasons, the random variables

X

and

Y

are sometimes called factors and the possible values of the variables categories.

Computational and Simulation Exercises

Computational Exercises

In each of the following exercises, specify the number of degrees of freedom of the chi-square statistic, give the value of the statistic and compute the

P

-value of the test.

A coin is tossed 100 times, resulting in 55 heads. Test the null hypothesis that the coin is fair.

Suppose that we have 3 coins. The coins are tossed, yielding the data in the following table:

	Heads	Tails
Coin 1	29	21
Coin 2	23	17
Coin 3	42	18

Test the null hypothesis that all 3 coin are fair.
Test the null hypothesis that the 3 coins have the same probability of heads.

A die is tossed 240 times, yielding the data in the following table:

Score	1	2	3	4	5	6
Frequency	57	39	28	28	36	52

Test the null hypothesis that the die is fair.
Test the null hypothesis that the die is an ace-six flat die (faces 1 and 6 have probability $1 4$ each while faces 2, 3, 4, and 5 have probability $1 8$ each).

Two dice are tossed, yielding the data in the following table:

Score	1	2	3	4	5	6
Die 1	22	17	22	13	22	24
Die 2	44	24	19	19	18	36

Test the null hypothesis that die 1 is fair and die 2 is an ace-six flat.
Tes the null hypothesis that all the dice have have the same probability distribuiton.

The Buffon trial data set gives the results of 104 repititions of Buffon's needle experiment. The number of crack crossings is 56. In theory, this data set should correspond to 104 Bernoulli trials with success probability $p 2$ . Test to see if this is reasonable.

The number of emissions from a piece of radioactive material was recorded for a sample of 100 one-second intervals. The frequency distribution is given in the following table:

Count	0	1	2	3	4	5	6	7	8	9	10	11
Frequency	3	1	2	12	23	13	20	9	6	6	4	1

Test to see if the data come from a Poisson distribution with parameter 5.
Test to see if the data come from a Poisson distribution generally.

A university classifies faculty by rank as instructors, assistant professors, associate professors, and full professors. The data, by faculty rank and gender, are given in the following contingency table. Test to see if faculty rank and gender are independent.

Faculty	Instructor	Assistant Professor	Associate Professor	Full Professor
Male	62	238	185	115
Female	118	122	123	37

Test to see if Michelson's velocity of light data come from a normal distribution.

Simulation Exercises

In the simulation exercises below, you will be able to explore the goodness of fit test empirically.

In the dice goodness of fit experiment, set the sampling distribution to fair, the sample size to 50, and the significance level to 0.1. Set the test distribution as indicated below and in each case, run the simulation 1000 times. In case (a), give the empirical estimate of the significance level of the test and compare with 0.1. In the other cases, give the empirical estimate of the power of the test. Rank the distributions in (b)-(d) in increasing order of apparent power. Do your results seem reasonable?

fair
ace-six flats
the symmetric, unimodal distribution
the distribution skewed right

In the dice goodness of fit experiment, set the sampling distribution to ace-six flats, the sample size to 50, and the significance level to 0.1. Set the test distribution as indicated below and in each case, run the simulation 1000 times. In case (a), give the empirical estimate of the significance level of the test and compare with 0.1. In the other cases, give the empirical estimate of the power of the test. Rank the distributions in (b)-(d) in increasing order of apparent power. Do your results seem reasonable?

fair
ace-six flats
the symmetric, unimodal distribution
the distribution skewed right

In the dice goodness of fit experiment, set the sampling distribution to the symmetric, unimodal distribution, the sample size to 50, and the significance level to 0.1. Set the test distribution as indicated below and in each case, run the simulation 1000 times. In case (a), give the empirical estimate of the significance level of the test and compare with 0.1. In the other cases, give the empirical estimate of the power of the test. Rank the distributions in (b)-(d) in increasing order of apparent power. Do your results seem reasonable?

the symmetric, unimodal distribution
fair
ace-six flats
the distribution skewed right

In the dice goodness of fit experiment, set the sampling distribution to the distribution skewed right, the sample size to 50, and the significance level to 0.1. Set the test distribution as indicated below and in each case, run the simulation 1000 times. In case (a), give the empirical estimate of the significance level of the test and compare with 0.1. In the other cases, give the empirical estimate of the power of the test. Rank the distributions in (b)-(d) in increasing order of apparent power. Do your results seem reasonable?

the distribution skewed right
fair
ace-six flats
the symmetric, unimodal distribution

Suppose that $D 1$ and $D 2$ are different distributions. Is the power of the test with sampling distribution $D 1$ and test distribution $D 2$ the same as the power of the test with sampling distribution $D 2$ and test distribution $D 1$ ? Make a conjecture based on your results in Exercises 24-26.

In the dice goodness of fit experiment, set the sampling and test distributions to fair and the significance level to 0.05. Run the experiment 1000 times for each of the following sample sizes. In each case, give the empirical estimate of the significance level and compare with 0.05.

$n 10$
$n 20$
$n 40$
$n 100$

In the dice goodness of fit experiment, set the sampling distribution to fair, the test distributions to ace-six flats, and the significance level to 0.05. Run the experiment 1000 times for each of the following sample sizes. In each case, give the empirical estimate of the power of the test. Do the powers seem to be converging?

$n 10$
$n 20$
$n 40$
$n 100$

6. Chi-Square Tests

Topics

The One-Sample Bernoulli Model

The Multi-Sample Bernoulli Model

The Completely Specified Case

The Equal Probability Case

The One-Sample Multinomial Model

The Multi-Sample Multinomial Model

The Completely Specified Case

The Equal PDF Case

A Goodness of Fit Test

The Completely Specified Case

The Partially Specified Case

A Test of Independence

Computational and Simulation Exercises

Computational Exercises

Simulation Exercises