Why Score Matches Can Still Be Poor Covariate Matches
Propensity-score matching compresses several pretreatment characteristics into one estimated probability of treatment. That reduction can be useful, but it can also obscure an important diagnostic question: after matching, do the treated and comparison groups actually look alike in the covariates that matter? A small difference in estimated scores is not itself evidence that every covariate is balanced.
Balance is a property of the matched sample, not a guarantee supplied by the matching algorithm. This matters because residual differences in measured covariates can undermine the credibility of a treatment-effect comparison. Balance checks can reveal such differences; they cannot establish that unmeasured confounding is absent or that the causal assumptions hold.
The Balancing Property and Its Limits
Let \(T\) indicate treatment, with \(T=1\) for treated units and \(T=0\) for comparison units, and let \(X\) denote the vector of pretreatment covariates. The propensity score is the conditional probability of treatment given those covariates. For a causal interpretation, a common set of assumptions includes consistency, conditional exchangeability, and positivity. In particular, conditional exchangeability requires that, given \(X\), treatment assignment is independent of the potential outcomes \(Y(1)\) and \(Y(0)\); positivity requires a nonzero probability of either treatment status at covariate values relevant to the target population.
The balancing result can be seen directly. Conditional on a value \(e\) of the true score, the treatment probability is the same for every covariate value in that score group. For any covariate set \(A\), the probability of a treated unit having covariates in \(A\), conditional on \(e\), therefore factors as follows:
The same argument applies to comparison units when \(0<e<1\), so the covariate distributions agree between treatment groups conditional on the true score. But matching usually uses an estimated score, \(\hat e(X)\), and selects a finite set of pairs or neighbors. The estimate need not preserve the true balancing property exactly. Even if the score model is well specified, estimation error and the discrete selection of matches can leave residual differences. And even when scores are nearly equal, a single scalar distance does not require each component of \(X\) to be close.
How to Assess Covariate Balance
A balance table compares each pretreatment covariate between the treated and comparison groups, both before and after matching. For a continuous covariate \(X_j\), a common summary is the standardized mean difference (SMD): the difference in group means divided by a standard deviation used as a reference. A useful convention is to use the pooled standard deviation from the original, unmatched sample for both the before- and after-matching values. Keeping the denominator fixed makes the two SMDs easier to compare.
Here the subscripts 1 and 0 indicate treated and comparison groups, and the reference standard deviation is calculated before matching. For a binary covariate, the analogous standardization uses the pooled proportion \(p_{\mathrm{pool}}\):
After matching with replacement, comparison units that are selected multiple times should contribute according to the matching weights when calculating post-match summaries. The same principle applies to any matching scheme that assigns unequal weights. The chosen convention for standardization should be stated, since different denominators can produce different numerical SMDs.
SMDs describe the size of a difference relative to a scale; they are not hypothesis tests. A small p-value can reflect a large sample even when a covariate difference is practically modest, while a nonsignificant test does not establish balance. Absolute SMDs below 0.1 are sometimes used as a rough screening convention, but this is not a universal pass/fail rule. Consider the covariates’ substantive importance, their distributions, and the matching design.
A table should be supplemented with plots when distributional features may be hidden by a mean. Useful displays include side-by-side histograms or empirical cumulative distribution functions for continuous variables, and dot plots of absolute SMDs across covariates (often called Love plots). Check not only the covariates entered as simple terms, but also substantively important nonlinear transformations and interactions. If matching was intended to retain treated units, also report how many treated units were retained and make clear that any resulting effect estimate pertains to the retained group rather than automatically to all treated units.
Worked Example: Calculating and Comparing Standardized Differences
Suppose a hypothetical study has 40 treated and 40 comparison units before matching. Before matching, mean age is 52 years in the treated group and 48 years in the comparison group; the sample standard deviations are 10 and 12 years. After matching, 32 units remain in each group, with mean ages of 51.2 and 50.0 years. Use the original pooled standard deviation as the denominator for both SMDs.
First calculate the reference standard deviation:
The absolute SMD before matching is
After matching, using the same reference denominator, it is
The age difference is smaller after matching, but its absolute SMD remains slightly above 0.1, a commonly used screening benchmark. This describes residual imbalance in age; it does not by itself determine whether the causal estimate is biased or whether the matched design is acceptable. Other covariates and distributional checks are also needed.
Why the Estimated Score Can Hide Imbalance
There are several distinct reasons that score matching can fail to balance observed covariates. First, the score model may omit nonlinear terms or interactions that help explain treatment assignment. Second, even a well-specified model produces estimates from finite data. Third, matching on one-dimensional score distance may not select the covariate combinations that best balance the full vector. Finally, ties, caliper choices, limited overlap, or the order in which neighbors are selected can affect which units are retained. These mechanisms are different: a table can reveal imbalance, but it does not by itself identify which mechanism caused it.
The following example shows the difference between a fitted-score match and balance on a covariate transformation. It is a hypothetical illustration of how a model using only a linear term can fail to represent a nonlinear treatment-assignment pattern.
Worked Example: A Linear Score Can Miss a Nonlinear Difference
Consider eight hypothetical units with a continuous covariate \(X\). The treated units have \(X=-3,-1,1,3\); the comparison units have \(X=-2,-1,1,2\). Suppose the estimated propensity score uses only an intercept and a linear term in \(X\).
The treated and comparison groups both have mean \(X=0\). Across all eight units, the sum of \(X\) is also 0. In a logistic regression with an intercept and \(X\), an intercept of 0 gives fitted probability \(0.5\) for every unit, and the score equation for the slope is satisfied: the sum of \(X(T-0.5)\) is \(0-0.5(0)=0\). Because the logistic likelihood is concave, these values give a fitted model with slope 0 and intercept 0. The estimated scores are therefore all \(0.5\), so matching on that score cannot distinguish the units.
Now check the nonlinear covariate \(X^2\). Its treated mean is \((9+1+1+9)/4=5\), while its comparison mean is \((4+1+1+4)/4=2.5\). The groups match exactly on the mean of \(X\) and have the same estimated score, yet differ by 2.5 units in mean \(X^2\). Using the original-sample pooled standard deviation for \(X^2\), the within-group sample variances are 21.33 for the treated group and 3.000 for the comparison group, giving
Thus, balance on the fitted score and on the linear covariate mean does not imply balance on a nonlinear feature of the covariates. In practice, the appropriate response is to diagnose important transformations and, if needed, revise the score specification or matching design and reassess balance.
Small Score Distances Do Not Guarantee Covariate Balance
Even without a nonlinear-model example, a score match can have tiny distances and still leave a large difference in an observed covariate. The next example separates two quantities that are easy to confuse: closeness of estimated scores and balance of covariates.
Worked Example: Close Scores but an Imbalanced Binary Covariate
Suppose four hypothetical treated units have estimated scores \(0.40, 0.45, 0.55, 0.60\), and four matched comparison units have scores \(0.41, 0.44, 0.56, 0.59\). Pair units in that order. Each pair differs by only 0.01 in score, and both groups’ average score is 0.50.
Let \(Z\) be a binary pretreatment covariate. Its values are \(0,0,1,1\) among the treated units and \(0,0,0,0\) among the comparison units. The group proportions are therefore \(\hat p_1=2/4=0.50\) and \(\hat p_0=0/4=0\). The pooled proportion across the eight matched units is \(p_{\mathrm{pool}}=2/8=0.25\). The absolute SMD is
The score means agree and every pair has a small score gap, yet \(Z\) remains substantially imbalanced. This is why a balance assessment must examine the original covariates, not only the estimated score.
A Practical Diagnostic Sequence
A disciplined diagnostic sequence makes it less likely that a single reassuring score summary will be mistaken for evidence of balance. It also keeps the design decision separate from the subsequent outcome comparison.
List pretreatment covariates relevant to treatment assignment and outcome prediction. Include important nonlinear terms or interactions when substantively justified.
Report group summaries and consistently defined SMDs. State the standardization denominator and account for matching weights or repeated use of comparison units.
Use plots to look for differences in spread, tails, or shape that means alone can miss. Report how many units were unmatched or excluded and describe the resulting target group.
Consider revising the score specification, caliper, distance measure, or matching design. Recalculate the diagnostics after each change; do not assume that changing the score automatically fixes balance.
Common Mistake: Treating a Good Score Match as a Balance Test
A tempting but incorrect conclusion is: “The matched scores are nearly equal, so the covariates must be balanced.” This confuses a theoretical property of the true propensity score in the population with a finite-sample result obtained by matching on an estimated score. Score closeness is useful design information, but it is not a substitute for checking each covariate and its important distributional features.
The correction is to report balance diagnostics directly and interpret them as diagnostics, not as proof of causal identification. Good measured-covariate balance supports the design’s comparability on those measured characteristics. It does not verify conditional exchangeability, guarantee overlap, or rule out bias from unmeasured factors. Conversely, a detected imbalance is a reason to investigate the design before interpreting an outcome contrast as a treatment effect.
Check Your Understanding
Answer these questions using the definitions and examples above.
- Why does the population balancing property of the true propensity score not guarantee balance after matching on an estimated score?
- A matched sample has equal treated and comparison means for a covariate. Name one distributional feature that could still differ.
- Why is a fixed, pre-match reference standard deviation useful when comparing SMDs before and after matching?
- In the binary-covariate example, the average estimated scores agreed. What calculation showed that the covariate remained imbalanced?
- What can good balance on measured covariates support, and what causal assumption does it not establish?