6.5 - Binomial and Poisson tests

6.5 - Binomial and Poisson tests

A very small count may suggest a reduced probability; a very large count may suggest an increased event rate. You will calculate the relevant null tail for binomial and Poisson data, find discrete critical boundaries, and use a normal approximation when its conditions support it. The observation must be compared with all results at least as far into the chosen tail.

Direct binomial tests

If a process historically has failure probability 0.40.4, does just one failure among 1010 independent trials show an improvement? The single probability P(X=1)\mathrm P(X=1) misses an even stronger result: no failures. A lower-tailed test uses P(X≤1)\mathrm P(X\le1).

You already know

For X∼B(n,p)X\sim\mathrm B(n,p), P(X=r)=(nr)pr(1−p)n−r\mathrm P(X=r)=\binom nrp^r(1-p)^{n-r}. The binomial probability formula is in MF19. In a test, substitute the null value of pp, with independent trials and a common success probability.

First state the hypotheses and the null distribution. Then calculate the chosen tail including the observed count, compare with the level, and conclude in context. A binomial “single observation” can be one observed count from a fixed number of trials; it does not mean that n=1n=1.

For a two-tailed equal-tail test at level α\alpha, compare the observed tail with α/2\alpha/2. Equivalently, double that tail and compare with α\alpha when using this convention. Do not switch to a one-tailed test just because the observed result lies on one side.

Worked example 1: a decrease in failure probability

Question. Historically, an independently tested component fails with probability 0.40.4. After a process change, 11 of 1010 components fails. Test at the 5%5\% level whether the failure probability has decreased, using the binomial distribution.

Plan. Use the lower tail under the historical probability, including zero failures.

Let pp be the failure probability after the change. Then

H0:p=0.4,H1:p<0.4.H_0:p=0.4,\qquad H_1:p<0.4.

Under H0H_0, X∼B(10,0.4)X\sim\mathrm B(10,0.4). Thus

P(X≤1)=0.610+(101)(0.4)(0.6)9=0.0463574016.\begin{aligned} \mathrm P(X\le1)&=0.6^{10}+\binom{10}1(0.4)(0.6)^9\\ &=0.0463574016. \end{aligned}

Since 0.0463574<0.050.0463574<0.05, reject H0H_0.

Answer. There is sufficient evidence at the 5%5\% level that the failure probability has decreased below 0.40.4.

Check. P(X=1)=0.040310784\mathrm P(X=1)=0.040310784 is only part of the tail; adding P(X=0)=0.0060466176\mathrm P(X=0)=0.0060466176 gives the required result.

Keep the same data but suppose the investigation was specified in advance as testing for any change. The tail allowance is halved.

Worked example 2: a two-tailed investigation

Question. For the same process and observed count, test H0:p=0.4H_0:p=0.4 against H1:p≠0.4H_1:p\ne0.4 at the 5%5\% level using equal tails.

Plan. Compare the lower-tail probability with 0.0250.025, since departures in both directions matter.

Under H0H_0, X∼B(10,0.4)X\sim\mathrm B(10,0.4), and

P(X≤1)=0.0463574016>0.025.\mathrm P(X\le1)=0.0463574016>0.025.

The observation is not in the lower rejection tail. The equal-tail critical region is actually X=0X=0 or X≥8X\ge8, since

P(X=0)=0.0060466176<0.025,P(X≥8)=0.0122945536<0.025,P(X≥7)=0.0547618816>0.025.\begin{aligned} \mathrm P(X=0)&=0.0060466176<0.025,\\ \mathrm P(X\ge8)&=0.0122945536<0.025,\\ \mathrm P(X\ge7)&=0.0547618816>0.025. \end{aligned}

The adjacent lower boundary fails because P(X≤1)>0.025\mathrm P(X\le1)>0.025.

Answer. Do not reject H0H_0. There is insufficient evidence at the 5%5\% level that the failure probability has changed.

