From Donor Pool to Counterfactual Prediction
Synthetic control constructs a comparison for one treated region as a weighted average of untreated regions. The calculation is concrete, but its credibility depends on several linked choices: which regions can serve as donors, which pre-intervention features to match, how the weights are estimated, and whether the weighted donor outcomes provide a plausible prediction after treatment begins.
This tutorial works through those choices using hypothetical regional data. The aim is to calculate a synthetic control and its post-intervention predictions, not to treat a good numerical match as automatic proof of a causal effect. The identifying assumptions and donor-pool logic remain essential.
The donor-list index \(i\) is not a region label. Keeping that distinction explicit avoids confusion: the \(i\)th matrix column, weight \(w_i\), and donor outcome \(Y_{d_i,t}\) all refer to the same region, \(d_i\). The region label can be a name or a number; it need not equal \(i\).
Assumptions That Give the Calculation a Causal Role
Before calculating weights, specify the setting in which the synthetic control is intended to represent the treated region’s untreated outcome. Let \(T_0\) be the last pre-intervention period and \(T_0+1\) the first period with treatment. Write \(Y_{r,t}(0)\) for region \(r\)’s outcome at time \(t\) without the intervention, and \(Y_{r,t}(1)\) for its outcome under the intervention. The treated region’s post-intervention \(Y_{1,t}(0)\) is unobserved.
A useful synthetic control analysis requires, at minimum, that the intervention is not received by donor regions and does not affect their outcomes through spillovers; the treated region is not affected before the stated intervention date; and the outcome and predictor definitions are comparable across regions and time. It also requires that the weighted donors provide a credible approximation to the treated region’s untreated post-intervention path. Close pre-intervention predictor and outcome fit supports that approximation but cannot establish it by itself: a new, treated-region-specific shock after \(T_0\), for example, could break the prediction.
The convex-weight restriction means the synthetic control is a weighted average of the selected donors, rather than an extrapolation outside their range. If the treated region is far outside the donors’ pre-intervention predictor combinations, good fit may be impossible. Donor selection and predictor construction therefore come before interpreting any post-intervention calculation.
Worked Example: Screening the Donor Pool
Worked Example: Screening the Donor Pool
Suppose a hypothetical regional program begins in 2021 in treated region A. Potential donors are regions B, C, D, and E. The analysis team plans to use the 2018–2020 period for pre-intervention predictors and annual outcomes. It verifies that B, C, and D did not receive the program, were not exposed to a spillover, and have comparable outcome definitions and complete data for all three pre-intervention years. Region E introduced a similar policy in 2020.
Region E is excluded: its 2020 outcome and predictor history may already reflect an intervention, so it is not an untreated comparison for the planned design. Regions B, C, and D satisfy the stated screening criteria and form the donor list. Thus \(J=3\), with \(d_1=\text{B}\), \(d_2=\text{C}\), and \(d_3=\text{D}\). Every subsequent donor column and weight will follow this ordering.
This screen does not prove that B, C, or D can reproduce A’s untreated path. It rules out a known source of contamination and establishes which regions enter the calculation. Eligibility should be judged using substantive knowledge as well as data availability; adding a region merely because its values improve the pre-period match can undermine the comparison.
Constructing Comparable Predictors
Predictors should be chosen using pre-intervention information and defined in the same way for the treated region and every donor. They may include pre-intervention outcomes, covariates, or summaries of outcome history. Averages over a specified window can reduce noise, while selected pre-intervention annual outcomes can preserve features of the trajectory. The window and definitions should be set so that they do not inadvertently include post-treatment information.
Predictors may be measured in different units. Since scale affects a weighted distance, the scaling rule and predictor-importance matrix should be stated rather than left implicit. In the example, each predictor is divided by a fixed, prespecified scale: ten outcome-index points for the pre-period mean, two percentage points for unemployment, and ten percentage points for manufacturing share. These are unit conversions, not standardizations estimated from the observed sample.
Worked Example: Building the Predictor Vectors
The hypothetical outcome is a regional index measured in points. For 2018–2020, region A’s annual outcomes are 48, 49, and 50, so its pre-period mean is
The comparable donor outcomes are B: 38, 40, and 42; C: 58, 60, and 62; and D: 53, 55, and 57. Their pre-period means are 40, 60, and 55 points, respectively. The other pre-intervention predictors are unemployment rate and manufacturing share. The raw values and the resulting rescaled predictor vectors are:
| Region | Mean outcome, points | Unemployment, % | Manufacturing share, % | Rescaled predictor vector |
|---|---|---|---|---|
| A (treated) | 49 | 7.4 | 27 | (4.9, 3.7, 2.7)' |
| B (donor 1) | 40 | 6.0 | 20 | (4.0, 3.0, 2.0)' |
| C (donor 2) | 60 | 10.0 | 40 | (6.0, 5.0, 4.0)' |
| D (donor 3) | 55 | 7.0 | 25 | (5.5, 3.5, 2.5)' |
For instance, A’s rescaled unemployment predictor is \(7.4/2=3.7\), and its rescaled manufacturing predictor is \(27/10=2.7\). In the donor-list ordering B, C, D, the treated vector and donor matrix are
The columns of \(X_0\) correspond, in order, to B, C, and D. This ordering will also be used for the weights and post-period outcomes.
Estimating the Weights
The standard synthetic control weight problem chooses a convex combination of donors to match the treated region’s predictors. For a specified positive semidefinite matrix \(V\), the optimization is
The matrix \(V\) determines the relative importance of predictor mismatches. It can, for example, place greater emphasis on selected predictors or account for their scales. Here the predictors have already been rescaled by the stated fixed divisors, so use \(V=I\), the identity matrix. This gives equal weight to squared mismatches on the rescaled predictor scales. In an applied analysis, the predictor scales and importance choices should be justified and should not be selected simply to manufacture a preferred post-intervention result.
Worked Example: Solving for the Donor Weights
With \(V=I\), a zero predictor mismatch is feasible if some convex combination of B, C, and D equals A’s rescaled predictor vector. Let \(w_1\), \(w_2\), and \(w_3\) be the weights for B, C, and D, respectively. The sum constraint gives \(w_3=1-w_1-w_2\). Substituting this into the predictor-matching equations and subtracting D’s vector yields
The last two equations are identical because the unemployment and manufacturing differences across these particular regions have the same pattern. From the first equation, \(-3w_1+w_2=-1.2\), so \(w_2=3w_1-1.2\). Substitute into the second:
Then \(w_2=3(0.5)-1.2=0.3\), and \(w_3=1-0.5-0.3=0.2\). All three weights are nonnegative and sum to one. Checking each predictor directly,
Thus \(X_0W=X_1\), and the objective equals zero. Since a squared-distance objective cannot be less than zero, these weights attain a global minimum. The equations above also determine \(w_1\) and \(w_2\) uniquely, and the sum constraint then determines \(w_3\), so this exact match has a unique feasible weight vector in this example.
Predicting the Treated Region’s Untreated Outcomes
Once weights are fixed using pre-intervention information, the synthetic control’s outcome in period \(t\) is the same weighted average of the donor outcomes in that period. With the donor-list notation used throughout, the prediction for the treated region’s outcome without the intervention is
This is a counterfactual prediction, not an observed value. Its causal interpretation relies on the assumptions described earlier: in particular, that the selected weighted donors continue to approximate the treated region’s untreated outcome after the intervention begins.
Worked Example: Calculating Post-Intervention Predictions
Suppose the hypothetical donor outcomes are as follows. The intervention begins in 2021, after the three pre-intervention years. The treated region’s observed post-intervention outcomes are included for context, but they are not used to estimate the weights.
| Year | A, treated region (observed) | B, donor 1 | C, donor 2 | D, donor 3 |
|---|---|---|---|---|
| 2021 | 55 | 43 | 63 | 58 |
| 2022 | 58 | 45 | 65 | 60 |
For 2021, apply the estimated weights to the donor outcomes in the same order, B, C, D:
For 2022, the calculation is
Under the synthetic control assumptions, 52.0 and 54.0 points are predictions of A’s outcomes in those years without the program. A’s observed outcomes, 55 and 58 points, are not those counterfactuals; comparing observed and predicted outcomes is a subsequent step in an analysis.
What the Worked Calculation Does—and Does Not—Establish
The example has exact matching on three constructed predictors, and the estimated weights are transparent. Those facts are useful, but they do not establish that the post-intervention predictions are correct. The result depends on the substantive validity of the donors, the timing of treatment, the absence of donor contamination and spillovers, the pre-period information selected, and the plausibility that the weighted donors track A’s untreated path after 2020.
A further limitation is that a synthetic control is constrained to the convex hull of the selected donors. If no nonnegative weights summing to one can reproduce the treated predictors, the optimization returns the closest feasible match rather than an exact one. In that case the residual mismatch should be reported and considered when judging the prediction. A numerical optimizer’s output is not itself evidence that the design is credible.
Exclude regions with treatment, spillovers, incompatible outcomes, or other substantive reasons they cannot represent an untreated comparison.
Use consistently defined pre-intervention information, state the windows and transformations, and consider how scale affects mismatch.
Minimize predictor mismatch subject to nonnegative weights that sum to one, keeping donor columns, weights, and region outcomes in one consistent order.
Apply the fixed weights to donor outcomes, then interpret the resulting counterfactual prediction in light of the design assumptions.
Check Your Understanding
Answer each question using the definitions and calculations in this tutorial.
- If donors are listed as \(d_1=\text{B}\), \(d_2=\text{C}\), and \(d_3=\text{D}\), which region’s outcome is multiplied by \(w_2\) in the counterfactual prediction?
- Why is region E excluded from the hypothetical donor pool, and what general problem does this exclusion address?
- Using the example weights, verify the predicted 2021 outcome by calculating the weighted average of B, C, and D’s outcomes.
- What does the nonnegative-weight and sum-to-one restriction imply about the range of predictor combinations the synthetic control can represent?
- Why does an exact pre-intervention predictor match not, by itself, guarantee a valid post-intervention counterfactual?