Statistical Inference About Means and Proportions for Two Populations

Learning Objectives

  • Element 1 (LO 1): Describe and explain the distribution form and parameters of the sampling distribution of the difference between the sample means of two different populations.
  • LO 2: Calculate and interpret the interval estimate at a given level of confidence for a difference between two population means when the standard deviations of the two populations are known.
  • LO 3: Conduct a hypothesis test about the difference between two population means when the standard deviations of the two populations are known.
  • LO 4: Calculate and interpret the interval estimate at a given level of confidence for a difference between two population means when the standard deviations of the two populations are unknown.
  • LO 5: Conduct a hypothesis test about the difference between two population means when the standard deviations of the two populations are unknown.
  • LO 6: Conduct a hypothesis test about the difference between two population means when using a matched sample design.
  • LO 7: Calculate and interpret the interval estimate at a given level of confidence for a difference between two population means when using a matched sample design.
  • LO 8: Describe and explain the distribution form and parameters of the sampling distribution of the difference between the sample proportions of two different populations.
  • LO 9: Calculate and interpret the interval estimate at a given level of confidence for a difference between two population proportions.
  • LO 10: Conduct a hypothesis test about the difference between two population proportions.

Overview of Inferences About Two Populations

  • Statistical inference extends beyond single-sample analyses of a single population mean or proportion to comparative studies involving two distinct populations.
  • The prime parameter of interest is either the difference between two population means, denoted as μ1−μ2\mu_1 - \mu_2, or the difference between two population proportions, denoted as p1−p2p_1 - p_2.
  • Practical scenarios include:
    • Estimating the difference between the mean starting salary for a population of males and the mean starting salary for a population of females.
    • Testing for differences in quality between two suppliers by evaluating the proportion of defective parts produced by Supplier A versus Supplier B.
  • Independent random samples are selected separately from each population: a random sample of size n1n_1 from Population 1 and an independent random sample of size n2n_2 from Population 2.

Inferences About Two Population Means: σ1\sigma_1 and σ2\sigma_2 Known

  • Parameters and Independent Samples:
    • Let μ1\mu_1 denote the mean of Population 1 and μ2\mu_2 denote the mean of Population 2.
    • The population standard deviations σ1\sigma_1 and σ2\sigma_2 are assumed to be known prior to collecting the sample data.
    • Independent random samples of sizes n1n_1 and n2n_2 yield sample means xˉ1\bar{x}_1 and xˉ2\bar{x}_2.
  • Point Estimator:
    • The point estimator of the difference between the two population means is the difference between the sample means:   xˉ1−xˉ2\bar{x}_1 - \bar{x}_2
  • Sampling Distribution of xˉ1−xˉ2\bar{x}_1 - \bar{x}_2:
    • Mean: The expected value of the sampling distribution is:   E(xˉ1−xˉ2)=μ1−μ2E(\bar{x}_1 - \bar{x}_2) = \mu_1 - \mu_2
    • Standard Error (Formula 10.2):σ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}}
    • Distribution Form: The sampling distribution of xˉ1−xˉ2\bar{x}_1 - \bar{x}_2 is normally distributed if both underlying populations follow a normal distribution, or if the sample sizes are sufficiently large (n1≥30n_1 \ge 30 and n2≥30n_2 \ge 30) that the Central Limit Theorem applies.
  • Interval Estimation Form:
    • General interval estimation formula: Point Estimate±Margin of Error\text{Point Estimate} \pm \text{Margin of Error}.
    • Margin of Error (Formula 10.3):Margin of Error=zα/2σxˉ1−xˉ2=zα/2σ12n1+σ22n2\text{Margin of Error} = z_{\alpha/2}\sigma_{\bar{x}_1 - \bar{x}_2} = z_{\alpha/2}\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}
    • Interval Estimate (Formula 10.4):(xˉ1−xˉ2)±zα/2σ12n1+σ22n2(\bar{x}_1 - \bar{x}_2) \pm z_{\alpha/2}\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}
  • Home Styles Furniture Store Application (Interval Estimation):
    • Home Styles operates two furniture retail stores in Buffalo, New York: one located in the inner city (Population 1) and one located in a suburban shopping center (Population 2).
    • Management seeks to analyze demographic differences in customer age.
    • Historical demographic data establish known standard deviations: σ1=9 years\sigma_1 = 9\,\text{years} and σ2=10 years\sigma_2 = 10\,\text{years}.
    • Sample data collected:
    • Inner city: n1=36n_1 = 36, xˉ1=40 years\bar{x}_1 = 40\,\text{years}
    • Suburban: n2=49n_2 = 49, xˉ2=35 years\bar{x}_2 = 35\,\text{years}
    • Point estimate computation:   xˉ1−xˉ2=40−35=5 years\bar{x}_1 - \bar{x}_2 = 40 - 35 = 5\,\text{years}
    • At a 95%95\% confidence level (α=0.05\alpha = 0.05), the critical value is zα/2=z0.025=1.96z_{\alpha/2} = z_{0.025} = 1.96.
    • Margin of error calculation:   1.969236+10249=1.968136+10049=4.06 years1.96\sqrt{\frac{9^2}{36} + \frac{10^2}{49}} = 1.96\sqrt{\frac{81}{36} + \frac{100}{49}} = 4.06\,\text{years}
    • The 95%95\% confidence interval is:   5±4.06 years5 \pm 4.06\,\text{years}
    • Lower limit: 5−4.06=0.94 years5 - 4.06 = 0.94\,\text{years}
    • Upper limit: 5+4.06=9.06 years5 + 4.06 = 9.06\,\text{years}
  • Excel Worksheet Implementation for Interval Estimation (homestyle_mod):
    • Open data file homestyle_mod. Column A stores the age data for 3636 inner city customers; Column B stores the age data for 4949 suburban customers.
    • Descriptive statistics: entered into cells E5:F6.
    • Known population standard deviations: entered into cells E8 and F8.
    • Cell E9 standard error formula:   =SQRT(E8^2/E5 + F8^2/F5)
    • Confidence coefficient of 0.950.95 is entered into cell E11; level of significance α=1−confidence coefficient\alpha = 1 - \text{confidence coefficient} is computed in cell E12.
    • Cell E13 computes the critical value using the function NORMSINV.
    • Cell E14 calculates the margin of error by multiplying cell E13 by cell E9.
    • Cell E16 calculates the point estimate (xˉ1−xˉ2\bar{x}_1 - \bar{x}_2).
    • Cell E17 computes the lower limit (0.940.94); cell E18 computes the upper limit (9.069.06).
    • Template generalization formulas: dynamic ranges in cells E5:F6 replace fixed ranges:
    • Cell E5: =COUNT(A:A)
    • Cell F5: =COUNT(B:B)
    • Cell E6: =AVERAGE(A:A)
    • Cell F6: =AVERAGE(B:B)

