EconBase
← Back to paper

Instrumental Variables with Time-Varying Exposure: New Estimates of Revascularization Effects on Quality of Life

The exact contents of citations.db main_text.text for this paper — one flattened LaTeX string, title through conclusion, appendix excluded, unmodified except for removing email addresses. This is what our citation measures are computed over.

56,865 characters

Instrumental Variables with Time-Varying Exposure: New Estimates of Revascularization Effects on Quality of Life


\setlength{\abovedisplayskip}{8pt}
\setlength{\belowdisplayskip}{8pt}

\begin{titlepage}
\title{\singlespacing Instrumental Variables with Time-Varying Exposure: \\ New Estimates of Revascularization Effects on Quality of Life\thanks{
Fatima Djalalova provided excellent research assistance. This manuscript was prepared using ISCHEMIA Research Materials obtained from the NHLBI Biologic Specimen and Data Repository Information Coordinating Center and does not necessarily reflect the opinions or views of ISCHEMIA or the NHLBI}}
\author{Joshua D. Angrist\thanks{ MIT Department of Economics and NBER. Email: [email removed]} \\   Bruno Ferman\thanks{Sao Paulo School of Economics, FGV. Email: [email removed]}\and Carol Gao\thanks{ MIT Operations Research Center. Email: [email removed]}  \\ Peter  Hull\thanks{ Brown University Department of Economics and NBER. Email: [email removed]} \and Otavio Tecchio\thanks{MIT Department of Economics. Email: [email removed]}
\\ Robert W. Yeh\thanks{Smith Center for Outcomes Research and Division of Cardiovascular Medicine, Department of Medicine, Beth Israel Deaconess Medical Center Email: [email removed]}
}
\date{January 2025}
\maketitle

\begin{abstract}
\begin{singlespace}
\noindent
\par

\noindent The ISCHEMIA Trial randomly assigned patients with ischemic heart disease to an invasive treatment strategy centered on revascularization with a control group assigned non-invasive medical therapy. As is common in such ``strategy trials,'' many participants assigned to treatment remained untreated while many assigned to control crossed over into treatment.  Intention-to-treat (ITT) analyses of strategy trials preserve randomization-based comparisons, but ITT effects are diluted by non-compliance.  Conventional per-protocol analyses that condition on treatment received are likely biased by discarding random assignment. In trials where compliance choices are made shortly after assignment, instrumental variables (IV) methods solve both problems---recovering an undiluted average causal effect of treatment for treated subjects who comply with trial protocol. In ISCHEMIA, however, some controls were revascularized as long as five years after random assignment. This paper extends the IV framework for strategy trials, allowing for such dynamic non-random compliance behavior. IV estimates of long-run revascularization effects on quality of life are markedly larger than previously reported ITT and per-protocol estimates. We also show how to estimate complier characteristics in a dynamic-treatment setting. These estimates reveal increasing selection bias in naive time-varying per-protocol estimates of revascularization effects. Compliers have baseline health similar to that of the study population, while control-group crossovers are far sicker.

\end{singlespace}
\end{abstract}

\onehalfspacing
\setcounter{page}{0}
\thispagestyle{empty}
\end{titlepage}

\section{Introduction}

Patients with chronic coronary artery disease may be treated by revascularization, an invasive strategy involving percutenaous coronary intervention (PCI) or coronary artery bypass grafting (CABG). Both are potentially traumatic and resource-intensive interventions. Alternatively, patients may be treated conservatively with a combination of lifestyle changes and drugs. The question of whether the benefits of revascularization for chronic coronary artery disease outweigh the associated risks and costs remains controversial. The recently concluded ISCHEMIA trial assesses the effects of an invasive strategy centered on revascularization on mortality and various measures of quality of life for patients with stable ischemic heart disease and angina. Intention-to-treat (ITT) analyses that compare ISCHEMIA trial participants based on assigned treatment arm show no statistically significant impact on mortality \citep{spertus2020ckd,spertus2020}. ITT estimates of treatment assignment effects on angina frequency and other Seattle Angina Questionaire (SAQ) domain scores suggest invasive treatment produces modest and fading effects on quality of life and angina frequency.

These findings need not be the last word from ISCHEMIA, however, as analyses of ISCHEMIA trial data face important methodological challenges. Problems arise from the fact that the trial saw high rates of treatment group non-adherence (i.e., some assigned to invasive treatment were not revascularized) as well as control group crossovers (i.e., some assigned to conservative treatment were revascularized).  Assignment to the invasive treatment group boosted revascularization rates substantially, but not deterministically. Conditional on treatment assigned, patients and their doctors chose whether to receive the revascularization treatment.  Such choices, of course, are determined in part by patient frailty and prognosis; they're not independent of potential outcomes.

The timing of control-group crossovers in ISCHEMIA also varied. Some crossovers were revascularized shortly after random assignment, but many waited years for invasive treatment.  Consequently, in each follow-up wave, control-group crossovers include a mix of recently-revascularized patients and patients who were revascularized several years ago. This has implications for the estimation of dynamic treatment effects, such as might be seen if revascularization  benefits fade.\footnote{The recent British Heart Foundation SENIOR-RITA  trial is likewise affected by imperfect compliance and time-varying crossovers \citep{kunadian2024}. In this trial, 6\% of controls, assigned to conservative treatment, crossed over to angiography within 7 days. As of the five-year follow-up, 24\% of controls underwent angiography, and 14\% were revascularized. Among patients assigned to invasive treatment, 90\% underwent angiography and 50\% received revascularization.}

In the face of self-selection into treatment in a randomized clinical trial, analyses of trial data typically feature intention-to-treat effects and ``as-treated'' analyses comparing participants who were treated as planned.  As we've explained elsewhere \citep{pnas2023} and review briefly in the next section, ITT effects are diluted by noncompliance with trial protocols---whether through nonadherence or crossovers. Equally important, estimates that group patients by treatment received effectively discard random assignment, opening the door to selection bias. The problem of selection bias can be especially acute for dynamic effects. Among ISCHEMIA participants, control group crossovers are far sicker than the study population (a result we show below). Since the share of crossovers grows over time, as-treated estimates are increasingly biased downwards over time.

This paper addresses these problems by applying the econometric method of instrumental variables (IV). Specifically, we use random assignment as an instrument for time-varying revascularization exposure, generating new estimates of the quality-of-life gains due to revascularization. On the methodological side, we extend IV methods to capture effects of a time-varying exposure under a set of identifying assumptions well-suited to strategy trials like ISCHEMIA. These results exploit novel \emph{wave ignorability} assumptions that leverage the repeated measurement of outcomes across multiple follow-up waves. Our framework yields average causal effects of revascularization exposure on quality of life for patients revascularized as a result of the trial. We also show how to estimate complier characteristics in a dynamic-treatment setting.

