Why a Point Prediction Is Not a Certainty
A predictive model may return a single number, such as an estimated claim cost of $900. That number is a point prediction: a concise summary of what the model predicts, not a promise about the outcome. Two distinct kinds of uncertainty matter when communicating it. The estimated average response may be uncertain because it was learned from a finite data set, while an individual outcome may also vary around that average.
There is a related but separate question: how uncertain is an estimate of average model performance, such as a test-set coverage rate? A performance estimate summarizes a population or a set of future cases; it does not tell us exactly what will happen for one particular case. Keeping these quantities separate prevents a common reporting error: presenting uncertainty about an average as if it described the range of individual outcomes.
A confidence interval for a conditional mean is not an interval containing most individual outcomes. Conversely, a prediction interval is not a direct measure of how precisely a model's average performance has been estimated. The intended quantity must be identified before selecting or interpreting an interval.
Mean Intervals and Individual Prediction Intervals
Suppose a model estimates a conditional mean response at predictor values \(x_0\), denoted by \(\hat{\mu}(x_0)\). Let \(SE_{\text{mean}}\) be the standard error of that estimated mean, and let \(\sigma\) be the standard deviation of individual errors around the conditional mean. Under an approximate normal model with independent mean-estimation error and a new individual error, the standard error for predicting one new outcome is the square root of the sum of their variances.
Using a standard normal approximation, an approximate 95% confidence interval for the mean is \(\hat{\mu}(x_0) \pm 1.96SE_{\text{mean}}\). An approximate 95% prediction interval for one new outcome is \(\hat{\mu}(x_0) \pm 1.96SE_{\text{new}}\). These are normal-approximation formulas, not universal formulas: they require the stated error structure and an appropriate estimate of each standard error. When parameters are estimated in a small sample, a suitable \(t\) critical value is generally used instead of 1.96.
Worked Example: Mean Estimate Versus One New Outcome
Worked Example: Mean Estimate Versus One New Outcome
A hypothetical model estimates a mean cost of $900 for cases with specified characteristics. The standard error of the estimated mean is $25, and the estimated standard deviation of individual outcomes around that mean is $300. Assume independent normal errors and use a 95% normal approximation. Find both intervals, in dollars, rounding to cents. The method is to calculate the mean interval and then combine the two independent variance contributions for an individual prediction.
Mean interval. Its standard error is $25, so the margin is \(1.96(25)=49.00\). The interval is \(900\pm49\), or [$851.00, $949.00].
Individual prediction interval. First combine the variances: \(25^2=625\), \(300^2=90{,}000\), and \(625+90{,}000=90{,}625\). Thus \(\sqrt{90{,}625}=301.0399\) to four decimal places. The margin is \(1.96(301.0399)=590.0382\), which rounds to $590.04. The interval is [$900-$590.04, $900+$590.04], or [$309.96, $1,490.04].
Check. Both intervals are centered at $900, and the individual interval is wider because it includes the $300 individual-outcome variation as well as uncertainty in the estimated mean. The common mistake is to use only \(1.96(25)\) for an individual outcome, which omits the much larger residual variation.
Why the Variances Add
The combined standard error follows from a simple variance argument. Write a new outcome as its conditional mean plus a new error. The difference between that outcome and the fitted mean is the estimation error in the fitted mean plus the new individual error. If these two errors are independent, the variance of their sum is the sum of their variances.
The second equality uses two conditions: the fitted-mean estimation error has variance \(SE_{\text{mean}}^2\), and the new outcome's error has variance \(\sigma^2\); their covariance is zero because they are assumed independent. Taking the square root gives \(SE_{\text{new}}\). If the errors are dependent, or if the individual error variance changes across predictor values, this particular calculation may not apply. The appropriate uncertainty method must reflect the model and data.
Worked Example: Using Regression Leverage
Worked Example: Using Regression Leverage
A hypothetical simple linear regression is fitted to \(n=25\) observations, estimating \(p=2\) coefficients (an intercept and a slope). At a specified predictor value, the fitted mean is 120 units, the residual standard deviation estimate is \(s=10\), and the leverage for that prediction point is \(h_0=0.04\). Assume independent normal errors with constant variance and use \(t_{0.975,23}=2.069\). Find the 95% intervals for the mean response and one new response. Use the regression standard-error formulas, then round interval endpoints to two decimals.
The residual degrees of freedom are \(n-p=25-2=23\). For a fitted mean, the standard error is \(s\sqrt{h_0}=10\sqrt{0.04}=10(0.2)=2\). Its margin is \(2.069(2)=4.138\), so the mean interval is \(120\pm4.138\), or [115.86, 124.14].
For one new response, the standard error includes the new-observation variance: \(s\sqrt{1+h_0}=10\sqrt{1.04}=10(1.019804)=10.19804\). The margin is \(2.069(10.19804)=21.10\) to two decimals. The prediction interval is \(120\pm21.10\), or [98.90, 141.10].
Check. The prediction interval is wider than the mean interval and both are centered at 120. The common mistake is to omit the \(1\) inside \(\sqrt{1+h_0}\), thereby calculating a mean interval while labeling it an individual prediction interval. These regression formulas rely on the stated linear-model assumptions; a different model may require a different interval method.
Uncertainty in Average Performance
An estimate of performance also varies from one evaluation sample to another. For example, if a procedure produces nominal 90% prediction intervals, the observed fraction of test outcomes inside those intervals is an estimate of its coverage. It is not guaranteed to equal 90% in every finite test set.
For a simple coverage calculation, assume the test cases are independent and each has the same probability \(p_c\) of being covered. The coverage indicator for each case is then a Bernoulli variable. If \(n\) cases are evaluated and \(\hat{p}_c\) are covered as a fraction, the estimated standard error is \(\sqrt{\hat{p}_c(1-\hat{p}_c)/n}\). A normal-approximation interval is \(\hat{p}_c\pm1.96SE\), when the approximation is reasonable. This interval concerns the underlying average coverage under these assumptions; it does not tell us whether a particular future outcome will fall inside its interval.
Worked Example: Uncertainty in Observed Test-Set Coverage
A hypothetical evaluation uses 100 independent cases to assess a procedure intended to produce 90% prediction intervals. Outcomes for 90 cases fall inside their corresponding intervals. Estimate the coverage rate and calculate a 95% normal-approximation interval. Then explain what the interval does and does not establish.
The observed coverage is \(\hat{p}_c=90/100=0.90\). Its estimated standard error is \(\sqrt{0.90(1-0.90)/100}=\sqrt{0.09/100}=\sqrt{0.0009}=0.03\). The margin is \(1.96(0.03)=0.0588\). The approximate interval is \(0.90\pm0.0588\), or [0.8412, 0.9588], to four decimal places.
Check. The interval describes sampling uncertainty in the average coverage rate under the independent, common-coverage Bernoulli assumption. It does not mean that any one future case has a 84.12% to 95.88% probability of being covered. The common mistake is to interpret uncertainty about a population rate as a probability statement about one specific case.
What Intervals Do—and Do Not—Include
An interval is only as credible as its assumptions and the data used to construct it. For example, the normal prediction formulas above do not automatically account for an incorrect model form, omitted variables, biased or poorly measured outcomes, or a change in the relationship between predictors and outcomes. Nor does a narrow interval prove that the model is appropriate for a particular decision. Earlier material on sensitivity to modeling and preprocessing choices addresses one source of variation in predictions; it should not be silently treated as if it were already included in a standard regression standard error.
A confidence interval is also not a guarantee that the fixed mean has a specified probability of lying in the particular interval after it has been calculated. In the frequentist interpretation, the confidence level describes the long-run coverage of the interval procedure under its assumptions. A prediction interval has a different target: a future individual outcome. Neither interval should be presented without naming that target and describing the important assumptions.
When reporting a point prediction, give context rather than false precision. For example, “The model estimates an average cost of about $900 for cases with these recorded characteristics; individual costs can vary substantially” is more informative than “The cost will be $900.” If an interval is available, state whether it describes uncertainty in the mean or a new individual outcome, its confidence level, and the assumptions that materially affect its interpretation. Avoid calling a range a guarantee.
Check Your Understanding
For each problem, identify the quantity being described before calculating. All money amounts are hypothetical.
- A model estimates a mean of $240, with standard error $12 and individual residual standard deviation $50. Under independent normal errors, calculate approximate 95% intervals for the mean and one new outcome using 1.96. (Answer: mean [$216.48, $263.52]; individual [$139.22, $340.78].)
- A regression uses \(n=16\) observations and estimates \(p=2\) coefficients. At a specified point, \(\hat{\mu}=75\), \(s=8\), and \(h_0=0.0625\). With \(t_{0.975,14}=2.145\), calculate the 95% mean and individual intervals. (Answer: mean [70.71, 79.29]; individual [57.31, 92.69].)
- Of 64 independent test cases, 56 outcomes fall inside their respective prediction intervals. Estimate the coverage rate and its approximate 95% Wald interval using the plug-in binomial standard error. (Answer: coverage 0.8750; interval [0.7940, 0.9560].)
- On 49 independent evaluation cases, the average signed residual (observed minus predicted) is \(-3\) units and the residual sample standard deviation is 14 units. Using a normal approximation, calculate a 95% interval for the mean signed residual and interpret whether it indicates clear average bias. (Answer: standard error 2; interval [\(-6.92\), \(0.92\)]; it includes zero, so it does not indicate clear average bias at this approximate level, but it does not establish individual accuracy.)
- A conditional mean estimate is 500, its standard error is 18, and the new-outcome residual standard deviation is 24. Under independent normal errors, calculate 90% intervals using \(z=1.645\), rounding to cents. (Answer: mean [470.39, 529.61]; individual [450.65, 549.35].)