UpSkill Sprint Consulting logoUpSkill Sprint Consulting
UpSkillSprint · Statistics · Beginner

Resampling in Minitab

Learn bootstrap confidence intervals and randomization tests with a ready-to-run practice dataset, the original 120 manufacturing records, or data you paste from your own work.

5 Minitab procedures5 live simulation labs120-row datasetNo formulas required to begin

1. Begin with the big idea

Resampling means repeatedly rebuilding samples with a computer. We use the patterns in those rebuilt samples to measure uncertainty or test a claim when a textbook distribution is not our best starting point.

Let us slow that sentence down

You normally have one real sample: perhaps 120 tensile-strength results, 40 customer wait times, or 25 fill weights from a packaging line. That sample gives one mean, one median, or one pass rate. The uncomfortable question is: if we collected another reasonable sample from the same process, how much might the answer move?

Collecting thousands of new physical samples would be slow, expensive, and—if your test destroys the specimen—a little hard on the budget. Resampling lets the computer create many rehearsal samples from the evidence already available. It calculates the statistic for each rehearsal, then studies how those answers vary.

Think of it as letting your dataset audition for the same role thousands of times. Some rows appear repeatedly; some miss a particular audition. The computer is fast, does not ask for overtime, and never complains about running sample number 9,999.

A factory example: “Is 51.5 MPa a stable answer?”

Imagine that these eight measurements are the only specimens tested today. Their mean is 51.5 MPa. That number is useful, but it is not carved into a steel billet. A different set of eight specimens would probably give a slightly different mean.

Observed sample — the measurements we actually paid to collect
4850505152525455
Observed mean = 51.5 MPa. These eight values are the evidence; resampling does not secretly manufacture a ninth measurement.
Step 1

Start with the observed values.

Step 2

Draw eight values with replacement.

Step 3

Calculate the new sample mean.

Step 4

Repeat thousands of times and compare the means.

One computer-built bootstrap sample
Select “Build another resample” to watch the idea.
Amber values are duplicates in this resample. Duplicates are expected because each selected value is returned before the next draw.

What do the repeated answers tell us?

How uncertain is our estimate?

If most resampled means stay close to 51.5 MPa, the estimated mean is fairly stable. If they spread widely, our estimate is less precise. Bootstrap confidence intervals use that spread.

Could chance explain our result?

A randomization test builds results under a stated null claim. If our real result sits deep in the simulated tail, ordinary chance under that claim becomes a weak explanation.

Why not always use a textbook curve?

Classical formulas are excellent when their assumptions fit. Resampling is especially helpful for awkward statistics, questionable normality, or when seeing the simulated logic makes the uncertainty easier to understand.

Important reality check: “Computer-generated sample” does not mean “new evidence.” If every thermometer was miscalibrated by 3°C, resampling will reproduce the mistake thousands of times—with impressive speed and absolutely no shame. The original sample still must be representative, measured correctly, and collected with the right sampling unit.

Three situations you can picture

Steel mill

You tested one coil per heat and want uncertainty in the average Charpy energy. Resampling individual coil rows may be wrong if coils from the same heat are dependent; the heat may need to be the sampling unit.

Coffee shop

You record 30 customer wait times. A few “large latte, six modifications” orders create a skewed distribution. Bootstrapping the median can estimate a typical wait without pretending the times form a perfect bell curve.

New process trial

You compare a Standard and New process. Shuffling group labels asks what mean differences would occur if the labels did not truly matter—like swapping team jerseys and checking whether the scoreboard still looks unusually one-sided.

Observed sample

The 120 rows you actually measured. It is your evidence and the source of every resample.

One resample

A computer-generated sample made according to a rule. Bootstrap draws rows with replacement; randomization imposes the null hypothesis.

Resampling distribution

The statistic from thousands of resamples. Its spread, percentiles, or extreme tails answer the question.

Bootstrap and randomization are cousins, not twins

Bootstrap → estimation
Ask: “What values are plausible for the population mean or mean difference?” Output: confidence interval.
Randomization → hypothesis testing
Ask: “If the null claim were true, how often would a result this extreme occur?” Output: p-value.

Data → apply resampling rule → calculate statistic → repeat many times → interpret distribution

Questions & Answers — 21 questions Open this section to explore common questions about bootstrapping and randomization tests.
What does “with replacement” mean?

After a row is drawn, it goes back into the pool before the next draw. A row can appear zero, one, or several times. The resample stays the same size as the original sample.

Does resampling create new real-world information?

No. It reveals what the observed data imply under a specific resampling rule. A biased, dependent, or unrepresentative sample remains problematic.

Does a confidence interval assume that the underlying distribution is normal?

Not always. It depends on what confidence interval you are calculating, the sample size, and the method used.

SituationDoes the raw data need to be normal?Why?
Small sample, for example n < 30Approximately yes for a classical t-based intervalThe method relies more heavily on the population being approximately normal.
Large sample, for example n ≥ 30Usually noThe Central Limit Theorem often makes the sampling distribution of the mean approximately normal.
Very skewed data or major outliersBe carefulEven n = 30 may not be enough.
Bootstrap confidence intervalNo normality assumption is requiredThe sampling distribution is estimated by resampling the observed data.

The key distinction is that a standard confidence interval for the mean does not necessarily require the raw data to be normally distributed. What matters is whether the sampling distribution of the mean is approximately normal.

For example, individual Charpy values could be right-skewed:

80, 95, 110, 115, 125, 140, 160, 220, 310, …

If you repeatedly took samples of 50 specimens and calculated the mean of each sample, those means would tend to form a more nearly normal distribution because of the Central Limit Theorem.

Keep the two distributions separate:
Population or data distribution → can be non-normal.
Sampling distribution of x̄ → becomes approximately normal as n increases.
Important: “n ≥ 30” is a rule of thumb, not a law. Severe skewness or outliers may require a much larger sample or a bootstrap interval.
What is the difference between randomization and bootstrapping?
Randomization asks: “Could this observed difference have happened by chance if there were really no effect?”

Bootstrapping asks: “How uncertain is my estimate?”

Suppose two rolling practices have an observed mean Charpy difference of about 31 J.

Randomization or permutation test: start with the null hypothesis that the two groups are equivalent. Combine the observations, shuffle the A/B labels, split them back into groups of the original sizes, calculate the mean difference, and repeat thousands of times. This creates a null distribution.