The paper is organized as follows. The next section introduces our econometric framework, beginning with a review of IV methods applied to clinical trials in a static setting. This is followed by our extension to time-varying exposure and dynamic effects. This section also shows how to identify the characteristics of latent groups defined by a dynamic response to treatment assignment.  Section \ref{sec:ests} applies the methods of Section \ref{sec:metrics} to ISCHEMIA.  The resulting estimates show substantially larger and more sustained revascularization effects on SAQ summary score and angina frequency scores than ITT or traditional per-protocol analyses would suggest.  We also document increasing selection bias in as-treated estimates over time, making these increasingly misleading guides to revascularization benefits. Section \ref{sec:concl} summarizes and discusses implications for other trials.

\section{Instrumental Variables for Clinical Exposure}\label{sec:metrics}

\subsection{Econometric Framework}

Our econometric framework is motivated by ISCHEMIA, which randomly assigned trial participants with stable ischemic heart disease to either invasive treatment (angiography as a prelude to intended revascularization) or conservative treatment involving medical therapy alone. Let $Z \in \{0,1\}$ be a variable indicating invasive assignment. Trial participants' quality of life is measured in each of $\bar{w}$ annual follow-up waves; $Y_w$ denotes quality of life measured in wave $w \in \{1,2,3,\hdots, \bar{w}\}$. Treatment is assigned some time in the year ahead of wave 1.\footnote{The trial collected data every few months in the first year and then semi-annually \citep{spertus2020}. Our analysis is limited to yearly follow-ups which were generally more extensive. \cite{spertus2020} uses a Bayesian framework to estimate high-frequency time-varying exposure effects on an ITT basis.}

Post-assignment exposure to revascularization is time-varying: many participants assigned conservative were revascularized later, while some assigned invasive were never revascularized or revascularized years after random assignment. Let $T_w \in \{0,1,2,\hdots, w\}$ denote \emph{revascularization exposure}, defined as the years a participant has lived since revascularization by wave $w$. In any wave, a trial participant is either never revascularized $(T_w=0)$, revascularized shortly after random assignment $(T_w=w)$, revascularized in the year of observation $(T_w=1)$, or revascularized at some other time in between $(1<T_w<w)$.

Potential outcomes capture heterogeneous causal effects of revascularization in each wave. Specifically, let $Y_w(t)$ denote a participant's outcome in wave $w$ when they've lived $T_w = t$ years since revascularization. Once a participant is revascularized, they're revascularized forever.  Formally, this means that for all $w\in\{1,...,\Bar{w}-1\}$, $T_w=t$ implies $T_{w+1}=t+1$. This treatment pattern mirrors a staggered adoption setup in differences-in-differences models.\footnote{See, e.g., \cite{callaway2021, ATHEY202262, borusyak2021revisiting}).} Measured in wave $w$, the incremental individual causal effect of revascularization exposure is $Y_w(t)-Y_w(t-1)$. This notation encompasses, for instance, the difference in quality of life that participants experience when revascularized this year versus last year and the difference in quality of life participants experience when revascularized this year versus never. The cumulative effect of revascularization for $t \le w$ years, observed in wave $w$, is $Y_w(t)-Y_w(0)$. The cumulative effect of revascularization
exposure for a participant revascularized shortly after random assignment is therefore $Y_w(w)-Y_w(0)$.

Just as potential outcomes are indexed against treatment received, it's useful to index potential treatment against treatment assigned.  To that end,
let $T_w(z)$ denote revascularization exposure in wave $w$ when $Z=z$. This is defined for all trial participants, regardless of assignment. The effect of invasive assignment on revascularization exposure in wave $w$ is therefore $T_w(1)-T_w(0)$.