Check. The two rejection tails have total probability 0.01834117120.0183411712, below 0.050.05. They cannot in general be made exactly 2.5%2.5\% each with whole-number counts.

These examples answer different research questions. Their different decisions are a reason to choose the alternative before collecting the data.

Exam tip

Show the probability expression, its numerical value and the comparison. A lone small probability or a single-term calculation does not justify a tail test.

Carry out an upper-tail test with new values.

Direct Poisson tests

A Poisson test follows the same logic, but its parameter is a mean count for a stated interval. If the historical rate is 44 events per hour, a two-hour observation has null mean 88, rather than 44. Use the same interval for the parameter, model and observed count.

You already know

For X∼Po⁡(λ)X\sim\operatorname{Po}(\lambda), P(X=r)=e−λλr/r!\mathrm P(X=r)=e^{-\lambda}\lambda^r/r!; this formula is in MF19. A suitable process has independent, singly occurring events at a constant rate. A longer observation interval scales the mean by its length.

To find a critical region, evaluate tails near the boundary and show that including the next count would exceed the level. In an upper tail,

P(X≥k)=1−P(X≤k−1).\mathrm P(X\ge k)=1-\mathrm P(X\le k-1).

Subtracting P(X≤k)\mathrm P(X\le k) instead omits the boundary value kk.

Worked example 3: a one-hour count

Question. Calls to a service are modelled by a Poisson process with historical rate 44 per hour. Investigate a suspected increase at the 5%5\% level. Find the one-hour critical region and conclude when 99 calls are received.

Plan. Use Po⁡(4)\operatorname{Po}(4) under the null and compare adjacent upper tails.

Let λ\lambda be the mean hourly call count. Use

H0:λ=4,H1:λ>4.H_0:\lambda=4,\qquad H_1:\lambda>4.

For a one-hour count XX, under H0H_0,

P(X≥8)=1−e−4∑r=074rr!=0.0511336,P(X≥9)=1−e−4∑r=084rr!=0.0213634.\begin{aligned} \mathrm P(X\ge8)&=1-e^{-4}\sum_{r=0}^{7}\frac{4^r}{r!}=0.0511336,\\ \mathrm P(X\ge9)&=1-e^{-4}\sum_{r=0}^{8}\frac{4^r}{r!}=0.0213634. \end{aligned}

The first exceeds 0.050.05; the second does not.

Answer. The critical region is X≥9X\ge9, with actual significance probability 0.02140.0214. Nine calls lead to rejection: there is sufficient evidence at the 5%5\% level that the hourly call rate has increased.

Check. The difference between the two tails is P(X=8)\mathrm P(X=8), about 0.02980.0298; including that one extra count would push the probability above 5%5\%.

Keep the hourly hypothesis and the observed count of nine, but double the observation time. The null mean of the observed count doubles too.

Worked example 4: a two-hour count

Question. For the same service, 99 calls are received in a two-hour observation. Carry out the upper-tailed 5%5\% test of the hourly rate.

Plan. Retain H0:λ=4H_0:\lambda=4, H1:λ>4H_1:\lambda>4, but model the two-hour total with mean 2λ2\lambda.

Under H0H_0, Y∼Po⁡(8)Y\sim\operatorname{Po}(8). The observed upper tail is

P(Y≥9)=1−e−8∑r=088rr!=0.4074527>0.05.\mathrm P(Y\ge9)=1-e^{-8}\sum_{r=0}^{8}\frac{8^r}{r!}=0.4074527>0.05.

Do not reject H0H_0.

Answer. There is insufficient evidence at the 5%5\% level that the hourly call rate exceeds 44.

Check. The observed rate is now 9/2=4.59/2=4.5 calls per hour, much closer to the historical rate of 44. A larger count is needed for rejection over two hours: P(Y≥13)=0.0637972>0.05\mathrm P(Y\ge13)=0.0637972>0.05 but P(Y≥14)=0.0341807<0.05\mathrm P(Y\ge14)=0.0341807<0.05.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account