Hypothesis Tests About μ1−μ2\mu_1 - \mu_2: σ1\sigma_1 and σ2\sigma_2 Known

  • Hypothesis Testing Formulations:
    • Let d0d_0 represent the hypothesized difference between μ1\mu_1 and μ2\mu_2.
    • Three structural forms are available:
    • Upper-tail test:     H0:μ1−μ2≤d0versusHa:μ1−μ2>d0H_0: \mu_1 - \mu_2 \le d_0 \quad \text{versus} \quad H_a: \mu_1 - \mu_2 > d_0
    • Lower-tail test:     H0:μ1−μ2≥d0versusHa:μ1−μ2<d0H_0: \mu_1 - \mu_2 \ge d_0 \quad \text{versus} \quad H_a: \mu_1 - \mu_2 < d_0
    • Two-tail test:     H0:μ1−μ2=d0versusHa:μ1−μ2≠d0H_0: \mu_1 - \mu_2 = d_0 \quad \text{versus} \quad H_a: \mu_1 - \mu_2 \ne d_0
    • When testing for equality, d0=0d_0 = 0, simplifying the null hypothesis to H0:μ1−μ2=0H_0: \mu_1 - \mu_2 = 0 (μ1=μ2\mu_1 = \mu_2).
  • Test Statistic (Formula 10.5):z=(xˉ1−xˉ2)−d0σ12n1+σ22n2z = \frac{(\bar{x}_1 - \bar{x}_2) - d_0}{\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}}
  • Training Center Quality Evaluation Application:
    • Individuals trained at Center A (Population 1) and Center B (Population 2) take a standardized exam to determine training quality differences.
    • Parameter definition: μ1\mu_1 is the mean exam score for Center A; μ2\mu_2 is the mean exam score for Center B.
    • Hypotheses:   H0:μ1−μ2=0versusHa:μ1−μ2≠0H_0: \mu_1 - \mu_2 = 0 \quad \text{versus} \quad H_a: \mu_1 - \mu_2 \ne 0
    • Historical standardized exams demonstrate score standard deviations close to 10 points10\,\text{points}, setting known parameters: σ1=10\sigma_1 = 10 and σ2=10\sigma_2 = 10.
    • Specified level of significance: α=0.05\alpha = 0.05.
    • Sample data:
    • Center A: n1=30n_1 = 30, xˉ1=82\bar{x}_1 = 82
    • Center B: n2=40n_2 = 40, xˉ2=78\bar{x}_2 = 78
    • Test statistic calculation:   z=(82−78)−010230+10240=410030+10040=1.66z = \frac{(82 - 78) - 0}{\sqrt{\frac{10^2}{30} + \frac{10^2}{40}}} = \frac{4}{\sqrt{\frac{100}{30} + \frac{100}{40}}} = 1.66
    • pp-Value Approach:
    • Standard normal distribution cumulative area for z=1.66z = 1.66 is 0.95150.9515.
    • Upper tail area: 1.0000−0.9515=0.04851.0000 - 0.9515 = 0.0485.
    • Two-tailed pp-value:     p-value=2(0.0485)=0.0970p\text{-value} = 2(0.0485) = 0.0970
    • Rejection rule: Reject H0H_0 if p-value≤αp\text{-value} \le \alpha.
    • Decision: Because p-value=0.0970>0.05p\text{-value} = 0.0970 > 0.05, do not reject H0H_0. There is insufficient evidence to conclude the training centers differ in quality.
    • Critical Value Approach:
    • For α=0.05\alpha = 0.05, zα/2=z0.025=1.96z_{\alpha/2} = z_{0.025} = 1.96.
    • Rejection rule: Reject H0H_0 if z≤−1.96z \le -1.96 or z≥1.96z \ge 1.96.
    • Decision: Since z=1.66z = 1.66 falls inside (−1.96,1.96)(-1.96, 1.96), do not reject H0H_0.
  • Excel Data Analysis Tool for Hypothesis Testing (examscores):
    • Tool: z-Test: Two-Sample for Means.
    • Data file: examscores. Column A contains 3030 scores for Center A; Column B contains 4040 scores for Center B.
    • Variance inputs: σ12=102=100\sigma_1^2 = 10^2 = 100 and σ22=102=100\sigma_2^2 = 10^2 = 100.
    • Procedure:
    1. Click Data tab on the ribbon.
    2. Click Data Analysis in the Analyze group.
    3. Select z-Test: Two-Sample for Means.
    4. Enter A1:A31 in Variable 1 Range.
    5. Enter B1:B41 in Variable 2 Range.
    6. Enter 0 in Hypothesized Mean Difference.
    7. Enter 100 in Variable 1 Variance (known).
    8. Enter 100 in Variable 2 Variance (known).
    9. Check Labels.
    10. Enter 0.05 in Alpha.
    11. Select Output Range and specify D4.
    12. Click OK.
    • Output interpretation:
    • Cell E11: Test statistic z=1.6562z = 1.6562.
    • Cell E12: One-tail pp-value labeled P(Z≤z) one-tailP(Z \le z)\text{ one-tail}.
    • Cell E14: Two-tail pp-value labeled P(Z≤z) two-tail=0.0977P(Z \le z)\text{ two-tail} = 0.0977.
    • Since 0.0977>0.050.0977 > 0.05, do not reject H0H_0.
  • Sample Size Considerations:
    • Random samples with n1≥30n_1 \ge 30 and n2≥30n_2 \ge 30 are adequate in most practical applications.
    • When either sample size is less than 3030, the underlying population distributions must be reasonably assumed to be approximately normal.