Given an \emph{exclusion restriction}, random assignment makes $Z$ independent of potential outcomes and potential treatments.  This \emph{independence assumption} is formalized as:
\begin{assumption}[Independence]\label{assump:independence}
For each wave $w\in\{1,...,\Bar{w}\}$, the random variables $Y_w(0),...,Y_w(w),\\ T_w(1), T_w(0)$ are jointly independent of $Z$.
\end{assumption}
\noindent The exclusion restriction in this case asserts that assignment to the invasive arm affects outcomes solely by increasing the likelihood of revascularization.\footnote{Exclusion is formalized by double indexing potential outcomes as in \cite{AIR}. Let $Y_w(t,z)$ denote a participant's wave-$w$ potential outcome given $t$ years of exposure and assignment $z$. The exclusion restriction says that $Y_w(t,z)= Y_w(t,z^\prime)$  for each $t \le w$, $z$, and $z^\prime\neq z$.} This assumption is plausible in the ISCHEMIA trial, since randomization to the invasive treatment likely had no direct effects on outcomes. Importantly, the exclusion restriction allows assignment to the invasive treatment to affect outcomes via the \emph{timing} of revascularization.

As in \cite{late94} and \cite{angrist_imbens95}, we assume that invasive assignment either induces revascularization, makes revascularization happen sooner, or leaves revascularization exposure unchanged. This is formalized as:
\begin{assumption}[Monotonicity]\label{assump:exposure_monotonicity}
    For each wave $w\in\{1,...,\Bar{w}\}$, $T_w(1) \ge T_w(0)$ almost surely.
\end{assumption}
\noindent Given monotonicity, trial participants can be classified in each wave as \emph{compliers} (with $T_w(1)>T_w(0)$), \emph{always-takers} (with $T_w(1)=T_w(0)>0$), or \emph{never-takers} (with $T_w(1)=T_w(0)=0$). For always- and never-takers, invasive assignment leaves revascularization exposure unchanged. For compliers, invasive assignment increases exposure.

Invasive assignment is also assumed to generate a \emph{first stage} in wave 1, making the instrument relevant for treatment received shortly after random assignment.  That is:
\begin{assumption}[Relevance]\label{assump:relevance}
    $E[T_1\mid Z=1]>E[T_1\mid Z=0]$.
\end{assumption}
\noindent Under Assumption \ref{assump:independence}, this implies a positive wave-1 first-stage effect on treatment:
$$E[T_1\mid Z=1]-E[T_1\mid Z=0]=E[T_1(1)-T_1(0)] > 0.$$
Because wave-$1$ compliers have $T_1(1)>T_1(0)$, the wave-1 first-stage is the probability of wave-1 compliance under Assumption \ref{assump:exposure_monotonicity}. Moreover, because treatment is irreversible, Assumption \ref{assump:relevance} implies the existence of compliers such that $T_w(1)= w>T_w(0)$ in each wave. In other words, a subset of compliers in each wave experiences the maximum possible $w$ years of revascularization exposure when assigned invasive.

Under Assumptions \ref{assump:independence}-\ref{assump:relevance}, a simple Wald-type IV estimand using wave-1 data identifies the average causal effect of one year of revascularization exposure for wave-1 compliers. Specifically, in wave 1, revascularization exposure, $T_1$, is a Bernoulli treatment that indicates participants revascularized shortly after random assignment. The \cite{late94} local average treatment effect (LATE) theorem applied to wave-1 data therefore implies that:
\begin{equation}
E[Y_1(1)-Y_1(0)\mid T_1(1)>T_1(0)] =
\frac{E[Y_1\mid Z=1]-E[Y_1\mid Z=0]}{E[T_1\mid Z=1]-E[T_1\mid Z=0]}.\label{eq:late}
\end{equation}
The left-hand side of this expression is the average effect of revascularization exposure, $Y_1(1)-Y_1(0)$, on wave-1 compliers. This is obtained by dividing the wave-1 \emph{reduced form} by the wave-1 first stage.  The former is defined as the difference in average quality of life between participants assigned invasive and participants assigned conservative, measured in the first wave. The first stage in the Wald denominator is given by the corresponding difference in wave-1 revascularization rates.

In clinical trials, the reduced form is an ITT effect that compares outcomes by treatment assigned, indicated here by $Z$. The share of the trial population for whom treatment status is changed by random assignment is given by the first stage.  Equation \eqref{eq:late} shows that in a static or short-run analysis of trial data, ITT is diluted relative to the effect of revascularization itself: because the first stage is necessarily between zero and one, the magnitude of LATE exceeds the corresponding ITT effect. Intuitively, because revascularization assignment is assumed to affect outcomes solely by inducing revascularization among compliers, ITT is diluted by averaging in causal effects of zero for always- and never-takers.

IV estimation of per-protocol treatment effects contrast with conventional as-treated analyses of clinical trials that ignore random assignment (see, for instance, the \cite{cabanaPP} analysis of the CABANA trial comparing treatments for atrial fibrillation). Our as-treated estimates come from ordinary least squares (OLS) regressions of $Y_w$ on $T_w$, with a few baseline controls. Non-random treatment take-up likely biases such estimates; in the LATE framework, this bias can be understood as stemming from differences in health between compliers, always-takers, and never-takers \citep{angrist2004,aronow2013}. We substantiate this view with a comparison of group characteristics in Section \ref{sec:popchar} below.

\subsection{Identification and Estimation of Dynamic Causal Effects}\label{sec:dynamicfx}

IV analysis of longer-run effects in the ISCHEMIA trial is complicated by time-varying exposure in a model with heterogeneous potential outcomes. As in \cite{angrist_imbens95} and \cite{rose2024}, the principle complication here arises from the fact that compliance occurs along an extensive margin in which participants who would never have revascularized are induced to revascularize by invasive assignment and an intensive margin in which assignment induces earlier revascularization among participants who would have revascularized anyway. Consequently, complier populations differs for each exposure level and change over time.  At the same time, the availability of repeated follow-ups (waves) gives us a handle on this problem that's not been fully exploited in previous applications of IV to models with dynamic effects.

We disentangle dynamic exposure effects with the aid of an assumption that limits effect heterogeneity across waves:
\begin{assumption}[Wave Ignorability]\label{assump:wave_ignore}
    Average incremental effects for compliers are common across waves: for each pair of waves $w$ and $v\le w$, and for each exposure time $t\le v$,
    \begin{equation}\label{eq:definelambda}
        E[Y_w(t)- Y_w(t-1) \mid T_w(1) \ge t > T_w(0)] = E[Y_{v}(t)- Y_{v}(t-1) \mid T_{v}(1) \ge t > T_{v}(0)] \equiv \lambda_t.
    \end{equation}
\end{assumption}
\noindent Wave ignorability allows incremental exposure effects, $\lambda_t$, to vary freely with exposure time while asserting that incremental effects are independent of the wave in which they're seen.  For instance, in a given wave, the incremental effect of revascularization last year $Y_w(1)-Y_w(0)$ may differ from the incremental effect of one year of early exposure, $Y_w(w)-Y_w(w-1)$. Assumption \ref{assump:wave_ignore} also allows effects to vary for compliers, always-takers, and never-takers. Wave ignorability mirrors restrictions implicit in event study regression models for dynamic causal effects, which typically index causal effects by event time rather than calendar time or treatment cohort (see, e.g., \cite{borusyak2021revisiting}).\footnote{\cite{ferman_tecchio_2023} explore the identifying power of assumptions like wave ignorability for a variety of dynamic treatment effect models.}

The following theorem shows that in the LATE framework outlined in the previous section, wave ignorability is sufficient to identify the set of incremental causal effects $\lambda_t$.
\begin{theorem}[Dynamic Incremental Causal Effects]\label{thm:lambda}
    For each $t\in\{1,...,\bar{w}\}$, let $D_{wt}=\mathbf{1}[T_w\ge t]$ indicate exposure of at least $t$ periods as of wave $w$. Stack the data across waves and consider a two-stage least squares (2SLS) estimator that uses $Z$ and $Z\times\mathbf{1}[w=t]$ with $t\in\{2,...,\bar{w}\}$ to instrument the set of $D_{wt}$ in the linear model
    \begin{equation}\label{eq:lambda_IV}
        Y_w = \mu + \sum_{t=2}^{\bar{w}}\alpha_t\mathbf{1}[w=t] + \sum_{t=1}^{\bar{w}} \lambda_t D_{wt} + \eta_w,
    \end{equation}
   where $\mu$ and $\alpha_t$ are constants and the 2SLS estimand is defined by $E[\eta_w\mid Z]=0$ for all $w$. Given Assumptions \ref{assump:independence}-\ref{assump:wave_ignore}, coefficients on $D_{wt}$ in this 2SLS estimand equal the vector of $\lambda_t$ defined in \eqref{eq:definelambda}.
\end{theorem}

\begin{proof}
See Appendix \ref{appendix:lambda_proof}.
\end{proof}

\noindent Theorem \ref{thm:lambda} is proved using the average causal response (ACR) theorem \citep{angrist_imbens95}. The ACR theorem shows that the wave-specific IV estimand using $Z$ to instrument $T_w$ with outcome $Y_{w}$ can be written as a weighted average of incremental complier revascularization effects for each period up to $w$. Under wave ignorability, these incremental effects are common across $w$ and can therefore be obtained by solving a system of linear equations linking wave-specific reduced forms and first stages. The 2SLS estimator described in the theorem implements this solution.

Beyond incremental effects of an additional year of exposure, it's useful to know the average causal effect of each level of exposure relative to a common reference outcome with no exposure---which, for ISCHEMIA, means no revascularization. Conceptually, this cumulates incremental exposure effects.  Identification of cumulative revascularization exposure effects can be obtained under a somewhat stronger version of Assumption \ref{assump:wave_ignore}:

\begin{assumption}[Strong Wave Ignorability]\label{assump:strong_wave_ignore}
    Average incremental effects for compliers at a given level of exposure are common across waves for compliers at any possible margin: for each $t,v,w,v^\prime,$ and $t^\prime$ such that $t\le v\le w$ and $t\le t^\prime\le v^\prime\le w$,
    \begin{equation}
        E[Y_w(t)- Y_w(t-1) \mid T_w(1) \ge t' > T_w(0)] = E[Y_{v}(t)- Y_{v}(t-1) \mid T_{v'}(1) \ge t' > T_{v'}(0)] \equiv \lambda_t.
    \end{equation}
