What Is Uncertain in a Balance Summary?
A balance summary describes the covariate differences in a particular treated and comparison sample after an adjustment procedure. For that realized sample and a fully specified procedure, a summary such as a standardized mean difference is a calculated quantity, not a random one. Uncertainty enters when the observed units are viewed as a sample from a wider population, or when the adjustment procedure itself can produce different matched samples.
These sources should not be conflated. Sampling uncertainty concerns how the summary might vary across samples of units from the target population. Procedure uncertainty can arise because propensity scores are estimated, matches depend on which units are observed, or the algorithm uses random tie-breaking. Specification sensitivity is variation in results across deliberate choices, such as different calipers or matching rules. A range across specifications is informative, but it is not automatically a confidence interval.
For a continuous covariate \(X\), one common summary is the standardized mean difference (SMD):
Here \(\bar{x}_T\) and \(\bar{x}_C\) are the treated and comparison-group sample means, \(s_T\) and \(s_C\) are their sample standard deviations, and \(n_T\) and \(n_C\) are their sample sizes. The SMD is unitless; its sign indicates the direction of the mean difference. For binary covariates, a difference in proportions can be standardized using a specified scale, but the scale should be reported because different conventions give different summaries.
Resampling and Its Interpretation
A bootstrap approximates repeated sampling by repeatedly drawing units with replacement from the observed data and recalculating a statistic. If \(\hat{\theta}\) is the observed balance summary and \(\hat{\theta}^{*(b)}\) is its value in bootstrap replicate \(b\), the bootstrap standard error is the sample standard deviation of the replicate values:
A percentile interval uses selected quantiles of the replicate summaries. For example, the 2.5th and 97.5th percentiles form a nominal 95% bootstrap percentile interval. The word “nominal” matters: an interval has a 95% coverage interpretation only when the resampling procedure adequately approximates the relevant sampling distribution and the method’s conditions support that coverage. The calculation alone does not guarantee it.
The resampling unit must match the independent sampling unit. If observations are clustered, resampling individual rows as if they were independent can understate variability; the clusters may need to be resampled as units. In a matched analysis, resampling individual treated and comparison observations separately can also break the matched structure. If the statistic is defined on matched pairs and the pairs are independent sampling units, resampling whole pairs preserves that structure. It does not, however, capture uncertainty from how those pairs were selected unless the matching procedure is repeated as part of the resampling.
Decide whether the goal is to describe sampling variation, variation from the matching procedure, or sensitivity to analytic choices. These are different questions.
Resample independent units, preserving clusters or pairs when those structures are part of the data.
To include propensity-score estimation and matching variation, refit and rematch in each replicate. Recomputing a summary on a fixed matched sample answers a narrower, conditional question.
Describe the resampling scheme and any instability. Do not interpret an interval as a balance guarantee or as evidence of causal validity.
Worked Balance-Uncertainty Examples
Worked Example: A Standardized Difference and Its Bootstrap Variation
Suppose a hypothetical matched sample has 40 treated and 60 comparison units. For a covariate measured in years, the treated mean is 52, the comparison mean is 49, and the sample standard deviations are 10 and 12 years. First calculate the SMD.
The pooled variance in the denominator is \( [(40-1)(10^2)+(60-1)(12^2)]/(40+60-2) \). Substituting the values gives \( [39(100)+59(144)]/98 = 12{,}396/98 = 126.4898 \). Thus the pooled standard deviation is \( \sqrt{126.4898}=11.247 \) years, and \( \operatorname{SMD}=(52-49)/11.247=0.2667 \). The observed treated mean is therefore about 0.267 pooled standard deviations higher than the comparison mean.
Now suppose, for illustration, that eight resampling replicates of the full specified calculation produce SMDs of 0.18, 0.22, 0.24, 0.27, 0.29, 0.31, 0.34, and 0.39. Their mean is \( \bar{\theta}^{*}=2.24/8=0.28 \). The sum of squared deviations from 0.28 is 0.032, so the bootstrap standard error is \( \sqrt{0.032/(8-1)}=\sqrt{0.004571}=0.06761 \) SMD units.
This small set of replicates illustrates the calculation, not a dependable interval estimate: eight replicates are far too few to characterize tail quantiles well. With an appropriately large number of replicates, percentile limits could be calculated, but their validity would still depend on what was resampled and which steps were repeated. The result describes variation in this balance summary under that resampling scheme; it does not say that the observed imbalance is acceptable.
Worked Example: Resampling Matched Pairs
In a hypothetical matched sample, suppose the covariate differences within four treated-comparison pairs are 1, 2, 3, and 6 units. The mean paired difference is \( (1+2+3+6)/4=3 \) units. If these pairs are treated as independent sampling units, a bootstrap replicate should draw four pairs with replacement, not draw treated and comparison rows independently.
Consider six illustrative bootstrap replicates of the pair differences: (1, 2, 3, 6), (1, 1, 3, 6), (2, 2, 3, 6), (1, 3, 3, 6), (1, 2, 6, 6), and (1, 2, 2, 3). Their means are 3, 2.75, 3.25, 3.25, 3.75, and 2.00 units. The mean of these replicate means is 3.00 units.
The sum of squared deviations of those six means from 3.00 is \(0^2+(-0.25)^2+0.25^2+0.25^2+0.75^2+(-1)^2=1.75\). Therefore the illustrative bootstrap standard error is \( \sqrt{1.75/(6-1)}=\sqrt{0.35}=0.5916 \) units.
The calculation preserves the pair structure, but six replicates do not support a reliable interval. More importantly, if the pairs were created by a matching algorithm applied to a larger cohort, resampling the already formed pairs conditions on the selected matches. To include variation from match selection, one would need to resample the appropriate original units and rerun the matching procedure, if that procedure and resampling method support the intended inference.
Worked Example: A Specification Range Is Not a Confidence Interval
Suppose an analyst compares three prespecified matching choices in a hypothetical data set. For one covariate, their absolute SMDs are 0.08, 0.11, and 0.24. The minimum is 0.08, the maximum is 0.24, and the range is \(0.24-0.08=0.16\) SMD units.
If the analyst had specified 0.10 as a practical balance tolerance before examining the results, only the first choice meets that tolerance for this covariate. The range from 0.08 to 0.24 shows that the conclusion is sensitive to the matching choice. It is not a 95% interval: the endpoints are results from different procedures, not quantiles from a sampling distribution. Nor does the first choice establish overall balance, since other covariates and distributional features may still differ.
What an Interval Can and Cannot Establish
An interval can communicate the precision of a summary under its stated resampling scheme. A narrow interval suggests little variation in that summary under that scheme; a wide one indicates greater variation. Neither width nor location replaces examination of the underlying balance diagnostics. For example, an interval for an SMD that includes zero does not establish balance: it may simply be too wide to distinguish meaningful differences from zero.
A balance threshold is a substantive or design criterion, not a hypothesis test. If an interval lies inside a prespecified tolerance range, that is evidence that the summary is consistently within that range under the stated resampling procedure. It still does not prove that every covariate is balanced, that distributions match beyond their means, or that confounding has been removed. Conversely, an interval that crosses the tolerance boundary signals uncertainty about whether this particular summary meets the criterion; it does not justify selecting whichever endpoint supports a preferred conclusion.
Resampling also has failure modes. Matching can be a non-smooth procedure: a small change in the sample may change which units are matched, and some bootstrap samples may produce poor overlap or few usable matches. Ordinary bootstrap intervals are not guaranteed to have accurate coverage for every matching method and balance statistic. Report failed or unstable replicates rather than silently discarding them, and do not present a nominal interval as calibrated unless the method’s validity is justified for the setting.
The uncertainty in a balance summary is separate from uncertainty in an estimated treatment effect. Balance diagnostics assess observed covariates after adjustment; they do not quantify unmeasured confounding or establish causal identification. In contrast with a single resampling interval, examining several diagnostics and prespecified design choices can reveal different kinds of imbalance or sensitivity. These approaches complement one another, but a specification range should remain labelled as sensitivity analysis, not sampling uncertainty.
Common Mistakes and Practical Reporting
- Calling a specification range a confidence interval. Values across matching choices describe sensitivity to those choices, not repeated-sampling variation.
- Resampling the wrong unit. Resampling rows independently can destroy cluster or pair dependence. State the sampling unit and preserve the relevant structure.
- Resampling only the final matched data without qualification. This may describe uncertainty conditional on the selected matches, not uncertainty from estimating scores or choosing matches.
- Treating an interval that includes zero as proof of balance. It may reflect limited precision; also assess magnitude, distributions, and the prespecified practical tolerance.
- Assuming a bootstrap automatically provides valid coverage. Matching can be unstable or non-smooth. The resampling method must be appropriate to the statistic and procedure.
- Interpreting balance uncertainty as causal-effect uncertainty. Balance summaries concern measured covariates, not the treatment effect or unmeasured confounding.
A clear report identifies the balance statistic, the resampling unit, the number of replicates, which estimation and matching steps were repeated, the interval method, and any failed or unstable replicates. If the purpose is to compare design choices, report the choices and their summary values separately from any sampling interval. This makes the scope of the uncertainty statement explicit.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A bootstrap resamples the final matched pairs but does not rerun propensity-score estimation or matching. Which uncertainty does it describe most directly?
- A 95% bootstrap interval for an SMD includes zero. Does this establish that the covariate is balanced? Explain.
- Three matching specifications produce absolute SMDs of 0.06, 0.14, and 0.20. What is their range, and why is it not a confidence interval?
- Data are sampled in intact clinics, with patients within each clinic potentially correlated. What should be resampled to respect the sampling structure?
Solutions. (1) It describes variation conditional on the already selected matched pairs; it does not directly include uncertainty from fitting scores or selecting matches. (2) No. The interval may be wide, and including zero does not show that relevant differences are small. Balance requires interpreting the magnitude and examining other diagnostics against an appropriate criterion. (3) The range is \(0.20-0.06=0.14\) SMD units. Its endpoints come from different specifications, not sampling quantiles. (4) Resample clinics as clusters, carrying their patients together, rather than treating every patient as an independent sampling unit.