---
title: Valid Prediction Time (VPT)
url: https://www.emergentmind.com/topics/valid-prediction-time-vpt
type: topic
---

# Valid Prediction Time (VPT)

Searching arXiv for recent papers on Valid Prediction Time and related conformal risk-control formulations.
Valid Prediction Time (VPT) denotes a prediction horizon or lower bound certified by an explicit validity criterion. In recent arXiv literature, the term appears in several technically distinct forms: an anytime-valid guarantee that a sequentially updated prediction set never exceeds a target risk level with probability at least \(1-\delta\) [2602.04364]; a vessel-trajectory reliability horizon \(a_{90/95}\), defined as the largest forecast horizon whose \(95\%\) lower-confidence bound on the probability of accurate prediction remains at least \(0.90\) [2508.14198]; the first time at which normalized forecast error in a chaotic system crosses a fixed threshold [2508.06730]; and a covariate-specific lower bound on a right-censored survival time satisfying training-set conditional validity [2505.02213]. This suggests that VPT is best viewed as a family of validity-qualified horizon concepts rather than a single universal metric.

## 1. Terminological scope and formal variants

The literature summarized here uses the same term for several non-equivalent mathematical objects. What unifies them is not a common formula but a common role: each VPT is a horizon up to which a prediction remains valid under a stated criterion.

| Setting | Formal object called VPT | Validity criterion |
|---|---|---|
| Anytime-valid conformal risk control | Implicitly, an infinite valid prediction time with high probability | \(\Pr[\forall n:\;E(\ell(C_n(X),Y)\mid \mathcal Z_n)\le \alpha]\ge 1-\delta\) |
| Vessel trajectory prediction | \(a_{90/95}\) | Largest \(h\) such that the \(95\%\) lower bound on \(\mathrm{POAP}(h)\) is at least \(0.90\) |
| Reservoir computing for Lorenz | Threshold-crossing time | \(\mathrm{VPT}=\min\{t>0\mid E(t)\ge 0.4\}\) |
| Right-censored survival | Lower prediction bound \(\hat L(x)\) | \(\Pr(T\ge \hat L(x)\mid \text{training data},X=x)\ge 1-\alpha\) |

In the sequential conformal setting, VPT is best understood as the duration over which a calibrated predictor continues to satisfy a risk constraint at every sample size. In vessel trajectory prediction, VPT is a reliability horizon attached to a detection-style confidence analysis. In chaotic forecasting, it is the first failure time of a normalized trajectory error. In right-censored survival analysis, VPT is not a time of future breakdown but a lower predictive bound for an event time, tailored to covariates and endowed with conditional validity.

## 2. Anytime-valid conformal risk control

In the sequential conformal formulation, calibration points \((X_1,Y_1),\dots,(X_n,Y_n)\) are observed i.i.d. from \(P\), possibly with an importance-weight \(\omega\) for test-time shift. One fixes a loss \(\ell(C(X),Y)\in[0,B]\) that is monotone in the size of the prediction set \(C\), a target risk level \(\alpha\), and a nested family of sets \(\{C_\lambda(X)\}_{\lambda\in\Lambda}\). Calibration produces a data-dependent \(\lambda_n\), and the output is \(C_n(X):=C_{\lambda_n}(X)\). The defining guarantee is
\[
\Pr\Bigl[\forall\,n\in\mathbb N:\;E\bigl[\ell(C_n(X),Y)\mid \mathcal Z_n\bigr]\le \alpha\Bigr]\ge 1-\delta,
\]
where \(\mathcal Z_n=\{(X_i,Y_i)\}_{i\le n}\). The associated empirical risk is
\[
R_n(\lambda)=\frac1n\sum_{i=1}^n \ell(C_\lambda(X_i),Y_i).
\]
The key distinction from standard conformal formulations is that error is not merely controlled on average over many possible calibration datasets of fixed size; it remains valid with high probability over a cumulatively growing calibration dataset at any time point [2602.04364].

