Statistical Process Control

Chapter 2
Statistical Foundations for Process Analysis

Post Graduate Program in Total Quality Management
Statistical Process Control   •   Post Graduate

Statistical Foundations for Process Analysis

Introduction

Statistical process control (SPC) is founded on the idea that process performance can be understood, monitored, and improved through data. Quality professionals rarely observe an entire population of outputs. Instead, they collect samples, summarize observations, model variation, and make decisions under uncertainty. Statistical foundations therefore determine whether a control chart, capability study, process comparison, or improvement project produces a trustworthy conclusion.

This chapter develops the statistical concepts required for process analysis: descriptive statistics, probability distributions, sampling distributions, estimation, confidence intervals, and hypothesis testing. Particular attention is given to assumptions, uncertainty, effect size, and the distinction between statistical significance and practical significance.

Learning Objectives

By the end of this chapter, the reader should be able to:

  1. Distinguish populations, samples, parameters, and statistics.
  2. Summarize process data using appropriate numerical and graphical methods.
  3. Explain common probability distributions used in quality analysis.
  4. Describe sampling distributions and the central limit theorem.
  5. Construct and interpret point estimates and confidence intervals.
  6. Formulate and evaluate statistical hypotheses.
  7. Recognize assumptions, errors, uncertainty, and sources of bias.
  8. Distinguish statistical significance from operational or economic significance.
  9. Apply statistical reasoning to an SPC-oriented process problem.

1. Statistical Thinking in Quality Management

A process is a set of interacting activities that transforms inputs into outputs. Variation is present because of materials, machines, methods, measurement systems, people, environment, and time. Statistical analysis does not eliminate variation; it helps determine whether observed variation is routine, unusual, consequential, or evidence of a process change.

Three principles are especially important.

All processes vary. A single measurement does not establish process performance. Decisions should be based on patterns and distributions.

Variation has causes. Common-cause variation is inherent in the current process design. Special-cause variation arises from identifiable, unusual factors. Statistical evidence can help separate these situations, but investigation and process knowledge remain essential.

Data are conditional on how they were produced. A precise calculation based on biased sampling, unstable measurement, or dependent observations can be misleading. Statistical technique cannot compensate for fundamentally poor data collection.

Populations, samples, parameters, and statistics

A population is the complete set of units or events of interest, such as all bottles produced during a month. A sample is a subset observed from that population. A parameter describes a population, whereas a statistic describes a sample.

QuantityPopulation parameterSample statistic
Meanμx̄
Standard deviationσs
Proportionpp̂
Varianceσ²s²

The objective of inference is to use sample statistics to learn about unknown population parameters while expressing the uncertainty of that learning.

2. Descriptive Statistics for Process Data

Measures of location

For observations x₁, x₂, …, xₙ, the arithmetic mean is:

x̄ = (∑ᵢ₌₁ⁿ xᵢ) / n

Here, n is the number of observations and xᵢ is observation i. The mean is appropriate for quantitative data when extreme values do not dominate interpretation. It has the same units as the observations.

The median is the middle ordered observation, or the average of the two middle observations when n is even. It is resistant to extreme values and is useful for skewed process data. The mode is the most frequent value and is particularly useful for categorical or discrete data.

Measures of dispersion

The range is:

R = xₘₐₓ − xₘᵢₙ

It is simple but depends only on two observations. The population variance is:

σ² = [∑ᵢ₌₁ᴺ (xᵢ − μ)²] / N

For a sample, the unbiased conventional estimator of variance is:

s² = [∑ᵢ₌₁ⁿ (xᵢ − x̄)²] / (n − 1)

The sample standard deviation is s = √s². The denominator n − 1, called the degrees-of-freedom adjustment, compensates for the fact that x̄ is estimated from the same data. Variance has squared units, while standard deviation has the original units.

The coefficient of variation is:

CV = (s / x̄) × 100%

It is useful for comparing relative variability across processes with different means, provided the mean is positive and meaningfully different from zero.

Percentiles divide ordered data into cumulative proportions. The interquartile range is:

IQR = Q₃ − Q₁

