Conformal prediction (CP) offers a principled framework for uncertainty quantification, but it fails to guarantee mask-conditional coverage when faced with missing covariates. In addressing the heterogeneity induced by various missing patterns, Mask-Conditional Valid (MCV) Coverage has emerged as a more desirable property than Marginal Coverage. In this work, we adapt split CP to handle missing values by proposing a preimpute-mask-then-correct framework that can offer valid coverage. We show that our method provides guaranteed Marginal Coverage and Mask-Conditional Validity for general missing data mechanisms. A key component of our approach is a reweighted conformal prediction procedure that corrects the prediction sets after distributional imputation (multiple imputation) of the calibration dataset, making our method compatible with standard imputation pipelines. We derive two algorithms and prove that they achieve both marginal validity and MCV. We evaluate them on synthetic and real-world datasets. It reduces significantly the width of prediction intervals w.r.t standard MCV methods, while maintaining the target guarantees.
Figures & tables
Figure 1: Illustration of the pre-imputation: the calibration data is imputed and then masked by the M of test point such that they share the same pattern. In this way, the split CP can give a prediction set conditional on M . But a distribution shift exists. Unlike CP-MDA, our method leverages all available calibration data without modifying the test point.
Figure 2: Synthetic dataset on MCAR: prediction interval width across 1023 masks (mean and 95% CI).
Figure 3: Synthetic dataset on MCAR: prediction coverage across 1023 masks (mean and 95% CI).
Mechanism
MDA-Exact
MDA-Nested
Weighted CP
ARC CP
MAR
87.8–89.7% / 5.84
87.9–89.8% / 5.67
89.0–91.2% / 5.79
88.3–90.2% / 5.64
MNAR
87.8–89.8% / 5.74
87.8–89.8% / 5.64
89.1–91.0% / 5.79
88.5–90.3% / 5.60
Table 1: Worst-case coverage (95% CI) and average prediction interval width over 31 masks.
Bin of Loracle(m)
ARC-CP cov
ARC-CP wid
WCP cov
WCP wid
Nested cov
Nested wid
[17.74,20.39]
0.911
27.92
0.884
27.14
0.926
29.58
[20.39,21.16]
0.910
27.99
0.885
27.20
0.927
29.76
[21.16,21.57]
0.910
28.04
0.885
27.16
0.928
29.90
[21.57,22.02]
0.911
28.08
0.881
27.08
0.927
29.91
[22.02,22.40]
0.911
28.20
0.882
27.17
0.928
30.04
[22.40,22.77]
0.911
28.14
0.884
27.21
0.927
30.17
Table 2: 100-dimensional synthetic setting: average mask-conditional coverage and interval width, grouped by oracle interval width Loracle(m) . Results of coverage and width are averaged over 1000 masks with 50 replications.
Appendix figures & tables20 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 4: Synthetic dataset: mean and 99% confidence interval of prediction interval width with respect to 1023 missingness patterns, results by using deterministic imputation.
Figure 5: Synthetic dataset: mean and 99% confidence interval of prediction coverage with respect to 1023 missingness patterns, results by using deterministic imputation.
Algorithm 4 Acceptance-rejection corrected CP with pre-imputation of Calibration
Figure 6: Coverage vs. Correlation of estimated and true density ratios. Each point represents one experiment. High correlation yields near-nominal coverage.
Figure 7: Miscoverage vs estimated bound. Points below y=x confirm the bound’s conservativeness.
Method
Synthetic dataset
Concrete dataset
MDA-Exact
17034s
2031s
MDA-Nested
510932s
43225s
Weighted CP
228638s
37983s
ARC CP
21320s
4766s
Appendix
Table 3: Comparison of running time of 4 methods.
Method
Avg. Cov
Avg.Width
MDA-Exact
92.9±0.7
∞
MDA-Nested
91.8±0.5
336.3±2.3
Weighted CP
91.9±0.6
306.4±2.2
ARC CP
93.0±0.6
309.9±3.1
Appendix
Table 4: Comparison on Bike sharing dataset for different methods. The average coverage and interval width of 15 different missingness patterns are listed.
Method
Coverage
Interval width
MDA-Exact
94.03−95.01
33.0
MDA-Nested
94.29−95.15
34.1
ARC CP
94.04−95.04
32.08
Appendix
Table 5: Comparison on Bike sharing dataset for different methods. The worst-case coverage and interval width of 32 different missingness patterns are listed.
Figure 8: Synthetic dataset: Comparison of prediction coverage (mean and 99% confidence interval) across 1023 missingness patterns, with and without the proposed correction method.
Figure 9: Synthetic dataset: Comparison of the width of prediction interval (mean and 99% confidence interval) across 1023 missingness patterns, with and without the proposed correction method.
Figure 10: Concrete dataset: mean and 99% confidence interval of prediction coverage with respect to 255 missingness patterns
Mech.
MDA-Exact
MDA-Nested
Weighted CP
ARC CP
MAR
89.1% ± 1.64% / 29.1
90.8% ± 1.43% / 31.4
89.1% ± 1.91% / 29.8
90.6% ± 1.42% / 28.9
Appendix
Table 6: Worst-case coverage (95% CI) and average prediction interval width over 15 masks, concrete dataset with 4 possible MAR missing covariates
Mech.
MDA-Exact
MDA-Nested
Weighted CP
ARC CP
MNAR
89.0% ± 1.78% / 30.1
88.7% ± 1.76% / 30.4
90.0% ± 1.67% / 30.8
89.0% ± 1.83% / 30.1
Appendix
Table 7: Worst-case coverage (95% CI) and average prediction interval width over 15 masks, concrete dataset with 4 possible missing covariates
Figure 11: Synthetic dataset with 20% missing values: mean and 99% confidence interval of prediction interval width with respect to 1023 missingness patterns, comparison for whether using correction method.
Figure 12: Synthetic dataset with 20% missing values: mean and 99% confidence interval of prediction coverage with respect to 1023 missingness patterns, comparison for whether using correction method.
Figure 13: Synthetic dataset with 20% missing values and independent features: mean and 99% confidence interval of prediction interval width with respect to 1023 missingness patterns, comparison for whether using correction method.
Figure 14: Synthetic dataset with 20% missing values and independent features: mean and 99% confidence interval of prediction coverage with respect to 1023 missingness patterns, comparison for whether using correction method.
Figure 15: Synthetic dataset: unsmoothed results for mean and 99% confidence interval of prediction width with respect to 1023 missingness patterns.
Figure 16: Synthetic dataset: unsmoothed results for mean and 99% confidence interval of prediction coverage with respect to 1023 missingness patterns.
Conformal prediction (CP) provides a distribution-free framework for uncertainty quantification with finite-sample marginal coverage guarantees. Yet, missing values introduce distributional heterogeneity, where marginal coverage can hide systematic undercoverage for certain missingness patterns. At the same time, existing methods based on distances or smoothness require a continuous representation of partially observed covariates. Although imputation restores continuity, it obscures the missingness patterns and introduces task-dependent overhead. We develop a new CP framework that is adaptive to missingness patterns, without imputing covariates. We rely on data augmentation by feature embedding that accounts for missingness. Our framework represents each input with missing values by a tree embedding together with its raw missingness mask, and then learns a regularized quantile score function. Varying the choices for the penalty function and the augmentation enables us to characterize three practically relevant validity levels: exact coordinate-wise validity, exact validity on observed masks, and approximate missingness-conditional guarantees for general mechanisms (MCAR, MAR, MNAR). Experiments on synthetic and real datasets support these guarantees and show that the proposed method reduces miscoverage risk while maintaining informative intervals, especially in non-MCAR settings. These results highlight the practical value of uncertainty quantification that adapts to missingness.
Conformal prediction provides rigorous distribution-free finite-sample guarantees for marginal coverage under the assumption of exchangeability, but may exhibit systematic undercoverage or overcoverage for specific subpopulations. Assessing conditional validity is challenging, as standard stratification methods suffer from the curse of dimensionality. We propose Conformal Prediction Assessment (CPA), a framework that reframes the evaluation of conditional coverage as a supervised learning task by training a reliability estimator that predicts instance-level coverage probabilities. Building on this estimator, we introduce the Conditional Validity Index (CVI), which decomposes reliability into safety (undercoverage risk) and efficiency (overcoverage cost). We establish convergence rates for the reliability estimator and prove the consistency of CVI-based model selection. Extensive experiments on synthetic and real-world datasets demonstrate that CPA effectively diagnoses local failure modes and that CC-Select, our CVI-based model selection algorithm, consistently identifies predictors with superior conditional coverage performance.
Zheng Zhou, Xiangfei Zhang, Chongguang Tao +1
Qiuzhen College, Tsinghua University · Yau Mathematical Sciences Center, Tsinghua University
A predictor with small overall bias can still exhibit systematic bias conditional on its predicted value, limiting its reliability for downstream uncertainty quantification and candidate screening. We propose Isotonic Conformal Inference, a unified framework for prediction and multiple testing that separates calibration from inference and exploits the level set structure of an isotonic predictor. This structure enables the computation of localized conformal quantiles and p-values while avoiding repeated isotonic refitting across candidate responses. For prediction, we introduce Split Isotonic Conformal Prediction (sICP) and Transductive Isotonic Conformal Prediction (tICP), establish prediction-conditional coverage guarantees for both methods, and characterize their finite-sample and asymptotic self-calibration properties. For selection, we introduce Isotonic Conformal Selection (ICS), which provides finite-sample stratum-conditional false discovery rate control. Experiments on simulated regression problems, a healthcare-utilization dataset, and a high-throughput ADMET drug discovery dataset demonstrate the practical performance of the framework and substantially lower computational costs relative to full transductive baselines.
Daniel Bensimon, Zhiquan Liu, Sean X. Yu +2
Department of Mathematics and Statistics, McGill University; Mila · Eli Lilly and Company · Corresponding author, Department of Mathematics and Statistics, McGill University; Mila