\end{assumption}

\noindent Assumption \ref{assump:strong_wave_ignore} implies  Assumption \ref{assump:wave_ignore} (to see this, set $v'=v$ and $t'=t$ in Assumption \ref{assump:strong_wave_ignore}). In general, however, Assumption \ref{assump:strong_wave_ignore} adds to Assumption \ref{assump:wave_ignore} the requirement that a given incremental causal effect be the same at each follow-up wave for complier groups defined over every relevant level of exposure (indexed by $t'$), not just the level at which the incremental effect in question is measured (indexed by $t$). Although stronger than Assumption \ref{assump:wave_ignore} in the sense of encompassing restrictions on additional latent types,  Assumption \ref{assump:strong_wave_ignore} likewise amounts to the assertion that the point in time at which effects are measured doesn't matter for incremental causal effects at a given level of exposure.


\begin{comment}
{\color{red}
For each $w$, $t\le w$, and $t^\prime\in\{t,\dots,w\}$,
    \begin{equation}
        E[Y_w(t)- Y_w(t-1) \mid T_w(1) \ge t' > T_w(0)] = E[Y_{t}(t)- Y_{t}(t-1) \mid T_{t'}(1) \ge t' > T_{t'}(0)] \equiv \lambda_t.
    \end{equation}
}
\end{comment}

Under strong wave ignorability, average cumulative causal effects are identified and equal to a 2SLS estimand described in the following theorem.
\begin{theorem}[Dynamic Cumulative Causal Effects]\label{corollary:Lambda}
Suppose Assumptions \ref{assump:independence}-\ref{assump:relevance} and \ref{assump:strong_wave_ignore} hold. Then:
\begin{itemize}
\item[i.] For all $w$ and $t\le w$,
 \begin{equation}
    \label{eq:corrol_def_Lambda}
       E[Y_w(t)- Y_w(0) \mid T_w(1) \ge t > T_w(0)]=\sum_{i=1}^t \lambda_i\equiv \Lambda_t.
    \end{equation}
\item[ii.] For each $t\in\{1,...,\bar{w}\}$, let $R_{wt}=\mathbf{1}[T_w=t]$. Stack the data across waves and consider a 2SLS estimator that uses $Z$ and $Z\times\mathbf{1}[w=t]$ with $t\in\{2,...,\bar{w}\}$ to instrument the set of $R_{wt}$ in the linear model
    \begin{align}\label{eq:Lambda_IV}
    Y_w = \phi + \sum_{t=2}^{\bar{w}}\delta_t\mathbf{1}[w=t] + \sum_{t=1}^{\bar{w}} \Lambda_t R_{wt} + \varepsilon_w,
    \end{align}
    where $\phi$ and $\delta_t$ are constants and the 2SLS estimand is defined by $E[\varepsilon_w\mid Z]=0$ for all $w$. The coefficients on $R_{wt}$ in this 2SLS estimand equal the vector of $\Lambda_t$ defined in \eqref{eq:corrol_def_Lambda}.
\end{itemize}
\end{theorem}
\begin{proof}
See Appendix \ref{appendix:Lambda_proof}.
\end{proof}

The first part of Theorem \ref{corollary:Lambda} says that cumulative revascularization effects are identified and equal the sum of incremental revascularization effects. Strong wave ignorability is key to this result: each $\lambda_t$ captures incremental causal effects for different complier groups. A scenario where trial participants have 0-2 years of exposure highlights the role played by Assumption \ref{assump:strong_wave_ignore} in this context. In this case, the relevant incremental effects at year 2 are $\lambda_1 = E[Y_{2}(1) - Y_{2}(0) \mid T_{2}(1) \ge 1 > T_{2}(0)]$ and $\lambda_{2} = E[Y_{2}(2) - Y_{2}(1) \mid T_{2}(1) \ge 2 > T_{2}(0)]$, which average over different complier groups. The appendix proof shows that under strong wave ignorability, effects for these groups are the same at a given increment: $E[Y_{2}(1) - Y_{2}(0) \mid T_{2}(1) \ge 1 > T_{2}(0)]=E[Y_{2}(1) - Y_{2}(0) \mid T_{2}(1) \ge 2 > T_{2}(0)]$. Hence, $\lambda_1+\lambda_2=\Lambda_2$. The second part of Theorem \ref{corollary:Lambda} shows that this summing can be automated by instrumenting indicators for the full set of exposure levels $R_{wt}$.

A useful benchmark for dynamic effects arises when incremental causal effects equal zero beyond the first year of revascularization exposure.  In this case, revascularization can be seen as a Bernoulli treatment that improves participants quality of life by a fixed and unchanging amount. This scenario justifies a conventional IV estimand that takes a dummy for any exposure as the mediator of causal effects in any wave. In the first wave, this estimand is the same as that for revascularization exposure. The following theorem justifies any-exposure IV for waves $w>1$.

\begin{theorem}[Any-Exposure IV]\label{thm:any-exposure_IV}
Let $V_w=\mathbf{1}[T_w>0]$ and suppose Assumptions \ref{assump:independence}-\ref{assump:relevance} hold. Moreover, for each $w\in\{2,...,\Bar{w}\}$, assume that:
\begin{itemize}
    \item[i.] The instrument is relevant: $E[V_w\mid Z=1]-E[V_w\mid Z=0]\neq0$.
    \item[ii.] Incremental effects are zero beyond one year of exposure:
    \begin{equation}\label{assumption:any}
     E[Y_w(t)-Y_w(t-1)\mid T_w(1)\ge t^\prime>T_w(0)]=0
    \end{equation}
    for all $t\in\{2,...,w\}$ and $t^\prime\in\{1,...,w\}$.
\end{itemize}
Consider the 2SLS estimator that, for a given $w$, uses $Z$ to instrument $V_w$ in
\begin{equation}
\label{eq:any-esxposure_linearmodel}
Y_w= \gamma_w + \tau_w V_w + \xi_w,
\end{equation}
where $\gamma_w$ and $\tau_w$ are constants and the 2SLS estimand is defined by $E[\xi_w\mid Z]=0$. Coefficient $\tau_w$ in this 2SLS estimand equals
\begin{equation}\label{eq:tau}
\tau_w = E[Y_w(T_w(1))-Y_w(0)\mid T_w(1)\ge 1 >T_w(0)].
\end{equation}
\end{theorem}
\begin{proof}
See Appendix \ref{appendix:any-exposure_IV}.
\end{proof}

