Two-Population Statistical Inference: Mean, Proportion, and Variance Comparison

Foundations of Two-Population Comparative Inference

  • Core Objective of Comparative Inference:

    • Statistical inference frequently extends beyond evaluating a single population parameter to comparing parameters across two distinct populations.

    • Comparative studies aim to determine whether an observed disparity between two sample metrics reflects a genuine underlying difference in the population parameters or simply sampling fluctuation.

  • Primary Parameters Subject to Comparison:

    • Population Means (μ1\mu_1 vs. μ2\mu_2): Evaluated to assess differences in central tendency, continuous measurements, or performance levels across two treatments, environments, or groups.

    • Population Variances (σ12\sigma_1^2 vs. σ22\sigma_2^2): Analyzed to compare dispersion, consistency, or variability between two processes or populations.

    • Population Proportions (P1P_1 vs. P2P_2 or π1\pi_1 vs. π2\pi_2): Evaluated for binary or categorical attributes to compare prevalence, preference, success rates, or defect rates.

  • Dual Methodological Pillars of Comparative Inference:

    • Confidence Intervals: Provide an estimated plausible range of numerical values for the true difference between two population parameters (e.g., μ1−μ2\mu_1 - \mu_2 or P1−P2P_1 - P_2) at a prescribed confidence level (1−α1 - \alpha).

    • Hypothesis Testing: Evaluates probabilistic evidence against a specific null claim regarding the parameters (typically that there is no difference, H0:μ1=μ2H_0: \mu_1 = \mu_2), facilitating formal decision-making under uncertainty.

  • Structural Transition from One-Sample to Two-Sample Design:

    • In a single population setting, a sample of size nn is extracted from a population with mean μ\mu and variance σ2\sigma^2, yielding the sample mean distribution:     Xˉ∼N(μ,σ2n)\bar{X} \sim N\left(\mu, \frac{\sigma^2}{n}\right)

    • In a two-population comparative framework, independent random samples are extracted from two separate populations:

    • Population I has parameter mean μ1\mu_1, variance σ12\sigma_1^2, sample size n1n_1, and sample mean Xˉ1\bar{X}_1.

    • Population II has parameter mean μ2\mu_2, variance σ22\sigma_2^2, sample size n2n_2, and sample mean Xˉ2\bar{X}_2.

    • The analytical focus shifts from individual point estimators to the distribution of their difference (Xˉ1−Xˉ2\bar{X}_1 - \bar{X}_2).

