NivaarExam PrepOfficial exam papers ↗

23-Ind-A6 Systems Simulation · December 2018

Question 2 of 8: Chi-Squared Goodness-of-Fit Test for the Hypothesized Poisson

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

National Exams — December 2018 — 17-Ind-A6 Systems Simulation. Three-hour, closed-book exam; one of two permitted calculators (Sharp or Casio), one 8.5″×11.0″ aid sheet (both sides). Format: three sections — Part A (Input Modelling): do 2 of 3 questions, 30 marks total; Part B (Modelling Concepts): do 1 of 2, 20 marks; Part C (Output Analysis): do 1 of 2, 20 marks — plus a 1-mark trivia bonus. All seven graded questions and the bonus are solved below for completeness. The exam's own front matter has two internal quirks, transcribed as printed: NOTES item "4." appears twice on page 1, and every page footer reads "17-Ind-A6/Dec. 2019" against a "December 2018" masthead.

Reference texts: Banks, Carson, Nelson & Nicol, Discrete-Event System Simulation (5th ed., Pearson) — input data analysis (ch. 9), random-variate generation, output analysis for a single system and comparing alternative systems (ch. 11–12), verification and validation (ch. 10).

Question 2 (Part A.2): Chi-Squared Goodness-of-Fit Test for the Hypothesized Poisson (10 marks)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

Given. The same 100-patrol survey of Question 1; hypothesized $H_0:$ data $\sim$ Poisson($\lambda=6$); $\alpha=0.05$.

Find. The $\chi^2$ test statistic, its critical value, and a reject/fail-to-reject decision.

Approach. Compute expected counts $E_x = n\cdot P(X=x)$ from the standard Poisson pmf (the exam's own printed formula, $f(x)=\lambda t e^{-\lambda t x}/x!$, is internally non-standard — treated as a transcription slip and replaced with $P(X=k)=\lambda^k e^{-\lambda}/k!$), combine adjacent cells so every expected count is $\ge 5$ (the standard chi-square validity rule), then compute $\chi^2=\sum (O-E)^2/E$ against $\chi^2_{0.95,df}$ with $df=(\text{bins})-1$ — not bins$-2$, because $\lambda=6$ is specified, not estimated from this sample.

  1. Expected counts and binning. $E_x = 100\cdot\dfrac{6^x e^{-6}}{x!}$ for $x=0,\ldots,9$, with the open tail $E_{\ge 10}=100\,P(X\ge 10)$. The raw cells for $x=0,1,2$ (0.25, 1.49, 4.46) and $x=9,10^+$ (6.88, 8.39) fall below 5 and must be pooled with neighbours:
Bin$\le 2$345678$\ge 9$
Observed $O$4512212419114
Expected $E$6.208.9213.3916.0616.0613.7710.3315.28
  1. Test statistic. Eight bins, $df=8-1=7$: $$\chi^2 = \sum\frac{(O-E)^2}{E} = \frac{(4-6.20)^2}{6.20}+\cdots+\frac{(4-15.28)^2}{15.28} = \boxed{18.44}.$$
  2. Critical value and decision. $$\chi^2_{0.95,7} = \boxed{14.07}.$$ Since $18.44 > 14.07$, reject $H_0$ — Robert's hypothesis is not correct: the data do not fit Poisson($\lambda=6$) at the 5% level. The dominant contributors to $\chi^2$ are the $\ge 9$ bin (observed 4 vs. expected 15.3, a large shortfall in the tail) and the $\le 2$ bin (observed 4 vs. expected 6.2) — the real distribution has noticeably less spread than a Poisson with the same rate, consistent with the mean/variance gap flagged as a check note in Question 1 (mean 5.74 well above sample variance 3.02, where a true Poisson requires equality).
QuantityResult
$\chi^2$ statistic18.44
$df$7
$\chi^2_{0.05,7}$ (critical)14.07
DecisionReject $H_0$ — Poisson(6) is not a good fit
Check: the exam's printed Poisson pmf, $f(x)=\lambda t e^{-\lambda tx}/x!$, does not match the standard form and is treated as a transcription slip; the standard pmf $P(X=k)=\lambda^k e^{-\lambda}/k!$ (with $t=1$) is used throughout, consistent with how $\lambda=6$ is used as a per-patrol rate in Question 1.