Bootstrapping: preserve the groups and repeatedly sample with replacement from each group. Calculate the statistic each time. This creates a bootstrap distribution that describes the precision and uncertainty of the estimate.

FeatureRandomizationBootstrapping
Main purposeHypothesis testingEstimation and uncertainty
Typical outputp-valueConfidence interval or standard error
Main questionCould this happen by chance under H₀?How variable is my estimate?
SamplingUsually without replacementWith replacement
Group labelsUsually shuffledUsually preserved
CreatesNull distributionSampling distribution
Memory trick: Randomization = shuffle; bootstrap = resample with replacement. Randomization asks about evidence of an effect. Bootstrap asks how precisely the effect was estimated.
What does “randomization rearranges labels without replacement” mean?

Suppose the observed data are:

Practice APractice B
200150
210160
220170

The observed difference in means is 50 J. The randomization test asks: if A and B really had no effect, could a 50 J difference arise just by randomly assigning these observations to two groups?

First remove the labels:

200, 210, 220, 150, 160, 170

Then randomly assign three observations to A and three to B. Every original observation appears exactly once in that randomization.

Valid: A = [200, 160, 170]    B = [210, 220, 150]

Not valid: A = [200, 200, 170]

The second assignment is invalid because the original dataset contained only one 200 J observation. Without replacement means that once an observation is assigned during one randomization, it cannot be selected again in that same randomization. Bootstrap sampling is different: repeated values are allowed.

How do we determine how often random shuffling produces a difference as large as the observed 31 J?

Begin with the observed difference:

Dobserved = 31 J

Assume the null hypothesis is true:

H₀: μA = μB

Shuffle the group labels thousands of times. After each shuffle, calculate:

D* = x̄A* − x̄B*

For a two-sided test, count randomized differences at least as extreme as the observed result in either direction:

|D*| ≥ |Dobserved|

Suppose 10,000 randomizations produce:

ResultNumber of randomizations
−31 < D* < +319,720
D* ≥ +31145
D* ≤ −31135
Total10,000

Extreme results = 145 + 135 = 280

p = 280 / 10,000 = 0.028

If A and B truly had no difference, about 2.8% of random label assignments produced a difference at least as extreme as the 31 J result observed.

Actual data → observed difference → assume H₀ → shuffle labels → calculate D* → count extreme values → p-value
What about a one-sample randomization test?

A one-sample randomization test follows the same broad idea, but there are no A/B group labels to shuffle. Suppose:

H₀: μ = 180 J

Sample: 185, 190, 175, 200, 195

The sample mean is 189 J, so the observed difference from the null value is 9 J. One possible one-sample method is a sign-randomization test. First calculate the deviations from 180:

+5, +10, −5, +20, +15

Keep the magnitudes but randomly flip their signs. For example:

+5, −10, +5, −20, +15

This randomized set has a mean deviation of −1 J, corresponding to a randomized sample mean of 179 J. Repeating the sign randomization creates a null distribution centered around the hypothesized value.

Technical caveat: sign-flipping is not a completely assumption-free test of the population mean. It relies on an appropriate symmetry or exchangeability assumption. Minitab’s one-sample randomization procedure uses the null-centering method explained below instead.
If we use the same one sample without replacement, wouldn’t the mean stay the same?

Yes. That objection is correct.

If we simply reorder:

185, 190, 175, 200, 195

into:

200, 175, 195, 185, 190

the mean does not change. Reordering one sample cannot create a useful null distribution. A one-sample randomization procedure must randomize something other than the observation order.

In a sign-randomization approach, the signs attached to the deviation magnitudes change. For example, the magnitudes 5, 10, 5, 20, and 15 could become:

+5, −10, +5, −20, +15

−5, +10, −5, +20, +15

Those different sign patterns produce different means.

Key insight: If you literally shuffle a single sample without replacement, nothing happens to its mean.
Is sign-flipping the method Minitab uses for a 1-sample randomization test?

No. Minitab uses a different approach.

Sample: 185, 190, 175, 200, 195    x̄ = 189    H₀: μ = 180

Minitab first shifts every observation by the same amount so the adjusted sample mean equals the hypothesized mean. The required shift is:

180 − 189 = −9

The null-adjusted data become:

176, 181, 166, 191, 186    adjusted mean = 180

Minitab then repeatedly resamples with replacement from this null-adjusted sample, keeping the original sample size. Each resample produces a new mean, building a null distribution centered around 180.

176, 176, 191, 181, 191 → mean = 183

166, 181, 166, 186, 176 → mean = 175

Bootstrap confidence intervalMinitab 1-sample randomization
Starts with the original dataStarts with the original data
Keeps the center at x̄Shifts the data so the center equals μ₀
Samples with replacementSamples with replacement
Builds a distribution around x̄Builds a distribution around μ₀
Estimates uncertainty or a confidence intervalTests H₀ and calculates a p-value

In short, Minitab uses bootstrap-style resampling after first forcing the data to represent the null hypothesis.

How does Minitab perform a 2-sample randomization test?

Minitab’s 2-sample randomization test follows the classic permutation logic. Suppose:

Practice APractice B
210175
195182
220168
205190
215177

x̄A = 209    x̄B = 178.4    Dobserved = 30.6 J

H₀: μA − μB = 0

Minitab pools all observations, randomly reallocates them into two groups while preserving the original group sizes, and calculates a new difference in means.

Randomized ARandomized B
210195
175220
182205
190215
177168

D* = x̄A* − x̄B*

Repeating this process thousands of times creates a distribution centered near zero, representing differences expected under the null hypothesis. For a two-sided test, Minitab counts the proportion satisfying:

|D*| ≥ 30.6

That proportion is the randomization p-value.

If my original bootstrap sample size is 6 and I resample 1,000 times, does every resample also contain 6 observations?

Yes. In ordinary nonparametric bootstrapping, if the original sample size is n = 6, every bootstrap resample also contains six observations. Those observations are selected with replacement.

Original sample: [100, 110, 120, 130, 140, 150]

SampleSix selected observationsMean
Original100, 110, 120, 130, 140, 150125.0
Bootstrap #1100, 100, 120, 140, 140, 150125.0
Bootstrap #2110, 110, 110, 130, 150, 150126.7
Bootstrap #3100, 120, 120, 120, 130, 140121.7

Every resample contains exactly six selections, but some original observations may appear several times while others may not appear at all. Do not confuse these numbers:

n = 6 → observations inside each bootstrap resample

B = 1,000 → number of separate bootstrap resamples

You create 1,000 separate samples of six observations each and usually calculate one statistic from each sample. Keeping the resample size at six approximates how the statistic would vary across repeated studies having the same sample size as the original study.

Original sample n = 6 → draw 6 values with replacement → calculate one statistic → repeat 1,000 times → 1,000 bootstrap statistics
What resampling methods should I practice in Minitab, and what questions do they answer?

The practice set was organized around five Minitab resampling methods:

MethodPractice Variable(s)Core Question
Bootstrapping for 1-Sample FunctionStrength_MPaWhat is the 95% confidence interval for the average strength?
Bootstrapping for 2-Sample MeansStrength_MPa + GroupIs the New process mean higher than Standard?
Randomization Test for 1-Sample MeanStrength_MPaIs the average strength different from 51 MPa?
Randomization Test for 1-Sample ProportionPassFailIs the pass rate greater than 90%?
Randomization Test for 2-Sample MeansStrength_MPa + GroupIs the New process statistically different from Standard?

Use 1,000 or 5,000 resamples to start. Higher numbers of resamples give smoother resampling distributions.

What does each resampling method mean, what is it used for, and what is a practical example?

These methods in Minitab fall under resampling statistics. They are powerful because they do not rely as heavily on normal-distribution assumptions as many traditional t-tests and z-tests. They use repeated random sampling or simulation to estimate uncertainty, significance, or confidence intervals.

1. Bootstrapping for 1-Sample Function

Bootstrapping repeatedly:

  1. Randomly samples your data.
  2. Samples with replacement.
  3. Calculates a statistic each time.
  4. Builds a distribution of that statistic.

“Function” simply means the statistic you care about: mean, median, standard deviation, percentile, capability metric, and so on.

It is used when you want confidence intervals, robust estimates, small-sample analysis, or non-normal-data analysis.

Real steel question: “What is the likely range of average Charpy toughness for this production run?”

For example, instead of assuming normality for 25 Charpy tests, bootstrap resampling takes repeated samples from those 25 observations, estimates the mean each time, and creates an empirical confidence interval.

2. Bootstrapping for 2-Sample Means

This uses the same concept but compares the mean of Group A with the mean of Group B. Both groups are repeatedly resampled.

It is used to estimate the difference between process means, the confidence interval for that difference, and the stability of the difference.

Example: “Did the new cooling practice improve yield strength?”

3. Randomization Test for 1-Sample Mean

This is a simulation-based hypothesis test. It asks whether a sample mean is significantly different from a target without relying on a theoretical t-distribution in the same way as the classical test.

Example: Target strength = 51 MPa; observed average = 53 MPa. “Could this difference happen by random chance?”

4. Randomization Test for 1-Sample Proportion

This tests whether a proportion differs from a target proportion. A proportion can represent pass/fail, good/bad, defect/no defect, yield, scrap rate, or acceptance rate.

Example: 93 of 100 pipes pass. Target = 90%. “Is the process truly better than 90%, or is this random variation?”

5. Randomization Test for 2-Sample Means

This tests whether two groups truly differ. It repeatedly shuffles group labels, recalculates the difference, and builds a null distribution.

Example: Compare coil exit temperatures for good coils versus failed coils. “Are the means statistically different?”

Key Difference: Bootstrap vs. Randomization

MethodMain PurposeMental ModelTypical Output
BootstrappingEstimate uncertainty / confidence intervals“What values are plausible?”Confidence intervals; estimate variability
Randomization TestTest statistical significance“Could this happen by chance?”p-values; significance evidence
How does bootstrapping for a 1-sample function work?

Bootstrapping answers this question:

“If I repeated this experiment many times, what range of answers would I likely get?”

But instead of collecting new real-world samples, we repeatedly resample from the data we already have.

The Big Concept

Suppose you collected only five yield-strength observations:

CoilYield Strength
C171
C269
C374
C473
C568

Bootstrapping says:

“Let’s simulate thousands of possible samples using the observed data itself.”

What “With Replacement” Means

This is the most important concept. From the original data:

71, 69, 74, 73, 68

a bootstrap sample might be:

71, 71, 74, 68, 73

Another might be:

69, 68, 68, 74, 71

Some values repeat and some disappear because sampling is with replacement. After selecting a value, it goes back into the pool.

Why This Works

Your original sample is treated as the best estimate of the real population. Bootstrapping simulates thousands of alternate realities based on your observed process variation.

What Is a “Function”?

FunctionMeaning
MeanAverage
MedianMiddle value
Standard deviationVariation
PercentileDistribution location
MaximumWorst case
CpkCapability
CorrelationRelationship strength

What Bootstrapping Produces

Usually, a confidence interval. For example:

StatisticResult
Mean YS71.2 ksi
95% CI69.8–72.5 ksi

Meaning:

“We are 95% confident the true process mean lies in this range.”

Practical Steel Example

Suppose coiling-temperature measurements are 610, 618, 605, 622, and 611. The question is:

“What is the likely true average coiling temperature?”

Bootstrap:

  1. Resamples the data thousands of times.
  2. Computes the mean each time.
  3. Builds a distribution of means.
  4. Calculates a confidence interval.

Step-by-Step in Minitab

Step 1 — Open the dataset. Use Strength_MPa.
Step 2 — Open the tool. Calc → Resampling → Bootstrapping for 1-Sample Function.
Step 3 — Fill the dialog. Samples = Strength_MPa; Function = Mean; Number of resamples = 1,000 or 5,000.
Step 4 — Confidence level. Choose 95%.
Step 5 — Run. Click OK.
Important insight: Bootstrap does not magically create more information. It estimates the uncertainty in your estimate.
Common mistake: “Bootstrapping creates fake data.” Not exactly. It creates simulated samples from the observed empirical distribution.
How does bootstrapping for 2-sample means work?
This directly addresses: “Are these two processes actually different?”

Instead of estimating one mean, we estimate the difference between two means.

Bootstrap repeatedly:

  1. Resamples Group A.
  2. Resamples Group B.
  3. Calculates both means.
  4. Computes their difference.
  5. Repeats thousands of times.

This builds a distribution of mean differences.

Example

GroupMean Strength
Standard51.2
New53.8
Observed difference = 53.8 − 51.2 = +2.6 MPa

