FAQ/Poisson - CBU statistics Wiki

Revision 9 as of 2008-07-23 13:11:41

Clear message
location: FAQ / Poisson

How do I compare observed numbers correct to those expected by chance in a multi-choice task?

Suppose a questionnaire has a set of k possible responses to N questions of which one of the k choices is the correct response. Assuming by chance each of the k responses to each question is equally likely the number expected by chance equals N/k.

Suppose we observe x $$\leq$$ N responses for a patient. We wish to see how likely x responses are given the patient responds at random to k choices on each of N questions.

To do this we can assume the number of correct responses, x, follows a Poisson distribution, Po(m) of general form:

$$P(X=x|\mbox{random responses}) = \frac{\mtext{x}}{x!}etext{-m}$$

where $$\mu$$ is the expected total number of correct responses for N questions. Since this equals N/k we can rewrite the above as:

$$P(X=x|\mbox{random responses}) = \frac{(N/k)text{x}}{x!}etext{-(N/k)}$$

one-tailed p-value = $$P(X \leq x) = \sum_text{1}^text{x}P(X=x|\mbox{random responses})$$ = 0.5(two-tailed p-value)

For large N the Poisson distribution is approximately Normally distributed with mean and variance both equal to N/k so we can analogously obtain a one-tailed p-value as:

$$P(X \leq x) = Probit(\frac{x-(N/k)}{\sqrt{N/k}})$$.

Example

Suppose we have 14 questions each with 3 possible responses of which only one is correct and a patient gets a total score of 3 correct responses. We wish to determine how likely it is we would observe 3 or fewer responses given the patient has responded at random (1 sided p-value).

The expected number of correct responses assuming the patient is responding at random is 14/3=4.67.

Using the Poisson distribution

P(X $$\leq$$ 3) = P(X = 0) + P(X = 1) + P(X = 2) + P(X = 3)

$$ = (1 + 4.67 + \frac{4.67text{2}}{2} + \frac{4.67text{3}}{6})e^text{-4.67}$$

= 0.31.

Since P(X $$\leq$$ 3) = P(X $$\leq$$ 14-3) because the Poisson distribution is symmetric the two-tailed p-value = 0.31x2= 0.628.

We conclude there is no evidence to suggest a score of 3 on a three-choice task of 14 questions differs from chance.

Using the Normal approximation:

$$P(X \leq 3) = Probit \frac{3-4.67}{\sqrt{4.67}} = Probit(-0.772) = 0.22 with a two-sided p-value of 0.44.

We again conclude there is no evidence to suggest a score of 3 on a three-choice task of 14 questions differs from chance.

For N > 30 the Poisson and Normal approximations should closely agree. For N < 30 the Poisson is preferable as it is a discrete distribution assuming, as in this example, only integer values can occur.