We propose ULFS-KDPE, a kernel debiased plug-in estimator based on the universal least favorable submodel, for estimating pathwise differentiable parameters in nonparametric models. The method constructs a data-adaptive debiasing flow in a reproducing kernel Hilbert space (RKHS), producing a plug-in estimator that achieves semiparametric efficiency without requiring explicit derivation or evaluation of efficient influence functions. We place ULFS-KDPE on a rigorous functional-analytic foundation by formulating the universal least favorable update as a nonlinear ordinary differential equation on probability densities. We establish existence, uniqueness, stability, and finite-time convergence of the empirical score along the induced flow. Under standard regularity conditions, the resulting estimator is regular, asymptotically linear, and attains the semiparametric efficiency bound simultaneously for a broad class of pathwise differentiable parameters. The method admits a computationally tractable implementation based on finite-dimensional kernel representations and principled stopping criteria. In finite samples, the combination of solving a rich collection of score equations with RKHS-based smoothing and avoidance of direct influence-function evaluation leads to improved numerical stability. Simulation studies illustrate the method and support the theoretical results.
Figures & tables
Parameter
Method
#Cov.
Bias( ×100 )
Var a
RMSE
ULFS-KDPE
475
-0.824
0.0561
0.0567
DGP1
ψATE
KDPE
343
-1.438
0.0546
0.0565
TMLE
500
-0.004
0.0618
0.0618
One-step TMLE
500
-0.247
0.0605
0.0606
ULFS-KDPE
475
-1.082
0.1724
0.1728
ψRR
KDPE
343
-3.213
0.1584
0.1616
Table 1. Bias, Convergence with 100 Steps, Variance and Root Mean Squared Error (RMSE) for ATE, RR and OR with different estimators (ULFS-KDPE, KDPE, TMLE, One-step TMLE).
Figure 1. Finite-sample behavior of the different ATE estimators (ULFS-KDPE, KDPE, TMLE, One-step TMLE). True asymptotic distribution is depicted in purple.
Parameter
Stopping Criterion
#Cov.
Bias ( ×100 )
Var a
RMSE
Density stabilization
500
-0.840
0.0553
0.0560
DGP1
ψATE
Empirical score stabilization
493
-0.806
0.0566
0.0572
Vanishing update direction
500
-1.357
0.0547
0.0564
Variance-dominated updates
500
-1.573
0.0550
0.0573
EIF approximately solved
447
-1.156
0.0561
0.0573
Density stabilization
500
-1.142
0.1700
0.1705
Table 2. Bias, Convergence with 200 Steps, Var, RMSE for ATE, RR and OR under different stopping criteria for ULFS-KDPE.
Figure 2. ULFS-KDPE under different stopping criteria. True asymptotic distribution is depicted in purple.
Iteration Limit
Stopping Criterion
#Cov.
ψATE
ψRR
ψOR
Density plateau
475
0.0567
0.1728
0.3975
DGP1
100
Score plateau
404
0.0573
0.1676
0.3792
Vanishing Direction
498
0.0563
0.1669
0.3784
Variance dominate
500
0.0567
0.1684
0.3803
EIF sovled
407
0.0547
0.1721
0.3943
Density plateau
500
0.0560
0.1702
0.3896
Table 3. Convergence and RMSE for ATE, RR and OR with ULFS-KDPE under different iteration limits (# of steps allowed).
Semiparametric inference quantifies uncertainty for low-dimensional targets when models also contain flexible, infinite-dimensional nuisance functions.
We study functional bilevel optimization with conditional moments, where plug-in hypergradients can inherit first-order errors from nonparametrically learned conditional expectations. We introduce OBiGrad, an orthogonalized hypergradient estimator whose bias is a product of nuisance errors. Under our regression and complexity conditions, projected gradient descent with the direct plug-in hypergradient yields expected stationarity error O(T−1/2+n−1/2+n−ν/b) for the best iterate, whereas OBiGrad achieves O(T−1/2+n−1/2+n−2ν/b), where b=ν(2+p)+γ combines the approximation and complexity exponents. Consequently, when ν(2−p)≥γ, OBiGrad attains the parametric statistical rate O(T−1/2+n−1/2), improving the sufficient sample-size bound relative to the direct plug-in bound for a fixed iteration budget. We illustrate these gains on instrumental-variable and fitted-Q-regression problems, examining stationary-point estimation, returned population stationarity, and regularization sensitivity.
Fares El Khoury, Houssam Zenati, Nathan Kallus +2
Université Grenoble Alpes, Inria, CNRS, Grenoble INP, LJK · Gatsby Computational Neuroscience Unit, University College London · Cornell University +1
Prior-data fitted networks (PFNs) have recently gained popularity as an end-to-end approach to amortized Bayesian inference. Yet, their statistical behavior for semi-parametric targets, such as means under missingness or causal effects, remains poorly understood. In this work, we develop theoretical foundations for debiased PFNs to perform inference on finite-dimensional targets in general semi-parametric settings. First, we show that, despite being Bayes-optimal with respect to the target, an idealized PFN can exhibit a posterior concentration rate that implicitly depends on the concentration rates of (latent) nuisance functions and can become slower as the prior becomes more diverse. Then, to speed up the target posterior concentration rates, we introduce two generic debiasing strategies based on efficient influence functions: (i) a label correction and (ii) a prior correction. Finally, we show that the debiased PFNs (i) & (ii) possess concentration rates that asymptotically match a convergence rate given by traditional frequentist bias-corrected estimators. In this way, we establish a novel connection between amortized Bayesian inference with PFNs and frequentist semi-parametric efficiency theory. Under additional assumptions, this allows us to guarantee local and global optimality of PFN-based estimators (namely, the Bernstein-von Mises property and minimax optimality, respectively). Therefore, our work provides guidelines for designing debiased PFNs for semi-parametric inference with principled and provable asymptotic guarantees.
Valentyn Melnychuk, Emil Javurek, Stefan Feuerriegel
LMU Munich & Munich Center for Machine Learning, Munich, Germany
This study investigates semiparametric efficient estimation of causal and structural parameters in a semi-supervised setting. In our setting, unlabeled auxiliary regressors are available in addition to labeled observations consisting of outcomes and regressors. Our goal is to construct estimators of causal and structural parameters whose asymptotic variances are smaller than those of estimators constructed using only labeled data. We refer to this framework as prediction-powered causal inference (PPCI). We first derive the efficient influence function and the efficiency bound, which imply that the use of auxiliary regressors can attain a smaller asymptotic variance than the efficiency bound attainable from labeled observations alone. Then, by combining the efficient influence function with the debiased machine learning (DML) framework, we propose methods that we call DML-PPCI. If we construct an estimating-equation estimator, we refer to the method as EE-DML-PPCI; if we construct a targeted-learning estimator, we refer to the method as TMLE-DML-PPCI. The asymptotic variances of both estimators match our derived efficiency bound. In the construction of the estimators, estimation of the efficient influence function plays an important role. In our study, the efficient influence function is also a Neyman orthogonal score, which depends on the Riesz representer and the regression function. For Riesz representer estimation, we develop semi-supervised generalized Riesz regression with convergence rate guarantees.
Masahiro Kato
Data Analytics Department, Mizuho-DL Financial Technology, Co., Ltd.