But is that real, or random variation? Bootstrap helps estimate uncertainty around that difference.

What Bootstrap Does Internally

Standard group: original 51, 52, 50, 53, 49 → possible bootstrap sample 51, 51, 53, 49, 50 → mean = 51.0.

New group: original 54, 55, 53, 52, 56 → possible bootstrap sample 55, 55, 53, 52, 54 → mean = 53.8.

Δ* = x̄New* − x̄Standard* = 53.8 − 51.0 = +2.8 MPa

Repeat thousands of times. Values such as 2.8, 2.1, 3.4, 1.9, and 2.6 form the bootstrap distribution of mean differences.

Key Interpretation Rule

Confidence IntervalInterpretation
Does not cross 0A meaningful difference is supported by the interval.
Crosses 0The uncertainty includes “no difference.”
Fully positiveNew process likely improves strength.
Fully negativeNew process likely worsened strength.

Why zero? Because:

Δ = μNew − μStandard = 0

means there is no difference between the group means.

Step-by-Step in Minitab

Step 1. Use Strength_MPa and Group.
Step 2. Calc → Resampling → Bootstrapping for 2-Sample Means.
Step 3. Samples = Strength_MPa; Corresponding groups = Group; First sample = Standard; Second sample = New.
Step 4. Use 5,000 resamples.
Step 5. Use a 95% confidence level.
Step 6. Click OK.
Important insight: Bootstrap gives magnitude + uncertainty, not just “significant/not significant.” Practical significance matters too.
How does a randomization test for a 1-sample mean work?

Now we move from estimating uncertainty with bootstrap to testing significance with randomization testing. This is a major conceptual shift.

“Could my observed sample mean reasonably occur by random chance if the true mean were actually the target value?”

This is fundamentally a hypothesis test, but performed using simulation instead of relying only on theoretical distributions.

Engineering Example

Suppose target yield strength = 51 MPa and the sample is 53, 52, 54, 51, 55. The observed mean is 53 MPa.

H0: μ = 51 MPa
H1: μ ≠ 51 MPa

Now the question becomes:

“Is 53 truly different from 51?” — or — “Could this difference happen just from random variation?”

Core Logic

Assume H0 is true → generate simulated samples under H0 → compute a mean each time → build a null distribution → compare the actual mean

The randomization distribution represents what random variation alone would produce if the null hypothesis were true.

The p-value

The p-value answers how unusual the observed result is under H0. A small p-value means the observed result is rare in the simulated null world; a large p-value means random variation could readily explain it.

Important engineering interpretation: The procedure does not prove H0 or H1. It evaluates whether the observed evidence is surprising under H0. Statistical significance is also not the same as practical significance.

Step-by-Step in Minitab

Step 1. Use Strength_MPa.
Step 2. Calc → Resampling → Randomization Test for 1-Sample Mean.
Step 3. Sample = Strength_MPa; Test mean = 51.
Step 4. Choose the alternative: ≠ for different, > for greater than, or < for less than. For practice, use ≠.
Step 5. Use 5,000 resamples.
Step 6. Click OK.
If we already have a target mean, why do we still need to resample?

The key is this:

We are NOT resampling to estimate the target. We are resampling to estimate what random variation would look like IF the target were true.

Suppose the target mean is:

μ0 = 51

and the observed sample mean is:

x̄ = 53

Now the question becomes:

“Is 53 unusually far from 51?” — or — “Could random chance naturally produce 53 even if the true mean were actually 51?”

Why We Cannot Just Compare 53 vs. 51 Directly

Because samples naturally vary. Even if the true process mean is 51, one sample may average 49, another 52, another 50.5, and another 53 because of random process and sampling variation.

So we need a way to determine how much variation is “normal.”

The Dartboard Analogy

Imagine a dartboard whose target center is 51. Darts naturally scatter around the center. If one dart lands at 53, is that normal scatter or evidence that the aim has shifted?

You need to know how wide the natural scatter normally is. That is what resampling estimates.

What Randomization Builds

It builds the distribution of expected random variation under the null hypothesis — the null distribution.

Suppose Minitab simulates 5,000 samples under:

H0: μ = 51

and generates means such as 50.8, 51.2, 50.5, 51.7, 49.9. Most are near 51 because under H0, 51 is assumed true.

Then compare the real mean, 53, and ask:

“How often do simulated means get as far as 53?”
OutcomeInterpretation
53 is rare in the null distribution53 is unusually far from 51 → evidence against H0.
53 is common in the null distributionNormal variation could readily explain it → insufficient evidence against H0.
Simplest possible mental model:
Target mean → “What are we testing against?”
Resampling → “How much random variation should we expect around that target?”

Randomization testing is basically:

“Let’s pretend the target is true and simulate reality thousands of times. Does our observed result fit naturally into that simulated reality?”
Are we resampling the original dataset, or are we simply simulating random samples?

We are resampling the original dataset, but in a way that simulates what random samples would look like under the null hypothesis. That subtle distinction is the whole trick.

MethodWhat Gets Resampled?
BootstrapOriginal dataset directly
Randomization test for 1-sample meanModified / recentered version of the original data

Bootstrap

“Assume the sample is reality.”

So bootstrap directly resamples the original data with replacement.

Randomization Test

“Assume the null hypothesis is reality.”

Before resampling, the data is shifted/recentered to align with the target mean. Then resampling/simulation occurs.

For example, suppose:

Data = 53, 52, 54, 51, 55    and    x̄ = 53
H0: μ = 51

The data are centered at 53, but under H0 they should be centered at 51. So Minitab first transforms the data.

OriginalRecentered
5351
5250
5452
5149
5553

Then it repeatedly samples from the recentered data, computes means, and builds the null distribution.

Bootstrap: Observed Data → Resample Directly → Estimate uncertainty
1-Sample Randomization: Observed Data → Recenter to H0 → Resample → Test significance
Extremely important insight: The randomization distribution represents what we would expect to see purely from random chance if the null hypothesis were true.
What is the purpose of the recentering step in a 1-sample randomization test?

The recentering step exists because:

We must force the data to behave as though the null hypothesis is true.

Without recentering, the simulation would still reflect your observed sample mean instead of the hypothesized mean. That would break the logic of the test.

The Core Goal

“If the null hypothesis were actually true, what kinds of sample means would naturally occur from random variation alone?”

Suppose:

x̄ = 53    but    H0: μ = μ0 = 51