The main theorem assumes that \(\ell(\cdot)\in[0,B]\) is monotone and right-continuous in \(\lambda\), and introduces the stitched correction term
\[
f_{B,m,\delta}(v)=1.44\,\sqrt{v\,\Bigl[2\ln\bigl(\ln_2(v/m)+1\bigr)+\ln\!\tfrac{\pi^2}{6\delta}\Bigr]}
+2.42\,B\;\Bigl[2\ln\bigl(\ln_2(v/m)+1\bigr)+\ln\!\tfrac{\pi^2}{6\delta}\Bigr].
\]
With
\[
m^*=\min\Bigl\{m\in\mathbb N:\; f_{B,m,\delta}\bigl(\alpha(B-\alpha)\,m\bigr)/m\le \alpha\Bigr\},
\]
the update rule is
\[
\lambda_n=\inf\Bigl\{\lambda\in\Lambda:\;R_n(\lambda)\le \alpha-\tfrac{f_{B,m^*,\delta}\bigl(\alpha(B-\alpha)\,n\bigr)}{n}\Bigr\}.
\]
The resulting sets satisfy
\[
\Pr\Bigl[\forall\,n:\;E\bigl[\ell(C_n(X),Y)\mid \mathcal Z_n\bigr]\le \alpha\Bigr]\ge 1-\delta.
\]

For split conformal miscoverage, the loss becomes the \(0\)–\(1\) loss \(\mathbf 1\{y\notin C\}\), so \(B=1\), and the guarantee specializes to
\[
\Pr\bigl[\forall n:\;\Pr(Y\notin C_n(X)\mid \mathcal Z_n)\le \alpha\bigr]\ge 1-\delta.
\]
The formulation therefore turns VPT into a statement about indefinite sequential validity: with probability at least \(1-\delta\), the valid prediction time is infinite.

## 3. Concentration, tightness, and online operation

The proof architecture is quantile-based and handles discontinuities by a randomized mixture at the threshold \(\hat\lambda(\alpha)\). For each \(\lambda\), the losses
\[
\xi_i(\lambda)=\ell(C_\lambda(X_i),Y_i)\in[0,B]
\]
are right-continuous and nonincreasing in \(\lambda\). Let \(\hat\lambda(\alpha)\) be the infimum of \(\lambda\) with \(R(\lambda)\le \alpha\), and define a randomized mixture at any discontinuity so that
\[
E\bigl[\pi\,\xi_i(\hat\lambda)+(1-\pi)\,\xi_i(\hat\lambda^-)\bigr]=\alpha.
\]
This yields the zero-mean process
\[
M_t(\alpha)=\sum_{i=1}^t\bigl(\alpha-\xi_i(\alpha)\bigr),
\]
whose increments lie in \([\alpha-B,\alpha]\). The variance proxy is \(V_t=\alpha(B-\alpha)\,t\) and the scale parameter is \(c=B\). A time-uniform stitched boundary gives
\[
\Pr\Bigl[\forall t:\;M_t(\alpha)<f_{B,m,\delta}\bigl(\max\{m,V_t\}\bigr)\Bigr]\ge 1-\delta.
\]
Rearranging the boundary and using the definition of \(\lambda_n\) shows \(\lambda_n\ge \hat\lambda(\alpha)\) on the same event, hence \(R_n(\lambda_n)\le \alpha\) for all \(n\) [2602.04364].

The same work proves a matching lower bound in the sense of asymptotic tightness. If the random leaps of \(R_n(\lambda)\) vanish,
\[
d_n=\sup_{\lambda\in\Lambda}\;\lim_{\epsilon\to0^+}\bigl[R_n(\lambda^-)-R_n(\lambda)\bigr]\xrightarrow{a.s.}0,
\]
then with probability \(1-2\delta\),
\[
\forall n:\quad \alpha-k_n\le E\bigl[\ell(C_n(X),Y)\mid \mathcal Z_n\bigr]\le \alpha,
\]
with
\[
k_n=d_{g(n)}+\gamma_{g(n)}'(\delta_{\lfloor\log n\rfloor}),\qquad
g(n)=2^{\lfloor\log_2 n\rfloor},\qquad
\gamma_n'(\delta)=O\Bigl(\sqrt{\tfrac{1}{n}\ln\tfrac1\delta}\Bigr).
\]
In particular,
\[
\alpha-E[\ell]\;=\;O\Bigl(\max\{d_n,\;B\,\sqrt{n^{-1}\ln\ln n}\}\Bigr).
\]
No faster time-uniform rate than \(\sqrt{\ln\ln n/n}\) can be achieved, so the rates are asymptotically optimal.

The online implementation is correspondingly simple. On each new calibration point, one updates the variance proxy \(V=\alpha(B-\alpha)n\), computes \(\Delta=f_{B,m^*,\delta}(\max\{m^*,V\})/n\), and selects the smallest \(\lambda\) on a discrete grid \(\Lambda\) such that the empirical risk is at most \(\alpha-\Delta\), with an additional monotonicity step \(\lambda_n=\min\{\lambda_n,\lambda_{\mathrm{prev}}\}\). For continuous \(\Lambda\), the infimum can be implemented by sorting scores, as in split conformal, or by binary search.

