Level 3 — ProductionStatistics & Probability — Intermediate

Statistics & Probability — Intermediate

45 minutes50 marksprintable — key stays hidden on paper

Level 3 Paper — Production (Derivations, Reasoning, Explanations)

Time limit: 45 minutes
Total marks: 50

Instructions: Show all working. Where a derivation is requested, build it from first principles — stating axioms/definitions used. Marks are indicated per part.


Question 1 — Derive Bayes' theorem and apply (9 marks)

(a) Starting from the definition of conditional probability, derive Bayes' theorem for two events AA and BB with P(B)>0P(B)>0. State every step and the reason for it. (3)

(b) A diagnostic test detects a disease that affects 2%2\% of a population. The test is 95%95\% sensitive (correctly positive if diseased) and 90%90\% specific (correctly negative if healthy). A randomly chosen person tests positive. Compute the probability they actually have the disease, giving the answer to 4 decimal places. (4)

(c) In one sentence, explain out loud why this posterior probability is so much lower than the test's sensitivity. (2)


Question 2 — Central tendency and dispersion from grouped data (11 marks)

The times (minutes) taken by 40 students to complete a task:

Time (min) 0–10 10–20 20–30 30–40 40–50
Frequency 4 10 14 8 4

(a) Estimate the mean. (2)

(b) Derive the grouped median using the interpolation formula   L+N2Ff×h\;L+\dfrac{\frac{N}{2}-F}{f}\times h. State what each symbol represents. (4)

(c) Estimate the modal value using the modal-class formula. (3)

(d) State whether the distribution is skewed and justify using your three averages. (2)


Question 3 — Variance from first principles (8 marks)

(a) Show algebraically that the population variance identity σ2=1nxi2xˉ2\sigma^2 = \frac{1}{n}\sum x_i^2 - \bar{x}^2 follows from the definition σ2=1n(xixˉ)2\sigma^2 = \frac{1}{n}\sum (x_i-\bar{x})^2. Justify each manipulation. (4)

(b) For the data set {2,4,4,6,9}\{2, 4, 4, 6, 9\}, compute the mean, variance, and standard deviation using your identity. (4)


Question 4 — Binomial distribution: mean & variance derivation (9 marks)

(a) Write the PMF of a binomial random variable XB(n,p)X\sim B(n,p) and state its two parameters. (2)

(b) Derive E[X]=npE[X]=np from first principles (you may use the identity r(nr)=n(n1r1)r\binom{n}{r}=n\binom{n-1}{r-1}). Explain the key step. (4)

(c) A fair six-sided die is rolled 12 times. Let XX be the number of sixes. State E[X]E[X], Var(X)\mathrm{Var}(X), and compute P(X=2)P(X=2) to 4 decimal places. (3)


Question 5 — Counting & binomial theorem (8 marks)

(a) Explain the difference between a permutation and a combination, and derive (nr)=nPrr!\binom{n}{r}=\dfrac{{}^nP_r}{r!}. (3)

(b) In how many ways can the letters of the word BALLOON be arranged so that the two L's are never adjacent? Show your reasoning. (3)

(c) Find the coefficient of x5x^5 in the expansion of (2x3x)9\left(2x-\dfrac{3}{x}\right)^{9} using the general term. (2)


Question 6 — Probability rules & explain-out-loud (5 marks)

(a) State the three Kolmogorov axioms of probability. (2)

(b) Two events AA, BB satisfy P(A)=0.5P(A)=0.5, P(B)=0.4P(B)=0.4, P(AB)=0.7P(A\cup B)=0.7. Determine whether AA and BB are independent, mutually exclusive, or neither. Show the test for each. (3)

Answer keyMark scheme & solutions

Question 1

(a) (3 marks)
By definition of conditional probability: P(AB)=P(AB)P(B)P(A\mid B)=\dfrac{P(A\cap B)}{P(B)} and P(BA)=P(AB)P(A)P(B\mid A)=\dfrac{P(A\cap B)}{P(A)}. (1)
From the second, P(AB)=P(BA)P(A)P(A\cap B)=P(B\mid A)P(A). (1)
Substitute into the first: P(AB)=P(BA)P(A)P(B)P(A\mid B)=\dfrac{P(B\mid A)P(A)}{P(B)}. (1) This is Bayes' theorem.