\noindent The causal effects defined in this theorem (the right-hand side of \eqref{eq:tau}) contrast potential outcomes at any positive exposure with potential outcomes under no exposure. Note that identification of this effect does not require wave ignorability (Assumption \ref{assump:wave_ignore}). Under wave ignorability and the conditions of the theorem, however, $\tau_w$ is constant across waves. Constant $\tau$ can be imposed by stacking specification \eqref{eq:any-esxposure_linearmodel} across waves. Under constant effects, the model is over-identified using $Z$ and $Z \times \mathbf{1}[w = t]$ with $t\in \{2, ..., \Bar{w}\}$ for $V_w$ in the stack. The associated over-identification test statistic is then a test of constant effects.

\subsection{Characterizing Dynamic Compliers and Always-Takers}\label{sec:character}

In a trial like ISCHEMIA with time-varying exposure, compliance is dynamic: assignment may induce treatment in any wave.  An \emph{immediate complier} is treated shortly after random assignment when assigned to be treated, but not otherwise. In later waves, compliance reflects behavior in earlier waves: compliers who would eventually have been revascularized are revascularized sooner when assigned treatment.

We aim to characterize the baseline health of different sorts of compliers, as well as the health of always-takers for whom revascularization timing is unchanged by assignment. In general, treated trial participants are a mix of compliers and always-takers. Complier characteristics are relevant for an assessment of the external validity of IV estimates. In particular, the case for clinical relevance of a set of IV estimates is stronger when the associated compliers have baseline health similar to that of the study population. Moreover, because the treated population includes always-takers, the contrast between complier and always-taker means sheds light on the biases we should expect in as-treated per-protocol estimates.

Identification of latent group characteristics is facilitated by a restriction on compliance timing that we call Immediate Compliers Only (IMCO).  In the ISCHEMIA trial, most participants assigned to invasive treatment were revascularized in the first year following assignment, with few induced to revascularize in later years. Turning this empirical observation into an assumption, we have:

\begin{assumption}[IMCO]\label{assump:imco}
    Assignment either shifts participants into treatment immediately or has no effect: in each wave $w\in\{1,...,\Bar{w}\}$, $P[T_w(1) = w \mid T_w(1) > T_w(0)] = 1$.
\end{assumption}

\noindent Assumption \ref{assump:imco} says that everyone for whom revascularization exposure is changed by random assignment is revascularized in the first wave.\footnote{This can be seen as a flipped version of the \cite{rose2024} restriction requiring extensive margin compliers only (EMCO). With an ordered treatment under EMCO, everyone whose treatment is increased by assignment moves from $T_w(0)=0$ to strictly positive $T_w(1)$.} Always-takers, by contrast, are revascularized in the same wave regardless of assignment, though not necessarily the first.

Covariate means for compliers and always-takers in our dynamic setting are identified by the following theorem. Appendix \ref{appendix:always_taker_mean_general} extends results in the theorem to identify marginal means of potential outcomes.

\begin{theorem}[Dynamic Characterizations]\label{thm:late_takers}
    Let $X$ be a baseline covariate, meaning that it's fixed across waves and unchanged by assignment. Suppose that $(X, T_w(1),T_w(0))$ is independent of $Z$ for all $w\in\{1,...,\Bar{w}\}$ and that Assumptions \ref{assump:exposure_monotonicity}-\ref{assump:relevance} hold. Then: \
\begin{itemize}
    \item[i.] \emph{Immediate complier means} are given by
    \begin{align}
    \label{eq:immediate_complier_X}
    E[X \mid T_1(1)>T_1(0)]
    =\frac{E[\mathbf{1}[T_1=1]\times X\mid Z=1] - E[\mathbf{1}[T_1=1]\times X\mid Z=0]}{E[\mathbf{1}[T_1=1]\mid Z=1] - E[\mathbf{1}[T_1=1]\mid Z=0]}.
    \end{align}
    \item[ii.] \emph{Immediate always-taker means} are given by
    \begin{equation}
    \label{eq:immediate_AT_X}
    E[X \mid T_1(1)=T_1(0)=1]=E[ X \mid T_1=1, Z=0].
    \end{equation}
\end{itemize}
     \noindent Moreover, if Assumption \ref{assump:imco} also holds, complier means are constant across waves and equal to immediate complier means:
    \begin{equation}
    \label{eq:latercomplier_ag_X}
    E[X \mid T_w(1)>T_w(0)]=E[X\mid T_1(1)>T_1(0)]
    \end{equation}
    for all $w\in\{1,...,\Bar{w}\}$. Also, given Assumption \ref{assump:imco},
\begin{itemize}
\setcounter{enumi}{2}
    \item[iii.] For each $w>1$ and $t\in\{0,\dots,w-1\}$ such that
    $$E[\mathbf{1}[T_w=t]\mid Z=1]-E[\mathbf{1}[T_w=t]\mid Z=0]\neq 0,$$ \emph{disaggregated complier means} are given by
    \begin{equation}
    \label{eq:latercomplier_dsg_X}
    E[X \mid T_w(1)=w, T_w(0)=t]=
    \frac{E[\mathbf{1}[T_w=t]\times X\mid Z=1]-E[\mathbf{1}[T_w=t]\times X\mid Z=0]}{E[\mathbf{1}[T_w=t]\mid Z=1]-E[\mathbf{1}[T_w=t]\mid Z=0]}.
    \end{equation}
    \item[iv.] For each $w>1$, \emph{later always-taker means} are given by
    \begin{equation}
    \label{eq:laterATs_X}
    E[X\mid w>T_w(1)=T_w(0)\ge 1]
    = E[X\mid w>T_w\ge1, Z=1].
    \end{equation}
   \emph{Marginal always-takers means}, which average immediate and later always-takers, can be obtained using
    \begin{equation}
    \label{eq:AT_X}
    \begin{split}
   & E[X\mid T_w(1)=T_w(0)\ge 1] \\
    &= \pi_w E[ X \mid T_1=1, Z=0]
    +
    (1-\pi_w)E[X\mid w>T_w\ge1, Z=1],
    \end{split}
    \end{equation}
    where $\pi_w=E[\mathbf{1}[T_w=w]\mid Z=0]/ \{E[\mathbf{1}[T_w=w]\mid Z=0] + E[\mathbf{1}[1\le T_w<w]\mid Z=1]\}$ is the share of immediate always-takers among all always-takers.