If we directly resample the original data, simulated samples will naturally center around 53, not 51. The simulation would represent a world where the true mean is 53, which contradicts the null hypothesis.

Then the test would effectively ask:

“Could 53 happen if the true mean were 53?”

which is meaningless. Of course 53 is common if reality is centered at 53.

The Recentering Equation

xi* = xi − x̄ + μ0

where:

  • xi = original observation
  • x̄ = observed sample mean
  • μ0 = hypothesized mean
  • xi* = recentered observation
OriginalRecentered
5351
5250
5452
5149
5553

What Recentering Preserves

Recentering changes the center, but preserves the spread, variability, shape, skewness, outliers, and process behavior.

Visually, it simply slides the whole distribution from a center of 53 to a center of 51 while leaving the shape the same.

Recentering exists to force the simulated data to obey the null hypothesis while preserving the original process variation.
How does a randomization test for 2-sample means work?
“Could the observed difference between two group means happen purely by random chance if the groups were actually equivalent?”

This is a significance test for the difference between two means.

Engineering Example

GroupMean Yield Strength
Standard MP51 MPa
New MP54 MPa
Δobserved = x̄New − x̄Standard = 54 − 51 = 3 MPa

The question is:

“Is this 3 MPa difference real?” — or — “Could random process variation easily create this difference?”

The Null Hypothesis

H0: μ1 = μ2
equivalently, H0: μ1 − μ2 = 0

Under H0, group labels should not matter. That is the entire philosophy of the test.

Why This Differs From 1-Sample Randomization

A 1-sample randomization test needed recentering because there was a target mean. A 2-sample randomization test usually does not need that same recentering step. Under H0, both groups are treated as coming from the same population, so instead we shuffle the group labels.

How the Randomization Works

Original data:

ValueGroup
51Standard
50Standard
52Standard
55New
54New
53New
Step 1 — Pool all data. Combine 51, 50, 52, 55, 54, 53 and temporarily ignore the original labels.
Step 2 — Randomly reassign labels. Keep the original group sizes, but randomly assign observations to Standard or New.
Step 3 — Compute the randomized difference. Calculate the two randomized group means and their difference.
Step 4 — Repeat thousands of times. Values such as 0.8, −1.2, 1.5, 0.3, −0.6 form the null distribution.

The null distribution represents:

“What kind of differences occur naturally if the groups are truly equivalent?”

Then compare the observed difference, +3 MPa, and ask:

“How often does random label shuffling create a difference this extreme?”

Why Label Shuffling Makes Sense

If H0 is true, group membership is meaningless. Therefore, randomly reassigning labels should produce equally plausible datasets.

If randomized differences resemble the observed difference, the observed difference is not unusual. If the observed difference is much larger than what randomization usually produces, that is evidence the groups differ.

Step-by-Step in Minitab

Step 1. Use Strength_MPa and Group.
Step 2. Calc → Resampling → Randomization Test for 2-Sample Means.
Step 3. Samples = Strength_MPa; Corresponding groups = Group; First sample = Standard; Second sample = New.
Step 4. Usually choose ≠ for a two-sided test.
Step 5. Use 5,000 resamples.
Step 6. Click OK.
Engineering insight: Statistical significance ≠ practical significance. A tiny difference can be statistically significant but operationally meaningless. Always evaluate engineering impact.
How does a randomization test for a 1-sample proportion work?

Now we move from means to proportions. This is extremely important in quality engineering because many industrial problems are naturally pass/fail, defect/no defect, cracked/not cracked, downgraded/not downgraded, or accepted/rejected.

“Could the observed proportion reasonably occur by random chance if the true process proportion were the target value?”

Quality Example

Suppose historical pass rate = 90%, so the target proportion is:

p0 = 0.90

Now suppose 100 pipes are tested and only 83 pass:

p̂ = 83/100 = 0.83

The question is:

“Is 83% significantly worse than 90%?” — or — “Could normal random variation explain this?”

Hypotheses

H0: p = 0.90

For a two-sided question:

H1: p ≠ 0.90

For a deterioration question:

H1: p < 0.90

The Key Concept

Under H0, the process still has a 90% pass probability. The randomization procedure simulates what pass/fail outcomes would naturally look like if that were true.

With proportions, the outcomes are binary. If:

p0 = 0.90

then each simulated observation has a 90% probability of Pass and a 10% probability of Fail.

For 100 pipes, one simulation might produce 91 Pass / 9 Fail; another 88 / 12; another 93 / 7. These results build the null distribution of pass rates.

Then compare the actual 83% and ask:

“How often would a true 90% process randomly produce 83% or worse?”
ResultMeaning
Small p-value, e.g. 0.002Only a very small proportion of null simulations are this poor → evidence the process degraded.
Large p-value, e.g. 0.28The observed result is not especially unusual under the target process.

Step-by-Step in Minitab

Step 1. Use PassFail, containing Pass and Fail.
Step 2. Calc → Resampling → Randomization Test for 1-Sample Proportion.
Step 3. Sample = PassFail; Event = Pass; Test proportion = 0.90.
Step 4. For a degradation study, use < to test whether the pass rate fell below target.
Step 5. Use 5,000 resamples.
Step 6. Click OK.
Final engineering insight: A statistically significant defect increase may still be operationally manageable, while a non-significant increase may still require monitoring. Statistics supports engineering judgment — it does not replace it.

2. Choose the correct method

Select what you want to do, the statistic, and the number of groups.

Choose one option on each row.
Method 1 of 5

Bootstrap for a 1-sample function

Estimation

First: what is this method?

A 1-sample bootstrap estimates the uncertainty around one statistic calculated from one sample—for example a mean, median, or another supported function. It answers: “If this sample represents the process, how much could my statistic change from sample to sample?”

How it works

  1. Calculate the statistic from the real sample.
  2. Draw a new sample of the same size with replacement.
  3. Calculate the same statistic from that resample.
  4. Repeat many times to form a bootstrap distribution.
  5. Use its spread for the standard error and its percentiles for a confidence interval.

When to use it

  • You have one representative sample and want an interval estimate.
  • The statistic has no simple textbook standard-error formula, such as a median.
  • The data may be skewed or normality is doubtful.
  • Do not resample individual rows when they are dependent; resample the independent unit instead.

Core equations

