Table of Contents
Fetching ...

Estimating causal effects of functional treatments with modified functional treatment policies

Ziren Jiang, Erjia Cui, Jared D. Huling

TL;DR

This work proposes a novel definition of the population average over a functional variable using a functional principal component analysis (FPCA) decomposition, and establishes the causal identifiability of the MFTP estimand.

Abstract

Functional data are increasingly prevalent in biomedical research. While functional data analysis has been established for decades, causal inference with functional treatments remains largely unexplored. Existing methods typically focus on estimating the causal average dose response functional (ADRF), which requires strong positivity assumptions and offers limited interpretability. In this work, we target a new causal estimand, the modified functional treatment policy (MFTP), which focuses on estimating the average potential outcome when each individual slightly modifies their treatment trajectory from the observed one. A major challenge for this new estimand is the need to define an average over an infinite-dimensional object with no density. By proposing a novel definition of the population average over a functional variable using a functional principal component analysis (FPCA) decomposition, we establish the causal identifiability of the MFTP estimand. We further derive outcome regression, inverse probability weighting, and doubly robust estimators for the MFTP, and provide theoretical guarantees under mild regularity conditions. The proposed estimators are validated through extensive simulation studies. Applying our MFTP framework to the National Health and Nutrition Examination Survey (NHANES) accelerometer data, we estimate the causal effects of reducing disruptive nighttime activity and low-activity duration on all-cause mortality.

Estimating causal effects of functional treatments with modified functional treatment policies

TL;DR

This work proposes a novel definition of the population average over a functional variable using a functional principal component analysis (FPCA) decomposition, and establishes the causal identifiability of the MFTP estimand.

Abstract

Functional data are increasingly prevalent in biomedical research. While functional data analysis has been established for decades, causal inference with functional treatments remains largely unexplored. Existing methods typically focus on estimating the causal average dose response functional (ADRF), which requires strong positivity assumptions and offers limited interpretability. In this work, we target a new causal estimand, the modified functional treatment policy (MFTP), which focuses on estimating the average potential outcome when each individual slightly modifies their treatment trajectory from the observed one. A major challenge for this new estimand is the need to define an average over an infinite-dimensional object with no density. By proposing a novel definition of the population average over a functional variable using a functional principal component analysis (FPCA) decomposition, we establish the causal identifiability of the MFTP estimand. We further derive outcome regression, inverse probability weighting, and doubly robust estimators for the MFTP, and provide theoretical guarantees under mild regularity conditions. The proposed estimators are validated through extensive simulation studies. Applying our MFTP framework to the National Health and Nutrition Examination Survey (NHANES) accelerometer data, we estimate the causal effects of reducing disruptive nighttime activity and low-activity duration on all-cause mortality.
Paper Structure (21 sections, 10 theorems, 35 equations, 9 figures)

This paper contains 21 sections, 10 theorems, 35 equations, 9 figures.

Key Result

Lemma 1

(Limit of $\mu_K$ exists) Under conditions C1 and C2, the sequence of $\{\mu_K\}_{K=1}^{\infty}$ is a Cauchy sequence. Therefore, the limit of $\mu_K$ exists.

Figures (9)

  • Figure 1: An illustration of minute-level physical activity (PA) records for 50 randomly sampled individuals in the NHANES study, where the PA trajectories of three individuals are highlighted in color and the remaining trajectories are shown in grey.
  • Figure 2: An illustration of the potential outcome for the MFTP. The left panel is the observed functional treatment, and the right panel is the modified treatment where we only modified part of the trajectory.
  • Figure 3: Simulation results with number of principal components under balancing equals 4. The y-axis is the log of mean squared error and the x-axis is the log of sample size. The scales of y-axis and x-axis are adjusted to be the same.
  • Figure 4: Coverage results for the estimated 95% confidence interval under various simulation scenarios. The number of principal components under balancing equals 4.
  • Figure 5: Log of mean squared error for the estimators under different number of principal components.
  • ...and 4 more figures

Theorems & Definitions (13)

  • Lemma 1
  • Lemma 2
  • Corollary 3.1
  • Theorem 3
  • Lemma 4
  • Theorem 5
  • Example 1
  • Example 2
  • Example 3
  • Lemma 6
  • ...and 3 more