The World-Model Evaluation Ladder › One-step prediction

Michael BrenndoerferJuly 31, 20267 min read
World Models HandbookAll chapters
Part 11 of 12Evaluation and Understanding
Chapter 61The World-Model Evaluation Ladder › One-step prediction
Reading section 2 of 9 · One-step prediction · 7 minSection 2 of 9 · 7 min

Book contents

Section 2 of 9 · 7 min read

The World-Model Evaluation Ladder › One-step predictionOne section · Sections in this chapter2 / 9 ↓
Reading tools Definitions & inline links

Choose your expertise level to adjust how many terms are explained. Beginners see more tooltips, experts see fewer to maintain reading flow. Hover over underlined terms for instant definitions.

Article links

Make inline references clickable

The one-step test scores a prediction at an observed input without feeding the model's previous prediction back into it. That makes it convenient for regression-style models, but the resulting score does not test recursive simulation or action selection. Teacher forcing also does not make transitions statistically independent; neighboring observations can remain correlated.

Let T(ds′∣s,a)T(ds'\mid s,a) be a controlled Markov transition kernel. It returns a probability law over successors, not a successor state vector. Let DD be a common evaluation distribution over (s,a)(s,a), and let m(s,a)=ET[s′∣s,a]m(s,a)=\mathbb{E}_T[s'\mid s,a] be the conditional mean. For a deterministic point predictor ff with finite squared risk, assume the successor and predictor have finite second moments under this joint law. Then

