Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA.
Department of Epidemiology and Biostatistics, McGill University, Montréal, Québec, Canada.
Stat Med. 2024 Jan 15;43(1):34-48. doi: 10.1002/sim.9940. Epub 2023 Nov 5.
Within the principal stratification framework in causal inference, the majority of the literature has focused on binary compliance with an intervention and modelling means. Yet in some research areas, compliance is partial, and research questions-and hence analyses-are concerned with causal effects on (possibly high) quantiles rather than on shifts in average outcomes. Modelling partial compliance is challenging because it can suffer from lack of identifiability. We develop an approach to estimate quantile causal effects within a principal stratification framework, where principal strata are defined by the bivariate vector of (partial) compliance to the two levels of a binary intervention. We propose a conditional copula approach to impute the missing potential compliance and estimate the principal quantile treatment effect surface at high quantiles, allowing the copula association parameter to vary with the covariates. A bootstrap procedure is used to estimate the parameter to account for inflation due to imputation of missing compliance. Moreover, we describe precise assumptions on which the proposed approach is based, and investigate the finite sample behavior of our method by a simulation study. The proposed approach is used to study the 90th principal quantile treatment effect of executive stay-at-home orders on mitigating the risk of COVID-19 transmission in the United States.
在因果推断的主要分层框架内,大多数文献都集中在二元干预和建模手段的合规性上。然而,在一些研究领域,合规性是部分的,研究问题——因此分析——关注的是对(可能较高的)分位数的因果效应,而不是对平均结果的偏移。对部分合规性进行建模具有挑战性,因为它可能存在可识别性问题。我们在主要分层框架内开发了一种估计分位数因果效应的方法,其中主要分层由二元干预的两个水平的(部分)合规性的双变量向量定义。我们提出了一种条件 Copula 方法来推断缺失的潜在合规性,并在较高分位数处估计主要分位数治疗效果表面,允许 Copula 关联参数随协变量而变化。使用 bootstrap 过程来估计参数,以弥补因缺失合规性推断而导致的膨胀。此外,我们描述了所提出方法的精确假设,并通过模拟研究调查了我们方法的有限样本行为。所提出的方法用于研究行政居家令对减轻美国 COVID-19 传播风险的第 90 个主要分位数治疗效果。