where Q₁ and Q₃ are the first and third quartiles. The IQR is resistant to extreme observations and is often used in box plots.

Shape and graphical analysis

A histogram displays the distribution of continuous measurements. A box plot displays the median, quartiles, and potential outliers. A run chart displays observations in time order and can reveal shifts, trends, cycles, or clustering that are hidden by an unordered summary. A scatterplot examines association between two variables.

Skewness describes asymmetry. In a positively skewed distribution, a long right tail often causes the mean to exceed the median. Kurtosis concerns tail weight and concentration relative to a reference distribution, although it should not be interpreted without graphical and process knowledge.

Worked example: filling volume

Eight measured volumes, in millilitres, are 498, 501, 500, 499, 502, 500, 497, and 503.

The sum is 4,000 ml, so:

x̄ = 4,000 / 8 = 500 ml

The deviations from the mean are −2, 1, 0, −1, 2, 0, −3, and 3 ml. Their squared deviations sum to 28 ml². Therefore:

s² = 28 / (8 − 1) = 4 ml²

and s = √4 = 2 ml. The sample indicates a centre near 500 ml and a typical dispersion of approximately 2 ml. It does not, by itself, establish long-term process stability or compliance with specifications.

3. Probability and Distributions

Probability quantifies uncertainty. A random variable assigns numerical values to outcomes of a chance mechanism. A probability distribution specifies the probabilities associated with those values.

Normal distribution

A normal random variable with mean μ and standard deviation σ has density:

f(x) = [1 / (σ√(2π))] exp[−(x − μ)² / (2σ²)], −∞ < x < ∞

The normal model is continuous, symmetric, and determined by μ and σ, with σ > 0. Approximately 68.27%, 95.45%, and 99.73% of observations lie within 1, 2, and 3 standard deviations of μ, respectively, when the normal model is appropriate.

Standardization uses the z-score:

z = (x − μ) / σ

The z-score is dimensionless and indicates how many standard deviations x lies from the mean. Normality should not be assumed merely because a method is convenient. Histograms, probability plots, subject-matter knowledge, and residual analysis should be considered.

Binomial distribution

The binomial distribution models the number X of successes in n independent trials, each with the same success probability p:

P(X = k) = C(n,k) pᵏ (1 − p)ⁿ⁻ᵏ, k = 0,1,…,n

Its mean is np and variance is np(1 − p). The model is suitable for binary outcomes such as conforming/nonconforming, provided trials are appropriately defined and independent.

Poisson distribution

The Poisson distribution models counts of events in a fixed interval when events occur independently at a stable average rate λ:

P(X = k) = e⁻λ λᵏ / k!, k = 0,1,2,…

Both its mean and variance equal λ. Overdispersion, clustering, or changing exposure can invalidate the simple Poisson model.

Exponential distribution

The exponential distribution is often used for waiting times between Poisson events. With rate λ > 0:

f(t) = λe⁻λᵗ, t ≥ 0

Its mean is 1/λ. It assumes a constant event rate and a memoryless waiting-time property, assumptions that may not hold for equipment failures subject to ageing.

Probability and quality decisions

A probability model is not a statement that nature must follow a formula exactly. It is an approximation used to quantify uncertainty under stated conditions. Model checking is therefore part of responsible quality analysis.

4. Sampling Distributions and the Central Limit Theorem

A sampling distribution is the probability distribution of a statistic over repeated samples of the same size from a population. It is not the distribution of individual measurements.

If observations are independent with population mean μ and finite variance σ², the sample mean has:

E(x̄) = μ

and

Var(x̄) = σ² / n

The standard deviation of the sampling distribution, called the standard error of the mean, is:

SE(x̄) = σ / √n

When σ is unknown, it is commonly estimated by s / √n. Increasing n reduces sampling error at a square-root rate: quadrupling sample size approximately halves the standard error.

The central limit theorem states that, under suitable independence and finite-variance conditions, the distribution of x̄ approaches normality as n increases, even when individual observations are not normal. The theorem does not guarantee independence, remove bias, or justify ignoring severe dependence or extreme distributions.

