Cited by 497Advanced Causal Inference TechniquesSpatial and Panel Data AnalysisStatistical Methods and Inference

The DOI matches its Crossref record; 79 of the 83 references it cites were found in a registry or an index; none of the 29 sources its claims cite is checked for support yet.

Matrix Completion Methods for Causal Panel Data Models

Open access
Athey, Susan, Bayati, Mohsen, Doudchenko, Nikolay, Imbens, Guido, Khosravi, Khashayar·2018·Journal of the American Statistical Association

Abstract

In this article, we study methods for estimating causal effects in settings with panel data, where some units are exposed to a treatment during some periods and the goal is estimating counterfactual (untreated) outcomes for the treated unit/period combinations. We propose a class of matrix completion estimators that uses the observed elements of the matrix of control outcomes corresponding to untreated unit/periods to impute the “missing” elements of the control outcome matrix, corresponding to treated units/periods. This leads to a matrix that well-approximates the original (incomplete) matrix, but has lower complexity according to the nuclear norm for matrices. We generalize results from the matrix completion literature by allowing the patterns of missing data to have a time series dependency structure that is common in social science applications. We present novel insights concerning the connections between the matrix completion literature, the literature on interactive fixed effects models and the literatures on program evaluation under unconfoundedness and synthetic control methods. We show that all these estimators can be viewed as focusing on the same objective function. They differ solely in the way they deal with identification, in some cases solely through regularization (our proposed nuclear norm matrix completion estimator) and in other cases primarily through imposing hard restrictions (the unconfoundedness and synthetic control approaches). The proposed method outperforms unconfoundedness-based or synthetic control estimators in simulations based on real data.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The paper proposes a matrix-completion estimator with nuclear norm minimization (MC-NNM) for causal panel data, imputing untreated potential outcomes for treated unit-periods. It shows that unconfoundedness, synthetic control, difference-in-differences, and matrix completion can be viewed as minimizing the same objective function but differing in restrictions or regularization, and finds MC-NNM performs well across thin, fat, and square matrices in simulations.

Key findings

  • MC-NNM uses nuclear norm regularization plus unregularized unit and time fixed effects; the authors argue this is computationally attractive and avoids the fixed-rank computational difficulties of interactive fixed effects methods.
  • The paper proves that horizontal regression (unconfoundedness), vertical regression (synthetic control), elastic net vertical regression, difference-in-differences, and nuclear norm matrix completion all minimize the same objective \(Q(Y; R, A, B, \gamma, \delta)\), differing only in hard restrictions or regularization (Theorem 1).
  • In the California smoking data (\(N=38\), \(T=31\)), MC-NNM performs well across \(T_0/T\) ratios; in staggered adoption, VT-EN performs poorly and MC-NNM is the superior approach.
  • In stock market data (daily returns for 2453 stocks over 3082 days; subsamples with \(N \times T = 4900\), from \((490,10)\) to \((10,490)\), 25% missing in a block design), MC-NNM adapts to both thin and fat regimes, does as well as the best estimator in both, and better than both in approximately square settings; it approximates data with rank 4 to 12.

Methods

  • Panel data setting with \(N\) units over \(T\) periods, binary treatment \(W_{it}\), and potential outcomes \(Y_{it}(0)\), \(Y_{it}(1)\); the target is average treatment effect on the treated, imputing missing \(Y(0)\) entries for treated unit-periods.
  • Estimator: \(\min_{L,\Gamma,\Delta} \frac{1}{|O|} \|P_O(Y - L - \Gamma 1_T^\top - 1_N \Delta^\top)\|_F^2 + \lambda \|L\|_*\), where \(O\) indexes observed control outcomes and \(\|\cdot\|_*\) is the nuclear norm; solved via soft-impute with cross-validation for \(\lambda\).
  • Simulations use real data matrices where no unit is treated, artificially designate some units/periods as treated, and compare RMSE of DID, VT-EN, HR-EN, SC-ADH, and MC-NNM.

Limitations

  • The paper states that studying the asymptotic distribution of \(L^* - L\) to construct confidence intervals is beyond scope and is future research; it suggests resampling methods instead.
  • The paper states MC-NNM does not take into account the time series nature of observations and may have autocorrelated errors; it proposes an autoregressive modification but does not implement it in the main simulations.
  • The paper states MC-NNM puts equal weight on all observed elements of \(Y-L\), and proposes propensity-score weighting as an extension; the theoretical bound depends on \(p_c = \min_i \pi_T^{(i)}\), which can be relaxed only under restrictions on \(L^*\).
  • Evident from the text: the simulations use artificial missingness in real datasets rather than actual treatment assignments, so they assess imputation accuracy under constructed treatment patterns, not real causal effects.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.3386/w25132

Cites

79 of 79 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Verification

Claims the paper makes

0 of 29 cited sources checked · 0 supported · 0 flagged

Show all 49 claims

  • It may state two separate claims as one. 18
  • Probably not a substantial claim. 18
  • Borderline: one of the tests for a substantial claim is unclear. 1