Confidence Intervals for Comparing Two Population Means

  • Theoretical Sampling Distribution of the Difference Between Means:

    • By the reproductive property of independent normal variables (and central limit theorem principles for sufficiently large samples):     Xˉ1∼N(μ1,σ12n1)\bar{X}_1 \sim N\left(\mu_1, \frac{\sigma_1^2}{n_1}\right)     Xˉ2∼N(μ2,σ22n2)\bar{X}_2 \sim N\left(\mu_2, \frac{\sigma_2^2}{n_2}\right)

    • The difference between the two independent sample estimators has the expectation:     E(Xˉ1−Xˉ2)=μ1−μ2E(\bar{X}_1 - \bar{X}_2) = \mu_1 - \mu_2

    • The variance of the difference is the sum of the individual sampling variances:     Var(Xˉ1−Xˉ2)=σ12n1+σ22n2\text{Var}(\bar{X}_1 - \bar{X}_2) = \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}

    • Therefore, the sampling distribution of the difference is:     Xˉ1−Xˉ2∼N(μ1−μ2,σ12n1+σ22n2)\bar{X}_1 - \bar{X}_2 \sim N\left(\mu_1 - \mu_2, \frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}\right)

  • Case 1: Common Known Population Variance (σ12=σ22=σ2\sigma_1^2 = \sigma_2^2 = \sigma^2 Known):

    • Distributional Model: When both populations share a common, known variance σ2\sigma^2, the standard error simplifies to:     σXˉ1−Xˉ2=σ1n1+1n2\sigma_{\bar{X}_1 - \bar{X}_2} = \sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}

    • Standardized Pivot:     Z=(Xˉ1−Xˉ2)−(μ1−μ2)σ1n1+1n2∼N(0,1)Z = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{\sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}} \sim N(0, 1)

    • Two-Sided 95%95\% Confidence Interval Derivation:     P(−1.96≤Z≤1.96)=0.95P\left(-1.96 \le Z \le 1.96\right) = 0.95

    • Interval Bounds Formulation:     [(Xˉ1−Xˉ2)−1.96×σ1n1+1n2,  (Xˉ1−Xˉ2)+1.96×σ1n1+1n2]\left[ (\bar{X}_1 - \bar{X}_2) - 1.96 \times \sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}} ,\; (\bar{X}_1 - \bar{X}_2) + 1.96 \times \sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}} \right]

    • Comprehensive Worked Example (Skill Development Programs):

    • Scenario: Two groups of employees in an organization are sent to two skill development programs. The average performance levels reached at completion are 32.532.5 and 25.725.7. The standard deviation of performance across similar programs is known to be 3.63.6. Group sizes are n1=27n_1 = 27 and n2=32n_2 = 32.

    • Identified Parameters: Xˉ1=32.5\bar{X}_1 = 32.5, Xˉ2=25.7\bar{X}_2 = 25.7, σ=3.6\sigma = 3.6, n1=27n_1 = 27, n2=32n_2 = 32, confidence level = 95%95\% (Z0.025=1.96Z_{0.025} = 1.96).

    • Point Estimate of Mean Difference:       Xˉ1−Xˉ2=32.5−25.7=6.8\bar{X}_1 - \bar{X}_2 = 32.5 - 25.7 = 6.8

    • Margin of Error Calculation:       E=1.96×3.6×127+132=1.96×3.6×0.037037+0.031250=1.96×3.6×0.068287=1.96×3.6×0.261318=1.844E = 1.96 \times 3.6 \times \sqrt{\frac{1}{27} + \frac{1}{32}} = 1.96 \times 3.6 \times \sqrt{0.037037 + 0.031250} = 1.96 \times 3.6 \times \sqrt{0.068287} = 1.96 \times 3.6 \times 0.261318 = 1.844

    • Lower Limit:       32.5−25.7−1.84=4.9632.5 - 25.7 - 1.84 = 4.96

    • Upper Limit:       32.5−25.7+1.84=8.6432.5 - 25.7 + 1.84 = 8.64

    • Result: The 95%95\% confidence interval for the true difference in mean performance levels (μ1−μ2\mu_1 - \mu_2) is [4.96,8.64][4.96, 8.64].

  • Case 2: Common Unknown Population Variance (σ12=σ22=σ2\sigma_1^2 = \sigma_2^2 = \sigma^2 Unknown):

    • Methodological Principle: When population variances are assumed equal but their true value σ2\sigma^2 is unknown, sample variances S12S_1^2 and S22S_2^2 are combined into an unbiased estimator called the pooled sample variance (Sp2S_p^2).

    • Pooled Variance Formulation:     Sp2=(n1−1)S12+(n2−1)S22n1+n2−2S_p^2 = \frac{(n_1 - 1)S_1^2 + (n_2 - 1)S_2^2}{n_1 + n_2 - 2}

    • Pooled Standard Deviation:     Sp=Sp2S_p = \sqrt{S_p^2}

    • Degrees of Freedom: The pooled variance pools squared deviations from two separate samples, resulting in degrees of freedom given by:     df=n1+n2−2\text{df} = n_1 + n_2 - 2

    • Student's tt-Distribution Pivot:     t=(Xˉ1−Xˉ2)−(μ1−μ2)Sp1n1+1n2∼tn1+n2−2t = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{S_p \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}} \sim t_{n_1 + n_2 - 2}

    • General (1−α)%(1 - \alpha)\% Confidence Interval Formulation:     [(Xˉ1−Xˉ2)−tn1+n2−2, α/2×Sp1n1+1n2,  (Xˉ1−Xˉ2)+tn1+n2−2, α/2×Sp1n1+1n2]\left[ (\bar{X}_1 - \bar{X}_2) - t_{n_1 + n_2 - 2,\, \alpha/2} \times S_p \sqrt{\frac{1}{n_1} + \frac{1}{n_2}} ,\; (\bar{X}_1 - \bar{X}_2) + t_{n_1 + n_2 - 2,\, \alpha/2} \times S_p \sqrt{\frac{1}{n_1} + \frac{1}{n_2}} \right]

    • Comprehensive Worked Example (Television Picture Tube Molding Process):

    • Scenario: An executive developed a low-cost molding process for color television picture tube glass stems. To evaluate whether tube lifetime (in hours) decreases under the developed process compared to the existing process, two samples were assessed under the assumption of equal variance:

      • Process 1 (Existing): n1=22n_1 = 22, Xˉ1=3200\bar{X}_1 = 3200, S12=120S_1^2 = 120

      • Process 2 (Developed): n2=20n_2 = 20, Xˉ2=2950\bar{X}_2 = 2950, S22=100S_2^2 = 100

    • Degrees of Freedom & Critical Value: For a 90%90\% confidence level (α=0.10\alpha = 0.10), the degrees of freedom are:       df=22+20−2=40\text{df} = 22 + 20 - 2 = 40       t40, 0.05=1.684t_{40,\, 0.05} = 1.684

    • Pooled Variance Calculation:       Sp2=(22−1)×120+(20−1)×10022+20−2=21×120+19×10040=2520+190040=442040=110.5S_p^2 = \frac{(22 - 1) \times 120 + (20 - 1) \times 100}{22 + 20 - 2} = \frac{21 \times 120 + 19 \times 100}{40} = \frac{2520 + 1900}{40} = \frac{4420}{40} = 110.5       Sp=110.5≈10.51S_p = \sqrt{110.5} \approx 10.51

    • Standard Error of Difference:       SXˉ1−Xˉ2=Sp122+120=10.51×0.045455+0.05=10.51×0.095455=10.51×0.308957≈3.247S_{\bar{X}_1 - \bar{X}_2} = S_p \sqrt{\frac{1}{22} + \frac{1}{20}} = 10.51 \times \sqrt{0.045455 + 0.05} = 10.51 \times \sqrt{0.095455} = 10.51 \times 0.308957 \approx 3.247

    • Margin of Error:       E=1.684×3.247≈5.47E = 1.684 \times 3.247 \approx 5.47

    • Lower Limit:       3200−2950−1.684×10.51×122+120=250−5.47=244.533200 - 2950 - 1.684 \times 10.51 \times \sqrt{\frac{1}{22} + \frac{1}{20}} = 250 - 5.47 = 244.53

    • Upper Limit:       3200−2950+1.684×10.51×122+120=250+5.47=255.473200 - 2950 + 1.684 \times 10.51 \times \sqrt{\frac{1}{22} + \frac{1}{20}} = 250 + 5.47 = 255.47

    • Result: The 90%90\% confidence interval for the decrease in mean tube lifetime (μ1−μ2\mu_1 - \mu_2) is (244.53,255.47)(244.53, 255.47).

  • Case 3: Different Known Population Variances (σ12≠σ22\sigma_1^2 \ne \sigma_2^2, Both Known):

    • Methodological Principle: When population variances are unequal but known independently, no pooling is conducted. Standard errors use the exact known variances.

    • Standard Error:     σXˉ1−Xˉ2=σ12n1+σ22n2\sigma_{\bar{X}_1 - \bar{X}_2} = \sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}

    • General (1−α)%(1 - \alpha)\% Confidence Interval Formulation:     [(Xˉ1−Xˉ2)−Zα/2×σ12n1+σ22n2,  (Xˉ1−Xˉ2)+Zα/2×σ12n1+σ22n2]\left[ (\bar{X}_1 - \bar{X}_2) - Z_{\alpha/2} \times \sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}} ,\; (\bar{X}_1 - \bar{X}_2) + Z_{\alpha/2} \times \sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}} \right]

    • Comprehensive Worked Example (Task Efficiency):

    • Scenario: In a large organization, 1515 employees completed a skills development program, while 1010 did not. Efficiency scores are recorded below with known population variances σ12=3.6\sigma_1^2 = 3.6 and σ22=2.8\sigma_2^2 = 2.8.

      • Group 1 (Trained, n1=15n_1 = 15): 8,5,7,9,5,6,7,6,8,9,6,8,5,7,98, 5, 7, 9, 5, 6, 7, 6, 8, 9, 6, 8, 5, 7, 9

      • Group 2 (Untrained, n2=10n_2 = 10): 7,6,5,8,5,4,6,8,3,57, 6, 5, 8, 5, 4, 6, 8, 3, 5

    • Sample Means Calculation:       Xˉ1=8+5+7+9+5+6+7+6+8+9+6+8+5+7+915=10515=7.0\bar{X}_1 = \frac{8 + 5 + 7 + 9 + 5 + 6 + 7 + 6 + 8 + 9 + 6 + 8 + 5 + 7 + 9}{15} = \frac{105}{15} = 7.0       Xˉ2=7+6+5+8+5+4+6+8+3+510=5710=5.7\bar{X}_2 = \frac{7 + 6 + 5 + 8 + 5 + 4 + 6 + 8 + 3 + 5}{10} = \frac{57}{10} = 5.7

    • Standard Error Calculation:       σXˉ1−Xˉ2=3.615+2.810=0.24+0.28=0.52≈0.7211\sigma_{\bar{X}_1 - \bar{X}_2} = \sqrt{\frac{3.6}{15} + \frac{2.8}{10}} = \sqrt{0.24 + 0.28} = \sqrt{0.52} \approx 0.7211

    • Confidence Level Critical Value: For 90%90\% confidence (α=0.10\alpha = 0.10), Z0.05=1.64Z_{0.05} = 1.64.

    • Lower Limit:       7.0−5.7−1.64×3.615+2.810=1.3−1.64×0.7211=1.3−1.1826=0.127.0 - 5.7 - 1.64 \times \sqrt{\frac{3.6}{15} + \frac{2.8}{10}} = 1.3 - 1.64 \times 0.7211 = 1.3 - 1.1826 = 0.12

    • Upper Limit:       7.0−5.7+1.64×3.615+2.810=1.3+1.1826=2.487.0 - 5.7 + 1.64 \times \sqrt{\frac{3.6}{15} + \frac{2.8}{10}} = 1.3 + 1.1826 = 2.48

    • Result: The 90%90\% confidence interval for the average efficiency improvement (μ1−μ2\mu_1 - \mu_2) is [0.12,2.48][0.12, 2.48].

  • Case 4: Different Unknown Population Variances (σ12≠σ22\sigma_1^2 \ne \sigma_2^2, Both Unknown):

    • Methodological Principle: When variances cannot be assumed equal and remain unknown, each population variance is estimated by its respective sample variance (S12S_1^2 and S22S_2^2). Pooling is strictly invalid.

    • Standard Error:     SXˉ1−Xˉ2=S12n1+S22n2S_{\bar{X}_1 - \bar{X}_2} = \sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}

    • Degrees of Freedom Formulation: Utilizes standard two-sample degrees of freedom approximation:     df=n1+n2−2\text{df} = n_1 + n_2 - 2

    • Test Statistic / Pivot:     t=(Xˉ1−Xˉ2)−(μ1−μ2)S12n1+S22n2∼tn1+n2−2t = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}} \sim t_{n_1 + n_2 - 2}

    • General (1−α)%(1 - \alpha)\% Confidence Interval Formulation:     [(Xˉ1−Xˉ2)−tn1+n2−2, α/2×S12n1+S22n2,  (Xˉ1−Xˉ2)+tn1+n2−2, α/2×S12n1+S22n2]\left[ (\bar{X}_1 - \bar{X}_2) - t_{n_1 + n_2 - 2,\, \alpha/2} \times \sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}} ,\; (\bar{X}_1 - \bar{X}_2) + t_{n_1 + n_2 - 2,\, \alpha/2} \times \sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}} \right]

    • Comprehensive Worked Example (Puppy Food Variety Weights):

    • Scenario: To compare the quality of two food varieties provided for puppies, two samples were fed the respective diets over one month. End-of-month puppy weights (in grams) yielded:

      • Food I: n1=18n_1 = 18, Xˉ1=1250 g\bar{X}_1 = 1250\,\text{g}, S12=240 gS_1^2 = 240\,\text{g}

      • Food II: n2=15n_2 = 15, Xˉ2=950 g\bar{X}_2 = 950\,\text{g}, S22=210 gS_2^2 = 210\,\text{g}

    • Degrees of Freedom & Critical Value: For a 95%95\% confidence interval:       df=18+15−2=31\text{df} = 18 + 15 - 2 = 31       t31, 0.025=2.08t_{31,\, 0.025} = 2.08

    • Standard Error Calculation:       S12n1+S22n2=24018+21015=13.333+14.0=27.333≈5.228\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}} = \sqrt{\frac{240}{18} + \frac{210}{15}} = \sqrt{13.333 + 14.0} = \sqrt{27.333} \approx 5.228

    • Margin of Error:       E=2.08×5.228≈10.87E = 2.08 \times 5.228 \approx 10.87

    • Lower Limit:       1250−950−2.08×24018+21015=300−10.87=289.11250 - 950 - 2.08 \times \sqrt{\frac{240}{18} + \frac{210}{15}} = 300 - 10.87 = 289.1

    • Upper Limit:       1250−950+2.08×24018+21015=300+10.87=310.91250 - 950 + 2.08 \times \sqrt{\frac{240}{18} + \frac{210}{15}} = 300 + 10.87 = 310.9

    • Result: The 95%95\% confidence interval for the weight difference between the two diet groups is [289.1,310.9][289.1, 310.9].

