High-Dimensional Feature Selection for Sample Efficient Treatment Effect Estimation
This addresses the challenge of sample inefficiency in causal inference for practitioners dealing with high-dimensional observational data, offering a method to reduce bias and improve estimation accuracy.
The paper tackles the problem of high-dimensional feature selection for treatment effect estimation by proposing a method that identifies a sparse subset of covariates sufficient for unbiased estimation, reducing sample complexity from scaling with the full set size to the sparse subset size and log of total covariates.
The estimation of causal treatment effects from observational data is a fundamental problem in causal inference. To avoid bias, the effect estimator must control for all confounders. Hence practitioners often collect data for as many covariates as possible to raise the chances of including the relevant confounders. While this addresses the bias, this has the side effect of significantly increasing the number of data samples required to accurately estimate the effect due to the increased dimensionality. In this work, we consider the setting where out of a large number of covariates $X$ that satisfy strong ignorability, an unknown sparse subset $S$ is sufficient to include to achieve zero bias, i.e. $c$-equivalent to $X$. We propose a common objective function involving outcomes across treatment cohorts with nonconvex joint sparsity regularization that is guaranteed to recover $S$ with high probability under a linear outcome model for $Y$ and subgaussian covariates for each of the treatment cohort. This improves the effect estimation sample complexity so that it scales with the cardinality of the sparse subset $S$ and $\log |X|$, as opposed to the cardinality of the full set $X$. We validate our approach with experiments on treatment effect estimation.