\end{itemize}
\end{theorem}
\begin{proof}
    See Appendix  \ref{appendix:always_taker_mean}.
\end{proof}

Parts (i) and (ii) of Theorem \ref{thm:late_takers} follow from \cite{imbensrubin97} and \cite{abadie2003}: identification of average baseline characteristics for immediate compliers and immediate always-takers is analogous to identification of complier characteristics in a static IV setup. These groups are, respectively, participants revascularized immediately when assigned invasive but not otherwise and participants revascularized immediately regardless of treatment assignment.

The rest of the Theorem \ref{thm:late_takers} uses IMCO (Assumption \ref{assump:imco}) to identify complier and always-takers means in a dynamic framework. Specifically, equation \eqref{eq:latercomplier_dsg_X} disaggregates compliers based on the exposure level they  attain when assigned conservative.  These complier means can be computed using an IV estimand that takes $\mathbf{1}[T_w=t]\times X$ as the outcome while instrumenting $\mathbf{1}[T_w=t]$. Equation \eqref{eq:laterATs_X} recovers average baseline characteristics of those revascularized after the first wave regardless of treatment assignment. Finally, equation \eqref{eq:AT_X} combines \eqref{eq:immediate_AT_X} and \eqref{eq:laterATs_X} to recover the average baseline characteristics of the full set of always-takers as of wave $w$.

\section{Revascularization Effects on Quality of Life}\label{sec:ests}

\subsection{The ISCHEMIA Trial}

The ISCHEMIA trial randomized 5,179 patients with moderate to severe cardiac ischemia to one of two care strategies. Patients assigned to the invasive treatment arm were meant to undergo diagnostic coronary angiography and subsequent revascularization when feasible (through PCI or CABG) while also receiving medical therapy. Conservative-arm patients  were meant to receive medical therapy alone, with possible invasive treatment when medical therapy was deemed inadequate \citep{maron2020, spertus2020}.

We analyze trial data from waves 1-5, that is, data collected from one to five years after random assignment. Sample sizes decrease over time since participants who enrolled in the trial later contribute fewer observations ahead of the last follow-up date in December 2018.  Roughly 4300 participants contribute to the analysis of quality of life in wave 1, a number that falls to 670 by wave 5. The proportion with followup data is similar by assignment group.\footnote{Earlier analyses show mortality is unchanged by treatment assignment, suggesting quality of life comparisons are not biased by differential survival.}

In practice, many ISCHEMIA participants received a treatment different from that assigned. Revascularization rates in the control group, reported in the first column of Table \ref{t:sample}, increase from 12\% in wave 1 to 29\% in wave 5. At the same time, only 80\% of those assigned invasive were revascularized initially.  Revascularization rates in the group assigned invasive increased from 80\% in wave 1 to 83\% in wave 5.  Importantly, the difference in the likelihood of revascularization by assignment status, reported in the third column of the table, is well below one. Moreover, this gap falls from a high of 68\% in the first wave to 54\% by wave 5.

Participants assigned to the invasive arm enjoyed a somewhat higher quality of life as a result. ITT effects, reported in columns 5 and 7 of Table \ref{t:sample}, show modest treatment-induced gains in SAQ summary scores ranging from $1.6$ to just under $3$ (these are estimated by differences in raw means, without adjustment for covariates). Angina gains fluctuate around similar magnitudes.  Both of these effects can be compared to means of 86-94, with a standard deviations of $13$-$15$. These positive statistically significant ITT effects are a headline finding from the trial. Yet, as explained above, high rates of nonadherence and crossover likely dilute the ITT estimand relative to the effect of revascularization itself.

The theoretical results in the previous section allow us to obtain average causal effects of revascularization exposure measured in years and ranging from 1-5. The assumptions underpinning our IV analysis are that treatment assigned is independent of potential outcomes and potential exposure and that some participants
assigned invasive are induced to earlier (or any) exposure as a result with no reversals. The exclusion restriction implicit in Assumption \ref{assump:independence} rules out scenarios in which assignment to the invasive strategy improves quality of life with no change in revascularization. It's hard to see a case for random assignment being such a revivifying morale-booster. Assignment to the invasive treatment surely facilitated (rather than inhibited) revascularization for all subjects,  ensuring Assumption \ref{assump:exposure_monotonicity}. Subjects assigned to the invasive arm were also much more likely to be revascularized in the first year, a condition sufficient for Assumption \ref{assump:relevance}.

 Theorem \ref{thm:any-exposure_IV} gives a benchmark for IV analyses of revasculazation exposure: when revascularization matters only to the extent that participants are ever exposed (not how long ago), an any-exposure dummy is the sole mediator of causal effects. Figure \ref{f:any exposure} plots 2SLS estimates generated by this IV strategy; given restriction \eqref{assumption:any}, these are estimates of $\tau_w$ in Theorem \ref{thm:any-exposure_IV}. The figure also plots the corresponding ITT esimates (which are the same regardless of how exposure is coded) along with OLS estimates of any-exposure effects. The latter are computed by regressing SAQ outcomes on an any-exposure dummy ($V_w$) in each wave with the same covariates used for IV.  In this context, OLS is a conventional as-treated per-protocol analysis strategy that discards random assignment.\footnote{Estimates in the figure (and tables below) control for baseline angina frequency scores (as in \cite{spertus2020}) and trial-enrollment region dummies. While not needed for unbiased estimation of revascularization effects under Theorem \ref{thm:any-exposure_IV} since the probability of treatment was constant at one-half, these controls boost the statistical precision of estimated effects.}


2SLS (IV) estimates of any-revascularization effects are larger and more stable over time than ITT effects. For both SAQ outcomes, the latter fall from around four in the first wave to under two in wave 5. This is consistent with the fact that the ITT estimand is diluted by a declining first stage (reported in column 3 of Table \ref{t:sample}). 2SLS estimates of SAQ summary score effects, by contrast, are consistently close to four, while 2SLS estimates of angina gains decline from around 5.5 in the initial follow-up year to a fairly stable estimate of three in later waves. OLS estimates of any revascularization effects decline steadily over time, becoming negative by year five. We show below that the increasing divergence between OLS and 2SLS estimates in this figure is a consequence of increasing selection bias in an as-treated analysis.

\subsection{Dynamic Revascularization Effects}

Estimates in Figure \ref{f:stacked} suggest that, with one exception, the any-exposure model is not too far off. These estimates of dynamic revascularization exposure effects (relative to no exposure) were computed using equation \eqref{eq:Lambda_IV} in Theorem \ref{corollary:Lambda}. Exposure effects on the SAQ summary score are remarkably stable around four, suggesting revascularization improves quality of life by four points immediately and thereafter. 2SLS estimates suggests revascularization improves angina by 5.5 points initially, with long-run gains constant at around 3 points. Dotted lines in the figure mark the mean of exposure effects for waves 1-5 for summary scores and the mean of exposure effects for waves 2-5 for angina. A test of the hypothesis of equal effects across all waves for summary scores yields a p-value of $0.14$. For angina frequency, a test of constant effects across waves 2-5 generates a p-value of $0.86$.\footnote{Testing equality of angina score estimates for all five exposure values generates $\chi^2(4)= 36$, a decisive rejection.}