Inferences About Two Population Means: σ1\sigma_1 and σ2\sigma_2 Unknown

  • Foundational Conditions:
    • When σ1\sigma_1 and σ2\sigma_2 are unknown, the sample standard deviations s1s_1 and s2s_2 serve as point estimators.
    • The interval estimation and hypothesis testing procedures are based on the tt-distribution rather than the standard normal distribution.
  • Interval Estimation Form (Formula 10.6):(xˉ1−xˉ2)±tα/2s12n1+s22n2(\bar{x}_1 - \bar{x}_2) \pm t_{\alpha/2}\sqrt{\frac{s_1^2}{n_1} + \frac{s_2^2}{n_2}}
  • Degrees of Freedom Formula (Formula 10.7):df=(s12n1+s22n2)21n1−1(s12n1)2+1n2−1(s22n2)2df = \frac{\left(\frac{s_1^2}{n_1} + \frac{s_2^2}{n_2}\right)^2}{\frac{1}{n_1 - 1}\left(\frac{s_1^2}{n_1}\right)^2 + \frac{1}{n_2 - 1}\left(\frac{s_2^2}{n_2}\right)^2}
    • When dfdf is non-integer, round down to the nearest integer to ensure a larger critical value and a more conservative estimate.
  • Clearwater National Bank Application (Interval Estimation):
    • The bank investigates checking account balances across two branches: Cherry Grove (Population 1) and Beechmont (Population 2).
    • Sample parameters:
    • Cherry Grove: n1=28n_1 = 28, xˉ1=$1025\bar{x}_1 = \$1025, s1=150s_1 = 150
    • Beechmont: n2=22n_2 = 22, xˉ2=$910\bar{x}_2 = \$910, s2=125s_2 = 125
    • Degrees of freedom calculation:   df=(150228+125222)2128−1(150228)2+122−1(125222)2=47.8df = \frac{\left(\frac{150^2}{28} + \frac{125^2}{22}\right)^2}{\frac{1}{28 - 1}\left(\frac{150^2}{28}\right)^2 + \frac{1}{22 - 1}\left(\frac{125^2}{22}\right)^2} = 47.8
    • Rounded down to integer: df=47df = 47.
    • At 95%95\% confidence, the critical value for 47 df47\,df is t0.025=2.012t_{0.025} = 2.012.
    • Point estimate:   xˉ1−xˉ2=1025−910=$115\bar{x}_1 - \bar{x}_2 = 1025 - 910 = \$115
    • Margin of error calculation:   2.012150228+125222=2.012(38.96)=$782.012\sqrt{\frac{150^2}{28} + \frac{125^2}{22}} = 2.012(38.96) = \$78
    • Interval bounds: 115±78115 \pm 78
    • Lower limit: 115−78=$37115 - 78 = \$37
    • Upper limit: 115+78=$193115 + 78 = \$193
  • Excel Worksheet Implementation (Clearwater National Bank / Figure 10.5):
    • Raw data entered in Column A (Cherry Grove) and Column B (Beechmont).
    • Dynamic formulas for summary statistics:
    • Sample sizes: E5 = COUNT(A:A), F5 = COUNT(B:B)
    • Sample means: E6 = AVERAGE(A:A), F6 = AVERAGE(B:B)
    • Sample standard deviations: E7 = STDEV.S(A:A), F7 = STDEV.S(B:B)
    • Confidence coefficient of 0.950.95 is entered into cell E12.
    • Point estimate in cell E18: 115115.
    • Lower limit in cell E19: 3737.
    • Upper limit in cell E20: 193193.
  • Hypothesis Testing Formulations (σ1,σ2\sigma_1, \sigma_2 Unknown):
    • Test Statistic (Formula 10.8):t=(xˉ1−xˉ2)−d0s12n1+s22n2t = \frac{(\bar{x}_1 - \bar{x}_2) - d_0}{\sqrt{\frac{s_1^2}{n_1} + \frac{s_2^2}{n_2}}}
    • Degrees of freedom determined via Formula 10.7.
  • Software Evaluation Study Application (Hypothesis Testing):
    • A study evaluates whether a new software package reduces information systems project completion time relative to current technology.
    • Populations: Systems analysts using current technology (Population 1) versus analysts using the new software package (Population 2).
    • Means: μ1\mu_1 is the mean completion time for current technology; μ2\mu_2 is the mean completion time for the new package.
    • Research hypothesis posits that the new software will yield a shorter completion time (μ2<μ1\mu_2 < \mu_1, which means μ1−μ2>0\mu_1 - \mu_2 > 0).
    • Stated hypotheses:   H0:μ1−μ2≤0versusHa:μ1−μ2>0H_0: \mu_1 - \mu_2 \le 0 \quad \text{versus} \quad H_a: \mu_1 - \mu_2 > 0
    • Significance level: α=0.05\alpha = 0.05.
    • Sample data (from 2424 analysts, 1212 per group):
    • Current technology: n1=12n_1 = 12, xˉ1=325 hours\bar{x}_1 = 325\,\text{hours}, s1=40s_1 = 40
    • New software: n2=12n_2 = 12, xˉ2=286 hours\bar{x}_2 = 286\,\text{hours}, s2=44s_2 = 44
    • Test statistic computation:   t=(325−286)−040212+44212=39160012+193612=2.27t = \frac{(325 - 286) - 0}{\sqrt{\frac{40^2}{12} + \frac{44^2}{12}}} = \frac{39}{\sqrt{\frac{1600}{12} + \frac{1936}{12}}} = 2.27
    • Degrees of freedom computation:   df=(40212+44212)2112−1(40212)2+112−1(44212)2=21.8df = \frac{\left(\frac{40^2}{12} + \frac{44^2}{12}\right)^2}{\frac{1}{12 - 1}\left(\frac{40^2}{12}\right)^2 + \frac{1}{12 - 1}\left(\frac{44^2}{12}\right)^2} = 21.8
    • Rounding down gives df=21df = 21.
    • From the tt-distribution table with 21 df21\,df, the value t=2.27t = 2.27 places the upper-tail area (pp-value) between 0.010.01 and 0.0250.025.
    • Because the p-value<0.05p\text{-value} < 0.05, reject H0H_0.
    • Conclusion: The research supports the finding that the new software package provides a shorter mean completion time (noted in text phrasing as supporting a slower completion time differential relative to the null baseline).
  • Excel Data Analysis Tool (softwareredist / Figures 10.6 and 10.7):
    • Tool: t-Test: Two-Sample Assuming Unequal Variances.
    • File: softwareredist.
    • Setup: Variable 1 Range = A1:A13, Variable 2 Range = B1:B13, Hypothesized Mean Difference = 0, Labels checked, Alpha = 0.05, Output Range = D1.
    • Tool output:
    • Cell E9: Test statistic t=2.2721t = 2.2721.
    • Cell E10: One-tail pp-value P(T≤t) one-tail=0.0166P(T \le t)\text{ one-tail} = 0.0166.
    • Cell E12: Two-tail pp-value P(T≤t) two-tailP(T \le t)\text{ two-tail}.
    • Decision: Because p-value=0.0166<0.05p\text{-value} = 0.0166 < 0.05, reject H0H_0.
  • Robustness and Sample Size Guidelines:
    • The tt-distribution procedure is robust to non-normality.
    • Total sample size n1+n2≥20n_1 + n_2 \ge 20 with equal or nearly equal sample sizes provides reliable results even with non-normal populations.
    • Larger sample sizes are required if data distributions exhibit substantial skewness or outliers.
    • Sample sizes smaller than 2020 require verified approximate normality of the parent populations.

