Predictive versus parameter uncertainty
Decompose uncertainty within a declared predictive model while separating posterior approximation, data variation, and model misspecification.
On this page
The question for today
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.
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.
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.
Work it through
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.
The weighted mean is (8 + 12)/2 = 10 milliseconds. This is a mean of predictions, not a confidence bound for the parameter.
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.
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 arithmetic identifies components within the specified mixture while leaving model adequacy as an empirical question.
The living notebook
Explore the model
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.
- Mean · np
- 5
- Variance · np(1−p)
- 2.5
- Total probability
- 1
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.
Pause & reason
Check your understanding
Your turn
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.
Keep the idea
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.