2SLS estimates of exposure effects are markedly larger---and statistically distinguishable from---the corresponding OLS estimates. 2SLS estimates, OLS estimates, and the difference between them along with associated standard errors appear in Table \ref{t:hausman} for both SAQ outcomes.\footnote{Unlike the specification test comparing IV and OLS in \cite{hausman1978}, these standard errors are computed allowing unrestricted correlation between estimators.}  As can be seen in columns 3 and 6 of the table, differences are statistically significant at conventional levels for most of the summary score estimates and for one of the angina score estimates. For both outcomes, joint tests of equality generate decisive rejections.  It's also noteworthy that all 2SLS estimates in the table are significantly different from zero while the wave-5 OLS estimates are not.

\subsection{Contrasting Treated Populations}\label{sec:popchar}

Comparisons of baseline health measures reveal important differences between participants who do and do not comply with their assigned treatment. Theorem \ref{thm:late_takers} characterizes compliers and always-takers in a dynamic setting by adding IMCO (Assumption \ref{assump:imco}) to our basic framework.

Histograms of revascularization exposure by assignment, plotted in Figure \ref{f:exposure histogram} for each wave, support IMCO:  few participants assigned to revascularize were revascularized after wave 1. In wave 2, for instance, participants assigned invasive had either been revascularized for two years and hence immediately after random assignment (i.e., $T_w(1)=2$) or not at all (i.e., $T_w(1)=0$). Likewise, in wave 5, participants assigned invasive had either been revascularized for five years or not at all. These patterns are consistent with the presumption that $T_w(1)=w$ for all compliers.

Complier baseline SAQ scores are virtually indistinguishable from those of the overall trial population, while always-takers are substantially sicker. This can be seen by comparing columns 1, 2, and 6 of Table \ref{t:compliers}, which report baseline scores for the full sample, for immediate compliers, and for always-takers, respectively. The latter quantity is a marginal mean that averages the baseline health of immediate and later always-takers. Means for immediate always-takers, shown in column 5, are similar to the column 6 average; disaggregated complier means and later always-taker means are omitted since---as indicated by Figure \ref{f:exposure histogram}---few compliers are moved from intermediate levels of exposure to $T_w(1)=w$. The gap between complier and always-taker health is on the order of 10 points; this can be compared with standard deviations of baseline scores of around 19 points when measured in wave 1.

Baseline scores for the population of any-exposure always-takers, moved from no exposure to revascularization in either the first wave or later, are also well below baseline scores for any-exposure compliers. These statistics appear in columns 3 and 8 of Table \ref{t:compliers}.\footnote{For the any-exposure treatment, complier and always-taker means are defined by potential any-exposure dummies $V_w(z) \equiv 1[T_w(z)>0]$ for$z=0,1$. Complier means for participants with $V_w(1)>V_w(0)$ and always-taker means for participants with $V_w(1)=V_w(0)$ are obtained as described for Bernoulli treatments in \cite{pnas2023}.} In contrast with latent groups defined on the basis of years of revascularization exposure, however, the share of any-exposure always-takers among the treated doubles between waves 1 and 5 (as can be seen in column 7 of the table). This reflects the fact that many categorized as any-exposure always-takers (revascularized at some point regardless of assignment) are \emph{exposure} compliers, revascularized immediately instead of later when assigned invasive.

Differences in baseline health between any-exposure compliers and any-exposure always-takers, along with the fact that share treated in the latter category doubles over time, is of substantive importance.  This constellation of health differences and changing always-taker shares explain the decline of OLS (as-treated) estimates in Figure \ref{f:any exposure}.  The composition of the treated population defined conventionally, which in our framework equals the population with any exposure, reflects major \emph{compositional} changes.  The ever-treated population includes an increasing share of sick patients. Notably, column 4 in Table \ref{t:compliers} shows that the share of always-takers among the treated in the exposure analysis is more stable---paralleling relatively stable OLS estimates in Table \ref{t:hausman}.

Finally, it's noteworthy that compliers' baseline health is reasonably constant over time. \cite{spertus2020} document differences in quality of life effects for groups with different baseline health.  Given this correlation between baseline health and revascularization impact, stable complier health supports wave ignorability, which asserts that incremental exposure effects are constant for compliers revascularized on different dates.

\section{Summary and Conclusions}\label{sec:concl}

This paper develops instrumental variables methods for randomized trials with dynamic treatment exposure and outcomes measured in repeated follow-up waves. Variation in exposure time and repeated follow-ups feature in countless randomized trials.  In strategy trials and other pragmatic trials, exposure variation is outside trialists' control once the assignment die has been cast. Initial random assignment notwithstanding, treatment exposure is not randomly assigned.

Dynamic effects in such settings have long posed a challenge for clinical research.  Treatment choices made after random assignment are likely correlated with post-trial outcomes.  Our dynamic potential outcomes framework identifies average causal effects as a function of exposure time without imposing conditional independence of exposure and potential outcomes or constant effects. Identification here exploits a novel identifying assumption that we've called wave ignorability. As in event-study models, wave ignorability says that exposure effects are unrelated to calendar time. Along with assumptions analogous to those in \cite{angrist_imbens95}, wave ignorability delivers identification of incremental and cumulative causal effects of time-varying exposure, both easily computed by 2SLS.

Our framework also enriches the original \cite{AIR} categorization of treated participants as either compliers or always-takers in two ways.  First, we disaggregate complier means into means for trial participants moved into treatment from different points of the untreated exposure distribution. Second, we distinguish means for immediate always-takers who are treated just after random assignment from later always-takers treated regardless of assignment--but not immediately after assignment.

Application of these tools to the ISCHEMIA trial reveal quality of life improvements substantially greater than those that previously-reported ITT and conventional per-protocol estimates would suggest.  2SLS estimates imply revascularization yields sustained gains in SAQ summary scores on the order of four points. ITT analysis, by contrast, generates effects that decline to around two points. For angina, 2SLS yields estimates stable at around three points after an initial bump up to 5.5. For angina too, ITT estimates fall to around two.

