Standard Distributional Synthetic Controls

• Standard Distributional Synthetic Controls (DSC) is a method that estimates entire counterfactual outcome distributions by optimally weighting donor quantile functions.
• The technique minimizes the integrated squared error between a treated unit's quantile function and a convex combination of donor quantile functions, enabling inference on quantiles, CDFs, and density estimates.
• Despite its computational efficiency, DSC may suffer from geometric failures like vanishing gradients under support mismatch and inability to capture multimodality, leading to potential bias.

Standard Distributional Synthetic Controls (DSC) extend the canonical synthetic control methodology to the estimation of counterfactual outcome distributions, rather than point-valued means or averages. By optimally weighting donor units' distributions to approximate the full pre-treatment law of a treated unit, standard DSC enables causal inference on entire distributional objects such as quantile functions, cumulative distribution functions (CDFs), or density estimates. The method's foundation is the minimization of integrated squared error between the quantile function of the treated unit and a convex combination of donor quantile functions, most commonly operationalized as a convex quadratic program under $L_2$ geometry. This quantile-averaging approach admits elegant computational and theoretical properties but exhibits fundamental geometric fragility, particularly in the face of support mismatch and multimodal target distributions.

1. Formalization of Standard Distributional Synthetic Controls

Suppose there are $J+1$ units observed over a common time horizon, with unit 0 subject to treatment after a specified intervention point, and units $1,\ldots,J$ serving as untreated donor controls. Consider the post-intervention outcome distribution $Y_0$ for the treated unit and $Y_j$ for each donor. The object of interest is the counterfactual ("no-treatment") distribution of the treated unit after intervention.

Let $Q_{Y_j}(\tau)$ represent the $\tau$-th quantile of unit $j$'s outcome distribution for $\tau \in [0,1]$. Standard DSC forms a synthetic control as a weighted mixture of donor quantile functions, with weights $w = (w_1,\ldots,w_J)$ in the simplex: $$\mathcal{W} = \left\{ w \in \mathbb{R}^J \mid w_j \ge 0,\, \sum_{j=1}^J w_j = 1 \right\}.$$ The optimal weights minimize the $L_2$-distance between the treated unit's quantile function and the convex combination of donor quantile functions: $$\widehat{w} = \arg\min_{w \in \mathcal{W}} \int_{0}^{1} \left[ Q_{Y_0}(\tau) - \sum_{j=1}^J w_j Q_{Y_j}(\tau) \right]^2\, d\tau.$$ This is equivalently viewed as finding the element in the convex hull of the control quantile functions that is closest in $L^2([0,1])$ norm to the treated unit's quantile function (Liu, 24 Jan 2026, Gunsilius, 2020, Gunsilius et al., 13 Jan 2025, Zhang et al., 2024).

2. Identification and Theoretical Properties

Point identification of the weight vector $w$ within DSC requires an affine independence condition on the donor quantile functions: $$\sum_{j=1}^J \alpha_j Q_{Y_j}(\tau) = 0\ \forall \tau \in [0,1],\quad \sum_{j=1}^J \alpha_j = 0 \implies \alpha_j = 0 \ \forall j.$$ This ensures that the mapping $w \mapsto \sum_j w_j Q_{Y_j}$ is injective, and thus the $L^2$-objective is strictly convex with a unique minimizer (Liu, 24 Jan 2026). Abstract identification further relies on the assumption that the pre-treatment outcome process for the treated unit lies in the Wasserstein (distributional) convex hull of the donors—a "distance-preserving" causal model where outcome distributions evolve via (scaled) isometries in 2–Wasserstein space (Gunsilius, 2020, Zhang et al., 2024). Under such conditions, the DSC estimator is asymptotically optimal among all convex-aggregation estimators for quantile-function targets and achieves the lowest possible squared prediction error in post-treatment distributional counterfactuals (Zhang et al., 2024).

3. Optimization and Computation

The DSC weight estimation problem is a convex quadratic program, solved efficiently using discretized quantile grids and standard QP routines (Gunsilius et al., 13 Jan 2025). The algorithmic steps are as follows:

1. Estimate empirical quantile functions $Q_{Y_j}(\tau_k)$ for all donors and the treated unit on a fine grid $\{\tau_k\}$.
2. Construct the design matrix and target vector from these quantiles.
3. Solve