L1(f)=E(s,a)∼D,  s′∼T(⋅∣s,a)∥f(s,a)−s′∥2=ED∥f(s,a)−m(s,a)∥2+σT,D2,σT,D2=EDET[∥s′−m(s,a)∥2∣s,a].\begin{aligned} L_1(f)&=\mathbb{E}_{(s,a)\sim D,\;s'\sim T(\cdot\mid s,a)} \|f(s,a)-s'\|^2\\ &=\mathbb{E}_{D}\|f(s,a)-m(s,a)\|^2+\sigma^2_{T,D},\\ \sigma^2_{T,D}&=\mathbb{E}_{D}\mathbb{E}_T[ \|s'-m(s,a)\|^2\mid s,a]. \end{aligned}

The first term measures conditional-mean prediction error. The second averages conditional squared noise over the same evaluation law. It depends on TT and DD, but not on the deterministic point predictor. To obtain the identity, write f−s′=(f−m)+(m−s′)f-s'=(f-m)+(m-s') and expand the squared norm. The conditional expectation of the cross term is zero because ET[s′−m∣s,a]=0\mathbb{E}_T[s'-m\mid s,a]=0. Noise is already part of L1L_1; it is not added again to an expected squared error against sampled successors.

Where:

  • ss and aa: the joint state-action input drawn from DD;
  • s′s': the sampled successor drawn from T(⋅∣s,a)T(\cdot\mid s,a);
  • m(s,a)m(s,a): the conditional mean successor;
  • f(s,a)f(s,a): a deterministic point prediction;
  • σT,D2\sigma^2_{T,D}: the average conditional squared noise norm.

Training data D={(si,ai,si′)}i=1N\mathcal D=\{(s_i,a_i,s'_i)\}_{i=1}^N may be collected by a behavior policy and a state-estimation pipeline. The evaluation law can be chosen to match that collection regime or to represent a different use. The conditional-mean decomposition applies to the selected law; it does not establish that this law covers deployment.

On MM held-out transitions the empirical score is

L^1(f)=1M∑j=1M∥f(sj,aj)−sj′∥2.\widehat L_1(f)=\frac1M\sum_{j=1}^M\|f(s_j,a_j)-s'_j\|^2.

Here jj indexes the evaluation transitions, MM is their count, and each sj′s'_j is the observed successor of input (sj,aj)(s_j,a_j). With independent evaluation draws from the specified joint law, conditional on a fixed fitted predictor, this sample average is unbiased for its population risk. Disjoint row indices alone do not establish independence from the fitted predictor or independence among evaluation rows. For trajectory data, the sampling and splitting unit needs to be stated.

Three quantities help interpret the score. Let F\mathcal F be a fixed hypothesis class, let fF∗f^*_{\mathcal F} minimize L1L_1 over that class when such a minimizer exists, and let f^∈F\hat f\in\mathcal F be the fitted model.

  • Approximation error is L1(fF∗)−σT,D2L_1(f^*_{\mathcal F})-\sigma^2_{T,D}, the best conditional-mean error available in the class.
  • Fitting excess is L1(f^)−L1(fF∗)L_1(\hat f)-L_1(f^*_{\mathcal F}). Finite-data estimation and imperfect optimization can both contribute to it.
  • Evaluation error is L^1(f^)−L1(f^)\widehat L_1(\hat f)-L_1(\hat f), a signed discrepancy due to the evaluation sample.

These definitions give an exact bookkeeping identity when the minimizer exists:

L^1(f^)=σT,D2+[L1(fF∗)−σT,D2]+[L1(f^)−L1(fF∗)]+[L^1(f^)−L1(f^)].\widehat L_1(\hat f)=\sigma^2_{T,D} +[L_1(f^*_{\mathcal F})-\sigma^2_{T,D}] +[L_1(\hat f)-L_1(f^*_{\mathcal F})] +[\widehat L_1(\hat f)-L_1(\hat f)].

The final bracket can be negative; the first two excess terms cannot. More training data does not remove approximation error in a fixed misspecified class. Under suitable sampling, identifiability, and consistency assumptions it can reduce estimation excess over the best member, even when the class is misspecified.

For example, let Y=XY=X with E[X]=0\mathbb E[X]=0 and Var⁡(X)=1\operatorname{Var}(X)=1, but allow only constant predictors fc(x)=cf_c(x)=c. The best constant is zero and its approximation risk is one. Fitting the sample mean on NN independent draws gives expected excess risk 1/N1/N. The approximation error remains while estimation excess vanishes. There is no universal rule that greater capacity always increases estimation error; that tradeoff depends on the class, algorithm, regularization, and data.

For a comparison on common evaluation transitions, uncertainty belongs to the paired loss differences. The marginal standard error of one model's MSE is not the standard error of the difference. A numeric gap such as 0.0050.005 cannot be judged from M=800M=800 alone.

Finally, averaging signed prediction errors is not the same as averaging squared errors. Positive and negative mistakes do not cancel in MSE. A conditional mean can be the optimal point predictor yet fail to represent multimodal futures needed by a particular planner. That is a reason to evaluate a predictive distribution or the downstream task separately, not to claim that smoothing is automatically low-error.

A small two-component oscillator makes the distinction concrete without copying the answer into the predictor. Its true state is zt=(xt,vt)z_t=(x_t,v_t), and its transition rotates that vector by ω=2π/20\omega=2\pi/20. The same misspecified model is used in both tests: it rotates by the correct angle but expands the state by a factor of 1.021.02. Teacher forcing resets each input to the true state; open-loop rollout feeds the model's previous prediction back into itself.

In[3]:
Code
import numpy as np

t = np.arange(41)
omega = 2 * np.pi / 20
rotation = np.array(
    [[np.cos(omega), np.sin(omega)], [-np.sin(omega), np.cos(omega)]]
)
true_z = np.column_stack((np.sin(omega * t), np.cos(omega * t)))
oscillator_model = 1.02 * rotation
teacher_z = (oscillator_model @ true_z[:-1].T).T
rollout_z = np.zeros_like(true_z)
rollout_z[0] = true_z[0]
for step in range(40):
    rollout_z[step + 1] = oscillator_model @ rollout_z[step]
teacher_mse = float(np.mean(np.sum((teacher_z - true_z[1:]) ** 2, axis=1)))
oscillator_rollout_rmse = float(
    np.sqrt(np.mean(np.sum((rollout_z[1:] - true_z[1:]) ** 2, axis=1)))
)
constant_mean_mse = float(np.mean((true_z[:, 0].mean() - true_z[:, 0]) ** 2))
Out[4]:
Console
Teacher-forced vector MSE: 0.000400
Open-loop vector RMSE: 0.643858
Constant-mean scalar MSE over all 41 states: 0.487805

The teacher-forced vector MSE is 0.022=0.00040.02^2=0.0004, because every true state has norm one. In open loop the norm becomes 1.02t1.02^t, so the error grows even though the rotation angle is correct. A constant-mean predictor does not track this sine wave under teacher forcing: its scalar MSE across the 41 plotted states is 20/41≈0.48780520/41\approx0.487805, not zero.

Out[5]:
Visualization
True sine wave, near-true teacher-forced predictions, and expanding model rollout.
The same slightly expansive rotation model is used for teacher-forced predictions and open-loop rollout. Resetting each input to the true two-component state limits one-step error, while recursive use expands the oscillation amplitude. The plot shows the first state component.

Three confounds deserve special attention.

  1. Distribution coverage. The evaluation law DD need not match the state-action occupancy induced by deployment. If data cover only near-upright robot states, predictions far outside that region are extrapolations unless other data or structural constraints inform them. A state absent from a finite sample is not necessarily outside the support of the population distribution.
  2. Stochasticity. Point-prediction MSE and distributional scores answer different questions. A probabilistic model can be scored by its mean for MSE, or by a suitable proper scoring rule for its distribution. Comparing a random sample from one model with a point estimate from another changes the scoring convention. Uncertainty-aware planning needs separate tests. [Perceptual and Generative Evaluation] is planned to examine observation-quality scoring, not certify a planning system.
  3. Observation error. A state estimated from an image or point cloud can be inaccurate before the dynamics predictor receives it. Separating perception and transition errors needs component-specific measurements, identifiable model/noise assumptions, or interventions. Aggregate prediction error alone does not distinguish those components. Their units, directions, dependence, and propagation matter; two nominal error magnitudes alone do not determine the combined prediction error.

A one-step result is therefore a distribution-specific check. It can be decisive if that is the intended task. It cannot, by itself, establish rollout quality, intervention accuracy, or useful control.