A synthetic regression experiment illustrates the operational meaning of VPT. The setup is \(Y=2X+\epsilon\), \(\epsilon\sim N(0,1)\), \(X\sim \mathrm{Uniform}[-3,3]\), with target miscoverage \(\alpha=5\%\), error probability \(\delta=10\%\), score \(s(x,y)=|y-f(x)|\), and prediction sets \(C_\lambda(x)=\{y:s(x,y)\le \lambda\}\). Over \(500\) independent runs, approximately \(90\%\) of runs remain green, matching the nominal \(1-\delta\), whereas standard split conformal has almost every run exceed \(\alpha\) at some \(n\), despite controlling marginal miscoverage on average. The realized miscoverage curves remain below \(\alpha\) and approach it from below as \(n\) grows.

## 4. Reliability-based VPT in vessel trajectory prediction

In vessel trajectory prediction, the paper does not introduce a separate symbol “VPT” but instead uses the quantity \(a_{90/95}\) as a time-horizon reliability metric. It is defined as the largest prediction horizon \(h\), in minutes, for which the \(95\%\) lower-confidence bound on the probability of making an error below a chosen threshold still exceeds \(90\%\). If \(\mathrm{POAP}(h)\) denotes the probability that the displacement error at horizon \(h\) stays below a decision threshold \(\hat a_{\mathrm{th}}\), the construction is
\[
\hat a=b+m\cdot h+\epsilon,\qquad \epsilon\sim N(0,\tau^2),
\]
\[
\mathrm{POAP}(h)=P[\hat a<\hat a_{\mathrm{th}}]
=\Phi\!\Bigl(\frac{\hat a_{\mathrm{th}}-(b+m\cdot h)}{\tau}\Bigr),
\]
and
\[
a_{90/95}:=\max\{h:\; (95\% \text{ lower-bound on }\mathrm{POAP}(h))\ge 0.90\}.
\]
The \(95\%\) lower-confidence bound is obtained by the Wald method [2508.14198].

The step-by-step procedure is explicit. Predicted trajectories are interpolated to \(3\,\mathrm{s}\) steps, displacement error is computed at each minute-level horizon \(h=1\,\mathrm{min},2\,\mathrm{min},\dots\), mean error versus horizon is fitted with the linear regression above by maximum likelihood, a decision threshold \(\hat a_{\mathrm{th}}=20\,\mathrm{m}\) is chosen, \(\mathrm{POAP}(h)\) and its \(95\%\) lower bound are computed, and \(a_{90/95}\) is read off as the largest admissible horizon. The reported VPT is therefore threshold-specific and confidence-specific.

The numerical results show clear scenario dependence. For the overall dataset of \(2{,}578\) samples, the VPTs are \(2.452\,\mathrm{min}\) for STT-R-CSCT, \(3.190\,\mathrm{min}\) for N-CSCT, and \(3.541\,\mathrm{min}\) for GMM-Trans-GRU. In Encounter-3, with \(24\) samples, the values are \(3.116\,\mathrm{min}\), \(4.767\,\mathrm{min}\), and \(>5.000\,\mathrm{min}\), respectively. In Encounter-2, with \(355\) samples, they are \(2.707\,\mathrm{min}\), \(3.583\,\mathrm{min}\), and \(4.134\,\mathrm{min}\). The only row in which GMM-Trans-GRU is not best is Encounter-1 & Overtaken-1, where STT-R-CSCT achieves \(1.490\,\mathrm{min}\), N-CSCT \(1.422\,\mathrm{min}\), and GMM-Trans-GRU \(1.375\,\mathrm{min}\).

These values are intended for operational use. A navigation-decision module can restrict itself to predictions with horizon \(h\le \mathrm{VPT}\), thereby guaranteeing, with \(95\%\) confidence, an error below the chosen threshold. Stratifying by traffic complexity permits adaptive horizon selection for collision avoidance, route planning, or pilot-assist systems.

## 5. Chaotic systems and reservoir computing

For chaotic dynamical systems, VPT is defined as a threshold-crossing time of a normalized trajectory error. Let \(u_{\rm true}(t)\in\mathbb R^D\) be the ground-truth state and \(u_{\rm pred}(t)\) the predicted state. The normalized mean-squared error is
\[
E(t)=\frac{\sum_{i=1}^D \bigl[u_{{\rm true},i}(t)-u_{{\rm pred},i}(t)\bigr]^2}{\mathrm{Var}(u_{\rm true})},
\]
with
\[
\mathrm{Var}(u_{\rm true})=
\frac{1}{N}\sum_{n=1}^N\sum_{i=1}^D
\Bigl[u_{{\rm true},i}(n\tau)-\overline{u_i}\Bigr]^2,\qquad
\overline{u_i}=\frac1N\sum_{n=1}^N u_{{\rm true},i}(n\tau).
\]
The VPT is
\[
\mathrm{VPT}=\min\{t>0\mid E(t)\ge 0.4\},
\]
measured in units of the largest Lyapunov time \(T_\lambda=1/\lambda_{\max}\). The threshold \(0.4\) is chosen because once \(E(t)\) exceeds it, the trajectory has typically jumped onto the wrong lobe of the Lorenz attractor [2508.06730].