Confidence Intervals for Comparing Two Population Proportions

  • Distributional Foundations of Proportions:

    • Let P1P_1 and P2P_2 denote true population proportions. Independent random samples of size n1n_1 and n2n_2 produce sample proportions π1\pi_1 (or p1p_1) and π2\pi_2 (or p2p_2).

    • By the Normal approximation to the binomial distribution (valid when niPi≥5n_i P_i \ge 5 and ni(1−Pi)≥5n_i (1 - P_i) \ge 5):     π1∼N(P1,P1(1−P1)n1)\pi_1 \sim N\left(P_1, \frac{P_1(1 - P_1)}{n_1}\right)     π2∼N(P2,P2(1−P2)n2)\pi_2 \sim N\left(P_2, \frac{P_2(1 - P_2)}{n_2}\right)

    • The distribution of the difference in sample proportions satisfies:     π1−π2∼N(P1−P2,P1(1−P1)n1+P2(1−P2)n2)\pi_1 - \pi_2 \sim N\left(P_1 - P_2, \frac{P_1(1 - P_1)}{n_1} + \frac{P_2(1 - P_2)}{n_2}\right)

    • The standardized ZZ-pivot is:     Z=(π1−π2)−(P1−P2)P1(1−P1)n1+P2(1−P2)n2∼N(0,1)Z = \frac{(\pi_1 - \pi_2) - (P_1 - P_2)}{\sqrt{\frac{P_1(1 - P_1)}{n_1} + \frac{P_2(1 - P_2)}{n_2}}} \sim N(0, 1)

  • Confidence Interval Formulation for P1−P2P_1 - P_2:

    • Because population parameters P1P_1 and P2P_2 are unknown, they are substituted by sample proportions π1\pi_1 and π2\pi_2 inside the standard error radical.

    • General (1−α)%(1 - \alpha)\% Confidence Interval:     [(π1−π2)−Zα/2×π1(1−π1)n1+π2(1−π2)n2,  (π1−π2)+Zα/2×π1(1−π1)n1+π2(1−π2)n2]\left[ (\pi_1 - \pi_2) - Z_{\alpha/2} \times \sqrt{\frac{\pi_1(1 - \pi_1)}{n_1} + \frac{\pi_2(1 - \pi_2)}{n_2}} ,\; (\pi_1 - \pi_2) + Z_{\alpha/2} \times \sqrt{\frac{\pi_1(1 - \pi_1)}{n_1} + \frac{\pi_2(1 - \pi_2)}{n_2}} \right]

  • Comprehensive Worked Example 1 (Displeasure with Government Policies):

    • Scenario: A survey in Verginia measured displeasure toward government policies by selecting independent random samples from two provinces:

    • Southern province: n1=450n_1 = 450, opposition proportion π1=0.68\pi_1 = 0.68

    • Central province: n2=650n_2 = 650, opposition proportion π2=0.55\pi_2 = 0.55

    • Confidence Level: 95%95\% confidence implies Z0.025=1.96Z_{0.025} = 1.96.

    • Point Estimate of Difference:     π1−π2=0.68−0.55=0.13\pi_1 - \pi_2 = 0.68 - 0.55 = 0.13

    • Standard Error Calculation:     SE=0.68(1−0.68)450+0.55(1−0.55)650=0.68×0.32450+0.55×0.45650=0.2176450+0.2475650\text{SE} = \sqrt{\frac{0.68(1 - 0.68)}{450} + \frac{0.55(1 - 0.55)}{650}} = \sqrt{\frac{0.68 \times 0.32}{450} + \frac{0.55 \times 0.45}{650}} = \sqrt{\frac{0.2176}{450} + \frac{0.2475}{650}}     SE=0.00048356+0.00038077=0.00086433≈0.0294\text{SE} = \sqrt{0.00048356 + 0.00038077} = \sqrt{0.00086433} \approx 0.0294

    • Margin of Error:     E=1.96×0.0294≈0.0576E = 1.96 \times 0.0294 \approx 0.0576

    • Lower Limit:     0.13−1.96×0.0294=0.13−0.058=0.0720.13 - 1.96 \times 0.0294 = 0.13 - 0.058 = 0.072

    • Upper Limit:     0.13+1.96×0.0294=0.13+0.058=0.1880.13 + 1.96 \times 0.0294 = 0.13 + 0.058 = 0.188

    • Result: The 95%95\% confidence interval for the difference between the two population proportions is [0.072,0.188][0.072, 0.188].

  • Comprehensive Worked Example 2 (Late Attendance Decline via Fingerprint Machine):

    • Scenario: Ashwee Private Limited introduced fingerprint machines to reduce late attendance. Employees unwilling to accept the system were permitted to leave with compensation.

    • Initial state (Sample 1): n1=256n_1 = 256, late attendance proportion π1=0.46\pi_1 = 0.46

    • Post-implementation state (Sample 2): 3939 employees left, leaving n2=256−39=217n_2 = 256 - 39 = 217 employees, among whom late attendance dropped to π2=0.27\pi_2 = 0.27

    • Confidence Level: 90%90\% confidence implies Z0.05=1.64Z_{0.05} = 1.64.

    • Point Estimate of Decline:     π1−π2=0.46−0.27=0.19\pi_1 - \pi_2 = 0.46 - 0.27 = 0.19

    • Standard Error Calculation:     SE=0.46(1−0.46)256+0.27(1−0.27)217=0.46×0.54256+0.27×0.73217=0.2484256+0.1971217\text{SE} = \sqrt{\frac{0.46(1 - 0.46)}{256} + \frac{0.27(1 - 0.27)}{217}} = \sqrt{\frac{0.46 \times 0.54}{256} + \frac{0.27 \times 0.73}{217}} = \sqrt{\frac{0.2484}{256} + \frac{0.1971}{217}}     SE=0.0009703+0.0009083=0.0018786≈0.04334\text{SE} = \sqrt{0.0009703 + 0.0009083} = \sqrt{0.0018786} \approx 0.04334

    • Margin of Error:     E=1.64×0.04334≈0.071E = 1.64 \times 0.04334 \approx 0.071

    • Lower Limit:     0.19−1.64×0.04334=0.19−0.071=0.1190.19 - 1.64 \times 0.04334 = 0.19 - 0.071 = 0.119

    • Upper Limit:     0.19+1.64×0.04334=0.19+0.071=0.2610.19 + 1.64 \times 0.04334 = 0.19 + 0.071 = 0.261

    • Result: The 90%90\% confidence interval for the percentage decline in late attendance is [0.119,0.261][0.119, 0.261] (or 11.9%11.9\% to 26.1%26.1\%).