(b) (4 marks)
Let DD = disease, ++ = positive test.
P(D)=0.02P(D)=0.02, P(+D)=0.95P(+\mid D)=0.95, P(+Dc)=10.90=0.10P(+\mid D^c)=1-0.90=0.10. (1)
Total: P(+)=0.95(0.02)+0.10(0.98)=0.019+0.098=0.117P(+)=0.95(0.02)+0.10(0.98)=0.019+0.098=0.117. (1)
P(D+)=0.0190.117=0.16239P(D\mid +)=\dfrac{0.019}{0.117}=0.16239\ldots (1) 0.1624\approx 0.1624. (1)

(c) (2 marks)
Because the disease is rare (2%), the large healthy population produces many false positives that swamp the smaller number of true positives, so a positive result is still probably a false alarm. (2)


Question 2

(a) (2 marks)
Midpoints: 5,15,25,35,45. fx=4(5)+10(15)+14(25)+8(35)+4(45)=20+150+350+280+180=980\sum fx = 4(5)+10(15)+14(25)+8(35)+4(45)=20+150+350+280+180=980. (1)
Mean =980/40=24.5=980/40=24.5 min. (1)

(b) (4 marks)
N/2=20N/2=20. Cumulative freq: 4,14,28,36,40 → median class 20–30. (1)
Symbols: L=20L=20 (lower boundary), F=14F=14 (cf before class), f=14f=14, h=10h=10. (1)
Median =20+201414×10=20+614(10)=20+\dfrac{20-14}{14}\times10=20+\dfrac{6}{14}(10). (1)
=20+4.2857=24.29=20+4.2857=24.29 min. (1)

(c) (3 marks)
Modal class 20–30: f1=14,f0=10,f2=8,L=20,h=10f_1=14, f_0=10, f_2=8, L=20, h=10. (1)
Mode =L+f1f02f1f0f2×h=20+42818×10=L+\dfrac{f_1-f_0}{2f_1-f_0-f_2}\times h=20+\dfrac{4}{28-18}\times10. (1)
=20+410(10)=24.0=20+\dfrac{4}{10}(10)=24.0 min. (1)

(d) (2 marks)
Mean 24.524.5 > Median 24.2924.29 > Mode 24.024.0: mean>>median>>mode indicates a slight positive (right) skew. (2)


Question 3

(a) (4 marks)
σ2=1n(xixˉ)2=1n(xi22xˉxi+xˉ2)\sigma^2=\frac1n\sum(x_i-\bar x)^2=\frac1n\sum(x_i^2-2\bar x x_i+\bar x^2). (1)
=1nxi22xˉnxi+1n(nxˉ2)=\frac1n\sum x_i^2 -\frac{2\bar x}{n}\sum x_i+\frac1n(n\bar x^2). (1)
Since 1nxi=xˉ\frac1n\sum x_i=\bar x, middle term =2xˉxˉ=2xˉ2=2\bar x\cdot\bar x=2\bar x^2. (1)
=1nxi22xˉ2+xˉ2=1nxi2xˉ2.=\frac1n\sum x_i^2-2\bar x^2+\bar x^2=\frac1n\sum x_i^2-\bar x^2. (1)

(b) (4 marks)
xˉ=(2+4+4+6+9)/5=25/5=5\bar x=(2+4+4+6+9)/5=25/5=5. (1)
xi2=4+16+16+36+81=153\sum x_i^2=4+16+16+36+81=153. (1)
σ2=153/525=30.625=5.6\sigma^2=153/5-25=30.6-25=5.6. (1)
σ=5.62.3664\sigma=\sqrt{5.6}\approx2.3664. (1)


Question 4

(a) (2 marks)
P(X=r)=(nr)pr(1p)nrP(X=r)=\binom{n}{r}p^r(1-p)^{n-r}, r=0,,nr=0,\dots,n. (1) Parameters: nn (trials), pp (success prob). (1)

(b) (4 marks)
E[X]=r=0nr(nr)prqnrE[X]=\sum_{r=0}^n r\binom{n}{r}p^r q^{n-r}; r=0r=0 term vanishes. (1)
Use r(nr)=n(n1r1)r\binom{n}{r}=n\binom{n-1}{r-1}: E[X]=r=1nn(n1r1)prqnrE[X]=\sum_{r=1}^n n\binom{n-1}{r-1}p^r q^{n-r}. (1)
Factor npnp and let k=r1k=r-1: =npk=0n1(n1k)pkqn1k=np\sum_{k=0}^{n-1}\binom{n-1}{k}p^k q^{n-1-k}. (1)
The sum =(p+q)n1=1=(p+q)^{n-1}=1, so E[X]=npE[X]=np. (1)