Inferences About Two Populations: Matched Samples

  • Study Design Comparison:
    • Independent Sample Design: Separate, unrelated random samples are assigned to each treatment/population.
    • Matched Sample Design: Each sampled element provides a paired pair of observations, receiving both treatments in a randomized order.
    • Primary Advantage: Testing treatments under identical conditions (using the exact same subjects) eliminates between-subject variability, reducing sampling error and increasing the precision of the estimator.
  • Hypothesis Formulation for Matched Samples:
    • Let did_i denote the observed difference for pair ii.
    • Let μd\mu_d denote the mean difference for the underlying population.
    • Null and alternative hypotheses:   H0:μd=0versusHa:μd≠0H_0: \mu_d = 0 \quad \text{versus} \quad H_a: \mu_d \ne 0
  • Manufacturing Task Completion Application:
    • A manufacturing company tests two production methods across a random sample of n=6n = 6 workers.
    • Observed paired differences (di=Method 1−Method 2d_i = \text{Method 1} - \text{Method 2}): 0.60.6, −0.2-0.2, 0.50.5, 0.30.3, 0.00.0, 0.60.6.
    • Sample mean difference:   dˉ=∑din=1.86=0.30\bar{d} = \frac{\sum d_i}{n} = \frac{1.8}{6} = 0.30
    • Sample standard deviation of differences:   sd=∑(di−dˉ)2n−1=0.565=0.335s_d = \sqrt{\frac{\sum (d_i - \bar{d})^2}{n - 1}} = \sqrt{\frac{0.56}{5}} = 0.335
    • Assumption: The population of differences is normally distributed.
    • Test Statistic (Formula 10.9):t=dˉ−μdsd/nt = \frac{\bar{d} - \mu_d}{s_d / \sqrt{n}}
    • Degrees of freedom: df=n−1=6−1=5df = n - 1 = 6 - 1 = 5.
    • Calculation:   t=0.30−00.335/6=0.300.1368=2.20t = \frac{0.30 - 0}{0.335 / \sqrt{6}} = \frac{0.30}{0.1368} = 2.20
    • Table lookup for 5 df5\,df shows upper-tail area between 0.0250.025 and 0.050.05.
    • For a two-tailed test, doubling the area yields a pp-value between 0.050.05 and 0.100.10.
    • Decision: Since p-value>0.05p\text{-value} > 0.05, do not reject H0H_0. Exact Excel p-value=0.0795p\text{-value} = 0.0795.
  • Interval Estimation for Matched Samples:
    • General form: dˉ±tα/2sdn\bar{d} \pm t_{\alpha/2}\frac{s_d}{\sqrt{n}}.
    • Critical value for 95%95\% confidence with 5 df5\,df is t0.025=2.571t_{0.025} = 2.571.
    • Calculation:   0.30±2.571(0.3356)=0.30±0.35 minutes0.30 \pm 2.571\left(\frac{0.335}{\sqrt{6}}\right) = 0.30 \pm 0.35\,\text{minutes}
    • Lower limit: 0.30−0.35=−0.05 minutes0.30 - 0.35 = -0.05\,\text{minutes}
    • Upper limit: 0.30+0.35=0.65 minutes0.30 + 0.35 = 0.65\,\text{minutes}
  • Excel Data Analysis Tool (match / Figures 10.8 and 10.9):
    • Tool: t-Test: Paired Two Sample for Means.
    • File: match. Column A contains worker IDs, Column B contains Method 1 times, Column C contains Method 2 times.
    • Setup: Variable 1 Range = B1:B7, Variable 2 Range = C1:C7, Hypothesized Mean Difference = 0, Labels checked, Alpha = 0.05, Output Range = E1.
    • Output interpretation:
    • Cell F11: One-tail pp-value P(T≤t) one-tailP(T \le t)\text{ one-tail}.
    • Cell F13: Two-tail pp-value P(T≤t) two-tail=0.0795P(T \le t)\text{ two-tail} = 0.0795.
    • Decision: Because p-value=0.0795>0.05p\text{-value} = 0.0795 > 0.05, do not reject H0H_0.