Hypothesis Testing Framework for Comparing Two Population Means

  • General Hypothesis Structures:

    • Two-Tailed Test (Testing for Any Significant Disparity):     H0:μ1=μ2  ⟺  H0:μ1−μ2=0H_0: \mu_1 = \mu_2 \quad \iff \quad H_0: \mu_1 - \mu_2 = 0     H1:μ1≠μ2  ⟺  H1:μ1−μ2≠0H_1: \mu_1 \ne \mu_2 \quad \iff \quad H_1: \mu_1 - \mu_2 \ne 0

    • One-Tailed Test (Right-Tailed / Greater Than):     H0:μ1≤μ2  ⟺  H0:μ1−μ2≤0H_0: \mu_1 \le \mu_2 \quad \iff \quad H_0: \mu_1 - \mu_2 \le 0     H1:μ1>μ2  ⟺  H1:μ1−μ2>0H_1: \mu_1 > \mu_2 \quad \iff \quad H_1: \mu_1 - \mu_2 > 0

    • One-Tailed Test (Left-Tailed / Less Than):     H0:μ1≥μ2  ⟺  H0:μ1−μ2≥0H_0: \mu_1 \ge \mu_2 \quad \iff \quad H_0: \mu_1 - \mu_2 \ge 0     H1:μ1<μ2  ⟺  H1:μ1−μ2<0H_1: \mu_1 < \mu_2 \quad \iff \quad H_1: \mu_1 - \mu_2 < 0

  • Case 1: Common Known Variance (σ2\sigma^2 Known):

    • Null Distribution & Test Statistic: Under H0:μ1=μ2H_0: \mu_1 = \mu_2, the test statistic simplifies to:     Z=Xˉ1−Xˉ2σ1n1+1n2∼N(0,1)Z = \frac{\bar{X}_1 - \bar{X}_2}{\sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}} \sim N(0, 1)

    • Comprehensive Worked Example (Durability of Road Crossing Paints):

    • Scenario: To evaluate the durability of two road-marking paints on passenger crossings, twelve 100 cm2100\,\text{cm}^2 plots were painted with each paint on a heavily traveled road. After one week of use, average durability readings were 10.410.4 and 9.69.6. Historical experience establishes a common variance of σ2=2.4\sigma^2 = 2.4. Test for equality of durability at α=0.10\alpha = 0.10.

    • Hypotheses:       H0:μ1=μ2H_0: \mu_1 = \mu_2       H1:μ1≠μ2H_1: \mu_1 \ne \mu_2

    • Given Values: Xˉ1=10.4\bar{X}_1 = 10.4, Xˉ2=9.6\bar{X}_2 = 9.6, σ2=2.4\sigma^2 = 2.4, n1=12n_1 = 12, n2=12n_2 = 12, α=0.10\alpha = 0.10.

    • Test Statistic Calculation:       Z=10.4−9.62.4(112+112)=0.82.4×212=0.82.4×0.16667=0.80.4=0.80.632456≈1.26Z = \frac{10.4 - 9.6}{\sqrt{2.4 \left(\frac{1}{12} + \frac{1}{12}\right)}} = \frac{0.8}{\sqrt{2.4 \times \frac{2}{12}}} = \frac{0.8}{\sqrt{2.4 \times 0.16667}} = \frac{0.8}{\sqrt{0.4}} = \frac{0.8}{0.632456} \approx 1.26

    • Decision Rule & Critical Value: For a two-tailed test at α=0.10\alpha = 0.10, the benchmark critical value for fail-to-reject bounds is verified against the rejection region (tabulated standard threshold Z0.05=1.645Z_{0.05} = 1.645, with an upper bound test threshold of 1.961.96). Since calculated Z=1.26Z = 1.26 falls well below the critical threshold (1.26<1.645<1.961.26 < 1.645 < 1.96), it does not fall into the critical region.

    • Conclusion: Do not reject H0H_0. There is insufficient statistical evidence to establish that the durability of the two road-marking paints differs significantly.

  • Case 2: Common Unknown Variance (σ2\sigma^2 Unknown, Equal Variance Assumed):

    • Test Statistic Formulation:     t=Xˉ1−Xˉ2Sp1n1+1n2∼tn1+n2−2t = \frac{\bar{X}_1 - \bar{X}_2}{S_p \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}} \sim t_{n_1 + n_2 - 2}     where:     Sp=(n1−1)S12+(n2−1)S22n1+n2−2S_p = \sqrt{\frac{(n_1 - 1)S_1^2 + (n_2 - 1)S_2^2}{n_1 + n_2 - 2}}

    • Comprehensive Worked Example (Mobile Phone Recharge Card Sales):

    • Scenario: To evaluate demand for Rs. 100 recharge cards across North and Eastern provinces during November 2011, sales per shop were recorded:

      • Province I (n1=10n_1 = 10): 12,15,8,21,18,10,5,17,12,1412, 15, 8, 21, 18, 10, 5, 17, 12, 14

      • Province II (n2=8n_2 = 8): 18,12,10,16,20,12,10,618, 12, 10, 16, 20, 12, 10, 6

    • Hypotheses:       H0:μ1=μ2H_0: \mu_1 = \mu_2       H1:μ1≠μ2H_1: \mu_1 \ne \mu_2

    • Sample Statistics:       Xˉ1=13.2,S12=4.83,n1=10\bar{X}_1 = 13.2, \quad S_1^2 = 4.83, \quad n_1 = 10       Xˉ2=13.0,S22=4.66,n2=8\bar{X}_2 = 13.0, \quad S_2^2 = 4.66, \quad n_2 = 8

    • Degrees of Freedom:       df=10+8−2=16\text{df} = 10 + 8 - 2 = 16

    • Pooled Variance Computation:       Sp2=(10−1)(4.83)+(8−1)(4.66)16=9(4.83)+7(4.66)16=43.47+32.6216=76.0916≈4.756S_p^2 = \frac{(10 - 1)(4.83) + (8 - 1)(4.66)}{16} = \frac{9(4.83) + 7(4.66)}{16} = \frac{43.47 + 32.62}{16} = \frac{76.09}{16} \approx 4.756       Sp=4.756≈2.1808S_p = \sqrt{4.756} \approx 2.1808

    • Standard Error Calculation:       SE=Sp110+18=2.1808×0.10+0.125=2.1808×0.225=2.1808×0.47434=1.0345\text{SE} = S_p \sqrt{\frac{1}{10} + \frac{1}{8}} = 2.1808 \times \sqrt{0.10 + 0.125} = 2.1808 \times \sqrt{0.225} = 2.1808 \times 0.47434 = 1.0345

    • Calculated Test Statistic:       t=13.2−13.01.0345=0.21.0345≈0.062t = \frac{13.2 - 13.0}{1.0345} = \frac{0.2}{1.0345} \approx 0.062

    • Decision Rule: Two-tailed critical value at α=0.05\alpha = 0.05 with 1616 degrees of freedom is t16, 0.025=±2.12t_{16,\, 0.025} = \pm 2.12.

    • Decision: The calculated test statistic t=0.062t = 0.062 falls within the non-rejection interval [−2.12,2.12][-2.12, 2.12].

    • Conclusion: Fail to reject H0H_0 at the 5%5\% significance level (95%95\% confidence). The sample data provide no evidence that average mobile card sales per shop differ between the two provinces.

  • Case 3: Different Known Variances (σ12≠σ22\sigma_1^2 \ne \sigma_2^2, Both Known):

    • Test Statistic Formulation:     Z=(Xˉ1−Xˉ2)−(μ1−μ2)σ12n1+σ22n2∼N(0,1)Z = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}} \sim N(0, 1)

    • Under H0:μ1−μ2=0H_0: \mu_1 - \mu_2 = 0, this becomes:       Z=Xˉ1−Xˉ2σ12n1+σ22n2Z = \frac{\bar{X}_1 - \bar{X}_2}{\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}}

    • Comprehensive Worked Example (Daily Production at Manufacturing Plants):

    • Scenario: A manufacturing firm operates two plants producing identical products. Management investigates whether Plant I has significantly higher daily production than Plant II. Operating data from the past month are summarized:

      • Plant I: n1=20n_1 = 20 days operated, Xˉ1=62\bar{X}_1 = 62 units, known σ12=24\sigma_1^2 = 24

      • Plant II: n2=22n_2 = 22 days operated, Xˉ2=60\bar{X}_2 = 60 units, known σ22=36\sigma_2^2 = 36

    • Hypotheses (One-Tailed Test):       H0:μ1≤μ2H_0: \mu_1 \le \mu_2       H1:μ1>μ2H_1: \mu_1 > \mu_2

    • Test Statistic Calculation:       Z=62−602420+3622=21.20+1.6364=22.8364=21.68415≈1.19Z = \frac{62 - 60}{\sqrt{\frac{24}{20} + \frac{36}{22}}} = \frac{2}{\sqrt{1.20 + 1.6364}} = \frac{2}{\sqrt{2.8364}} = \frac{2}{1.68415} \approx 1.19

    • Decision Rule: For a right-tailed test at α=0.10\alpha = 0.10, the critical standard normal value is Z0.10=1.28Z_{0.10} = 1.28.

    • Decision: The calculated Z=1.19Z = 1.19 is less than the critical value 1.281.28 and does not fall in the critical region (Z<1.28Z < 1.28).

    • Conclusion: Fail to reject H0H_0 at the 10%10\% level. Sample data provide no statistically significant evidence that Plant I produces more items daily than Plant II.

  • Case 4: Different Unknown Variances (σ12≠σ22\sigma_1^2 \ne \sigma_2^2, Both Unknown):

    • Test Statistic Formulation:     t=Xˉ1−Xˉ2S12n1+S22n2∼tn1+n2−2t = \frac{\bar{X}_1 - \bar{X}_2}{\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}} \sim t_{n_1 + n_2 - 2}

    • Large Sample Equivalence: When total degrees of freedom exceed 3030 (n1+n2−2>30n_1 + n_2 - 2 > 30), the Student's tt-distribution converges toward the standard normal distribution N(0,1)N(0, 1), allowing the use of normal critical thresholds.

    • Comprehensive Worked Example (Sri Lankan Share Market Equity Prices):

    • Scenario: An investor compares the closing share prices of two companies over the previous month to test whether Company II has a higher mean close price at α=0.05\alpha = 0.05:

      • Company I: n1=20n_1 = 20 trading days, Xˉ1=37.50\bar{X}_1 = 37.50, S12=12.5S_1^2 = 12.5

      • Company II: n2=22n_2 = 22 trading days, Xˉ2=42.50\bar{X}_2 = 42.50, S22=15.0S_2^2 = 15.0

    • Hypotheses (Left-Tailed Test for μ1−μ2\mu_1 - \mu_2):       H0:μ1≥μ2(or μ1−μ2≥0)H_0: \mu_1 \ge \mu_2 \quad (\text{or } \mu_1 - \mu_2 \ge 0)       H1:μ1<μ2(Company II closing price is higher)H_1: \mu_1 < \mu_2 \quad (\text{Company II closing price is higher})

    • Standard Error Calculation:       SE=S12n1+S22n2=12.520+15.022=0.625+0.68182=1.30682≈1.14316\text{SE} = \sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}} = \sqrt{\frac{12.5}{20} + \frac{15.0}{22}} = \sqrt{0.625 + 0.68182} = \sqrt{1.30682} \approx 1.14316

    • Test Statistic Calculation:       t=37.50−42.501.14316=−5.01.14316≈−4.37t = \frac{37.50 - 42.50}{1.14316} = \frac{-5.0}{1.14316} \approx -4.37

    • Degrees of Freedom & Critical Value:       df=20+22−2=40\text{df} = 20 + 22 - 2 = 40

      • Because df=40>30\text{df} = 40 > 30, the standard normal distribution provides the critical threshold. For a one-tailed test at α=0.05\alpha = 0.05, the left-tail critical value is −Z0.05=−1.64-Z_{0.05} = -1.64.

    • Decision: The calculated test statistic t=−4.37t = -4.37 is substantially less than the critical threshold −1.64-1.64 (−4.37<−1.64-4.37 < -1.64), lying deep within the rejection region.

    • Conclusion: Reject H0H_0 at the 5%5\% significance level. There is strong statistical evidence (95%95\% confidence) that the closing share price of Company II is significantly higher than that of Company I.