No less important, conventional as-treated per-protocol estimates---computed here using OLS regressions on a dummy for any exposure---decline steeply over time. Application of our theoretical results on dynamic latent-group characterization reveal this to be an artifact of the poor health of always-takers. Control group crossovers who are revascularized regardless of assignment are much sicker than compliers.  This and the fact that the size of this group grows over time pulls as-treated estimates down to the point where they are negative and not significantly different from zero by wave 5. In sum, conventional per-protocol estimates give an exceptionally misleading view of the payoff to revascularization from a patient's point of view. This offers an interesting contrast to a recent (static) IV analysis of a mammography screening trial. \cite{kowalski2023behaviour} shows mammography always-takers to be healthier than screening compliers, leading conventional per-protocol estimates to exaggerate the gains from screening.

Our analysis of dynamic complier populations stands as weighs against the concern that the IV complier population is likely to be narrow and unrepresentative of any clinical populations of interest (expressed, for instance, in \cite{hernan_per-protocol_2017}). ISCHEMIA compliers have baseline health much like that of the overall study population. This, of course, is an empirical finding and not a theorem. This finding is not unique to ISCHEMIA, however.  In a static IV analysis of colorectal cancer screening trials,  \cite{pnas2023} show complier characteristics much like those of the population invited for cancer screening.  Our work here illuminates a path likely leading to a revision of previously accepted findings in other trials. We've also provided tools to address the question of whether these new estimates are of clinical relevance.



\singlespacing
\setlength\bibsep{5pt}
\bibliographystyle{bib/aea}
\bibliography{bib/references.bib}

\clearpage


\pagebreak
\onehalfspacing
\section*{} \label{sec:exhibits}
\addcontentsline{toc}{section}{Figures and Tables}
\captionsetup[subfigure]{font=scriptsize,labelfont=scriptsize}

\begin{comment}
\begin{figure}[H]
    \caption{Revascularization Rate by Assignment Status}
    \label{f: revasc_rate}
    \centering
        {\centering{\includegraphics[width=\linewidth]{exhibits/revasc_rate.pdf}}\par}
        \footnotesize \textit{Notes:} This figure plots the the revascularization rate by Assignment Status over follow-up years.\par
\end{figure}
\end{comment}

\begin{figure}[H]
    \caption{Any-Exposure Revascularization Effects by Wave}
    \label{f:any exposure}
    \centering
    \begin{minipage}{16.5cm}
        {\centering{\includegraphics[width=\linewidth]{exhibits/f1-coefplot-by-wave-12-15-24.pdf}}\par}
        \footnotesize \textit{Notes:} 2SLS (IV) estimates use assignment to instrument any-exposure dummies in each wave (year) following assignment; ITT is the corresponding reduced form. OLS is an as-treated effect. All estimates control for baseline angina frequency scores (as in \cite{spertus2020}) and enrollment regions. \par
    \end{minipage}
\end{figure}

\begin{figure}[H]
    \caption{2SLS Estimates of Revascularization Exposure Effects}
    \label{f:stacked}
    \centering
    \begin{minipage}{16.5cm}
        {\centering{\includegraphics[width=\linewidth]{exhibits/kappa_coefplot_noOLS.pdf}}\par}
        \footnotesize \textit{Notes:} This figure plots 2SLS estimates of average causal effects of 1-5 years of revascularization exposure, relative to never have been revascularized. These estimates control for baseline angina frequency scores (as in \cite{spertus2020}) and enrollment regions.\par
    \end{minipage}
\end{figure}

\begin{figure}[H]
    \caption{Histogram of Revascularization Exposure by Treatment Assigned}
    \label{f:exposure histogram}
    \centering
    \begin{minipage}{16.5cm}
        {\centering{\includegraphics[width=\linewidth]{exhibits/f3-hist-by-wave.pdf}\par}}
        \footnotesize \textit{Notes:} This figure plots the histogram of revascularization exposure by follow-up wave.  Exposure $T_w=0$ for participants not revascularized as of wave $w$. \par
    \end{minipage}
\end{figure}


\begin{table}[H]
    \caption{Means, Treatment Rates, and ITT Effects}
    \label{t:sample}
        \centering
    \begin{minipage}{6.5in}
        {\centering\includegraphics[width=16cm]{exhibits/t1-sample.pdf}\par}
        \footnotesize \textit{Notes:} Column 1 reports the revascularization rate by wave for patients randomized to the conservative treatment group. Column 2 reports the corresponding revascularization rate for those assigned invasive. Column 3 reports the first-stage effect of treatment assignment on revascularization by wave: this is column 2 minus column 1. Columns 4 and 7 report sample means. Columns 5 and 7 report ITT estimates for the effect of treatment assigned. SAQ scores are measured at the time of follow up. Patients who were deceased or did not complete a follow-up questionnaire are omitted. Standard deviations appear in square brackets and standard errors in parentheses. \par
    \end{minipage}
\end{table}



\begin{table}[H]
    \caption{2SLS and OLS Estimates of Revascularization Exposure Effects}
    \label{t:hausman}
        \centering
    \begin{minipage}{5.5in}
        {\centering\includegraphics[width=13.5cm]{exhibits/t2-hausman-12-17-24.pdf}\par}
        \footnotesize \textit{Notes:} This table compares 2SLS and OLS estimates of the effects of 1-5 years of revascularization exposure computed using equation \eqref{eq:Lambda_IV}, stacking data from all waves. Columns 3 and 6 report \cite{hausman1978}-type t-tests for the difference between 2SLS and OLS estimates, where standard errors are computed using the variance of the difference of estimates. Chi-square statistics at the bottom of the table test 2SLS-OLS joint equality. This statistic has a $\chi^2(5)$ distribution under the null. Estimates were computed  with controls for baseline angina frequency scores and enrollment regions. Standard errors, clustered on person, are reported in parentheses. P-values for joint tests appear in brackets in the last row.\par
    \end{minipage}
\end{table}


\begin{table}[H]
\caption{Complier and Always-Taker Characteristics}
\label{t:compliers}
\centering
\begin{minipage}{6.2in}
{\centering\includegraphics[width=16cm]{exhibits/t3-compliers-12-17-24.pdf}\par}
\footnotesize
\textit{Notes:} This table reports complier and always-taker means obtained using Theorem \ref{thm:late_takers} for baseline summary scores (Panel A) and baseline angina frequency scores (Panel B). Column 1 shows overall sample means. Column 2 reports immediate complier means computed using equation \eqref{eq:immediate_complier_X} for each wave. Column 3 reports any-exposure complier means. The share of always-takers among the treated appears in column 4; this is the denominator of $\pi_w$ divided by sample share of treated. Column 5 reports immediate always-taker means computed using equation \eqref{eq:immediate_AT_X}; column 6 reports marginal always-taker means computed using equation \eqref{eq:AT_X}. Column 7 reports the share of any-exposure always-takers among the treated and column 8 reports means for this group. Standard deviations appear in square brackets and robust standard errors appear in parentheses. \par
\end{minipage}
\end{table}
\clearpage