For a sample proportion p̂ based on n independent Bernoulli trials:

p̂ = X / n

with approximate standard error:

SE(p̂) = √[p(1 − p) / n]

In practice, p is unknown and may be replaced by p̂. Normal approximations are more credible when expected successes and failures are sufficiently large; exact or alternative methods may be preferable for small samples.

5. Estimation and Confidence Intervals

A point estimate is a single numerical estimate of a parameter. An interval estimate provides a range of plausible values together with a confidence level.

Confidence interval for a mean

If the population standard deviation is unknown and observations are independent and approximately normal, a confidence interval for μ is:

x̄ ± t₁₋α/₂, n₋₁ × (s / √n)

Here, α is the significance level, t₁₋α/₂,n−1 is the relevant Student’s t critical value, and n − 1 is the degrees of freedom. The t distribution accounts for additional uncertainty caused by estimating σ with s. For large samples, the t critical value approaches the standard normal critical value.

For a known σ, or in situations where a normal approximation is justified and σ is treated as known, the corresponding interval is:

x̄ ± z₁₋α/₂ × (σ / √n)

Confidence interval for a proportion

A commonly used large-sample Wald interval is:

p̂ ± z₁₋α/₂ √[p̂(1 − p̂) / n]

However, this interval can perform poorly for small samples or proportions near 0 or 1. Wilson or exact intervals are often preferable.

Interpretation

A 95% confidence procedure captures the true parameter in 95% of repeated samples under the model and sampling method. It is not technically correct to say that there is a 95% probability that a fixed parameter lies in one already calculated interval under the frequentist interpretation. Practitioners may informally describe the interval as a plausible range, but the repeated-sampling meaning should be understood.

Worked example: mean cycle time

A random sample of 25 cycle times has x̄ = 42.0 seconds and s = 4.0 seconds. A 95% interval uses t₀.₉₇₅,₂₄ ≈ 2.064:

42.0 ± 2.064 × (4.0 / √25)

The margin of error is 1.651 seconds, giving approximately (40.35, 43.65) seconds. This conclusion assumes representative sampling, independence, and a sufficiently appropriate t-based model.

6. Hypothesis Testing

Hypothesis testing provides a structured method for assessing evidence against a baseline claim.

The null hypothesis, H₀, usually represents a status quo or specified value. The alternative hypothesis, H₁ or Hₐ, represents a departure of interest. A test statistic measures the discrepancy between the sample result and H₀ relative to expected sampling variation.

For a one-sample mean with unknown σ:

t = (x̄ − μ₀) / (s / √n)

where μ₀ is the mean specified by H₀. Under the usual assumptions, this statistic follows a t distribution with n − 1 degrees of freedom when H₀ is true.

A p-value is the probability, assuming H₀ is true and the model is correct, of obtaining a result at least as extreme as the observed result in the direction defined by the alternative. A small p-value indicates incompatibility with H₀; it does not measure the probability that H₀ is true, nor does it measure practical importance.

A decision rule may compare the p-value with α, commonly 0.05:

  • If p ≤ α, reject H₀.
  • If p > α, fail to reject H₀.

“Fail to reject” is not equivalent to proving H₀ true. The data may simply be insufficiently informative.

Errors and power

A Type I error rejects a true H₀. Its long-run probability is α when the test assumptions and decision rule are appropriate. A Type II error fails to reject a false H₀. Its probability is β. Statistical power is:

Power = 1 − β

Power depends on sample size, variability, effect size, and α. Before collecting data, analysts should define the smallest effect that matters operationally and plan sample size accordingly.

One-sided and two-sided tests

A two-sided alternative tests whether a parameter differs in either direction. A one-sided alternative tests a specified direction. Direction must be justified before examining the data; choosing it after observing results inflates misleading claims.

Multiple comparisons

Testing many hypotheses increases the chance of at least one false positive. Corrections such as Bonferroni, Holm, or false discovery rate procedures may be appropriate, depending on the analytical purpose. In SPC, repeated monitoring also creates a sequential decision problem; control-chart rules should be designed and interpreted with this context in mind.

