Nonparametric Bootstrap Inference for the Targeted Highly Adaptive LASSO Estimator

TL;DR

This paper introduces a nonparametric bootstrap method for the HAL-TMLE estimator, providing consistent inference and improved confidence interval coverage for complex, high-dimensional models with minimal assumptions.

Contribution

It establishes the bootstrap's consistency for HAL-TMLE and proposes a novel selection method for the sectional variation norm to optimize confidence interval coverage.

Findings

Bootstrap provides consistent normal limit distribution estimation.

Proposed method improves finite sample coverage of confidence intervals.

Simulation results demonstrate excellent coverage in practical examples.

Abstract

The Highly-Adaptive-LASSO Targeted Minimum Loss Estimator (HAL-TMLE) is an efficient plug-in estimator of a pathwise differentiable parameter in a statistical model that at minimal (and possibly only) assumes that the sectional variation norm of the true nuisance functional parameters (i.e., the relevant part of data distribution) are finite. It relies on an initial estimator (HAL-MLE) of the nuisance functional parameters by minimizing the empirical risk over the parameter space under the constraint that the sectional variation norm of the candidate functions are bounded by a constant, where this constant can be selected with cross-validation. In this article, we establish that the nonparametric bootstrap for the HAL-TMLE, fixing the value of the sectional variation norm at a value larger or equal than the cross-validation selector, provides a consistent method for estimating the normal limit distribution of the HAL-TMLE. In order to optimize the finite sample coverage of the nonparametric bootstrap confidence intervals, we propose a selection method for this sectional variation norm that is based on running the nonparametric bootstrap for all values of the sectional variation norm larger than the one selected by cross-validation, and subsequently determining a value at which the width of the resulting confidence intervals reaches a plateau. We demonstrate our method for 1) nonparametric estimation of the average treatment effect based on observing on each unit a covariate vector, binary treatment, and outcome, and for 2) nonparametric estimation of the integral of the square of the multivariate density of the data distribution. In addition, we also present simulation results for these two examples demonstrating the excellent finite sample coverage of bootstrap-based confidence intervals.