Inferences About Two Population Proportions: p1−p2p_1 - p_2

  • Parameters and Estimators:
    • Let p1p_1 denote the proportion in Population 1 and p2p_2 denote the proportion in Population 2.
    • Samples of size n1n_1 and n2n_2 yield sample proportions pˉ1\bar{p}_1 and pˉ2\bar{p}_2.
    • Point Estimator (Formula 10.10):pˉ1−pˉ2\bar{p}_1 - \bar{p}_2
    • Standard Error (Formula 10.11):σpˉ1−pˉ2=p1(1−p1)n1+p2(1−p2)n2\sigma_{\bar{p}_1 - \bar{p}_2} = \sqrt{\frac{p_1(1-p_1)}{n_1} + \frac{p_2(1-p_2)}{n_2}}
    • Normality Assumption: The sampling distribution of pˉ1−pˉ2\bar{p}_1 - \bar{p}_2 is approximately normal if all four sample condition thresholds are met:   n1p1≥5,n1(1−p1)≥5,n2p2≥5,n2(1−p2)≥5n_1 p_1 \ge 5, \quad n_1(1-p_1) \ge 5, \quad n_2 p_2 \ge 5, \quad n_2(1-p_2) \ge 5
  • Interval Estimation:
    • Population proportions are unknown; sample proportions are substituted into the standard error formula.
    • Margin of Error (Formula 10.12):Margin of Error=zα/2pˉ1(1−pˉ1)n1+pˉ2(1−pˉ2)n2\text{Margin of Error} = z_{\alpha/2}\sqrt{\frac{\bar{p}_1(1-\bar{p}_1)}{n_1} + \frac{\bar{p}_2(1-\bar{p}_2)}{n_2}}
    • Confidence Interval (Formula 10.13):(pˉ1−pˉ2)±zα/2pˉ1(1−pˉ1)n1+pˉ2(1−pˉ2)n2(\bar{p}_1 - \bar{p}_2) \pm z_{\alpha/2}\sqrt{\frac{\bar{p}_1(1-\bar{p}_1)}{n_1} + \frac{\bar{p}_2(1-\bar{p}_2)}{n_2}}
  • Tax Preparation Firm Application (Interval Estimation):
    • The firm evaluates erroneous returns across Office 1 (Population 1) and Office 2 (Population 2).
    • Sample data:
    • Office 1: n1=250n_1 = 250, erroneous returns = 3535pˉ1=35250=0.14\bar{p}_1 = \frac{35}{250} = 0.14
    • Office 2: n2=300n_2 = 300, erroneous returns = 2727pˉ2=27300=0.09\bar{p}_2 = \frac{27}{300} = 0.09
    • Point estimate:   pˉ1−pˉ2=0.14−0.09=0.05\bar{p}_1 - \bar{p}_2 = 0.14 - 0.09 = 0.05
    • At 90%90\% confidence (α=0.10\alpha = 0.10), critical value zα/2=z0.05=1.645z_{\alpha/2} = z_{0.05} = 1.645.
    • Margin of error calculation:   1.6450.14(1−0.14)250+0.09(1−0.09)300=1.6450.1204250+0.0819300=0.0451.645\sqrt{\frac{0.14(1 - 0.14)}{250} + \frac{0.09(1 - 0.09)}{300}} = 1.645\sqrt{\frac{0.1204}{250} + \frac{0.0819}{300}} = 0.045
    • Interval estimate bounds: 0.05±0.0450.05 \pm 0.045
    • Lower limit: 0.05−0.045=0.0050.05 - 0.045 = 0.005
    • Upper limit: 0.05+0.045=0.0950.05 + 0.045 = 0.095
  • Excel Worksheet Implementation for Interval Estimation (tax_mod / Figure 10.10):
    • File: tax_mod. Columns A and B contain Yes or No error indicators.
    • Sample size: Cells E5 and F5 use COUNT (250250 and 300300).
    • Response of interest: Type Yes into cells E6 and F6.
    • Error counts: Cells E7 and F7 use COUNTIF.
    • Sample proportions: Computed in cells E8 and F8.
    • Confidence coefficient of 0.900.90 entered into cell E10; significance level α=0.10\alpha = 0.10 computed in cell E11.
    • Critical value zz computed in cell E12 using NORMSINV.
    • Estimated standard error computed in cell E14.
    • Margin of error computed in cell E15.
    • Point estimate computed in cell E17 (0.050.05).
    • Limits: Lower limit in cell E18 (0.00480.0048); Upper limit in cell E19 (0.09520.0952).

