Conformal prediction provides a powerful framework for constructing distribution-free prediction regions with finite-sample coverage guarantees. While extensively studied in univariate settings, its extension to multi-output problems presents additional challenges, including complex output dependencies and high computational costs, and remains relatively underexplored. In this work, we present a unified comparative study of nine conformal methods with different multivariate base models for constructing multivariate prediction regions within the same framework. This study highlights their key properties while also exploring the connections between them. Additionally, we introduce two novel classes of conformity scores for multi-output regression that generalize their univariate counterparts. These scores ensure asymptotic conditional coverage while maintaining exact finite-sample marginal coverage. One class is compatible with any generative model, offering broad applicability, while the other is computationally efficient, leveraging the properties of invertible generative models. Finally, we conduct a comprehensive empirical evaluation across 13 tabular datasets, comparing all the multi-output conformal methods explored in this work. To ensure a fair and consistent comparison, all methods are implemented within a unified code base.
@inproceedings{Dheur2025-br,title={A unified comparative study with generalized conformity scores for multi-output conformal regression},author={Dheur, Victor and Fontana, Matteo and Estievenart, Yorick and Desobry, Naomi and Ben Taieb, Souhaib},booktitle={The 42th International Conference on Machine Learning},year={2025},}
We present a new method for generating confidence sets within the split conformal prediction framework. Our method performs a trainable transformation of any given conformity score to improve conditional coverage while ensuring exact marginal coverage. The transformation is based on an estimate of the conditional quantile of conformity scores. The resulting method is particularly beneficial for constructing adaptive confidence sets in multi-output problems where standard conformal quantile regression approaches have limited applicability. We develop a theoretical bound that captures the influence of the accuracy of the quantile estimate on the approximate conditional validity, unlike classical bounds for conformal prediction methods that only offer marginal coverage. We experimentally show that our method is highly adaptive to the local data structure and outperforms existing methods in terms of conditional coverage, improving the reliability of statistical inference in various applications.
@inproceedings{Plassier2025-dh,title={Rectifying conformity scores for better conditional coverage},author={Plassier*, Vincent and Fishkov*, Alexander and Dheur*, Victor and Guizani, Mohsen and Ben Taieb, Souhaib and Panov, Maxim and Moulines, Eric},booktitle={The 42th International Conference on Machine Learning},year={2025},note={(<em>* denotes equal contribution.</em>)}}
With the growing utilization of machine learning models in real-world applications, improving neural network calibration has become a primary concern. Various methods have been proposed, including post-hoc methods that adjust predictions after training and regularization methods that act during training. In regression, post-hoc methods have shown superior calibration performance compared to regularization methods. However, the base model is trained without any direct connection to the subsequent post-hoc step. To address this limitation, we introduce a novel end-to-end method called Quantile Recalibration Training, integrating post-hoc calibration directly into the training process. We also propose an algorithm unifying our method with other post-hoc and regularization methods. We demonstrate the performance of Quantile Recalibration Training in a large-scale experiment involving 57 tabular regression datasets, showcasing improved predictive accuracy. Additionally, we conduct an ablation study, revealing the significance of different elements in our method. Aiming for reproducibility and fair comparisons, we have implemented our experiments in a common code base.
@inproceedings{Dheur2024-AISTATS,title={Probabilistic Calibration by Design for Neural Network Regression},booktitle={The 27th International Conference on Artificial Intelligence and Statistics},author={Dheur, Victor and Ben Taieb, Souhaib},year={2024},}
Sequences of labeled events observed at irregular intervals
in continuous time are ubiquitous across various fields.
Temporal Point Processes (TPPs) provide a mathematical
framework for modeling these sequences, enabling inferences
such as predicting the arrival time of future events and
their associated label, called mark. However, due to model
misspecification or lack of training data, these
probabilistic models may provide a poor approximation of the
true, unknown underlying process, with prediction regions
extracted from them being unreliable estimates of the
underlying uncertainty. This paper develops more reliable
methods for uncertainty quantification in neural TPP models
via the framework of conformal prediction. A primary
objective is to generate a distribution-free joint
prediction region for the arrival time and mark, with a
finite-sample marginal coverage guarantee. A key challenge
is to handle both a strictly positive, continuous response
and a categorical response, without distributional
assumptions. We first consider a simple but overly
conservative approach that combines individual prediction
regions for the event arrival time and mark. Then, we
introduce a more effective method based on bivariate highest
density regions derived from the joint predictive density of
event arrival time and mark. By leveraging the dependencies
between these two variables, this method exclude unlikely
combinations of the two, resulting in sharper prediction
regions while still attaining the pre-specified coverage
level. We also explore the generation of individual
univariate prediction regions for arrival times and marks
through conformal regression and classification techniques.
Moreover, we investigate the stronger notion of conditional
coverage. Finally, through extensive experimentation on both
simulated and real-world datasets, we assess the validity
and efficiency of these methods.
@inproceedings{Dheur2024-ML,title={Distribution-Free Conformal Joint Prediction Regions for Neural Marked Temporal Point Processes},booktitle={Machine Learning, Call for Papers: Conformal Prediction and Distribution-Free Uncertainty Quantification},author={Dheur*, Victor and Bosser*, Tanguy and Izbicki, Rafael and Ben Taieb, Souhaib},year={2024},note={(<em>* denotes equal contribution.</em>)}}
Accurate estimation of predictive uncertainty is essential for
optimal decision making. However, recent works have shown that
current neural networks tend to be miscalibrated, sparking
interest in different approaches to calibration. In this paper,
we conduct a large-scale empirical study of the probabilistic
calibration of neural networks on 57 tabular regression
datasets. We consider recalibration, conformal and
regularization approaches, and investigate the trade-offs they
induce on calibration and sharpness of the predictions. Based on
kernel density estimation, we design new differentiable
recalibration and regularization methods, yielding new insights
into the performance of these approaches. Furthermore, we find
conditions under which recalibration and conformal prediction
are equivalent. Our study is fully reproducible and implemented
in a common code base for fair comparison.
@inproceedings{Dheur2023-ICML,title={A Large-Scale Study of Probabilistic Calibration in Neural Network Regression},booktitle={The 40th International Conference on Machine Learning},author={Dheur, Victor and Ben Taieb, Souhaib},year={2023},}