STAT 101
M05 · L03
Module 5 · Lesson 3
Bootstrap Methods
What if you don’t know the population distribution? Pull yourself up by your own data. Resample from what you have and let computation do the math.
01 / 12
STAT 101
M05 · L03
Bradley Efron, 1979
The Audacious Idea
Classical inference needs a formula for every statistic. Bootstrap’s insight: treat the sample as a stand-in for the population. Simulate the sampling process using data you already have.
The substitution
Unknown F → Empirical F̂ₙ
Derive → Simulate B times
Derive → Simulate B times
02 / 12
STAT 101
M05 · L03
Four simple steps
The Algorithm
- Step 1: Draw n observations with replacement from your data → bootstrap sample x*
- Step 2: Compute the statistic T* = T(x*)
- Step 3: Repeat B times (B = 1,000–10,000)
- Step 4: Analyze the distribution of {T*₁, …, T*_B}
03 / 12
STAT 101
M05 · L03
No formula derivation needed
Bootstrap Standard Error
Bootstrap SE
\widehat{\mathrm{SE}}_{\mathrm{boot}} = \sqrt{\dfrac{1}{B-1}\sum_{b=1}^{B}(T^*_b - \bar{T}^*)^2}
Standard deviation of the B bootstrap statistics. Works for any estimator — mean, median, correlation, ratio.
04 / 12
STAT 101
M05 · L03
Simplest bootstrap CI
Percentile Interval
95% Percentile CI
\bigl[\,T^*_{(\alpha/2)},\; T^*_{(1-\alpha/2)}\,\bigr]
Take the 2.5th and 97.5th percentiles of the bootstrap distribution. Intuitive and widely applicable when n is large.
05 / 12
STAT 101
M05 · L03
Bias-corrected version
Basic Bootstrap CI
Pivot / Reflection CI
\bigl[\,2\hat{T} - T^*_{(1-\alpha/2)},\; 2\hat{T} - T^*_{(\alpha/2)}\,\bigr]
Gold standard
BCa (bias-corrected & accelerated) — adjusts for bias AND skewness. Best coverage for small n.
06 / 12
STAT 101
M05 · L03
The power of universality
Works for Any Statistic
- Median: no closed-form SE exists — bootstrap provides it
- Correlation: more robust than Fisher z-transform
- Ratio of means: no parametric formula
- AUC / F1-score: standard for ML metrics
- Regression coefficients: robust to non-normal residuals
07 / 12
STAT 101
M05 · L03
Two flavors
Parametric vs. Nonparametric
Nonparametric
Resample the data. No assumption. Universal.
Parametric
Simulate from fitted model. Efficient if model correct.
08 / 12
STAT 101
M05 · L03
Theoretical foundation
Why It Works
- F̂ₙ converges to F (Glivenko–Cantelli theorem)
- Bootstrap distribution of T*−T̂ converges to the sampling distribution of T̂−θ
- Bootstrap CIs are asymptotically valid
- Quality depends on how well F̂ₙ approximates F
- Larger n → better approximation
09 / 12
STAT 101
M05 · L03
Known failure modes
When Bootstrap Fails
- Sample max/min: bootstrap max can’t exceed observed max
- Heavy tails: infinite variance → bootstrap may not converge
- Tiny samples (n < 15): F̂ₙ too sparse; poor coverage
- Dependent data: use block bootstrap for time series
- Model selection: LASSO variable selection is not bootstrap-stable
10 / 12
STAT 101
Summary
Recap
What you learned
- Resample with replacement B times; compute the statistic each time
- Bootstrap SE = std dev of the B replicates — no formula needed
- Percentile CI = [T*_(α/2), T*_(1−α/2)] — simplest method
- BCa CI — best coverage when bootstrap is skewed or biased
- Nonparametric: resample data. Parametric: simulate from fitted model
- Fails for extremes, heavy tails, small n, and dependent data
Module 5 complete
Up next — Module 6: Hypothesis Testing →
12 / 12