Graph-Regularized Synthetic Control for Stable Counterfactual Estimation under Donor Redundancy
Organizations: University of Illinois Urbana-Champaign
Abstract
Synthetic control methods estimate counterfactual outcomes by weighting untreated donor units to match a treated unit before treatment. In recent applications, however, the donor pool can be large, highly correlated, and short in pre-treatment history, yielding many near-equivalent solutions with unstable and difficult-to-interpret donor weights. We propose Graph-Regularized Synthetic Control (Graph-SC), a convex estimator that augments the classical pre-treatment matching objective with a graph Laplacian penalty. The graph encodes similarity among donor trajectories, encouraging coherent weights over structurally related donors while preserving the transparency of synthetic controls. We provide theoretical results characterizing the bias-variance tradeoff induced by graph regularization and establish efficient convergence guarantees with an accelerated proximal gradient method. Across a controlled simulation under donor redundancy, the California Proposition 99 benchmark, and observational counterfactual analyses on the Waymo Open Motion Dataset, Graph-SC improves weight stability, reduces sensitivity to redundant donor pools, and maintains competitive pre- and post-treatment accuracy relative to classical and regularized synthetic control baselines. Our results suggest that graph regularization is most valuable not as a universal accuracy booster, but as a principled mechanism for selecting stable and interpretable counterfactuals from near-equivalent synthetic controls with comparable pre-treatment fit.
Figures & tables
| Method | Pre-RMSPE | Fit gain vs. SC | Post-RMSPE gap size | gap magnitude | Gap retained vs. SC | Active donors | Weight entropy | Criteria met |
| Sparse-only (Ablation) | 0.0058 | 99.6% | 13.0945 | 12.3503 | 63.3% | 20 | 2.5791 | 1/3 |
| Graph-only (Ablation) | 0.0663 | 96.0% | 18.4235 | 17.3714 | 89.0% | 38 | 3.3940 | 2/3 |
| Graph-SC (Ours) | 0.0893 | 94.6% | 19.8918 | 18.6245 | 95.4% | 38 | 3.4112 | 3/3 |
| SC Abadie and Gardeazabal (2003) | 1.6564 | 0.0% | 20.6056 | 19.5136 | 100.0% | 6 | 1.4752 | 1/3 |
| AugSC Ben-Michael et al. (2021) | 0.0001 | 100.0% | 13.1322 | 12.3747 | 63.4% | 6 | 1.4752 | 1/3 |
| Synthetic-DiD Arkhangelsky et al. (2021) | 1.6564 | 0.0% | 18.8482 | 17.6478 | 90.4% | 6 | 1.4752 | 0/3 |
| Outcome | Method | Pre-RMSPE | Fit gain vs. SC | main gap | Gap retained vs. best-fit | placebo | Criteria met |
| headway std | SC Abadie and Gardeazabal (2003) | 1.1738 | 0.0% | 5.3654 | 82.0% | 1.2637 | 1/3 |
| AugSC Ben-Michael et al. (2021) | 0.0148 | 98.7% | 6.0753 | 92.9% | 1.2366 | 2/3 | |
| Synthetic-DiD Arkhangelsky et al. (2021) | 1.1738 | 0.0% | 5.3654 | 82.0% | 1.2637 | 1/3 | |
| Generalized-SC Xu (2017) | 1.6611 | 12.2329 | 187.0% | 2.3937 | 1/3 | ||
| Robust-SC Amjad et al. (2018) | 2.0516 | 34.0756 | 520.9% | 3.1861 | 1/3 | ||
| Penalized-SCM Abadie and L’hour (2021) | 1.2156 | 5.5137 | 84.3% | 1.1103 | 1/3 |
Appendix figures & tables1 asset
Supplementary material from the paper’s appendix.
Appendix
| Estimator / graph setting | ATT MSE | Bias 2 | Variance | Mean Pre-RMSPE | Mean post-RMSPE vs. noiseless CF |
| Sparse-only | 0.1074 | 0.0936 | 0.0276 | 0.0057 | 1.4927 |
| SC | 0.1438 | 0.0449 | 0.1978 | 0.9794 | 1.0407 |
| Graph-only oracle | 0.7024 | 0.6122 | 0.1803 | 0.0677 | 1.6139 |
| Graph-SC oracle | 0.5208 | 0.4747 | 0.0921 | 0.0736 | 1.5328 |
| Graph-SC learned: Euclidean | 0.5360 | 0.5204 | 0.0312 | 0.0615 | 1.5816 |
| Graph-SC learned: Euclidean | 0.5337 | 0.4952 | 0.0770 | 0.0658 | 1.5374 |