Systematic Summary of Estimators, Standard Errors, and Distributions

  • Comparison Summary Across Variance Conditions for Two Population Means:

    • Common Known Variance (σ1=σ2=σ\sigma_1 = \sigma_2 = \sigma):

    • Standard Error Formulation:       σ1n1+1n2\sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}

    • Test Statistic / Pivot:       Z=(Xˉ1−Xˉ2)−(μ1−μ2)σ1n1+1n2Z = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{\sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}}

    • Sampling Distribution:       N(0,1)N(0, 1)

    • Confidence Interval:       (Xˉ1−Xˉ2)±Zα/2×σ1n1+1n2(\bar{X}_1 - \bar{X}_2) \pm Z_{\alpha/2} \times \sigma \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}

    • Common Unknown Variance (σ1=σ2=σ\sigma_1 = \sigma_2 = \sigma Unknown):

    • Standard Error Formulation:       Sp1n1+1n2S_p \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}       where:       Sp2=(n1−1)S12+(n2−1)S22n1+n2−2S_p^2 = \frac{(n_1 - 1)S_1^2 + (n_2 - 1)S_2^2}{n_1 + n_2 - 2}

    • Test Statistic / Pivot:       t=(Xˉ1−Xˉ2)−(μ1−μ2)Sp1n1+1n2t = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{S_p \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}}

    • Sampling Distribution:       tn1+n2−2t_{n_1 + n_2 - 2}

    • Confidence Interval:       (Xˉ1−Xˉ2)±tn1+n2−2, α/2×Sp1n1+1n2(\bar{X}_1 - \bar{X}_2) \pm t_{n_1 + n_2 - 2,\, \alpha/2} \times S_p \sqrt{\frac{1}{n_1} + \frac{1}{n_2}}

    • Different Known Variances (σ12≠σ22\sigma_1^2 \ne \sigma_2^2 Known):

    • Standard Error Formulation:       σ12n1+σ22n2\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}

    • Test Statistic / Pivot:       Z=(Xˉ1−Xˉ2)−(μ1−μ2)σ12n1+σ22n2Z = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}}

    • Sampling Distribution:       N(0,1)N(0, 1)

    • Confidence Interval:       (Xˉ1−Xˉ2)±Zα/2×σ12n1+σ22n2(\bar{X}_1 - \bar{X}_2) \pm Z_{\alpha/2} \times \sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}

    • Different Unknown Variances (σ12≠σ22\sigma_1^2 \ne \sigma_2^2 Unknown):

    • Standard Error Formulation:       S12n1+S22n2\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}

    • Test Statistic / Pivot:       t=(Xˉ1−Xˉ2)−(μ1−μ2)S12n1+S22n2t = \frac{(\bar{X}_1 - \bar{X}_2) - (\mu_1 - \mu_2)}{\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}}

    • Sampling Distribution:       tn1+n2−2t_{n_1 + n_2 - 2} (converges to N(0,1)N(0, 1) when n1+n2−2>30n_1 + n_2 - 2 > 30)

    • Confidence Interval:       (Xˉ1−Xˉ2)±tn1+n2−2, α/2×S12n1+S22n2(\bar{X}_1 - \bar{X}_2) \pm t_{n_1 + n_2 - 2,\, \alpha/2} \times \sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}

  • Inference Summary for Difference Between Two Population Proportions (P1−P2P_1 - P_2):

    • Standard Error Formulation:     π1(1−π1)n1+π2(1−π2)n2\sqrt{\frac{\pi_1(1 - \pi_1)}{n_1} + \frac{\pi_2(1 - \pi_2)}{n_2}}

    • Test Statistic / Pivot:     Z=(π1−π2)−(P1−P2)π1(1−π1)n1+π2(1−π2)n2Z = \frac{(\pi_1 - \pi_2) - (P_1 - P_2)}{\sqrt{\frac{\pi_1(1 - \pi_1)}{n_1} + \frac{\pi_2(1 - \pi_2)}{n_2}}}

    • Sampling Distribution:     N(0,1)N(0, 1)

    • Confidence Interval:     (π1−π2)±Zα/2×π1(1−π1)n1+π2(1−π2)n2(\pi_1 - \pi_2) \pm Z_{\alpha/2} \times \sqrt{\frac{\pi_1(1 - \pi_1)}{n_1} + \frac{\pi_2(1 - \pi_2)}{n_2}}