Observed statistic: T̂ = t(x₁, x₂, …, xₙ)
Bootstrap statistic: T*ᵦ = t(x*₁ᵦ, x*₂ᵦ, …, x*ₙᵦ)
SE_boot = √[ Σᵦ(T*ᵦ − T̄*)² / (B − 1) ]
Percentile CI = [ Qα/2(T*), Q1−α/2(T*) ]
t(·) = chosen function, such as meann = original sample sizeB = number of resamplesQ = percentile of bootstrap statisticsα = 1 − confidence level* = value created by resampling

Relatable example

A lab tests 30 tensile specimens and obtains an average strength. Testing another 5,000 physical sets would destroy a heroic amount of steel. The bootstrap reuses the 30 observed results in 5,000 computer-built samples to estimate how stable the mean is. It saves specimens—not standards of evidence.

Now apply it to the lesson dataset

Practice question: What is the 95% bootstrap CI for overall mean Strength_MPa?

Each bootstrap sample draws 120 strength values with replacement. Minitab calculates the mean. Repeating this creates a distribution of plausible sample means.

Calc›Resampling›Bootstrapping for 1-Sample Function
  1. Choose One or more samples, each in a column; enter Strength_MPa.
  2. Choose Mean as the function. Under Options, select the confidence level and number of bootstrap samples.
  3. Use at least 5,000 resamples for practice; a larger count makes simulation noise smaller.
  4. Click OK. Read the observed statistic, standard error, and confidence interval—not just the graph.
Run the simulation to create and interpret the interval.
Method 2 of 5

Bootstrap for 2-sample means

Estimation

First: what is this method?

A 2-sample bootstrap estimates uncertainty in the difference between two independent population means. It answers: “How large might the true New − Standard mean difference be, given the sampling variation in both groups?”

How it works

  1. Calculate the observed difference in group means.
  2. Resample with replacement within New and separately within Standard.
  3. Calculate New resample mean − Standard resample mean.
  4. Repeat many times and take percentiles of the differences.

The groups remain separate because we are estimating their difference, not pretending they are one population.

When to use it

  • You want a confidence interval for the difference between two independent group means.
  • Normal-theory assumptions are doubtful or simulation aids explanation.
  • Both samples should represent their populations and be independent within and between groups.
  • For matched before/after data, use a paired analysis; this independent-groups method is not the right wrench.

Core equations

Observed difference: Δ̂ = x̄New − x̄Standard
Resample b: Δ*ᵦ = x̄*New,b − x̄*Standard,b
Percentile CI = [ Qα/2(Δ*), Q1−α/2(Δ*) ]
x̄ = sample meanΔ̂ = observed mean differenceΔ*ᵦ = difference in resample bB = number of two-group resamples0 in CI = no difference remains plausibleOrder matters = reversing groups reverses the sign

Relatable example

A rolling mill compares average yield strength under a Standard practice and a New practice. A CI of 0.6 to 3.0 MPa for New − Standard says the plausible mean improvement is positive but modest. It does not say every New coil beats every Standard coil—the coils did not sign that contract.

Now apply it to the lesson dataset

Practice question: What is the 95% bootstrap CI for New mean − Standard mean?

Bootstrap within each group, keeping 60 New and 60 Standard observations. Compute New − Standard each time. If the confidence interval excludes zero, a zero population difference is not supported at the matching two-sided level.

Calc›Resampling›Bootstrapping for 2-Sample Means
  1. Select Samples in one column; Samples = Strength_MPa; Sample IDs = Group.
  2. Check the displayed order. This lesson defines the difference as New − Standard; reversing order reverses the sign.
  3. Choose confidence level and resamples, then click OK.
  4. Translate the interval into engineering units and direction.
Run the simulation to estimate the difference.
Method 3 of 5

Randomization test for a 1-sample mean

Hypothesis test

First: what is this method?

This test evaluates whether the mean of one population is reasonably consistent with a claimed value, μ₀. Unlike a confidence interval, it starts by temporarily assuming the null claim is true and asks: “Would our observed mean be unusual in that null world?”

How it works

  1. Calculate the observed mean and its distance from μ₀.
  2. Shift every value so the sample is centered at μ₀ while preserving its shape and spread.
  3. Resample from these null-centered values.
  4. Count how often a simulated mean is at least as extreme as the observed mean.

When to use it

  • You have one quantitative sample and a specific benchmark or target mean.
  • You want a p-value without relying solely on the one-sample t model.
  • Observations must be independent and representative.
  • Choose greater, less, or different before inspecting the answer; the p-value is not a buffet.

Core equations

H₀: μ = μ₀    H₁: μ ≠ μ₀, μ > μ₀, or μ < μ₀
Null-centered value: xᵢ⁰ = xᵢ − x̄ + μ₀
Observed distance: d_obs = x̄ − μ₀
p ≈ [1 + number of null results at least as extreme] / (B + 1)
μ = true population meanμ₀ = mean claimed by H₀x̄ = observed sample meanxᵢ⁰ = observation shifted into the null worldB = number of randomization samplesp-value = tail proportion under H₀

Relatable example

A filler is supposed to average 500 mL. Your sample averages 503 mL. The test asks whether a 3 mL departure is surprisingly large if the true mean really is 500 mL. It does not ask whether 503 mL is economically painful—that decision belongs to engineering, finance, and perhaps the person buying all that extra product.

Now apply it to the lesson dataset

Practice question: Is overall mean Strength_MPa different from 51 MPa?

Under H₀, the population mean is μ₀. A shift-based randomization test recenters every observation by subtracting the observed mean and adding μ₀. Resamples from that null-centered data show means expected when H₀ is true.

H₀: μ = 51 MPa    H₁: μ ≠ 51 MPa

Calc›Resampling›Randomization Test for 1-Sample Mean
  1. Sample = Strength_MPa; Hypothesized mean = 51.
  2. Choose the alternative: not equal, less than, or greater than. The wording must be set before seeing results.
  3. Set randomization samples (5,000 for practice) and run.
  4. The p-value is the fraction of null resamples at least as extreme as the observed result.
Run the test to compare the observed mean with the null distribution.
Method 4 of 5

Randomization test for a 1-sample proportion

Hypothesis test

First: what is this method?

This test is for a binary outcome—Pass/Fail, Yes/No, Defect/No Defect—and checks whether the population event proportion agrees with a claimed value p₀. It compares the observed event rate with rates generated when p₀ is assumed true.