Hypothesis Tests About p1−p2p_1 - p_2

  • Formulation of Hypotheses (Testing No Difference):
    • Upper-tail test:   H0:p1−p2≤0versusHa:p1−p2>0H_0: p_1 - p_2 \le 0 \quad \text{versus} \quad H_a: p_1 - p_2 > 0
    • Lower-tail test:   H0:p1−p2≥0versusHa:p1−p2<0H_0: p_1 - p_2 \ge 0 \quad \text{versus} \quad H_a: p_1 - p_2 < 0
    • Two-tail test:   H0:p1−p2=0versusHa:p1−p2≠0H_0: p_1 - p_2 = 0 \quad \text{versus} \quad H_a: p_1 - p_2 \ne 0
  • Pooled Estimator of pp:
    • Under the assumption that H0H_0 is true as an equality, p1=p2=pp_1 = p_2 = p.
    • Standard error under H0H_0 (Formula 10.14):   σpˉ1−pˉ2=p(1−p)(1n1+1n2)\sigma_{\bar{p}_1 - \bar{p}_2} = \sqrt{p(1-p)\left(\frac{1}{n_1} + \frac{1}{n_2}\right)}
    • Because pp is unknown, sample proportions are combined into a pooled estimator (Formula 10.15):   pˉ=n1pˉ1+n2pˉ2n1+n2\bar{p} = \frac{n_1\bar{p}_1 + n_2\bar{p}_2}{n_1 + n_2}
  • Test Statistic (Formula 10.16):z=pˉ1−pˉ2pˉ(1−pˉ)(1n1+1n2)z = \frac{\bar{p}_1 - \bar{p}_2}{\sqrt{\bar{p}(1-\bar{p})\left(\frac{1}{n_1} + \frac{1}{n_2}\right)}}
    • Applicable when sample sizes satisfy: n1p1≥5n_1 p_1 \ge 5, n1(1−p1)≥5n_1(1-p_1) \ge 5, n2p2≥5n_2 p_2 \ge 5, and n2(1−p2)≥5n_2(1-p_2) \ge 5.
  • Tax Preparation Firm Application (Hypothesis Testing):
    • Test whether error rates differ between Office 1 and Office 2 at α=0.10\alpha = 0.10.
    • Hypotheses:   H0:p1−p2=0versusHa:p1−p2≠0H_0: p_1 - p_2 = 0 \quad \text{versus} \quad H_a: p_1 - p_2 \ne 0
    • Data: n1=250n_1 = 250, pˉ1=0.14\bar{p}_1 = 0.14; n2=300n_2 = 300, pˉ2=0.09\bar{p}_2 = 0.09.
    • Pooled estimator computation:   pˉ=250(0.14)+300(0.09)250+300=35+27550=62550=0.1127\bar{p} = \frac{250(0.14) + 300(0.09)}{250 + 300} = \frac{35 + 27}{550} = \frac{62}{550} = 0.1127
    • Test statistic calculation:   z=0.14−0.090.1127(1−0.1127)(1250+1300)=0.050.1000(0.004+0.003333)=1.85z = \frac{0.14 - 0.09}{\sqrt{0.1127(1 - 0.1127)\left(\frac{1}{250} + \frac{1}{300}\right)}} = \frac{0.05}{\sqrt{0.1000(0.004 + 0.003333)}} = 1.85
    • pp-Value Calculation:
    • Cumulative normal area for z=1.85z = 1.85 is 0.96780.9678.
    • Upper tail area: 1.0000−0.9678=0.03221.0000 - 0.9678 = 0.0322.
    • Two-tail pp-value:     p-value=2(0.0322)=0.0644p\text{-value} = 2(0.0322) = 0.0644
    • Decision: Because p-value=0.0644<0.10p\text{-value} = 0.0644 < 0.10, reject H0H_0. Error rates differ significantly between the two offices.
  • Excel Worksheet Implementation for Hypothesis Testing (tax_mod / Figure 10.11):
    • Hypothesized difference entered in cell E10 (00).
    • Point estimate in cell E11 (0.050.05).
    • Pooled estimate pˉ\bar{p} computed in cell E13 (0.11270.1127).
    • Standard error estimate computed in cell E14 via Formula 10.14.
    • Test statistic computed in cell E15 via formula:   =(E11 - E10)/E14   (yields z=1.8462z = 1.8462).
    • Cumulative probabilities via NORM.S.DIST:
    • Cell E17: Lower-tail pp-value.
    • Cell E18: Upper-tail pp-value.
    • Cell E19: Two-tail pp-value computed as twice the minimum of the one-tailed values (0.06490.0649).
    • Decision: Reject H0H_0 since 0.0649<0.100.0649 < 0.10.