The relation to chaos is explicit. If \(\lambda_{\max}\) is the maximal Lyapunov exponent, nearby trajectories diverge on average like \(e^{\lambda_{\max} t}\), so \(E(t)\) scales approximately as \(E(0)e^{2\lambda_{\max} t}\). In practice, this means that VPT is controlled by the initial prediction error and the known \(\lambda_{\max}\approx 0.9056\) of the Lorenz system. The same argument motivates an efficient shortcut: VPT can be estimated from the error after only a small number of autonomous prediction steps, avoiding full rollout during hyperparameter search.

The paper also introduces Valid Ground Truth Time (VGTT), defined as the time up to which independent ODE solvers agree to within the same error threshold. Using the 5th-order Adams–Bashforth–Moulton method with \(10^{-12}\) tolerances and \(\tau=10^{-3}\), VGTT is approximately \(30\mbox{–}40\) Lyapunov times for the Lorenz system. A reported VPT exceeding VGTT is not meaningful, because beyond VGTT the reference trajectory itself becomes solver-dependent.

The reservoir computer is an echo-state architecture with state update
\[
\mathbf r_{n+1}=f\!\bigl(A\,\mathbf r_n+W_{\rm in}\,z_n\bigr),
\]
readout
\[
\hat z_{n+1}=W_{\rm out}\,\mathbf r_{n+1},
\]
and ridge objective
\[
\sum_{n=1}^{T/\tau}\|\hat z_n-z_n\|^2+\lambda\,\mathrm{Tr}(W_{\rm out}W_{\rm out}^T),
\]
whose closed-form solution is
\[
W_{\rm out}=Z\,R^T\,(R\,R^T+\lambda I)^{-1}.
\]
As \(\lambda\to 0\), both the reservoir computer’s VPT and the benchmark VPT increase and asymptote to a limit; for \(0<\lambda<10^{-6}\), the reservoir computer achieves approximately \(70\%\) of the benchmark VPT. Maximum VPT values exceed \(30\) Lyapunov times.

A two-regime structure appears in the spectral radius \(\rho\). One regime is near zero, \(\rho\approx 10^{-4}\), essentially independent of reservoir size \(N\), where the reservoir behaves like an extreme learning machine: stable, with little temporal memory, but recoverable by sufficiently large \(N\). The other is an edge-of-chaos regime at \(\rho\approx O(1)\), which retains more memory, slightly outperforms the near-zero regime, and requires more delicate tuning in \((N,\rho)\).

## 6. Right-censored survival outcomes

In right-censored survival analysis, VPT is a covariate-specific lower prediction bound for the event time. The observed data are i.i.d. copies of
\[
O_i=(W_i,\Delta_i,Y_i),\qquad Y_i=\min(T_i,C_i),\qquad \Delta_i=\mathbf 1\{T_i\le C_i\},
\]
where \(W\in\mathbb R^p\) are covariates, \(T\ge 0\) is the event time, and \(C\ge 0\) is the censoring time. A Valid Prediction Time is a function \(\hat L(x)\) such that, for a new subject with covariates \(X=x\),
\[
\Pr\bigl(T\ge \hat L(x)\mid \text{training data},X=x\bigr)\ge 1-\alpha.
\]
The paper studies asymptotic training-set conditional validity, also called APAC, after splitting the sample into a training set of size \(m\) and a calibration set of size \(n\) [2505.02213].

The target coverage functional is
\[
\Psi(P;L)=\Pr_{(W,T)\sim P}\bigl(T>L(W)\bigr)
=E_P\{S(L(W)\mid W)\},
\]
identified under two conditions: \(T\perp C\mid W\) on \([0,t_0]\), and \(\Pr(C\ge t_0\mid W)>0\). The method uses a semiparametric one-step estimator. First, nuisance estimators \((S_n,G_n)\) are fitted on the training data. Next, a candidate VPT is defined by
\[
L_{n,\tau}(w)=\min\bigl\{S_n^{-1}(1-\tau\mid w),\;G_n^{-1}(\eta_2\mid w)\bigr\},
\]
for a small \(\eta_2\). On the calibration set, one computes the plug-in coverage and corrects it by the average efficient influence function:
\[
\hat\psi_{n,\tau}
=\Psi_\tau(\hat P_n;L_{n,\tau})
+\frac1n\sum_{i=1}^n D(\hat P_n,G_n,S_n;L_{n,\tau})(O_i).
\]