How it works

  1. Define the event being counted, such as Pass.
  2. Calculate the observed proportion, p̂ = number of events ÷ sample size.
  3. Under H₀, generate many samples of n independent outcomes, each having event probability p₀.
  4. Count the simulated proportions at least as extreme as p̂ in the chosen direction.

When to use it

  • You have one sample of independent binary outcomes.
  • You want to compare the true event rate with a target or historical rate.
  • The event label must be explicit: testing Pass is the opposite direction from testing Fail.
  • If several coils come from one heat, independence may fail; 120 rows are not automatically 120 independent opportunities.

Core equations

Observed proportion: p̂ = X / n
H₀: p = p₀    Under H₀: X* ~ Binomial(n, p₀)
Simulated proportion: p̂* = X* / n
Right-tail p ≈ [1 + count(p̂* ≥ p̂)] / (B + 1)
X = observed event countn = number of independent trialsp̂ = observed sample proportionp₀ = proportion claimed under H₀Binomial(n,p₀) = null count modelB = number of simulated samples

Relatable example

A supplier promises at least 98% on-time delivery. In 100 deliveries, 99 were on time. That sounds excellent, but the test asks whether 99% is convincingly above 98% or merely one lucky extra truck. A nice-looking percentage and strong statistical evidence are not always the same guest.

Now apply it to the lesson dataset

Practice question: Is the overall Pass proportion greater than 90%?

Pass/Fail is binary. Under H₀: p = 0.90, the computer generates samples of 120 independent Bernoulli outcomes with a 90% chance of Pass. The right-tail p-value is how often the simulated pass proportion is at least the observed 92.5%.

H₀: p = 0.90    H₁: p > 0.90

Calc›Resampling›Randomization Test for 1-Sample Proportion
  1. Sample = PassFail; define the event as Pass.
  2. Hypothesized proportion = 0.90; alternative = greater than.
  3. Run 5,000 or more randomizations.
  4. Distinguish practical performance (92.5% observed) from statistical evidence (the p-value).
Run the test to see whether 92.5% is convincingly above 90%.
Method 5 of 5

Randomization test for 2-sample means

Hypothesis test

First: what is this method?

This test evaluates whether the means of two independent populations differ. It asks: “If the group label truly had no relationship with the response, how often would random label assignments produce a difference at least this large?”

How it works

  1. Calculate the observed mean difference.
  2. Pool all measured responses but keep the two group sizes fixed.
  3. Randomly shuffle the New and Standard labels.
  4. Calculate the shuffled mean difference and repeat many times.
  5. The p-value is the proportion of shuffled differences at least as extreme as the observed one.

When to use it

  • You want to test a difference between two independent group means.
  • Group labels are exchangeable under H₀—often supported by random assignment.
  • The groups should not differ systematically in product mix, time period, grade, operator, or other confounders.
  • Use a paired randomization method for matched pairs; do not shuffle away the pairing.

Core equations

H₀: μNew − μStandard = 0
Observed statistic: Δ_obs = x̄New − x̄Standard
Permutation b: Δ*ᵦ = x̄*New,b − x̄*Standard,b
Right-tail p ≈ [1 + count(Δ*ᵦ ≥ Δ_obs)] / (B + 1)
μ = population meanx̄ = sample meanΔ_obs = observed mean differenceΔ*ᵦ = difference after shuffle bExchangeable = labels can be swapped under H₀B = number of label shuffles

Relatable example

Suppose operators were randomly assigned either a New setup or Standard setup, and cycle time was recorded. Shuffling the setup labels creates the “setup does not matter” world. If the real difference is far into the tail, the labels appear to matter. If all New trials happened on night shift, however, the shuffle cannot separate setup from shift—statistics owns no time machine.

Now apply it to the lesson dataset

Practice question: Is New mean Strength_MPa higher than Standard mean?

Under H₀, the group labels are exchangeable because both groups come from populations with the same mean. Shuffle the 120 New/Standard labels, keep 60 observations in each group, and calculate New − Standard. This breaks the observed group relationship while preserving all measured strengths.

H₀: μNew − μStandard = 0    H₁: μNew − μStandard > 0

Calc›Resampling›Randomization Test for 2-Sample Means
  1. Samples = Strength_MPa; Sample IDs = Group.
  2. Set the direction to test New greater than Standard and verify the subtraction order.
  3. Run 5,000 shuffled-label samples.
  4. Use the right-tail p-value; then state the size of the observed 1.813 MPa difference separately.
Run the test to measure evidence for the New process.

3. How to read Minitab output

For a bootstrap interval

  • Original sample estimate: your observed statistic.
  • Bootstrap mean: centre of resampled statistics.
  • Standard error: resampling variability.
  • CI: plausible population values under the chosen method.

For a randomization test

  • Observed statistic: distance/difference in real data.
  • Null distribution: results expected if H₀ were true.
  • p-value: proportion at least as extreme.
  • Decision: reject H₀ only when p ≤ α.
Do not say: “p is the probability H₀ is true.” A p-value assumes H₀ is true and measures how unusual the observed statistic would be.

Assumptions resampling does not erase

  • The original observations should represent the target process or population.
  • Ordinary row resampling assumes independence. Time series, batches, heats, coils, repeated measures, and clusters may need block or cluster resampling.
  • Two-group label permutation requires exchangeability under H₀; major design differences or confounding can invalidate the shuffle.
  • Very small samples cannot reveal unseen tails or modes; resampling only reuses observed information.

4. Practice dataset

This is the same fictitious dataset supplied with the original lesson: 60 Standard and 60 New records. Pass means Strength_MPa ≥ 47. Defect is a separate fictitious quality field and is not used by the five exercises.

Rows120
Overall mean51.527
New − Standard1.813
Pass rate92.5%
SampleID Group ProductionDate Strength_MPa PassFail Defect

5. Statistics Implementation

Minitab is the main tool for this lesson, but Excel is useful for checking the mechanics. Excel does not have a single built-in “bootstrap” or “randomization test” command, so the formulas below assemble one resample at a time. That is excellent for learning and small demonstrations; for formal work, use Minitab’s dedicated procedures and record the resample count and seed.

Excel Functions