34 of 34
  • The authors state this as an assumption.

    The stable-unit-treatment-value assumption holds and there are no dynamic effects: potential outcomes are indexed only by the contemporaneous treatment for that unit, not by past treatments or treatments for other units.

    “Implicit in our set up is that we rule out dynamic effects and make the stable-unit-treatment-value assumption (Rubin (2006); Imbens and Rubin (2015)): the potential outcomes are indexed only by the contemporaneous treatment for that unit and not by past treatments or treatments for other units.”
  • The authors state this as an assumption.

    The complete outcome matrix is the sum of a low-rank matrix L* and mean-zero noise ε.

    “In the absence of covariates we model the N × T matrix of complete outcomes data matrix Y as Y = L * + ε, where E[ε|L * ] = 0 .”
  • The authors state this as an assumption.

    The noise ε is independent of L*, and its elements are σ-sub-Gaussian and independent of each other.

    “Assumption 1. ε is independent of L * , and the elements of ε are σ-sub-Gaussian and independent of each other.”
  • The authors state this as an assumption.

    L* is a deterministic matrix.

    “We also assume that L * is a deterministic matrix.”
  • The authors state this as an assumption.

    Conditional on L*, the adoption dates t_i are independent of each other and of ε.

    “Assumption 2. Conditional on L * , the adoption dates t i are independent of each other and of ε.”
  • The authors state this as an assumption.

    The noise variables {ε_it} are independent of each other and are independent of {t_i}.

    “We also assume that the noise variables {ε it } it∈[N ]×[T ] are independent of each other and are independent of {t i } N i=1 .”
  • The authors state this as an assumption.

    Lmax is a positive constant such that ||L*||max ≤ Lmax.

    “Let L max be a positive constant such that L * max ≤ L max (recall that L * max = max i,t |L * it |).”
  • The authors state this as an assumption.

    The treatment effect is not constant but is such that the matrix Y(1) has a low-rank expectation.

    “Assume the treatment effect is not constant but is such that the matrix Y(1) has a low-rank expectation.”
  • The authors state this as an assumption.

    In §6.1, all entries (i, t) of Y for t ≤ t_i are observed.

    “In §6.1 we assumed that all entries (i, t) of Y for t ≤ t i are observed.”
  • The authors state this as an assumption.

    In the extension, each pre-treatment entry (i, t) can be missing independently with probability p, for p not too large.

    “where each (i, t) ∈ ∪ N i=1 {(i, 1), (i, 2), . . ., (i, t i )} can be in O miss , independently, with probability p for p that is not too large.”
  • The paper relies on this without support.

    The error terms ε_it are mean zero, σ-sub-Gaussian, and independent across i and t.

    “random variables εit -E[ε it ] are centered, τ 2 -sub-Gaussian, and independent of each other (as i or t vary)”
11 of 11

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

  • Rank is incorrectly defined as a limit of Schatten norms

    Minor

    The stated Schatten norm formula does not converge to rank as p approaches zero. For a matrix with two singular values equal to 1, it equals 2^(1/p), which diverges, while the rank is 2.

  • Conditioning data switch from rows to columns

    Minor

    The remark first identifies the conditioning data as the first T -1 rows of Y, then identifies them as the first T -1 columns. These are different parts of the matrix.

  • Nuclear-norm factorization penalizes A twice and omits B

    Minor

    The claimed factorization identity has two terms involving A and none involving B, whereas Theorem 1 penalizes both factors. As written, rescaling A toward zero and B upward while keeping AB fixed makes the right-hand infimum zero even for a matrix with positive nuclear norm.

  • Vertical prediction range includes a pre-treatment period

    Minor

    Treatment is stated to begin in period T 0 + 1, but the vertical-regression description includes T 0 in its range of treated periods.

4 of 4

Citations

Cites

Flagged first, then by citations

Show all 83 references

4 of the 83 references couldn’t be matched to a known paper.

Research works printed but not linked: 4

  1. Alberto Abadie. Using synthetic controls: Feasibility, data requirements, and methodological aspects. Journal of Economic Literature, 2019.
  2. Joshua Angrist and Steve Pischke. Mostly Harmless Econometrics: An Empiricists' Com- panion. Princeton University Press, 2008.
  3. Guido W Imbens and Donald B Rubin. Causal Inference in Statistics, Social, and Biomedical Sciences. Cambridge University Press, 2015.
  4. Donald B Rubin. Matched sampling for causal effects. Cambridge University Press, 2006.

Cited by

Sorted by citations

Items 301–400 of 493

Metadata

Authors
Athey, S., Bayati, M., Doudchenko, N., Imbens, G., & Khosravi, K.
Affiliations
  • Stanford Institute for Economic Policy Research (US)
  • Stanford University (US)
  • Google (United States) (US)
Corresponding author
Imbens, Guido
Venue
Journal of the American Statistical Association
Published
2018-10
Volume
116
Issue
536
Pages
1716-1730
ISSN
0162-1459 (ISSN-L)
OpenAlex
W2766301804
Cited by
497
References
83
Citations by year
2026: 65 · 2025: 71 · 2024: 87 · 2023: 77 · 2022: 50 · 2021: 68 · 2020: 25 · 2019: 40 · 2018: 10
FWCI
29.67
Type
journal article
Language
en
Open access
Locations
  • Journal of the American Statistical Association (publishedVersion)
  • National Bureau of Economic Research (publishedVersion)
  • RePEc: Research Papers in Economics (submittedVersion)
  • arXiv (Cornell University) (submittedVersion)
  • SSRN Electronic Journal (submittedVersion)
  • National Bureau of Economic Research (submittedVersion)
  • arXiv (Cornell University)
Funders
  • National Science Foundation (1554140)
  • Office of Naval Research (N00014-17-1-2131), N00014-17-1, N00014)
  • Division of Civil, Mechanical and Manufacturing Innovation
Topics
  • Physical Sciences › Mathematics › Statistics and Probability › Advanced Causal Inference Techniques
  • Social Sciences › Economics, Econometrics and Finance › Economics and Econometrics › Spatial and Panel Data Analysis
  • Physical Sciences › Mathematics › Statistics and Probability › Statistical Methods and Inference
Keywords
matrix completion, potential outcomes, nuclear norm regularization, missing data patterns, synthetic control method, unconfoundedness, program evaluation, treatment effect estimation

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-09, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-09 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review