7. Statistical Versus Practical Significance

Statistical significance depends partly on sample size. With a very large sample, a trivial difference may produce a small p-value. With a small sample, an operationally important difference may fail to reach significance.

Practical significance asks whether the magnitude of an effect affects safety, customer requirements, capability, cost, delivery, or risk. Useful measures include:

  • absolute mean difference;
  • percentage difference;
  • confidence interval for the effect;
  • standardized effect size;
  • change in defect probability;
  • expected financial or operational impact.

For two independent means, a standardized effect may be expressed as:

d = (x̄₁ − x̄₂) / sₚ

where sₚ is a pooled standard deviation under the assumption of comparable variances. Standardization supports comparison across scales but does not replace process-specific importance criteria.

A strong quality conclusion therefore reports the estimated effect, confidence interval, assumptions, p-value where relevant, and operational meaning.

8. Assumptions, Data Quality, and Measurement Systems

Independence and autocorrelation

Many elementary formulas assume independent observations. Production data collected sequentially may be autocorrelated: nearby observations can be more similar than distant observations. Ignoring autocorrelation can underestimate uncertainty and create excessive false alarms. Time-series plots, autocorrelation analysis, rational subgrouping, and appropriate time-series methods may be required.

Representative sampling

Convenience samples can exclude certain shifts, machines, operators, products, or environmental conditions. Stratified or systematic sampling may improve representation. Sampling plans should identify the population, sampling frame, selection method, frequency, and handling of missing observations.

Normality and robustness

The t procedure for a mean is reasonably robust to moderate nonnormality in many settings, particularly with adequate sample size and no severe outliers. It is not robust to every departure. Transformations, robust methods, bootstrap intervals, or distribution-specific models may be appropriate.

Measurement error

Observed variation combines process variation and measurement variation. If Y is an observed measurement, a simple model is:

Y = X + E

where X is the true process value and E is measurement error. A gauge repeatability and reproducibility study, calibration, resolution review, and bias assessment should precede important process conclusions. A highly sensitive statistical test applied to an unstable measurement system can produce precise but invalid findings.

Data integrity and ethics

Quality analysts should preserve raw data, document exclusions, distinguish planned from exploratory analyses, and avoid selectively reporting favourable results. Traceability, access control, version history, and appropriate retention are governance requirements in regulated industries. Personal or employee data should be minimized and protected. Statistical evidence must not be manipulated to conceal safety or compliance risks.

9. Applied Case: Reducing Seal-Leak Defects

The following is an illustrative case designed for learning purposes.

Background

A packaging line produces sealed pouches. The historical defect rate is approximately 4%. Engineering introduces a revised sealing-temperature setting and samples 500 pouches after stabilization. Twelve are defective, giving:

p̂ = 12 / 500 = 0.024

The team wants to know whether the observed reduction is sufficiently credible and operationally meaningful.

Problem and analysis

The baseline claim is H₀: p = 0.04. The alternative is H₁: p < 0.04. Using the large-sample approximation under H₀:

SE₀ = √[0.04(1 − 0.04) / 500] ≈ 0.00876

The test statistic is:

z = (0.024 − 0.04) / 0.00876 ≈ −1.83

For a one-sided normal test, this corresponds to a p-value of approximately 0.034. At α = 0.05, the team would reject H₀ under this approximation.

The estimated absolute reduction is 1.6 percentage points, and the relative reduction is 40%. If each defect costs inspection, rework, and potential customer disruption, the effect may be economically important. However, the conclusion requires qualification.

Approach and solution

The quality team verifies that the 500 pouches represent all shifts and relevant materials, confirms that the inspection system detects leaks consistently, and checks whether observations are independent. It also examines the process over time rather than relying on one aggregate proportion. A confidence interval for the new rate should be calculated using a method suitable for proportions; a Wilson interval is preferable to an uncritical Wald interval.

The team then conducts a follow-up run across several production days. If the lower defect rate persists and no adverse effects occur in seal strength, throughput, or energy consumption, the new setting can be incorporated into standardized work and monitored with an appropriate attribute control chart.

