---
title: Changes-in-Changes (CiC)
url: https://www.emergentmind.com/topics/changes-in-changes-cic
type: topic
---

# Changes-in-Changes (CiC)

Searching arXiv for recent and foundational CiC-related papers to ground the article.
The Changes-in-Changes (CiC) estimator is a distributional causal-inference framework for two-group, two-period designs. In its classical form, it provides fully nonparametric identification of quantile treatment effects in a two-group, two-period setting under a rank-invariance and support assumption, rather than the usual parallel-trends condition of difference-in-differences [2211.14870]. Across subsequent work, CiC has been generalized to ordered discrete outcomes with underreporting, extreme quantiles, mediation, targeted policies via triple differences, multi-category discrete treatments, endogenous sample selection, attrition, and settings with group-level heterogeneity [2401.00618] [2402.12583] [2411.01617] [2502.08614] [2203.12740] [2307.15313]. The common theme is recovery of a missing counterfactual distribution by monotone quantile mapping rather than by mean-trend extrapolation.

## 1. Classical formulation and identification logic

In the classical CiC model one supposes a latent “no-treatment” outcome
\[
Y_N = h(U,T),
\]
where \(U\) is an unobserved scalar with distribution that may differ by group \(G\in\{0,1\}\) but satisfies \(U\perp T\mid G\), and \(h(\cdot,t)\) is strictly increasing [2211.14870]. If treatment \(I=G\cdot T\), the observed outcome is
\[
Y = Y_N(1-I) + Y_I \cdot I.
\]
Under the support condition that the support of \(U\mid G=1\) is contained in that of \(U\mid G=0\), it can be shown that for all \(y\),
\[
F_{Y_{11}^N}(y) = F_{Y_{10}\Bigl(F_{Y_{00}^{-1}(F_{Y_{01}(y))\Bigr),
\]
and hence the \(q\)-th quantile treatment effect is identified by a composition of observed empirical distribution functions and inverses [2211.14870].

A closely related formulation writes untreated potential outcomes as
\[
Y_{it}(0)\mid (G_i=g)=m\bigl(\mathcal U_{it},t\bigr),
\quad
\mathcal U_{i1}\mid G_i=g \stackrel{d}{=} \mathcal U_{i2}\mid G_i=g,
\]
with \(m(\cdot,t)\) strictly increasing and support\((\mathcal U\mid G_i=1)\subseteq\)support\((\mathcal U\mid G_i=0)\) [2502.08614]. This restates the same identification logic: untreated outcomes evolve through a monotone time transformation of a latent rank whose within-group distribution is stable over time.

The central contrast with difference-in-differences is explicit in the literature. Parallel trends imposes
\[
E[Y^0(t_1)-Y^0(t_0)\mid D=1]=E[Y^0(t_1)-Y^0(t_0)\mid D=0],
\]
which is only a first-moment condition, whereas CiC imposes a higher-order and distributional restriction through monotone quantile-to-quantile evolution [2402.12583]. This suggests that CiC is most naturally interpreted as a model for the full untreated distribution rather than as a correction to mean comparisons.

## 2. Structural assumptions and quantile mapping

The key assumptions recur in multiple variants. One common statement is that there exists a scalar latent \(U\) with common support such that
\[
Y^0(t)=h(U;t),
\]
with \(u\mapsto h(u;t)\) strictly increasing for each \(t\), the distribution of \(U\) the same at \(t_0\) and \(t_1\), and support overlap sufficient for inversion and matching [2402.12583]. In another formulation for binary treatment,
\[
Y_{td}=g_t^{(d)}(U), \qquad U\sim \mathrm{Uniform}(0,1),
\]
with a strictly increasing nonparametric structural model and rank-invariance, also called “copula stability” [2411.01617].

These assumptions imply that untreated evolution can be represented by a quantile–quantile map. In the simplest CiC notation,
\[
T_d(y)=F^{-1}_{Y^0(t_1)\mid D=d}\circ F_{Y^0(t_0)\mid D=d}(y),
\]
and, under the classical no-drift-across-groups implication, the same increasing map transports the period-\(t_0\) untreated distribution into the period-\(t_1\) untreated distribution across groups [2402.12583]. In the debiased semiparametric efficient extension, the corresponding “distributional bridge” assumption states that there exists a function \(\gamma(y,x)\), nondecreasing in \(y\) for each \(x\), such that
\[
\gamma(Y_0^0,X)\mid (X,U)\overset{d}{=}Y_1^0\mid (X,U),
\]
or, for continuous outcomes, that the \(0\to1\) quantile–quantile map \(Q_{Y_1^0\mid X,U}\circ F_{Y_0^0\mid X,U}\) does not depend on \(U\) [2507.07228].

Under these assumptions, the average treatment effect on the treated can be written as
\[
\theta
=
\mathbb E\Bigl[
Y_1 -
Q_{Y_1\mid A=0,X}\!\bigl(F_{Y_0\mid A=0,X}(Y_0)\bigr)
\Bigm|
A=1
\Bigr],
\]
and the counterfactual distribution on the treated satisfies
\[
F_{Y_1^0\mid A=1}(y)
=
\Pr\{\gamma(Y_0,X)<y\mid A=1\}
\]
[2507.07228]. The defining mechanism is therefore composition of an observed pre-treatment distribution for the treated with an observed untreated time map estimated from controls.

## 3. Ordered outcomes, false zeros, and other nonstandard outcomes

A major recent extension develops a Difference-in-Differences model for discrete, ordered outcomes, building upon elements from a continuous Changes-in-Changes model, with a focus on outcomes derived from self-reported survey data eliciting socially undesirable, illegal, or stigmatized behaviors like tax evasion or substance abuse, where too many “false zeros”, or more broadly, underreporting are likely [2401.00618]. The ordered-outcome framework introduces a latent continuous variable observed only through ordered categories \(0,1,\dots,J\):
\[
Y^*_{it}(d)=\eta_{gt}(X_{it})+\lambda_{gt}(X_{it})\cdot \epsilon_{it},
\qquad
Y_{it}(d)=j \ \text{iff}\ \tau_{j-1}<Y^*_{it}(d)\le \tau_j.
\]
For the three-category case \((J=2)\), one may normalize \(\tau_0=0,\tau_1=1\), giving categories \(0\), \(1\), and \(2\) according to whether the latent variable lies below \(0\), in \((0,1]\), or above \(1\) [2401.00618].

Within this threshold-crossing model, the CiC assumption becomes a conditional quantile-mapping restriction on latent untreated outcomes:
\[
F_{Y^*_{10}(0)\mid X}^{-1}\bigl(F_{Y^*_{11}(0)\mid X}(v\mid x)\bigm|x\bigr)
=
F_{Y^*_{00}(0)\mid X}^{-1}\bigl(F_{Y^*_{01}(0)\mid X}(v\mid x)\bigm|x\bigr).
\]
Under smoothness and invertibility, this is equivalent to restrictions on the location-scale parameters \(\eta_{gt}(x)\) and \(\lambda_{gt}(x)\) [2401.00618].

When reported outcomes are contaminated by underreporting, two identification strategies are described. One is partial identification via nonparametric bounds under one-sided misreporting and an upper bound \(\alpha\) on underreporting probability. The other is point identification via a semiparametric consumption–reporting model in which
\[
C(d)=\min\{Y(d),R(d)\},
\]
with a latent reporting bound \(R^*(d)\) and parametric latent indexes
\[
Y^*_{gt}(0)=X'\beta_{gt}+\sigma_{gt}\epsilon_{gt},
\qquad
R^*_{gt}(0)=Z'\gamma_{gt}+\tau_{gt}\nu_{gt},
\]
where \((\epsilon,\nu)\) are jointly real-analytic copula-distributed, independent of \((X,Z)\), and the CiC condition is imposed on each margin separately [2401.00618].

Another nonstandard-outcome direction concerns extreme tails. Existing changes-in-changes estimators are tailored to middle quantiles and do not work well for subpopulations with extreme outcomes, such as infants with extremely low birth weights [2211.14870]. “Extreme Changes in Changes” proposes a new CIC estimator for extreme quantiles by combining the usual CiC structure with extreme-value tail extrapolation under regular variation. The paper recommends use of the extreme CIC estimator for extreme, such as below \(5\%\) and above \(95\%\), quantiles, while the conventional CIC estimator should be used for intermediate quantiles [2211.14870].

## 4. Major extensions of the framework

CiC has been extended well beyond the canonical binary-treatment, two-group design. One line of work generalizes the model to causal mediation with a binary mediator. Under strict monotonicity,
\[
Y_t(d,m)=h(d,m,t,U),
\]
no anticipation, distributional invariance \(U\perp T\mid(D=d,M=m)\), and common support, within-cell quantile–quantile transforms
\[
Q_{d,m}(y)=F^{-1}_{Y_1\mid D=d,M=m}\circ F_{Y_0\mid D=d,M=m}(y)
\]
identify average and quantile direct and indirect effects for various subgroups [1909.04981]. With random assignment and mediator monotonicity, the paper further identifies direct and indirect effects on principal strata such as never-takers, always-takers, and compliers [1909.04981].

A second line extends CiC to targeted policies through a triple-difference analogue. In the triple-changes estimator, there are two states \(s\in\{s_0,s_1\}\), two eligibility groups \(d\in\{d_0,d_1\}\), and two times \(t_0<t_1\). Defining
\[
T_{s,d}(y)=F^{-1}_{Y(t_1)\mid s,d}\circ F_{Y(t_0)\mid s,d}(y),
\]
and the within-state drift
\[
T^*_s(y)=T_{s,d_1}\circ T_{s,d_0}^{-1}(y),
\]
the key new assumption is state-independent drift,
\[
T^*_{s_0}\equiv T^*_{s_1}\equiv T^*.
\]
Under this condition, the missing counterfactual distribution for the treated eligible subgroup in the treated state is point-identified as
\[
F_{Y^0(t_1)\mid s_1,d_1}(y)
=
F_{Y(t_0)\mid s_1,d_1}\circ T^{-1}_{s_0,d_1}\circ T_{s_0,d_0}\circ T^{-1}_{s_1,d_0}(y)
\]
[2402.12583].

A third line generalizes CiC to discrete treatments with more than two categories. Let \(D\in\{0,1,\dots,K\}\) and
\[
Y_t=\sum_{d=0}^K \mathbf 1\{D=d\}Y_{t,d}.
\]
The paper assumes a fully nonparametric rank-invariant representation for each treatment arm,
\[
Y_{t,d}=g_{t,d}(U), \quad U\sim\mathrm{Uniform}(0,1),
\]
and distinguishes weak rank stability, which pertains to untreated potential outcomes within each group, from strong rank stability, which requires invariance of the rank structure of each potential outcome over time even across groups [2411.01617]. Under strong rank stability, the counterfactual distribution for arm \(d\) in group \(d'\) is
\[
F_{Y_{1,d}\mid D=d'}(y)
=
F_{Y_0\mid D=d'}\!\Bigl(
Q_{Y_0\mid D=d}\!\bigl(F_{Y_1\mid D=d}(y)\bigr)
\Bigr),
\]
from which quantile treatment effects, attended quantile effects, ATEs, and ATTs are recovered [2411.01617].

A fourth line introduces group-heterogeneous CiC. In that setting, untreated potential outcomes are
\[
Y^N_{igt}=h(U_{igt},V_{gt},t),
\]
where \(U_{igt}\) is an individual-level unobservable and \(V_{gt}\) is a group-level unobservable, with \(h\) strictly increasing in both arguments [2307.15313]. Identification proceeds by matching both levels of latent heterogeneity through a two-stage quantile system over many control groups. There exists a pair \((\tau'_U,\tau'_V)\) satisfying
\[
Q_{Y_{N0}(\tau'_U)}(\tau'_V)=Q_{Y_{I0}(\tau_U^*)}(\tau_V^*),
\]
and then
\[
Y^N_{I1}(\tau_U^*,\tau_V^*)=Q_{Y_{N1}(\tau'_U)}(\tau'_V)
\]
[2307.15313]. This suggests a distributional matching interpretation in which control subgroups are selected to align both within-group and across-group quantiles with the treated groups’ pre-treatment position.

## 5. Sample selection, attrition, and semiparametric efficiency

A substantial literature studies cases in which outcomes are not always observed. One contribution shows that sample selection arises endogenously when treatment affects whether certain units are observed, and that the conventional ATT estimand may not be well defined, while the DiD estimand cannot be interpreted causally without additional assumptions [2502.08614]. Using principal stratification, it targets treatment effects for the Always-Observed subgroup:
\[
ATT_{AO}
=
E\bigl[Y_{i2}(1)-Y_{i2}(0)\mid G_i=1,V_i=AO\bigr].
\]
Combining CiC counterfactual quantile identification with Lee-style trimming yields sharp lower and upper bounds for \(QTT_{AO}(q)\) and, by integration, for \(ATT_{AO}\) [2502.08614]. The paper also develops a CiC selection model to identify the trimming proportions \(\pi_1,\pi_0\) under selection monotonicity and a latent-variable structure for selection [2502.08614].

A related paper corrects attrition bias using Changes-in-Changes in two-period panels where baseline outcomes are always observed and follow-up outcomes are observed only if \(R=1\) [2203.12740]. The model is
\[
Y_0(0)=\mu_0(0,U_0), \qquad Y_1(d)=\mu_1(d,U_1),
\]
with time-invariance of \(U\) within \((G,R)\) cells and strict monotonicity of the structural functions. Under these assumptions there exists a strictly increasing untreated map
\[
T_0(y)=F^{-1}_{Y_0\mid 0,1}\bigl(F_{Y_1\mid 0,1}(y)\bigr)
\]
such that
\[
F_{Y_1(0)\mid G=g,R=r}(y)=F_{Y_0\mid G=g,R=r}\bigl(T_0(y)\bigr),
\]
and similarly a treated map
\[
T_1(y)=F^{-1}_{Y_0\mid 1,1}\bigl(F_{Y_1\mid 1,1}(y)\bigr)
\]
[2203.12740]. These transformations identify ATT-R, ATE-R, and, under random assignment, ATE for the entire study population. The paper emphasizes that CiC requires no exclusion or “missing-at-random” restriction on response; instead it imposes structural restrictions on the outcome model, while \(R\) may depend arbitrarily on \(G\) and other unobservables \(V\) provided the joint law of \((U_t,V)\) is time-invariant [2203.12740].

A further development addresses semiparametric efficiency and inference with high-dimensional covariates and unmeasured confounding. “Debiased Semiparametric Efficient Changes-in-Changes Estimation” introduces a novel extension of CiC that permits high-dimensional unmeasured confounders and non-monotonic relationships between confounders and outcomes, and constructs efficient estimators that are Neyman orthogonal to infinite-dimensional nuisance parameters [2507.07228]. The efficient influence function for the ATT is
\[
\varphi(W;\theta,\eta)
=
\frac{A}{\pi}\{Y_1-\gamma(Y_0,X)-\theta\}
+
\frac{1-A}{\pi}\int_{Y_1}^{\gamma(Y_0,X)}\nu(u,X)\,du,
\]
with \(\eta=(\gamma,\nu,\pi)\), and it satisfies both \(\mathbb E[\varphi(W;\theta,\eta)]=0\) and the Neyman-orthogonality condition \(\mathbb E[\partial_\eta \varphi(W;\theta,\eta)]=0\) [2507.07228]. Estimation uses \(K\)-fold cross-fitting and arbitrary machine learning methods for nuisance functions, with \(\sqrt n\)-consistency and asymptotic normality under \(L^2\) convergence rates \(o_p(n^{-1/4})\) and standard complexity bounds [2507.07228].

## 6. Estimation, inference, applications, and limitations

Most CiC estimators are nonparametric plug-in procedures. In the basic multi-arm formulation, estimation proceeds by empirical CDFs \(\widehat F_{Y_t\mid D=d}\), empirical quantile functions \(\widehat Q_{Y_t\mid D=d}\), construction of counterfactual distributions by composition, inversion to obtain counterfactual quantiles, and optional rearrangement to enforce monotonicity and prevent crossing quantile curves [2411.01617]. In the triple-changes estimator, empirical maps \(\hat T_{s,d}(y)=\hat F^{-1}_{Y(t_1)\mid s,d}\circ \hat F_{Y(t_0)\mid s,d}(y)\) are composed exactly as in the identification formula, and nonparametric bootstrap within cells is used for confidence intervals [2402.12583]. In the extreme-tail setting, Hill estimators and tail extrapolation are plugged into the CiC identifying formula, yielding asymptotic normality and either plug-in or bootstrap confidence intervals [2211.14870].

Applications in the literature span several substantive domains. For ordered discrete outcomes with false zeros, recreational marijuana legalization for adults in several U.S. states is studied using “Monitoring the Future” repeated cross-sections of U.S. 8th-graders, 2015–2018, with past-30-day marijuana occasions coded \(0=\) none, \(1=\) 1–2 times, \(2=\) 3+ times [2401.00618]. Adding student- and state-level covariates uncovers \(\tau^c(0)\approx -0.02\), \(\tau^c(1)\approx +0.009\), and \(\tau^c(2)\approx +0.011\), all statistically significant at \(5\%\); accounting for underreporting via the semiparametric CiC model further amplifies the estimated effects by roughly \(50\%\) at each level, while there is no statistically significant treatment effect on reporting intentions [2401.00618]. In the extreme-quantile application, the 1993 EITC reform is associated with strictly positive and significant effects across low birth-weight quantiles \(q\in(0,0.20]\), including the most extreme quantiles [2211.14870]. In the mediation application, the JOBS II programme yields \(\hat\theta_1^n\simeq -0.04\) with \(p=0.40\), \(\hat\theta_1^c(0)\simeq +0.06\) with \(p=0.26\), \(\hat\delta_1^c(1)\simeq -0.17\) with \(p=0.04\), and \(\hat\Delta_1^c\simeq -0.11\) with \(p=0.06\) [1909.04981]. In the triple-changes application to Medicaid expansion, classical triple-difference is approximately \(0.170\) preventive-care visits with \(90\%\) CI \([-0.004,0.344]\), while triple-changes is approximately \(0.145\) with \(90\%\) CI \([-0.010,0.331]\) [2402.12583]. In the selection application to a Colombian job-training program, estimated \(\hat\pi_1=0.93\) and \(\hat\pi_0=0.96\), DiD-Lee bounds are \([-0.013,0.322]\), and CiC-Lee bounds are \([-0.109,0.429]\); the naïve complete-case CiC point estimate \(0.158\) has \(95\%\) CI \((0.07,0.246)\), but the selection-corrected bounds include zero [2502.08614]. In an early policy application to COVID-19 deaths, CiC estimates are heterogeneous across countries, with Germany and the United States showing negative effects relative to Sweden, while several other countries exhibit positive estimated effects [2006.12251].

A recurrent misconception is that CiC is merely a nonlinear DiD. The papers instead present it as a model based on monotone latent-rank evolution, quantile mapping, and support conditions, with DiD sometimes appearing only as a special or limiting case [2211.14870] [2502.08614]. Another misconception is that CiC automatically weakens assumptions relative to DiD. Several extensions explicitly require stronger distributional assumptions, support overlap, strict monotonicity, or parametric structure when confronting discrete outcomes, underreporting, or sample selection [2401.00618] [2402.12583] [2502.08614]. The literature therefore treats CiC not as assumption-free, but as a framework that trades mean-trend restrictions for structural restrictions on the evolution of outcome distributions.

Source: https://www.emergentmind.com/topics/changes-in-changes-cic