Case Problem: PAR Incorporated

  • Business Context:
    • PAR Incorporated is a manufacturer of golf equipment seeking to increase market share by developing a cut-resistant, durable golf ball coating.
    • Research concerns center on whether the new coating adversely affects driving distance.
    • Driving distance for the new cut-resistant ball must remain comparable to the current ball model.
  • Testing Methodology:
    • A mechanical hitting machine was used to hit balls of both models, eliminating human swinging variation so that distance differences reflect true product differences.
    • Driving distances were measured to the nearest yard and cataloged in file LO four five.
  • Managerial Report Requirements:
    1. Formulate and justify a hypothesis test comparing the driving distances of current and new balls.
    2. Analyze data, report the test statistic, compute the pp-value, state the statistical conclusion, and provide managerial recommendations.
    3. Provide descriptive statistical summaries of driving distance data for each ball model.
    4. Calculate the 95%95\% confidence interval for the population mean driving distance of each model, as well as the 95%95\% confidence interval for the difference between the two population means.
    5. Evaluate and discuss whether larger sample sizes and additional testing are required.

Appendix: Inferences About Two Populations with R

  • Difference Between Two Population Means: σ1\sigma_1 and σ2\sigma_2 Known:
    • Script: Dipmaid's nodes.R.
    • Context: Evaluating Center A (μ1\mu_1) versus Center B (μ2\mu_2) standardized exam scores with known parameters σ1=10\sigma_1 = 10 and σ2=10\sigma_2 = 10.
    • Working file: examscores_r.csv.
    • Procedural steps in RStudio console pane:
    • Step 1: Set the working directory to the folder containing examscores_r.csv.
    • Step 2: Read the CSV file into a data frame.
    • Step 3: Use subset to create a data frame containing only the variable Center.A.
    • Step 4: Use na.omit to remove missing values from center_a_df.
    • Step 5: Use nrow to calculate the sample size for Center A.
    • Step 6: Create a numeric vector of scores for Center A.
    • Step 7: Use mean to compute the sample mean for Center A.
    • Step 8: Repeat Steps 3 through 7 for Center B.
    • Step 9: Calculate the standard error of the difference between sample means.
    • Step 10: Calculate the test statistic zz.
    • Step 11: Use pnorm to calculate the pp-value.
    • Statistical results: Test statistic z=1.656z = 1.656, two-tailed p-value=0.097p\text{-value} = 0.097.
    • One-tailed modifications:
    • Lower-tail (Ha:μ1−μ2<0H_a: \mu_1 - \mu_2 < 0): adjust Step 11 to calculate lower-tail area directly.
    • Upper-tail (Ha:μ1−μ2>0H_a: \mu_1 - \mu_2 > 0): adjust Step 11 to calculate upper-tail area directly.
  • Difference Between Two Population Means: σ1\sigma_1 and σ2\sigma_2 Unknown:
    • Script: Difmeid's unknowns.R.
    • Context: Clearwater National Bank checking balances comparing Cherry Grove (n1=28n_1 = 28) and Beechmont (n2=22n_2 = 22).
    • Working file: checkscached_r.csv.
    • Procedural steps:
    • Step 1: Set working directory to the location of checkscached_r.csv.
    • Step 2: Read the CSV file into a data frame.
    • Step 3: Create a vector of balances for Beechmont.
    • Step 4: Create a vector of balances for Cherry Grove.
    • Step 5: Use the R function t.test to execute the two-sample hypothesis test.
    • Statistical results: Test statistic t=2.9546t = 2.9546, two-tailed p-value=0.0048p\text{-value} = 0.0048.
    • Alternative directions (lower-tail or upper-tail) are implemented by modifying the alternative argument in t.test in Step 5.
  • Difference Between Two Population Means: Matched Samples:
    • Script: DIPMANES match.R.
    • Context: Production completion times across n=6n = 6 workers for Method 1 and Method 2.
    • Working file: matched_r.csv.
    • Procedural steps:
    • Step 1: Set working directory to the location of matched_r.csv.
    • Step 2: Read the CSV file into a data frame.
    • Step 3: Create a vector of Method 1 completion times.
    • Step 4: Create a vector of Method 2 completion times.
    • Step 5: Use the R function t.test with the argument paired = TRUE to perform the matched-sample hypothesis test.