Outcome

The initial evidence supports a reduction, but it does not prove permanent improvement. The decision is based on statistical evidence, operational impact, measurement validity, process stability, and risk assessment together.

Learning points

  1. A sample proportion is an estimate, not the exact process rate.
  2. A p-value depends on a model and sampling design.
  3. Practical impact must be quantified separately from statistical significance.
  4. Time-order monitoring is essential after a process change.
  5. A process adjustment should not be standardized until its broader effects are verified.

10. Technology and Contemporary Practice

Statistical software makes calculation, visualization, resampling, and model checking accessible. Automated data collection from sensors and manufacturing execution systems enables high-frequency analysis, but volume does not guarantee validity. Algorithms may detect patterns while failing to distinguish causation from correlation.

AI and machine learning can support anomaly detection, predictive maintenance, and nonlinear modelling. Their use should be controlled through data validation, drift monitoring, explainability appropriate to the risk, human review, and protection against feedback loops. For routine SPC, transparent and auditable methods are often preferable to opaque models when the decision consequences are high.

Contemporary quality practice increasingly combines classical inference with bootstrap methods, Bayesian updating, hierarchical models, mixed-effects models, and time-series analysis. These methods are valuable when data are grouped, sequential, sparse, or otherwise inconsistent with elementary assumptions. The underlying principles remain unchanged: define the process, understand the data-generating mechanism, quantify uncertainty, and connect evidence to decisions.

Best-Practice Checklist

Before interpreting a statistical process analysis, ask:

  1. What population and process period do the data represent?
  2. Was the sampling method representative and documented?
  3. Is the measurement system adequate?
  4. Are observations independent, or is there time dependence?
  5. Are the selected distributional assumptions credible?
  6. Have outliers been investigated rather than automatically removed?
  7. What is the estimated effect and its uncertainty?
  8. Is the effect practically important?
  9. Were multiple comparisons or repeated monitoring considered?
  10. Can another analyst reproduce the analysis from retained records?

Chapter Summary

Statistical process analysis begins with a clear distinction between populations, samples, parameters, and statistics. Descriptive measures such as the mean, median, standard deviation, range, percentiles, and graphical displays summarize process behaviour, but time order and data quality must not be ignored. Probability distributions provide models for continuous measurements, binary outcomes, counts, and waiting times; each model has assumptions and limitations.

Sampling distributions explain why estimates vary from sample to sample. The standard error of a mean decreases with the square root of sample size, and the central limit theorem supports approximate inference under suitable conditions. Confidence intervals quantify estimation uncertainty, while hypothesis tests assess evidence against a defined baseline. Correct interpretation requires attention to Type I and Type II errors, power, p-values, and multiple comparisons.

Statistical significance is not the same as practical significance. Quality decisions should consider effect size, confidence intervals, process risk, cost, customer requirements, and measurement capability. Independence, representative sampling, normality, autocorrelation, and measurement error are central assumptions. Ethical statistical practice requires transparent records, honest reporting, protection of sensitive data, and reproducible analysis. These foundations prepare the analyst to design and interpret control charts and other SPC methods responsibly.

Key Terms

Population: Complete set of units or events of interest.

Sample: Subset of a population that is observed.

Parameter: Numerical characteristic of a population.

Statistic: Numerical characteristic calculated from a sample.

Mean: Arithmetic average of quantitative observations.

Median: Middle value of ordered observations.

Variance: Average squared deviation from a mean, with the specified population or sample denominator.

Standard deviation: Square root of variance, expressed in original units.

Standard error: Standard deviation of a sampling distribution.

Probability distribution: Model describing probabilities of possible values.

Normal distribution: Symmetric continuous distribution determined by mean and standard deviation.

Binomial distribution: Distribution for successes in a fixed number of independent binary trials.

Poisson distribution: Distribution for counts occurring at a stable average rate.

Central limit theorem: Result describing approximate normality of sample means under suitable conditions.

Confidence interval: Interval produced by a procedure with a stated long-run coverage probability.

Null hypothesis: Baseline claim evaluated in a hypothesis test.