$$\min_{w \ge 0,\, \sum w_j = 1} \sum_{k} \left( Q_{Y_0}(\tau_k) - \sum_j w_j Q_{Y_j}(\tau_k) \right)^2.$$

1. Apply the resulting estimated weights to the donor post-treatment distributions to construct the counterfactual distribution as the weighted quantile-function average.

Empirical densities and CDFs can also be combined via appropriate $L^2$ (Cramér–von Mises) or $L^1$ distances, enabling CDF-based or density-based variants of DSC (Kato et al., 2023, Gunsilius et al., 13 Jan 2025).

4. Geometric Pathologies and Limitations

Despite its computational tractability, standard DSC suffers from two salient geometric failures (Liu, 24 Jan 2026):

• Vanishing gradient under support mismatch: If the treated unit's outcome distribution has mass outside the convex hull of the donors' supports, the $L_2$-objective can become locally flat with respect to $w$, precluding informative weight updates. This is evident both in quantile-function space (flat regions in the loss) and in density space, where gradients $\nabla_w \mathcal{L}(w)$ vanish almost everywhere in the case of disjoint supports.
• Structural artifacts in presence of multimodality: Averaging quantile functions of unimodal donors cannot generate genuinely multimodal synthetic distributions, as convex combinations of unimodal quantiles remain unimodal. As a result, true multimodal targets (e.g., bimodal mixtures) are projected into unimodal synthetic distributions, producing collapsed or artifact–laden estimates that cannot reflect the structural complexity of the target.

These failure modes reflect a fundamental restriction of quantile-function averaging: the operation does not generally correspond to averaging outcome distributions unless all donors and target belong to the same location–scale family.

5. Statistical Inference, Bootstrap, and Extensions

Standard DSC admits vector-valued (quantile-wise) treatment effect estimators and supports comprehensive inference procedures:

• Confidence bands via nonparametric bootstrap of quantile differences.
• Permutation ("placebo") tests by re-estimating counterfactuals pretending each donor is treated and benchmarking observed treated–synthetic discrepancies.
• Formal hypothesis testing for goodness-of-fit (Wasserstein distance) and (first- or second-order) stochastic dominance between observed and synthetic distributions (Gunsilius et al., 13 Jan 2025, Gunsilius, 2020).

Extensions include:

• Adapting DSC to panel and cross-sectional settings, including scenarios with as few as two observed periods.
• Accommodating continuous, discrete, or mixed data types, noting convergence rate differences for discrete distributions.
• Allowing for regularization and group–heterogeneous constructions, in which donor groups or times are matched not just by period but also by group–level unobserved heterogeneity (Chen et al., 2023).

6. Comparison to Alternatives and Robustification

A growing literature highlights the geometric brittleness of standard DSC. Liu (2026) demonstrates that replacing the $L_2$ quantile loss with an optimal transport ($W_1$) Wasserstein GAN objective remedies both pathologies: informative gradients are preserved even under support mismatch, and synthetic mixtures can recover genuinely multimodal targets (Liu, 24 Jan 2026). In addition, moment-matching and density-matching SCM variants retain consistency and asymptotic unbiasedness under the mixture model but do not fundamentally address geometric artifacts induced by quantile aggregation (Kato et al., 2023).

Methodological advances such as Functional Synthetic Controls (Okano et al., 12 Jan 2026) provide isometric Hilbert space embeddings to sidestep non–Euclidean difficulties, and distributionally robust synthetic control (DRoSC) estimators (Koo et al., 4 Nov 2025) hedge against instability induced by high donor correlations or post-intervention weight shifts.

7. Practical Guidance and Empirical Performance

Standard DSC achieves substantial gains in estimating the full outcome distribution and quantile treatment effect, particularly when the true counterfactual lies in the distributional convex hull of the donors (Zhang et al., 2024, Gunsilius, 2020). Empirical studies confirm sharp fits to observed distributions with appropriate tuning and demonstrate lower MSEs for treatment effect estimation compared to mean-based SCM. Nevertheless, in applications where support alignment or multimodality are suspected, care must be taken, as standard quantile-averaging can yield biased or uninformative counterfactuals. The disco R package provides reference implementation, visualization, and inference tools for DSC analyses in contemporary empirical work (Gunsilius et al., 13 Jan 2025).