The final lower bound is selected by a confidence-based thresholding rule. For each \(\tau\) in a grid \(\mathcal T_n\subset[0,1]\), an estimated standard error \(\hat\sigma_{n,\tau}\) is formed, and the Wald lower confidence bound is
\[
\mathrm{CLB}_{n,\tau}
=\hat\psi_{n,\tau}-z_{1-\beta}\,\frac{\hat\sigma_{n,\tau}}{\sqrt n}.
\]
One then chooses the largest \(\hat\tau\) such that, for all \(\tau'\le \hat\tau\), \(\mathrm{CLB}_{n,\tau'}\ge 1-\alpha\), and outputs \(\hat L_n(x)=L_{n,\hat\tau}(x)\). The main efficiency theorem gives the asymptotic linear expansion
\[
\hat\psi_{n,\tau}
=\Psi(P_0;L_{n,\tau})
+\frac1n\sum_{i=1}^n D(P_0,G_0,S_0;L_{n,\tau})(O_i)
+o_p(n^{-1/2}),
\]
uniformly over \(\tau\in[0,1]\), implying asymptotic normality. The threshold selection then yields
\[
\Pr_{\mathcal D}\{\Pr_{(W,T)}(T\ge \hat L_n(W)\mid \mathcal D)\ge 1-\alpha\}\to 1-\beta.
\]

The empirical study reports six synthetic settings spanning univariate versus multivariate covariates, independent versus covariate-dependent censoring, and moderate versus high dimension. In Setting 3 with \(\alpha=0.1\), TCsurv has coverage \(0.902\), average VPT \(5.23\), and standard deviation \(0.45\); DFT-fixed has \(0.938\), \(5.95\), \(0.60\); DFT-adaptive-T has \(0.884\), \(5.10\), \(0.50\); and DFT-adaptive-CT has \(0.890\), \(5.15\), \(0.48\). The real-data application tracks users’ active times on a mobile application, with \(N\approx 2500\), covariates including gender, age, and number of children, and outcome defined as time to the ninth active day with artificial right-censoring at days \(11\)–\(12\). The implementation cost is dominated by nuisance fitting on the training set, while calibration corrections scale as \(O(n|\mathcal T_n|)\) with typically \(|\mathcal T_n|\approx 50\). The method is available as an R package and uses survSuperLearner for \((S_n,G_n)\).

## 7. Comparative interpretation and limitations

The four formulations reveal that VPT is domain-specific both in its mathematical definition and in its operational meaning. In sequential conformal risk control, it is a high-probability statement that a data-conditional risk constraint holds for every calibration size; in vessel trajectory prediction, it is the largest horizon supported by a lower confidence bound on \(\mathrm{POAP}(h)\); in chaotic prediction, it is the first threshold-crossing time of normalized error; and in right-censored survival analysis, it is a predictive lower bound on the event time itself [2602.04364] [2508.14198] [2508.06730] [2505.02213].

This suggests several common misconceptions should be avoided. A longer VPT is not directly comparable across these literatures, because each construction is tied to a distinct threshold or validity convention: \(\alpha\) and \(\delta\) in anytime-valid conformal risk control, \(\hat a_{\mathrm{th}}=20\,\mathrm{m}\) together with the \(90/95\) rule in vessel prediction, the error threshold \(0.4\) and Lyapunov-time normalization in chaotic forecasting, and the pair \((\alpha,\beta)\) together with censoring assumptions in survival analysis. Likewise, “infinite” VPT in the conformal setting does not mean zero error; it means that the data-conditional risk remains below a target level at all times with probability at least \(1-\delta\). In chaotic systems, a VPT larger than VGTT is explicitly not meaningful. In survival analysis, APAC validity depends on \(T\perp C\mid W\), positivity, and sufficiently fast nuisance-rate conditions, whereas the marginal guarantee follows from a simpler threshold rule.

Taken together, these works show that VPT functions as a calibrated horizon concept. Its substantive interpretation depends on what is being guaranteed: bounded risk, bounded displacement error, bounded trajectory divergence, or lower-bounded survival time.

Source: https://www.emergentmind.com/topics/valid-prediction-time-vpt