Table of Contents
Fetching ...

Randomization Inference When N Equals One

Tengyuan Liang, Benjamin Recht

TL;DR

The paper addresses inference for N-of-1 trials where treatment effects propagate over time, violating SUTVA due to temporal interference. It develops a linear time-invariant (LTI) framework with impulse response g and uses a generalized estimator $\widehat{\tau}^{\mathcal{L}}(q)$ based on random treatment schedules to recover linear estimands $\tau^{\mathcal{L}}(q) = \langle D_T q, g \rangle$, along with a circular-convolution analogue $\widehat{\tau}(q)$. The authors prove asymptotic normality under decay conditions on $g$ and $q$, derive higher-moment formulas, and propose a plug-in variance estimator for valid confidence intervals, showing that circular and linear models are asymptotically equivalent. They also demonstrate a rapid interleaving design that improves signal-to-noise and enables estimation of a broad class of estimands beyond standard N-of-1 designs. The work connects dynamic causal inference, system identification, and N-of-1 experimentation, and outlines directions for optimal design and extensions to nonstationary or nonlinear dynamics.$

Abstract

N-of-1 experiments, where a unit serves as its own control and treatment in different time windows, have been used in certain medical contexts for decades. However, due to effects that accumulate over long time windows and interventions that have complex evolution, a lack of robust inference tools has limited the widespread applicability of such N-of-1 designs. This work combines techniques from experiment design in causal inference and system identification from control theory to provide such an inference framework. We derive a model of the dynamic interference effect that arises in linear time-invariant dynamical systems. We show that a family of causal estimands analogous to those studied in potential outcomes are estimable via a standard estimator derived from the method of moments. We derive formulae for higher moments of this estimator and describe conditions under which N-of-1 designs may provide faster ways to estimate the effects of interventions in dynamical systems. We also provide conditions under which our estimator is asymptotically normal and derive valid confidence intervals for this setting.

Randomization Inference When N Equals One

TL;DR

The paper addresses inference for N-of-1 trials where treatment effects propagate over time, violating SUTVA due to temporal interference. It develops a linear time-invariant (LTI) framework with impulse response g and uses a generalized estimator based on random treatment schedules to recover linear estimands , along with a circular-convolution analogue . The authors prove asymptotic normality under decay conditions on and , derive higher-moment formulas, and propose a plug-in variance estimator for valid confidence intervals, showing that circular and linear models are asymptotically equivalent. They also demonstrate a rapid interleaving design that improves signal-to-noise and enables estimation of a broad class of estimands beyond standard N-of-1 designs. The work connects dynamic causal inference, system identification, and N-of-1 experimentation, and outlines directions for optimal design and extensions to nonstationary or nonlinear dynamics.$

Abstract

N-of-1 experiments, where a unit serves as its own control and treatment in different time windows, have been used in certain medical contexts for decades. However, due to effects that accumulate over long time windows and interventions that have complex evolution, a lack of robust inference tools has limited the widespread applicability of such N-of-1 designs. This work combines techniques from experiment design in causal inference and system identification from control theory to provide such an inference framework. We derive a model of the dynamic interference effect that arises in linear time-invariant dynamical systems. We show that a family of causal estimands analogous to those studied in potential outcomes are estimable via a standard estimator derived from the method of moments. We derive formulae for higher moments of this estimator and describe conditions under which N-of-1 designs may provide faster ways to estimate the effects of interventions in dynamical systems. We also provide conditions under which our estimator is asymptotically normal and derive valid confidence intervals for this setting.
Paper Structure (33 sections, 12 theorems, 141 equations, 3 figures, 1 table)

This paper contains 33 sections, 12 theorems, 141 equations, 3 figures, 1 table.

Key Result

Lemma 1

Recall $H_{ij}$'s defined in eqn:H-and-L, and assume that $q, g$ are supported only on the first $K$-entries with each entry bounded. Then

Figures (3)

  • Figure 1: From top to bottom row: Standard Treatment imd, Standard Treatment cum, and Our Rapid Interleaving.
  • Figure 2: Circular convolution model $\mathbf{y} = \mathbf{x} \oast g$: we plot the histogram of $\widehat{\tau}(\mathbf{1}_{<K})$ defined in \ref{['eqn:estimator-lag-K']} based on 20000 Monte Carlo simulations, the red solid vertical line corresponds to $\tau(\mathbf{1}_{<K})$ defined in \ref{['eqn:lag-K']}, the red dashed lines are two standard deviations from the expectation using the variance formula in Proposition \ref{['prop:second-moment']}.
  • Figure 3: Linear convolution model $\mathbf{y} = \mathbf{x} * g$: we plot the histogram of $\widehat{\tau}^{\mathcal{L}}(\mathbf{1}_{<K})$ defined in \ref{['eqn:estimator-lin']} based on 20000 Monte Carlo simulations, the red solid vertical line corresponds to $\tau^{\mathcal{L}}(\mathbf{1}_{<K})$ defined in \ref{['eqn:estimand-lin']}, the red dashed lines are two standard deviations from the expectation using the variance formula in Proposition \ref{['prop:second-moment']}.

Theorems & Definitions (21)

  • Lemma 1
  • proof
  • Proposition 1: Second Moment
  • Theorem 1: Quadratic Term
  • Lemma 2
  • Lemma 3: Fourth Moment Estimate
  • proof : Proof of Theorem \ref{['thm:Quadratic-Term']}
  • Theorem 2: Linear Term
  • proof : Proof of Theorem \ref{['lem:Linear-Term']}
  • Theorem 3: Uniform Consistency
  • ...and 11 more