research29 min studyLesson 861 in Engineering

Predictive versus parameter uncertainty

Decompose uncertainty within a declared predictive model while separating posterior approximation, data variation, and model misspecification.

On this page

What you’ll understand

  • Apply total variance to a specified posterior predictive mixture.
  • Distinguish parameter disagreement from future-outcome variation.
  • Explain why an uncertainty decomposition depends on the model and its approximation.

A regression service can return the same predictive variance for very different reasons: noisy observations within each plausible model, or disagreement between plausible models. These reasons suggest different experiments, but neither can be read directly from a single confidence score. Begin with an explicit joint predictive model for parameters and outcomes. The resulting decomposition describes that model’s uncertainty; it does not certify that the model family, posterior approximation, or deployment distribution is correct.

Name the random quantities before naming the uncertainty

Condition on a dataset D and an input x. In a Bayesian description, a parameter Θ is distributed according to a posterior, and the future outcome Y is distributed conditionally on Θ. Integrating over Θ produces a predictive distribution. An interval for a parameter and an interval for a future observation answer different questions, even when both use posterior samples. For scalar regression, total variance separates expected conditional outcome variance from variance of conditional means, provided the second moments exist.

The first component is often called aleatoric uncertainty and the second epistemic uncertainty. Those labels require context. A sensor effect treated as noise in one model could become predictable after adding a measured covariate in another. Conversely, a posterior can concentrate tightly around the wrong model. A small parameter component then reflects agreement within a restricted family, while systematic prediction errors remain large. Record the variables observed, likelihood, prior, fitting procedure, and units before interpreting either component as reducible or irreducible.

Var⁡(Y∣x,D)=EΘ∣D[Var⁡(Y∣x,Θ)]+Var⁡Θ∣D[E(Y∣x,Θ)]\operatorname{Var}(Y\mid x,D)=\mathbb E_{\Theta\mid D}[\operatorname{Var}(Y\mid x,\Theta)]+\operatorname{Var}_{\Theta\mid D}[\mathbb E(Y\mid x,\Theta)]
Total variance for a scalar outcome with finite second moment under the declared conditional model. Both terms have squared outcome units. The likelihood is assumed conditionally independent of D given x and Θ; the identity is not a claim that this model equals the deployment process.↔ Scroll wide equations horizontally.

A mixture can hide the source of its uncertainty

For binary outcomes, suppose each model assigns probability pΘ to label one. The marginal predictive probability is the posterior mean of pΘ. Its Bernoulli variance can be decomposed into expected within-model Bernoulli variance and variance across those probabilities. Models that all say one half and models split between zero and one have identical single-outcome marginal predictions, yet opposite allocations between the two terms. The decomposition therefore needs access to the chosen distribution over models, not just the averaged output.

Deep ensembles and repeated dropout passes can expose useful disagreement, but disagreement is not automatically an exact posterior variance or a reliable detector of every unfamiliar input. Shared architecture, training data, and optimization choices can make all members confidently wrong together. Approximate Bayesian interpretations also have assumptions about priors, objectives, and variational families. Evaluate whether a proposed uncertainty score predicts the errors relevant to the application, and compare it against simple baselines rather than promoting the method name into a guarantee.

pˉ=E[pΘ],pˉ(1−pˉ)=E[pΘ(1−pΘ)]+Var⁡(pΘ)\bar p=\mathbb E[p_\Theta],\qquad \bar p(1-\bar p)=\mathbb E[p_\Theta(1-p_\Theta)]+\operatorname{Var}(p_\Theta)
Conditional on the same x and D, pΘ is a dimensionless probability in [0,1]. This identity describes one binary outcome under a specified mixing distribution. Repeated outcomes sharing a latent parameter need not be independent Bernoulli trials.↔ Scroll wide equations horizontally.

Choose an intervention that can distinguish explanations

When disagreement is high, collecting labels near the disputed inputs may help distinguish candidate models. When every model predicts a broad outcome distribution, better measurements or new features may matter more than adding nearly identical examples. Test these hypotheses by tracking held-out prediction quality after targeted interventions. Preserve the before-and-after data definitions: a narrower reported interval caused by a different posterior approximation is not necessarily an improvement in accuracy, coverage, or useful decision quality.

The probability instrument exposes only a marginal law. Reset, set p to 0.5 and n to 10, predict mean five and variance 2.5, then Step. Compare a fresh Reset with both p at 0.5 and n at 20 before Step: mean ten, variance five. These are exact independent-binomial counts. The lab has no posterior, model ensemble, or shared latent parameter and cannot identify either uncertainty component from the displayed probability.

Separate two kinds of variance in an explicit mixture

Two hypothetical regression models have posterior weights one half each. Their conditional means are 8 and 12 milliseconds, and their conditional variances are 4 and 9 square milliseconds. Compute the predictive mean and both variance components without assuming the mixture is Gaussian.

  1. The weighted mean is (8 + 12)/2 = 10 milliseconds. This is a mean of predictions, not a confidence bound for the parameter.

  2. Expected within-model variance is (4 + 9)/2 = 6.5 square milliseconds. Variance of the means is ((8 − 10)² + (12 − 10)²)/2 = 4 square milliseconds.

  3. Total predictive variance is 6.5 + 4 = 10.5 square milliseconds. Its standard deviation is approximately 3.2404 milliseconds, but a Gaussian interval needs an additional distributional assumption.

The result

The arithmetic identifies components within the specified mixture while leaving model adequacy as an empirical question.

Explore the model

Interactive lab

Interactive lab · probability

A marginal probability does not reveal a posterior

Reset, set p = 0.5 and n = 10, predict mean 5 and variance 2.5, then Step. Compare a fresh Reset with both p = 0.5 and n = 20 before Step: mean 10, variance 5. Exact independent-binomial masses contain no posterior or ensemble decomposition.

A marginal probability does not reveal a posteriorThese are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Current readings: Mean · np: 5; Variance · np(1−p): 2.5; Total probability: 1. All parameters are adjustable with the labeled controls below.03581000.070.150.220.3Successes kP(X = k)k = 0 · P = 9.77e-4
Exact probabilityHighlighted outcomeMean np
Mean · np
5
Variance · np(1−p)
2.5
Total probability
1
t = 0 s

These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forIdentical marginal probabilities can arise from different model-uncertainty explanations.

Check your understanding

All ensemble members assign probability 0.5 to one binary label. Which statement follows within that ensemble?

Try it yourself

Two equally weighted scalar models predict means 1 and 5 meters, each with conditional variance 4 square meters. Compute the predictive mean and total variance. What changes if all models are biased upward by two meters?

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Average the conditional variances and add the variance of the two conditional means. A common shift changes the mean without changing either variance component.

Compare with the solution

The predictive mean is 3 meters. Expected conditional variance is 4 and variance of the means is ((1 − 3)² + (5 − 3)²)/2 = 4, so total variance is 8 square meters. A shared upward bias of two meters moves the mean to 5 meters while both variance terms remain unchanged. Consequently these uncertainty summaries alone cannot diagnose that systematic bias; compare predictions with independent outcomes.

Take it with you

  • Predictive and parameter uncertainty refer to different random quantities.
  • Variance components depend on the declared model and posterior approximation.
  • Agreement and narrow uncertainty cannot rule out shared misspecification.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.