Tutorials › Propensity Score Matching and Instrumental Variables › Full Matching For A Comparative Effectiveness Study

Worked Examples · Tutorial 676 of 699

Full Matching For A Comparative Effectiveness Study

Learn to construct ATT weights from full matched sets, evaluate balance, and calculate and interpret a weighted outcome contrast with an uncertainty estimate.

Intermediate 12 min read

What You'll Learn

  • Define full matching sets and calculate their ATT-oriented weights.
  • Express the matched ATT as a weighted difference in outcome means.
  • Assess covariate balance using weighted summaries and standardized differences.
  • Estimate uncertainty for a stratified matched-set contrast and interpret an interval.
  • Recognize failures of overlap, balance, and small-sample inference.

Why Full Matching Can Help

A comparative-effectiveness cohort may contain many treated individuals but relatively few comparable controls in some parts of the covariate distribution, and the reverse in other parts. A fixed 1:1 match can discard usable observations or force a uniform ratio where the data do not support one. Full matching instead groups retained treated and control individuals into matched sets with variable treatment-to-control ratios. For an average treatment effect on the treated (ATT), the sets can be weighted so that each treated individual contributes equally to the treated outcome mean and the controls collectively represent the treated group.

The method does not make treatment assignment random. Its causal interpretation still depends on measured covariates adequately controlling confounding, adequate overlap, and a correctly defined treatment, outcome, and follow-up period. The introductory material on propensity score matching provides the broader context; here the emphasis is on the mechanics and checks in a worked comparative-effectiveness study.

Definition: Let \(A_i=1\) denote treatment and \(A_i=0\) denote control, \(X_i\) the measured baseline covariates, and \(Y_i\) the outcome. A full matching set contains at least one treated and at least one control individual. Sets may have different numbers of treated and control individuals. If set \(j\) contains \(n_{Tj}\) treated individuals and \(n_{Cj}\) controls, ATT-oriented weights are 1 for each treated individual and \(n_{Tj}/n_{Cj}\) for each control. Thus the total control weight in set \(j\) equals \(n_{Tj}\), the number of treated individuals in that set.

The weights make the control group’s weighted size equal the treated group’s size, both overall and within each matched set. If the retained treated group has \(n_T=\sum_j n_{Tj}\) individuals, the weighted ATT contrast is the treated mean outcome minus the weighted control mean outcome. It targets the ATT for the treated individuals retained in the analysis, not automatically for every treated individual in the original cohort.

$$ \widehat{\mathrm{ATT}} = \frac{1}{n_T}\sum_j \left[ \sum_{i\in T_j}Y_i - \frac{n_{Tj}}{n_{Cj}}\sum_{i\in C_j}Y_i \right] $$

Equivalently, define the within-set contrast \(\widehat{\Delta}_j=\bar{Y}_{Tj}-\bar{Y}_{Cj}\). Then \(\widehat{\mathrm{ATT}}=\sum_j(n_{Tj}/n_T)\widehat{\Delta}_j\). Sets contribute in proportion to their number of treated individuals: this is important when set sizes differ.

Assumptions Before Interpreting the Contrast

The matching calculation is descriptive until causal identification assumptions are made. Let \(Y_i(1)\) and \(Y_i(0)\) be individual \(i\)’s potential outcomes under treatment and control. The ATT is \(E[Y(1)-Y(0)\mid A=1]\) in a specified target population. A causal interpretation for the matched estimate requires, at minimum, the following conditions:

  • Consistency and a well-defined treatment: the observed outcome equals the potential outcome for the treatment actually received, and treatment versions and follow-up are defined clearly.
  • Conditional exchangeability: given the measured baseline covariates \(X\), treatment is independent of the potential outcomes. In particular, there are no unmeasured confounders remaining after conditioning on \(X\).
  • Positivity for the treated target: treated individuals have a nonzero probability of receiving control at their covariate values, so suitable controls exist across the covariate profiles represented by the target treated group.
  • No interference: one individual’s treatment does not change another individual’s outcome, unless the estimand and design explicitly account for such interference.

Full matching may use an estimated propensity score \(e(X)=P(A=1\mid X)\), a distance based on covariates, or a specified combination. A good-looking match cannot establish exchangeability, and weighting cannot manufacture control information where there is no overlap. The target population must also be reported accurately if some treated individuals are excluded for lack of comparable controls.

Worked Example: Constructing Full-Matching Weights

Suppose a hypothetical comparative-effectiveness cohort yields four matched sets after matching on baseline covariates. Sets A, C, and D each contain two treated individuals and two controls. Set B contains two treated individuals and three controls. All eight treated individuals are retained.

SetTreated \(n_{Tj}\)Controls \(n_{Cj}\)Weight per control \(n_{Tj}/n_{Cj}\)Total control weight
A2212
B232/32
C2212
D2212