Summary of Formulas

  • Formula 10.1 — Point Estimator of the Difference Between Two Population Means:xˉ1−xˉ2\bar{x}_1 - \bar{x}_2
  • Formula 10.2 — Standard Error of xˉ1−xˉ2\bar{x}_1 - \bar{x}_2 (σ1,σ2\sigma_1, \sigma_2 Known):σ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}}
  • Formula 10.3 — Margin of Error (σ1,σ2\sigma_1, \sigma_2 Known):zα/2σ12n1+σ22n2z_{\alpha/2}\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}
  • Formula 10.4 — Interval Estimate of μ1−μ2\mu_1 - \mu_2 (σ1,σ2\sigma_1, \sigma_2 Known):(xˉ1−xˉ2)±zα/2σ12n1+σ22n2(\bar{x}_1 - \bar{x}_2) \pm z_{\alpha/2}\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}
  • Formula 10.5 — Test Statistic for μ1−μ2\mu_1 - \mu_2 (σ1,σ2\sigma_1, \sigma_2 Known):z=(xˉ1−xˉ2)−d0σ12n1+σ22n2z = \frac{(\bar{x}_1 - \bar{x}_2) - d_0}{\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}}
  • Formula 10.6 — Interval Estimate of μ1−μ2\mu_1 - \mu_2 (σ1,σ2\sigma_1, \sigma_2 Unknown):(xˉ1−xˉ2)±tα/2s12n1+s22n2(\bar{x}_1 - \bar{x}_2) \pm t_{\alpha/2}\sqrt{\frac{s_1^2}{n_1} + \frac{s_2^2}{n_2}}
  • Formula 10.7 — Degrees of Freedom for Two Independent Random Samples (tt-Distribution):df=(s12n1+s22n2)21n1−1(s12n1)2+1n2−1(s22n2)2df = \frac{\left(\frac{s_1^2}{n_1} + \frac{s_2^2}{n_2}\right)^2}{\frac{1}{n_1 - 1}\left(\frac{s_1^2}{n_1}\right)^2 + \frac{1}{n_2 - 1}\left(\frac{s_2^2}{n_2}\right)^2}
  • Formula 10.8 — Test Statistic for μ1−μ2\mu_1 - \mu_2 (σ1,σ2\sigma_1, \sigma_2 Unknown):t=(xˉ1−xˉ2)−d0s12n1+s22n2t = \frac{(\bar{x}_1 - \bar{x}_2) - d_0}{\sqrt{\frac{s_1^2}{n_1} + \frac{s_2^2}{n_2}}}
  • Formula 10.9 — Test Statistic for Matched Samples:t=dˉ−μdsd/nt = \frac{\bar{d} - \mu_d}{s_d / \sqrt{n}}
  • Formula 10.10 — Point Estimator of the Difference Between Two Population Proportions:pˉ1−pˉ2\bar{p}_1 - \bar{p}_2
  • Formula 10.11 — Standard Error of pˉ1−pˉ2\bar{p}_1 - \bar{p}_2:σpˉ1−pˉ2=p1(1−p1)n1+p2(1−p2)n2\sigma_{\bar{p}_1 - \bar{p}_2} = \sqrt{\frac{p_1(1-p_1)}{n_1} + \frac{p_2(1-p_2)}{n_2}}
  • Formula 10.12 — Margin of Error for p1−p2p_1 - p_2:zα/2pˉ1(1−pˉ1)n1+pˉ2(1−pˉ2)n2z_{\alpha/2}\sqrt{\frac{\bar{p}_1(1-\bar{p}_1)}{n_1} + \frac{\bar{p}_2(1-\bar{p}_2)}{n_2}}
  • Formula 10.13 — Interval Estimate of p1−p2p_1 - p_2:(pˉ1−pˉ2)±zα/2pˉ1(1−pˉ1)n1+pˉ2(1−pˉ2)n2(\bar{p}_1 - \bar{p}_2) \pm z_{\alpha/2}\sqrt{\frac{\bar{p}_1(1-\bar{p}_1)}{n_1} + \frac{\bar{p}_2(1-\bar{p}_2)}{n_2}}
  • Formula 10.14 — Standard Error of pˉ1−pˉ2\bar{p}_1 - \bar{p}_2 When p1=p2=pp_1 = p_2 = p:σpˉ1−pˉ2=p(1−p)(1n1+1n2)\sigma_{\bar{p}_1 - \bar{p}_2} = \sqrt{p(1-p)\left(\frac{1}{n_1} + \frac{1}{n_2}\right)}
  • Formula 10.15 — Pooled Estimator of pp:pˉ=n1pˉ1+n2pˉ2n1+n2\bar{p} = \frac{n_1\bar{p}_1 + n_2\bar{p}_2}{n_1 + n_2}
  • Formula 10.16 — Test Statistic for p1−p2p_1 - p_2:z=pˉ1−pˉ2pˉ(1−pˉ)(1n1+1n2)z = \frac{\bar{p}_1 - \bar{p}_2}{\sqrt{\bar{p}(1-\bar{p})\left(\frac{1}{n_1} + \frac{1}{n_2}\right)}}