Prediction-Powered Conditional Inference
Organizations: University of California, Los Angeles
Abstract
We study prediction-powered conditional inference in the setting where labeled data are scarce, unlabeled covariates are abundant, and a black-box machine-learning predictor is available. The goal is to perform statistical inference on conditional functionals evaluated at a fixed target point, such as conditional means, without imposing a parametric model for the conditional relationship. Our approach combines localization with prediction-based variance reduction. First, we introduce an RKHS localization method that learns a data-adaptive weight from covariates and reformulates the target conditional moment at the target point as a weighted unconditional moment. Second, we incorporate machine-learning predictions through a correction-based decomposition of this localized moment, yielding a prediction-powered estimator and confidence interval that reduce variance when the predictor is informative while preserving validity regardless of predictor accuracy. We establish nonasymptotic error bounds and, in the abundant-unlabeled regime, minimax-optimal convergence rates for the resulting estimator, prove pointwise asymptotic normality with consistent variance estimation, and provide an explicit variance decomposition that characterizes how machine-learning predictions and unlabeled covariates improve statistical efficiency. Numerical experiments on simulated and real datasets demonstrate valid conditional coverage and substantially sharper confidence intervals than alternative methods.
Figures & tables
Appendix figures & tables7 assets
Supplementary material from the paper’s appendix.
Appendix
| Experiment | Screen constants | |||
|---|---|---|---|---|
| Simulation | ||||
| Predictor quality | ||||
| Income (sex ) | ||||
| Income (sex ) | ||||
| BlogFeedback |
| Implementation | Coverage | Bias | RMSE | Width | Wall time (s) | Peak memory (MiB) |
|---|---|---|---|---|---|---|
| Two-fold | 0.964 | 0.0467 | 0.3465 | 1.3506 | 11.35 (1.00) | 2,680 (1.00) |
| No-split | 0.964 | 0.0472 | 0.3469 | 1.3524 | 37.15 (3.27) | 7,855 (2.93) |
| Component | Configuration | Coverage | Bias | RMSE | Width | Fallback |
|---|---|---|---|---|---|---|
| 0.952 (0.930) | 0.0124 | 0.3876 | 1.4510 | 1.000 | ||
| 0.952 (0.931) | 0.0126 | 0.3862 | 1.4465 | 0.897 | ||
| 0.950 (0.925) | 0.0222 | 0.3629 | 1.3640 | 0.000 | ||
| 0.945 (0.909) | 0.0402 | 0.3392 | 1.2721 | 0.000 | ||
| 0.939 (0.891) | 0.0577 | 0.3245 | 1.2076 | 0.000 | ||
| Bandwidth cap | 0.950 (0.925) | 0.0222 | 0.3629 | 1.3640 | 0.000 |
| Data | Reference | Shift | Coverage | RMSE | Width | |||
|---|---|---|---|---|---|---|---|---|
| PPCI | LO | PPCI | LO | PPCI | LO | |||
| Income | Matérn, | 0.1976 | 0.9329 | 0.9329 | 0.5181 | 0.5782 | 1.9079 | 2.1528 |
| Matérn, | 0 | 0.9429 | 0.9375 | 0.5571 | 0.6134 | 1.9079 | 2.1528 | |
| Matérn, | 0.1681 | 0.9233 | 0.9242 | 0.6404 | 0.6898 | 1.9079 | 2.1528 | |
| Gaussian, ESS matched | 0.0630 | 0.9400 | 0.9408 | 0.5364 | 0.5946 | 1.9079 | 2.1528 | |
| BlogFeedback | Matérn, | 0.0200 | 0.9175 | 0.9404 | 0.0425 | 0.0603 | 0.1442 | 0.2227 |
| Method | Coverage | Bias | RMSE | Width | |
|---|---|---|---|---|---|
| 0 | NW-CV | 0.964 | 0.0302 | 0.1818 | 0.7110 |
| 0 | NW-US | 0.970 | 0.0145 | 0.2027 | 0.7958 |
| 0 | RKHS localization | 0.970 | -0.0005 | 0.1674 | 0.6808 |
| 2 | NW-CV | 0.960 | 0.0288 | 0.1833 | 0.7171 |
| 2 | NW-US | 0.970 | 0.0135 | 0.2044 | 0.8014 |
| 2 | RKHS localization | 0.970 | -0.0005 | 0.1674 | 0.6813 |
| Method | Coverage | Bias | RMSE | ||
|---|---|---|---|---|---|
| 0.08 | NW | 0.956 | 0.0260 | 0.1766 | 0.007690 |
| 0.08 | RKHS localization | 0.980 | 0.0023 | 0.3806 | 0.000047 |
| 0.20 | NW | 0.652 | 0.1851 | 0.2175 | 0.047563 |
| 0.20 | RKHS localization | 0.968 | 0.0019 | 0.2601 | 0.000127 |
| 0.40 | NW | 0.000 | 0.5879 | 0.5943 | 0.147705 |
| 0.40 | RKHS localization | 0.952 | -0.0049 | 0.1960 | 0.000259 |