Bootstrap resampling and dependent data
Identify the resampling unit, calculate a complete small bootstrap distribution, and preserve cluster or time dependence.
On this page
The question for today
What you’ll understand
- Explain conditional resampling from an empirical distribution.
- Calculate a bootstrap variance without confusing it with an exact population variance.
- Choose a resampling structure that respects independent clusters or serial dependence.
Repeated measurements would reveal how an estimator fluctuates, but repeating an entire experiment may be impractical. A bootstrap approximates that sampling behavior using a distribution estimated from the observed data. Its success depends on what is resampled and which features the approximation preserves. It cannot recover unobserved systematic variation by duplicating rows. We will enumerate every resample of a two-point data set, then contrast resampling independent clusters with incorrectly treating repeated readings inside each cluster as independent.
Replace an unknown population with an empirical one
For independent identically distributed observations, the nonparametric bootstrap samples the original number of observations with replacement from their empirical distribution. Conditional on the observed sample, each original row has equal selection probability, including repeated equal-valued rows as separate observations. Recompute the full estimator for each resample. A parametric bootstrap instead draws from a fitted probability model. These are approximation procedures for sampling distributions; neither changes the actual observed sample size or creates additional independent experimental information.
Distinguish two sources of approximation. Even an infinite number of bootstrap replicates uses an estimated population distribution rather than the true one. A finite number of replicates adds Monte Carlo error to summaries of that conditional distribution. More replicates reduce the second problem while leaving the first. For smooth estimators under suitable regularity the approximation can improve with sample size, but sparse tails, boundary estimates, nonsmooth extrema and weak identification can defeat familiar bootstrap intervals. The sampling assumptions must therefore be checked rather than assumed from the method's name.
Dependence determines the independent resampling unit
Measurements within one detector run can share alignment, temperature or source fluctuations. If different runs are independent clusters, resampling whole runs preserves the within-run pattern better than resampling individual readings. Equal-size cluster examples are straightforward; unequal sizes require deciding whether the target weights runs or readings. Paired measurements should normally travel together when their pairing is part of the sampling design. Resampling entries independently can erase exactly the covariance that dominates a fitted parameter's uncertainty.
For a stationary time series, nearby observations can be correlated even without natural independent clusters. The variance of its mean contains lagged covariances. An iid bootstrap ignores those terms. Block bootstraps resample stretches of consecutive observations to retain local dependence, while choosing block length balances preserved correlations against the number of available blocks. Consistency needs assumptions such as sufficiently weak dependence, appropriate moments and a suitable growing block length; no one block size or theorem covers arbitrary nonstationary or long-memory records.
Resample the procedure relevant to the target
A fit based on estimated calibration, selection or tuning can inherit uncertainty from those operations. A bootstrap that treats them as fixed answers a conditional question. Repeating them within each appropriate resampling scheme addresses a broader procedure, but only if the relevant data and dependence structure are represented. A fixed unknown calibration offset absent from the data cannot be learned by an ordinary row bootstrap. External calibration must enter separately or be resampled at its genuinely independent experimental level.
Report the resampling unit, scheme, number of replicates, statistic and interval rule, together with any block-length or model sensitivity. A percentile interval is not automatically a confidence interval with its nominal coverage. Exact enumeration of a tiny bootstrap distribution is valuable for understanding mechanics, but two observations or two clusters cannot establish an asymptotic coverage guarantee. Independent validation and physics-based uncertainty budgets remain necessary even when a bootstrap histogram appears smooth and numerically stable.
Work it through
Enumerate a two-reading bootstrap
The iid observed readings are 0 and 2 mV. Enumerate all ordered bootstrap samples of size two, find the conditional distribution of their means, and calculate its exact variance.
Each draw chooses either observed value with probability one half. The four equally probable ordered samples are (0,0), (0,2), (2,0) and (2,2), with means 0, 1, 1 and 2 mV.
The conditional mean is one mV. The exact bootstrap variance is ; divide by four because the full distribution is being enumerated.
The original unbiased sample variance is two mV squared, so its usual plug-in mean-variance estimate is one mV squared. The bootstrap's half is a finite empirical-distribution effect, not proof that the true mean is known more accurately.
The bootstrap means are 0, 1, 2 mV with probabilities 1/4, 1/2, 1/4, and exact conditional variance 0.5 mV².
The living notebook
Explore the model
Interactive lab · sampling
Increasing nominal sample size assumes independence
Reset and set BOTH Success probability p = 0.5 and Independent samples n = 4; predict variance 0.0625 for the sample mean. For a fresh comparison, Reset and set p = 0.5 and n = 8; predict variance 0.03125.
- Mean sample proportion
- 0.5
- Sample-mean variance
- 0.0125
- Uniform CDF error bound
- ≤ 0.1062
Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forThe lab calculates an exact iid binomial sampling law and deterministic highlighted outcomes. It does not bootstrap data, resample clusters, estimate autocorrelation or justify treating dependent readings as independent.
Pause & reason
Check your understanding
Your turn
Try it yourself
Two independent equal-size clusters contain readings (0,0) and (2,2) mV. Resample two whole clusters with replacement and average all four readings. Find the exact conditional variance. Compare with an iid row bootstrap of the four original readings.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The cluster means are zero and two. A whole-cluster resample has the same mean distribution as the worked example.
Compare with the solution
Whole-cluster resampling gives mean variance . An iid row bootstrap draws four values from an empirical distribution with variance one mV squared, giving mean variance . It underestimates the cluster-bootstrap variance by a factor two because it destroys the perfect within-cluster association. With only two clusters neither calculation establishes reliable interval coverage for a general population of runs; the example diagnoses the consequence of the resampling unit.
Keep the idea
Take it with you
- Bootstrap replication approximates sampling variation without creating new information.
- Dependence and experimental grouping determine the resampling unit.
- Finite-replicate precision and statistical validity are different questions.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.