For example, in set B each control receives weight \(2/3\), so the three controls together contribute \(3(2/3)=2\), equal to the two treated individuals’ total weight. Across the four sets the controls’ total weight is \(2+2+2+2=8\), the same as the treated sample size. Giving each control weight 1 instead would give set B three control contributions for two treated contributions and would no longer implement this ATT weighting scheme.

Outcome Contrast and a Justification for the Weights

Consider hospital length of stay, measured in days, as a hypothetical outcome. Suppose the treated outcomes in sets A through D are respectively \(10,12\); \(8,10\); \(7,9\); and \(11,13\). Suppose the control outcomes are respectively \(12,14\); \(10,11,12\); \(9,11\); and \(13,15\). Each set’s treated mean is two days below its control mean.

Worked Example: Estimating the ATT in Days

The treated outcome total is \(22+18+16+24=80\), so the treated mean is \(80/8=10\) days. The control contributions, applying the weights from the previous calculation, are \(26\) in set A, \((2/3)(10+11+12)=22\) in set B, \(20\) in set C, and \(28\) in set D. Their total is \(96\); their total weight is 8. The weighted control mean is therefore \(96/8=12\) days.

Thus the weighted contrast is

$$ \widehat{\mathrm{ATT}}=10-12=-2.00\text{ days}. $$

In this hypothetical matched cohort, the treated group had an estimated mean stay two days shorter than the weighted control group. This is an estimated causal effect for the retained treated population only if the causal assumptions are credible; the calculation alone establishes a weighted association, not causation.

The weighting logic follows directly from the set structure. In set \(j\), the weighted control mean is \(\bar{Y}_{Cj}\), because multiplying every control by \(n_{Tj}/n_{Cj}\) changes the total weight but not the mean. The set’s weighted control total is \(n_{Tj}\bar{Y}_{Cj}\). Summing across sets and dividing by \(n_T\) gives a control mean in which each set is represented in proportion to its treated count. Subtracting that mean from the treated mean gives the stated stratified ATT contrast.

Assessing Balance Before Reading the Outcome Result

Matching should be evaluated using baseline covariates, not by asking whether treatment and control groups have statistically significant differences in those covariates. Standardized mean differences (SMDs), weighted means and proportions, and distributional summaries describe the magnitude of imbalance without making balance depend on sample size. Examine covariates used in the design and important prognostic variables, including nonlinear terms or interactions when scientifically relevant. Similar means alone can hide differences in spread or distribution.

For a continuous covariate, one common diagnostic is the difference between the treated mean and weighted control mean divided by a prespecified reference standard deviation, often the pooled standard deviation before matching. For a binary covariate, the difference in proportions can be divided by the standard deviation based on the average of the two proportions. The chosen convention should be stated and used consistently; an SMD is a diagnostic, not a hypothesis test or a universal pass/fail guarantee.

Worked Example: Reading Two Balance Diagnostics

In the hypothetical cohort, ages for the treated individuals are \(60,64; 50,54; 70,74; 55,59\) across sets A through D. Their total is \(486\), giving a treated mean age of \(486/8=60.75\) years. Let control ages be \(60,64; 50,52,54; 70,74; 55,59\). After applying the set weights, the control total is \(124+(2/3)(156)+144+114=486\), with total weight 8. The weighted control mean is also \(60.75\) years. If the prespecified reference standard deviation is 10.0 years, the age SMD is \((60.75-60.75)/10.0=0.000\).

Now suppose five of the eight treated individuals have a high baseline severity indicator. Suppose the weighted control count with high severity is 4, out of total control weight 8. The proportions are \(5/8=0.625\) and \(4/8=0.500\). Using the average proportion \(0.5625\), the binary standardizer is \(\sqrt{0.5625(1-0.5625)}=0.4961\), so the SMD is \((0.625-0.500)/0.4961=0.252\). Age is well aligned by this diagnostic, while severity still shows a meaningful imbalance that calls for investigation. Excellent balance on age does not cancel imbalance on severity.

Balance diagnostics should be computed after applying the actual analysis weights. Also inspect propensity-score overlap and the distribution of weights. A full-matching set with many treated individuals but very few controls gives each control a larger weight; such a set may reveal limited local overlap and can make the estimate sensitive to a small number of controls. Report set sizes, exclusions, weight ranges, and important weighted balance results so readers can judge what population and comparisons support the estimate.

Uncertainty for a Matched-Set Contrast

One model-based variance calculation treats matched sets as independent and estimates uncertainty within each set. Write \(s^2_{Tj}\) and \(s^2_{Cj}\) for the sample variances of the treated and control outcomes in set \(j\). Under independent sampling within treatment groups and sets, and conditional on the matched-set structure, the estimated variance of the set contrast is \(s^2_{Tj}/n_{Tj}+s^2_{Cj}/n_{Cj}\). This sample-variance estimator requires \(n_{Tj}\ge 2\) and \(n_{Cj}\ge 2\) in every set; if either group is a singleton, this formula is undefined for that group, and an alternative variance model or design-appropriate method is needed. Since the overall ATT weights set \(j\) by \(n_{Tj}/n_T\), independence across sets gives the following estimated variance:

$$ \widehat{\mathrm{Var}}(\widehat{\mathrm{ATT}}) = \sum_j \left(\frac{n_{Tj}}{n_T}\right)^2 \left( \frac{s^2_{Tj}}{n_{Tj}}+ \frac{s^2_{Cj}}{n_{Cj}} \right). $$

This is a conditional, model-based calculation, not a universal variance formula for every matching procedure. In particular, it does not automatically account for uncertainty from estimating the propensity score or selecting the matched sets. A design-appropriate robust or resampling procedure may be preferable, and resampling should repeat the relevant design steps when the aim is to include their uncertainty.

Worked Example: Standard Error and Interval

For the length-of-stay data, each treated pair has sample variance \(2\) days squared. The control sample variances in sets A through D are \(2,1,2,2\), respectively. Each set has two treated individuals, so each set’s weight in the ATT is \(2/8=1/4\). The estimated variance is

$$ \begin{aligned} \widehat{\mathrm{Var}}(\widehat{\mathrm{ATT}}) &=\left(\frac14\right)^2 \left[ \left(\frac22+\frac22\right) +\left(\frac22+\frac13\right) +\left(\frac22+\frac22\right) +\left(\frac22+\frac22\right) \right]\\ &=\frac1{16}\left(2+\frac43+2+2\right) =\frac{11}{24} \approx 0.4583\text{ days}^2. \end{aligned} $$

Therefore \(SE=\sqrt{0.4583}=0.6770\) days. A normal-approximation 95% interval is \(-2.00\pm1.96(0.6770)\), which is approximately \([-3.327,-0.673]\) days. Under the stated model and approximation, the interval describes a range of values compatible with the data and procedure; it is not a 95% probability statement about the fixed ATT. There are only four matched sets here, so this normal approximation is for illustration and may be unreliable in a real analysis. A suitable small-sample method and enough independent sets are important for dependable inference.

Comparisons, Failure Modes, and Limits

Full matching is useful when the data support different local treated-to-control ratios and the analyst wants to retain a broad set of comparable observations. Fixed-ratio matching instead imposes the same number of controls per treated individual and may discard otherwise useful observations. Direct outcome regression can adjust for measured covariates without forming matched sets, but depends on its outcome-model specification. Weighting without matching is another alternative, but can be especially sensitive to extreme weights. None of these approaches resolves unmeasured confounding merely by changing the calculation.

Common warning signs include poor weighted balance, treated individuals without comparable controls, a few controls carrying a large share of the total weight, and a small number of matched sets driving the outcome contrast. If matching discards treated individuals, the estimand changes to the ATT among those retained unless a broader target is justified. If the outcome is missing differentially or follow-up differs between groups, the simple weighted contrast may not answer the intended comparative-effectiveness question.

  • Mistake: “The matching algorithm formed sets, so baseline confounding is controlled.” Correction: Matching only uses the variables and distance structure specified; inspect weighted balance and explain why measured covariates are adequate, while acknowledging possible unmeasured confounding.
  • Mistake: “Every control should have weight 1 because every person is one observation.” Correction: For this ATT construction, each control in set \(j\) has weight \(n_{Tj}/n_{Cj}\), ensuring the set’s total control weight equals its treated count.
  • Mistake: “A small p-value for a baseline difference proves the match failed,” or “a large p-value proves balance.” Correction: Use weighted effect-size and distributional diagnostics; balance tests are strongly affected by sample size.
Key takeaway: Full matching gives a transparent ATT weighting scheme when every matched set contains treated and control individuals. A defensible result requires more than the weighted contrast: state the retained target population, verify weighted covariate balance and overlap, and use uncertainty calculations suited to the matched design.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A matched set has three treated individuals and six controls. What ATT-oriented weight does each control receive, and what is their combined weight?
  2. Two matched sets have within-set contrasts of 1 day and 4 days. The first contains one treated individual and the second contains three. What is the overall ATT contrast?
  3. After weighting, a continuous covariate has treated mean 40, control mean 43, and reference standard deviation 10. What is its SMD, and what does its sign indicate?
  4. Why does a balanced matched cohort not by itself establish a causal effect?

Solutions. 1. Each control receives \(3/6=0.5\); the six controls together have weight \(6(0.5)=3\). 2. The treated-count weights are \(1/4\) and \(3/4\), so the contrast is \((1/4)(1)+(3/4)(4)=3.25\) days. 3. The SMD is \((40-43)/10=-0.300\); the treated mean is 0.300 reference standard deviations below the weighted control mean. 4. Balance only concerns measured covariates included in the assessment. Causal interpretation also requires assumptions such as conditional exchangeability, positivity, consistency, and no interference; matching cannot demonstrate that unmeasured confounding is absent.

Clinical Biostats