p-value: Probability, under the null model, of a result at least as extreme as observed.

Type I error: Rejection of a true null hypothesis.

Statistical power: Probability of rejecting a false null hypothesis.

Practical significance: Operational, economic, technical, or safety importance of an effect.

Review Questions

Short Answer Questions

  1. Distinguish a parameter from a statistic and give one example of each.
  2. Why does the sample variance use n − 1 rather than n in its denominator?
  3. What does the standard error of the mean represent?
  4. State two assumptions required for a binomial model.
  5. Explain why failing to reject H₀ does not prove that H₀ is true.

Long Answer Questions

  1. Explain how descriptive statistics and graphical analysis should be combined when investigating process variation.
  2. Compare the normal, binomial, and Poisson distributions, including their assumptions and quality applications.
  3. Derive the logic of a confidence interval for a population mean when σ is unknown.
  4. Explain Type I error, Type II error, power, sample size, and practical significance in process improvement.
  5. Discuss how measurement error, autocorrelation, and nonrepresentative sampling can invalidate statistical conclusions.

Analytical / Application Questions

  1. A process sample of 16 observations has x̄ = 75 units and s = 8 units. Explain which confidence-interval method is appropriate and identify the assumptions that should be checked before calculating the interval.
  2. A manager reports p = 0.01 after analysing 20,000 units and claims that the improvement is important. What additional information should a quality professional request before accepting the claim?
  3. A defect-count series shows increasing variability and clustering over time. Evaluate whether a simple Poisson model is appropriate and recommend a statistical investigation plan.

Multiple Choice Questions with Answers

Question 1

Which quantity is a population parameter?

A. The sample mean x̄

B. The sample standard deviation s

C. The population mean μ

D. The observed sample proportion p̂

Correct Answer: C. The population mean μ

Question 2

For a sample with n observations, the conventional unbiased sample variance uses which denominator?

A. n − 1

B. n + 1

C. n²

D. 2n

Correct Answer: A. n − 1

Question 3

What does the standard error of the sample mean primarily measure?

A. The spread of individual observations

B. The variability of the mean across repeated samples

C. The probability that the null hypothesis is true

D. The measurement resolution of an instrument

Correct Answer: B. The variability of the mean across repeated samples

Question 4

Which distribution is most directly suited to modelling the number of defects in a fixed opportunity interval under a stable event rate?

A. Normal

B. Binomial

C. Uniform

D. Poisson

Correct Answer: D. Poisson

Question 5

If sample size is increased by a factor of four, the standard error of a mean approximately:

A. Doubles

B. Is halved

C. Is quartered

D. Does not change

Correct Answer: B. Is halved

Question 6

A p-value is best defined as the probability of:

A. The null hypothesis being true

B. Obtaining a result at least as extreme as observed, assuming the null model is true

C. Making a Type II error

D. The practical effect being important

Correct Answer: B. Obtaining a result at least as extreme as observed, assuming the null model is true

Question 7

Which statement about a 95% frequentist confidence interval is correct?

A. There is a 95% probability that the fixed parameter lies in this particular interval

B. 95% of individual observations must lie inside the interval

C. The procedure captures the true parameter in 95% of repeated samples under its assumptions

D. The interval proves that the estimate is unbiased

Correct Answer: C. The procedure captures the true parameter in 95% of repeated samples under its assumptions

Question 8

A statistically significant result may still lack practical significance when:

A. The sample size is large and the effect is operationally trivial

B. The p-value is greater than one

C. The measurement uses original units

D. The sample is randomly selected

Correct Answer: A. The sample size is large and the effect is operationally trivial

Question 9

Which issue can make ordinary statistical tests appear more certain than they should be when production data are collected sequentially?

A. Using a median

B. Autocorrelation among observations

C. Reporting units

D. Using a histogram

Correct Answer: B. Autocorrelation among observations

Question 10

If a test has Type II error probability β, its statistical power is:

A. α + β

B. αβ

C. 1 − β

D. 1 − α − β

Correct Answer: C. 1 − β

Statistical Process Control
Chapter 2 of 15