FunctionSyntaxPurposeWhen to use it
AVERAGE=AVERAGE(number1,[number2],...)Calculates a sample or resample mean.Use for the one-sample mean and both two-sample mean procedures.
MEDIAN=MEDIAN(number1,[number2],...)Calculates the middle value.Use when bootstrapping a median or another 1-sample location statistic.
STDEV.S=STDEV.S(number1,[number2],...)Calculates sample standard deviation.Use to describe the original data and compare spread before interpreting a resampling result.
COUNTIF=COUNTIF(range,criteria)Counts outcomes meeting a condition.Use to calculate Pass counts and observed proportions.
ROWS=ROWS(array)Returns the sample size.Use so a resample automatically contains the same number of observations as the source sample.
INDEX=INDEX(array,row_num,[column_num])Returns observations at randomly generated row positions.Use with RANDARRAY to draw bootstrap samples with replacement.
RANDARRAY=RANDARRAY([rows],[columns],[min],[max],[whole_number])Generates random row positions or random sorting keys.Use to build bootstrap draws, Bernoulli trials, or shuffled labels. Results recalculate unless values are frozen.
SORTBY=SORTBY(array,by_array1,[sort_order1],...)Randomly reorders labels when sorted by random keys.Use to create one permutation for a two-sample randomization test.
FILTER=FILTER(array,include,[if_empty])Returns values assigned to a selected group.Use after shuffling labels to calculate each permuted group mean.
LET=LET(name1,name_value1,calculation)Names intermediate arrays inside one readable formula.Use to keep resampling formulas understandable and auditable.
PERCENTILE.INC=PERCENTILE.INC(array,k)Returns a percentile from simulated statistics.Use the 0.025 and 0.975 percentiles for a two-sided 95% percentile bootstrap interval.

Excel Use Cases

One bootstrap mean

If Strength_MPa is in D2:D121, enter this in H2:

=LET(x,$D$2:$D$121,n,ROWS(x),draw,INDEX(x,RANDARRAY(n,1,1,n,TRUE)),AVERAGE(draw))

Fill the formula down to H5001. Each row is one bootstrap mean. Then calculate the 95% percentile interval with =PERCENTILE.INC(H2:H5001,0.025) and =PERCENTILE.INC(H2:H5001,0.975). Pressing Calculate again rebuilds the random results, so paste values if you need to preserve a run.

One shuffled mean difference

With group labels in B2:B121 and strengths in D2:D121, enter:

=LET(y,$D$2:$D$121,g,$B$2:$B$121,gp,SORTBY(g,RANDARRAY(ROWS(g))),AVERAGE(FILTER(y,gp="New"))-AVERAGE(FILTER(y,gp="Standard")))

This shuffles only the group labels and returns one New − Standard difference under the null. Fill down for many permutations, then count results at least as extreme as the observed difference to estimate the p-value.

One null-centred mean resample

To test H₀: μ = 51 MPa, use:

=LET(x,$D$2:$D$121,mu0,51,shifted,x-AVERAGE(x)+mu0,n,ROWS(x),AVERAGE(INDEX(shifted,RANDARRAY(n,1,1,n,TRUE))))

The data are shifted so their mean equals the null mean before sampling with replacement. Repeat the formula to build the null distribution.

One null pass-rate simulation

For H₀: p = 0.90, use:

=AVERAGE(--(RANDARRAY(ROWS($E$2:$E$121),1)<0.90))

Each random number below 0.90 acts like a Pass under the null. Repeating the formula builds the reference distribution for the pass proportion.

Minitab Navigation

QuestionExact Minitab pathChoose it when…
Estimate one mean, median, or other statisticCalc → Resampling → Bootstrapping for 1-Sample FunctionYou want a confidence interval for one population function.
Estimate a difference between independent meansCalc → Resampling → Bootstrapping for 2-Sample MeansYou want a confidence interval for Group 1 − Group 2.
Test one population meanCalc → Resampling → Randomization Test for 1-Sample MeanYou have one continuous sample and a claimed null mean.
Test one population proportionCalc → Resampling → Randomization Test for 1-Sample ProportionYou have binary outcomes and a claimed event proportion.
Test a difference between independent meansCalc → Resampling → Randomization Test for 2-Sample MeansYou want to test whether two group means differ by shuffling labels under the null.

Which option should you choose? Start with the decision, not the menu: use Bootstrapping to estimate uncertainty with a confidence interval; use a Randomization Test to test a null claim with a p-value. Then identify one sample versus two independent groups and continuous versus binary data.

Exam Tips

  • Know the resampling unit. ASQ CSSBB/CQE questions may hide dependence inside batches, heats, patients, shifts, or repeated readings. Resample the independent unit, not automatically every row.
  • Separate estimation from testing. A bootstrap confidence interval estimates a plausible range for a parameter; a randomization test evaluates how unusual the observed statistic is under H₀.
  • Remember “with replacement.” Bootstrap samples can contain duplicates and omit original observations. That is expected, not a computer having a small crisis.
  • Read the subtraction order. A positive New − Standard difference means New is higher. Reversing the order reverses the sign, although a two-sided significance conclusion is unchanged.
  • Do not confuse significance with importance. Compare the effect with specifications, risk, cost, and the minimum practically meaningful difference before recommending action.
  • Report reproducibility details. In real analysis, record the data source, statistic, resample count, confidence level or alpha, alternative hypothesis, and random seed.

6. Final workflow checklist

  1. Write the engineering question and define the population.
  2. Identify response type: continuous mean/function or binary proportion.
  3. Decide: estimate with a CI or test a null claim.
  4. Identify one sample or two groups; verify subtraction direction and event label.
  5. Check independence, sampling unit, representativeness, and confounding.
  6. Choose the matching Minitab procedure and alternative hypothesis before examining the p-value.
  7. Run many resamples; record the seed for reproducibility.
  8. Report the observed effect, interval or p-value, sample sizes, resample count, method, and practical meaning.
Engineering reporting model: “The New process averaged 52.434 MPa versus 50.621 MPa for Standard (difference 1.813 MPa). A 5,000-resample analysis [insert CI or p-value]… This applies to the sampled process conditions; independence and group comparability were assumed.”

7. Knowledge Check

Comprehension quiz

Pick the strongest answer, then submit. Explanations appear after you submit.

1. You need a 95% interval for a median with no normality assumption. What is the best choice?
2. What changes during a two-sample randomization test?
3. A 95% confidence interval for New − Standard is entirely above zero. What is defensible?
4. Why can ordinary row-by-row resampling be unsafe for multiple coils from the same heat?
5. A randomization test gives p = 0.18 at α = 0.05. What should you conclude?
← Back to Statistics lessons