Modern conformal forecasting systems often maintain several adaptive pipelines that differ in base forecasters, conformity scores, calibration windows, and update rules. Comparing them is difficult because coverage is a hard constraint, whereas efficiency should be optimized only among feasible pipelines. We formulate this problem as sequential inference for a stochastic constrained argmin. At each time, the target is the set of minimum-cost pipelines satisfying multiple prefix-average conditional miscoverage constraints. We introduce Coverage-Constrained Sequential Model Confidence Sets (CC-SMCS), which separate certifiably feasible, possibly feasible, and possibly constrained-optimal pipelines. Using simultaneous martingale confidence sequences, CC-SMCS projects a rectangular confidence region onto the constrained argmin and admits an exact closed-form rule. With probability at least 1−δ, it contains every constrained-optimal pipeline simultaneously over all times. This finite-sample guarantee requires no stationarity or mixing assumptions and remains valid under data-dependent stopping. We also establish an impossibility result for safe certification at the coverage boundary and extend the construction to delayed multi-horizon feedback and outcome-dependent efficiency objectives.
Figures & tables
Simulations
δ
PTU
radˉt0.25
radˉt0.50
radˉt0.75
pτ,T
Simulation 1
0.10
1.000
0.0361
0.0234
0.0186
0.001
0.25
1.000
0.0326
0.0211
0.0168
0.001
Simulation 2
0.10
1.000
0.0402
0.0283
0.0231
–
0.25
1.000
0.0363
0.0255
0.0209
–
Simulation 3
0.10
0.999
0.0372
0.0252
0.0209
–
0.25
0.999
0.0336
0.0227
0.0188
–
Table 1: Simulation results
Figure 1: The performance of simulations with δ=0.10 . The vertical dashed line in Figure 1(b) indicates the structural break.
Figure 2: Empirical results for DLinear at 95% nominal coverage with δ=0.10 .
Appendix figures & tables8 assets
Supplementary material from the paper’s appendix.
Appendix
Symbol
Meaning
Cmrt
Prediction set from method m under constraint r at round t
Bmrt
Miscoverage indicator 1{Yt∈/Cmrt}
art
Predictable exposure or subgroup weight
Γmr,t
Prefix-average conditional excess miscoverage
Qm,t
Observed prefix-average set cost
Vt⋆
True feasible method set
Appendix
Table 2: Main notation.
The selection of ρ
δ
Simulation Data Set
Data Set 1
Data Set 2
Data Set 3
Data Set 4
ρ⋆
0.10
17.616
20.494
29.584
30.946
8.151
0.25
21.752
25.305
36.529
38.211
10.064
Appendix
Table 3: Selected values of ρ
Figure 3: The performance of simulations with δ=0.25 . The vertical dashed line in Figure 3(b) indicates the structural break.
Figure 4: Empirical results for DLinear at 95% nominal coverage with δ=0.25 .
Figure 5: Empirical results for PatchTST at 95% nominal coverage with δ=0.10 .
Figure 6: Empirical results for PatchTST at 95% nominal coverage with δ=0.25 .
Figure 7: Empirical results for TimesNet at 95% nominal coverage with δ=0.10 .
Figure 8: Empirical results for TimesNet at 95% nominal coverage with δ=0.25 .
Adaptive conformal inference (ACI) of Gibbs and Cand{è}s and its variants are the standard approach to online conformal prediction under distribution shift, but they suffer from three fundamental limitations. First, their guarantees control only the signed long-run coverage error: persistent miscoverage in one direction can be masked by compensating errors later, so a method can satisfy the theoretical guarantee while being badly wrong for extended periods. Second, existing guarantees say nothing about prediction-set size, so validity can be achieved trivially at the cost of unduly wide prediction sets. Third, the efficiency guarantees that do exist compare against a fixed predictor chosen in hindsight, a benchmark that becomes increasingly less meaningful once the data-generating distribution shifts, since the very notion of an optimal threshold then changes over time. We consider a unified online learning framework that simultaneously controls absolute, non-cancelling coverage violation and prediction-set efficiency against a dynamically evolving benchmark for three important models. In the fully adversarial setting, exploiting the fact that the standard ACI update is exactly projected online gradient descent on the pinball loss, we derive simultaneous coverage and efficiency guarantees for arbitrary monotone Lipschitz efficiency objectives, with no distributional or {\it convexity} assumptions. In the stochastic setting with full-score feedback, we propose a sliding-window quantile tracker and establish a matching minimax lower bound showing our algorithm is rate-optimal. In the covariate-dependent stochastic setting, we develop a partitioned ACI algorithm that tracks a function-valued oracle threshold, and derive simultaneous coverage and efficiency guarantees.
This article considers an online version of conformal inference, called adaptive conformal inference [ACI] and introduced by Gibbs and Candès (2021): prediction sets are issued sequentially, after observing features and before the outcomes are revealed. These sets are evaluated both in terms of validity (the fraction of rounds where the outcome was lying in the prediction set) and efficiency (the average lengths of the prediction sets). The two criteria point to different directions (validity favors larger sets). We also target a wide range of scenarios, with exchangeable data and arbitrary data (lack of any stochastic guarantees) as two extremes. A series of existing strategies for ACI typically guarantee that empirical coverage converges to the desired level for arbitrary sequences, but they generally lack simultaneous efficiency guarantees. To provide a unified study, we first formulate ACI as a repeated two-player game with finite action sets and vector-valued payoffs encoding validity and efficiency. Building on this reformulation, we introduce a strategy based on Blackwell approachability and on its opportunistic extension by Bernstein et al. (2014) that ensures validity while adapting the efficiency of the prediction intervals to the underlying degree of stochasticity of the opponent player. The resulting guarantee is "best of many worlds": it recovers the relevant efficiency guarantees in exchangeable and adversarial settings, and provides guarantees in intermediate settings that arise in typical applications such as the forecasting of time series.
Guillaume Principato, Gilles Stoltz
Universit´e Paris-Saclay, CNRS, Inria, Laboratoire de math´ematiques d’Orsay, 91405, Orsay, France · EDF R&D, Palaiseau, France · HEC Paris, Jouy-en-Josas, France
Conformal prediction provides a principled framework for uncertainty quantification with finite-sample coverage guarantees. While recent work has extended conformal prediction to online and sequential settings, existing methods typically focus on a single coverage level and do not ensure consistency across multiple confidence levels. In many real-world applications, such as weather forecasting, macroeconomic prediction, and risk management, different users operate under heterogeneous risk tolerances and require calibrated uncertainty estimates across a range of coverage levels. In such settings, it is desirable to produce prediction sets corresponding to different coverage levels that are nested and valid simultaneously. In this paper, we propose two novel online conformal prediction methods that output nested prediction sets across a range of coverage levels, enabling simultaneous uncertainty quantification across the entire risk spectrum. Beyond interpretability, jointly estimating multiple coverage levels is known to improve statistical efficiency in classical quantile regression by enforcing non-crossing constraints and sharing information across quantiles. Our approaches leverage an online optimization perspective with small regret that translates to quantile estimation error control while enforcing nestedness of prediction sets. Empirical results on synthetic and real-world datasets, including applications in forecasting tasks with heterogeneous risk requirements, demonstrate that our method achieves stable coverage across all levels, strictly nested prediction sets, and improved efficiency compared to existing online conformal baselines.