(c) (3 marks)
n=12,p=1/6n=12,p=1/6: E[X]=12/6=2E[X]=12/6=2. Var=np(1p)=121656=1061.6667\mathrm{Var}=np(1-p)=12\cdot\frac16\cdot\frac56=\frac{10}{6}\approx1.6667. (1)
P(X=2)=(122)(1/6)2(5/6)10=66136(5/6)10P(X=2)=\binom{12}{2}(1/6)^2(5/6)^{10}=66\cdot\frac{1}{36}\cdot(5/6)^{10}. (1)
=0.2961=0.2961. (1)


Question 5

(a) (3 marks)
Permutation = ordered selection; combination = unordered. (1)
Each combination of rr objects can be ordered in r!r! ways, and these together give all permutations: nPr=(nr)r!{}^nP_r=\binom{n}{r}\cdot r!. (1)
Hence (nr)=nPrr!\binom{n}{r}=\dfrac{{}^nP_r}{r!}. (1)

(b) (3 marks)
BALLOON: 7 letters with L×2, O×2. Total arrangements =7!2!2!=1260=\dfrac{7!}{2!2!}=1260. (1)
LL together (glue): 6 units, O×2: 6!2!=360\dfrac{6!}{2!}=360. (1)
Never adjacent =1260360=900=1260-360=900. (1)

(c) (2 marks)
General term Tr+1=(9r)(2x)9r(3/x)r=(9r)29r(3)rx92rT_{r+1}=\binom{9}{r}(2x)^{9-r}(-3/x)^r=\binom{9}{r}2^{9-r}(-3)^r x^{9-2r}. (0.5)
92r=5r=29-2r=5\Rightarrow r=2. (0.5)
Coeff =(92)27(3)2=361289=41472=\binom{9}{2}2^{7}(-3)^2=36\cdot128\cdot9=41472. (1)


Question 6

(a) (2 marks)
(1) P(A)0P(A)\ge0 for all events; (2) P(Ω)=1P(\Omega)=1; (3) for mutually exclusive AiA_i, P(Ai)=P(Ai)P(\bigcup A_i)=\sum P(A_i). (2)

(b) (3 marks)
P(AB)=P(A)+P(B)P(AB)=0.5+0.40.7=0.2P(A\cap B)=P(A)+P(B)-P(A\cup B)=0.5+0.4-0.7=0.2. (1)
Independent test: P(A)P(B)=0.5×0.4=0.2=P(AB)P(A)P(B)=0.5\times0.4=0.2=P(A\cap B) ✓ → independent. (1)
Mutually exclusive would need P(AB)=0P(A\cap B)=0, false. So independent, not mutually exclusive. (1)


[
  {"claim":"Bayes posterior P(D|+)=0.1624","code":"num=Rational(95,100)*Rational(2,100); den=num+Rational(10,100)*Rational(98,100); result=round(float(num/den),4)==0.1624"},
  {"claim":"Grouped mean=24.5, median=24.2857","code":"mean=Rational(980,40); med=20+Rational(20-14,14)*10; result=(mean==Rational(49,2)) and (round(float(med),4)==24.2857)"},
  {"claim":"Variance of {2,4,4,6,9} is 5.6","code":"xs=[2,4,4,6,9]; m=Rational(sum(xs),5); v=Rational(sum(x*x for x in xs),5)-m**2; result=v==Rational(28,5)"},
  {"claim":"Binomial P(X=2), n=12,p=1/6 = 0.2961","code":"p=binomial(12,2)*Rational(1,6)**2*Rational(5,6)**10; result=round(float(p),4)==0.2961"},
  {"claim":"BALLOON non-adjacent L arrangements=900","code":"total=factorial(7)/(factorial(2)*factorial(2)); glue=factorial(6)/factorial(2); result=(total-glue)==900"},
  {"claim":"Coefficient of x^5 in (2x-3/x)^9 is 41472","code":"c=binomial(9,2)*2**7*(-3)**2; result=c==41472"}
]