Tutorials › Causal Inference › Solving the Synthetic Control Weight Problem

Synthetic Control and Interrupted Series · Tutorial 659 of 972

Solving the Synthetic Control Weight Problem

Learn to formulate, solve, and interpret the constrained least-squares problem that chooses synthetic control donor weights.

Advanced 11 min read

What You'll Learn

  • Express predictor matching as a quadratic objective in donor weights
  • Explain why nonnegative weights summing to one restrict the synthetic control to a convex combination of donors
  • Derive the quadratic-program form and establish its convexity
  • Use first-order conditions to understand which donors receive positive weights
  • Work through interior, boundary, and predictor-weighted solutions
  • Interpret predictor balance without treating it as proof of causal identification

From Predictor Matching to Donor Weights

Once a donor pool and a set of pre-intervention predictors have been chosen, the synthetic control is defined by weights on the donor units. The goal of the weight problem is to make a weighted combination of donor predictors resemble the treated unit’s predictors. The constraints matter: they determine which combinations are allowed and therefore what “best match” means.

Let there be \(J\) predictors and \(K\) donor units. Write \(X_1\) for the \(J \times 1\) vector of predictors for the treated unit, and \(X_0\) for the \(J \times K\) matrix whose \(k\)th column contains the predictors for donor \(k\). The vector \(W=(w_1,\ldots,w_K)'\) contains the donor weights. For a fixed \(J \times J\) predictor-weighting matrix \(V\), the discrepancy between the treated predictors and the synthetic predictors is measured by a weighted quadratic distance.

Definition: The synthetic control weight problem is to choose \(W\) to minimize \((X_1-X_0W)'V(X_1-X_0W)\), subject to \(w_k \geq 0\) for every donor \(k\) and \(\sum_{k=1}^K w_k=1\). Here \(V\) is symmetric positive semidefinite: \(V=V'\) and \(a'Va \geq 0\) for every \(J \times 1\) vector \(a\).

The synthetic predictor vector is \(X_0W\). The constraints make it a convex combination of the donor predictor vectors: no donor receives a negative weight, and the weights add to one. Thus, the synthetic predictors lie in the convex hull of the donor predictors. This is a substantive restriction, not just a computational convenience. It prevents the solution from extrapolating by subtracting one donor’s predictors from another’s.

The matrix \(V\) defines how discrepancies across predictors contribute to the objective. For example, a diagonal entry may give greater importance to one predictor than another. In this tutorial \(V\) is treated as fixed; the central task is to solve the weight problem for that chosen matrix. The objective measures predictor mismatch under that metric. It does not guarantee exact balance, since the treated unit may lie outside the donors’ convex hull or the selected predictors may not all be balanced by one set of weights.

Quadratic-Program Form and Convexity

The original objective can be written as a quadratic function of \(W\). Define \(Q=X_0'VX_0\), \(b=X_0'VX_1\), and \(c=X_1'VX_1\). Because \(V\) is symmetric, the two cross-products in the expansion are equal:

$$ \begin{aligned} (X_1-X_0W)'V(X_1-X_0W) &=X_1'VX_1-X_1'VX_0W-W'X_0'VX_1+W'X_0'VX_0W\\ &=c-2b'W+W'QW. \end{aligned} $$

Therefore the weight problem is a quadratic program: minimize \(W'QW-2b'W+c\) over the simplex of nonnegative weights summing to one. The constant \(c\) does not depend on \(W\), so it can be omitted when computing the minimizer. It is retained in the displayed expansion because it is part of the original mismatch.

The objective is convex. For any \(K \times 1\) vector \(u\),

$$ u'Qu=u'X_0'VX_0u=(X_0u)'V(X_0u)\geq 0. $$

This establishes that \(Q\) is symmetric positive semidefinite: symmetry follows from \(V=V'\), and the displayed inequality establishes nonnegative quadratic forms. Hence the quadratic objective is convex, and any feasible point satisfying the minimization conditions is a global minimizer. The simplex is also closed and bounded, so a minimizer exists. It need not be unique: different weight vectors can produce the same predictor combination, or the objective can be flat along a feasible direction.

Worked Example: Two Donors and an Interior Solution

Suppose the treated unit has two predictors \(X_1=(5,6)'\), and two donors have predictor vectors \(X_A=(2,7)'\) and \(X_B=(8,3)'\). Let \(V\) be the identity matrix, so the objective is the sum of squared predictor discrepancies. Denote donor B’s weight by \(t\). The constraints require donor A’s weight to be \(1-t\) and \(0 \leq t \leq 1\).

The synthetic predictors are \((1-t)(2,7)'+t(8,3)'=(2+6t,7-4t)'\). Thus the objective as a function of \(t\) is

$$ \begin{aligned} f(t)&=(5-(2+6t))^2+(6-(7-4t))^2\\ &=(3-6t)^2+(-1+4t)^2\\ &=10-44t+52t^2. \end{aligned} $$

The derivative is \(f'(t)=-44+104t\), so the unconstrained minimum occurs at \(t=44/104=11/26\), which is between zero and one and is therefore feasible. The minimizing weights are \(w_A=15/26\) and \(w_B=11/26\). The synthetic predictors are \((59/13,69/13)'\), approximately \((4.538,5.308)'\). The residuals, treated minus synthetic, are \((6/13,9/13)'\), and the minimized mismatch is

$$ \left(\frac{6}{13}\right)^2+\left(\frac{9}{13}\right)^2 =\frac{117}{169}\approx 0.6923. $$

Both weights are positive and sum to one. The solution balances the two predictors as well as this constrained, equally weighted squared-distance criterion allows; it does not make either predictor match exactly.

How the Constraints Shape the Solution

For a general number of donors, the nonnegativity constraints mean the optimum can lie on a face of the simplex. A donor with weight zero is excluded from the synthetic combination. The first-order conditions clarify how the active donors and excluded donors differ. For the objective \(f(W)=c-2b'W+W'QW\), its gradient is \(2(QW-b)\). Introduce a multiplier \(\nu\) for the constraint \(\mathbf{1}'W=1\), where \(\mathbf{1}\) is a \(K \times 1\) vector of ones, and nonnegative multipliers \(\mu_k\) for \(w_k \geq 0\). The Karush–Kuhn–Tucker conditions are:

$$ 2(QW-b)+\nu\mathbf{1}-\mu=0,\qquad \mathbf{1}'W=1,\qquad W\geq 0,\qquad \mu\geq 0,\qquad \mu_k w_k=0\ \text{for each }k. $$

For any donor with positive weight, complementary slackness gives \(\mu_k=0\), so its gradient component is balanced by the common equality-constraint multiplier. For a zero-weight donor, the corresponding multiplier can be positive; the condition prevents a feasible move that adds that donor from improving the objective. Since the objective is convex and the constraints are linear, these conditions are sufficient for a global optimum as well as necessary.

1
Fix the predictors and \(V\).
Specify the predictor matrix \(X_0\), treated vector \(X_1\), and symmetric positive semidefinite weighting matrix \(V\).
2
Form the quadratic objective.
Compute \(Q=X_0'VX_0\) and \(b=X_0'VX_1\); minimize \(W'QW-2b'W+c\), subject to the simplex constraints.
3
Check the active constraints.
Verify that the weights sum to one and are nonnegative. For a proposed solution, check the first-order conditions for its positive-weight donors and its zero-weight donors.
4
Interpret the fitted predictors.
Calculate \(X_0W\) and compare it with \(X_1\) predictor by predictor. The objective gives a weighted summary of these discrepancies, not a guarantee of exact balance.

Worked Example: A Boundary Solution

Suppose there are three donors with two predictors each: \(X_A=(0,0)'\), \(X_B=(2,0)'\), and \(X_C=(0,2)'\). The treated predictor vector is \(X_1=(2.5,1.5)'\), and \(V\) is the identity matrix. Write the weights as \(w_A,w_B,w_C\). The synthetic predictors are \((2w_B,2w_C)'\), with \(w_A+w_B+w_C=1\) and all weights nonnegative.

The treated vector lies outside the donors’ convex hull because its coordinates sum to \(4\), while every feasible synthetic vector has coordinates summing to at most \(2\). To find the closest feasible point, first consider whether an interior solution is possible. If the sum constraint were temporarily ignored, matching each coordinate separately would require \(w_B=2.5/2=1.25\) and \(w_C=1.5/2=0.75\). These values already sum to \(2\), leaving no nonnegative weight for donor A, so that unconstrained match is infeasible.

On the edge with \(w_A=0\), write \(w_B=t\) and \(w_C=1-t\). The objective is

$$ \begin{aligned} f(t)&=(2.5-2t)^2+(1.5-2(1-t))^2\\ &=(2.5-2t)^2+(-0.5+2t)^2\\ &=6.5-12t+8t^2. \end{aligned} $$

Its derivative is \(-12+16t\), so the edge minimum is \(t=12/16=0.75\). This is feasible, giving \(w_A=0,w_B=0.75,w_C=0.25\). The synthetic predictors are \((1.5,0.5)'\), and the residual is \((1,1)'\), so the squared mismatch is \(1^2+1^2=2\).

To check that this edge solution is also optimal over the full simplex, the gradient of the objective with respect to the synthetic predictor vector is \(2(X_0W-X_1)=(-2,-2)'\). The directional derivatives for moving weight from either positive-weight donor B or C to donor A are both \(4>0\): donor A’s predictor vector is \((0,0)'\), so such a move increases the objective. The edge minimum is therefore not improved by entering the interior. The zero weight is a consequence of the nonnegativity and sum constraints, rather than an arbitrary decision to discard donor A.

Worked Example: Changing Predictor Importance Changes the Weights

Suppose the treated unit’s predictors are \(X_1=(3,0)'\), and two donors have \(X_A=(0,0)'\) and \(X_B=(4,2)'\). Let \(t\) be donor B’s weight, so donor A’s weight is \(1-t\) and \(0 \leq t \leq 1\). Now use the fixed weighting matrix \(V=\begin{pmatrix}1&0\\0&4\end{pmatrix}\). This gives the second predictor four times the contribution of the first to the squared discrepancy.

The synthetic predictor vector is \((4t,2t)'\). Substitution into the weighted objective gives

$$ \begin{aligned} f(t)&=(3-4t)^2+4(0-2t)^2\\ &=9-24t+32t^2. \end{aligned} $$

The derivative is \(-24+64t\), which equals zero at \(t=24/64=0.375\). This lies in the feasible interval, so the solution is \(w_A=0.625,w_B=0.375\). The synthetic predictors are \((1.5,0.75)'\), with residual \((1.5,-0.75)'\). The weighted mismatch is

$$ (1.5)^2+4(-0.75)^2=2.25+2.25=4.5. $$

For comparison, if \(V\) were the identity matrix, the objective would be \((3-4t)^2+(0-2t)^2=9-24t+20t^2\), minimized at \(t=24/40=0.600\). Giving the second predictor more importance reduces donor B’s weight because that donor differs from the treated unit on that predictor. The weights are selected to balance predictors according to the specified metric, not to minimize every predictor discrepancy separately.

What the Weights Do—and Do Not—Establish

The optimization makes the predictor-matching rule explicit. If a predictor has a large discrepancy and receives substantial weight through \(V\), the solution has an incentive to reduce that discrepancy, subject to the available donors and the simplex constraints. But the solution cannot overcome a lack of suitable donors: it may leave appreciable residual imbalance when the treated predictors are outside the donor convex hull or when matching one predictor worsens another.

A small minimized objective means close fit under the chosen predictors and \(V\). It is not, by itself, evidence that the donor combination reproduces the treated unit’s untreated post-intervention trajectory. Causal interpretation still depends on the design’s identification assumptions and on defensible choices of donor pool, predictors, and pre-intervention window. The weights solve a specified optimization problem; they do not validate those choices.

  • Mistake: allowing negative weights to improve the fit. An unconstrained least-squares solution may subtract one donor from another. That is not the synthetic control problem defined here, because it permits extrapolation beyond convex combinations of donors.
  • Mistake: assuming every selected predictor must match exactly. The minimizer balances the aggregate weighted discrepancy. Exact matching is possible only when the treated predictor vector can be represented by a feasible donor combination.
  • Mistake: reading a positive donor weight as a causal claim. A positive weight identifies a donor’s contribution to the constructed predictor match. Causal interpretation requires more than the optimization result.
Key takeaway: With \(V\) fixed and symmetric positive semidefinite, synthetic control weights minimize a convex quadratic mismatch over the simplex. The resulting weights form a convex combination of donors and balance the selected pre-intervention predictors according to the specified metric, within the limits of the donor pool.

Check Your Understanding

Try these questions without looking back at the worked calculations.

  1. Why does requiring \(w_k \geq 0\) and \(\sum_k w_k=1\) restrict synthetic predictors to the donors’ convex hull?
  2. In the quadratic expansion, where is symmetry of \(V\) needed to combine the two cross-product terms?
  3. What property of \(V\) ensures \(Q=X_0'VX_0\) is positive semidefinite?
  4. If a donor has zero weight at the optimum, what do the constraints and first-order conditions imply about adding it through a feasible change in weights?
  5. Can a small minimized predictor discrepancy alone establish that a synthetic control identifies a causal effect? Explain.
Clinical Biostats