Tutorials › Introduction to Statistics › Why Prediction Intervals Are Wider Than Mean Response Intervals

Inference for the Regression Slope · Tutorial 913 of 996

Why Prediction Intervals Are Wider Than Mean Response Intervals

Learn how residual variation enters individual predictions and how it changes the standard error and width of a regression interval.

Intermediate 10 min read

What You'll Learn

  • Distinguish uncertainty about a mean response from variation among individual responses.
  • Derive the variance and standard error for predicting an individual response.
  • Compare the formulas for mean response confidence intervals and individual prediction intervals.
  • Calculate and interpret both intervals for a hypothetical linear regression.
  • Explain why prediction intervals are wider and identify a common misinterpretation.

Two Different Questions at the Same Predictor Value

A fitted regression line describes the mean response associated with each predictor value. But a particular individual response usually differs from that mean. This distinction leads to two related intervals at a predictor value \(x_0\): a confidence interval for the mean response at \(x_0\), and a prediction interval for one individual response at \(x_0\).

Both intervals account for uncertainty in estimating the regression line. A prediction interval must also account for the variation of individual responses around that line. The extra source of variation is why an individual prediction interval is wider, even though the intervals share the same center.

Definition: The mean response at \(x_0\) is the population average \(E(Y\mid x_0)\). A mean response confidence interval estimates this average. An individual prediction interval estimates a range for one future response \(Y_{\mathrm{new}}\) at \(x_0\), not for the average of responses at \(x_0\).

The Regression Model and the Two Standard Errors

Use the simple linear regression model \(Y_i=\beta_0+\beta_1x_i+\varepsilon_i\). For the usual \(t\)-based intervals, assume the mean response is linear in \(x\), the errors are independent with mean zero and constant variance \(\sigma^2\), and the errors are normally distributed. The predictor values are treated as fixed, or inference is understood conditionally on the observed predictor values. The usual residual standard deviation \(s\) estimates \(\sigma\), and the residual degrees of freedom are \(n-2\).

Write \(\bar{x}\) for the sample mean of the predictor and \(S_{xx}=\sum_{i=1}^{n}(x_i-\bar{x})^2\). At a specified \(x_0\), the fitted mean response is \(\hat{y}_0=b_0+b_1x_0\). Define

$$ h_0=\frac{1}{n}+\frac{(x_0-\bar{x})^2}{S_{xx}}. $$

This quantity describes the sampling uncertainty in the fitted mean at \(x_0\), relative to \(\sigma^2\). The standard error for estimating the mean response is \(s\sqrt{h_0}\). For predicting a new individual response, the standard error is \(s\sqrt{1+h_0}\). The added \(1\) is the residual-variation term.

Using \(t^*=t_{1-\alpha/2,n-2}\), the two-sided \(100(1-\alpha)\%\) intervals are:

$$ \text{Mean response CI: }\quad \hat{y}_0\pm t^*s\sqrt{h_0} $$
$$ \text{Individual prediction interval: }\quad \hat{y}_0\pm t^*s\sqrt{1+h_0}. $$

The formulas use the same fitted center and the same \(t\) critical value when they are calculated at the same \(x_0\), confidence level, and fitted model. Their standard errors—and therefore their margins of error—differ.

Worked Example: Comparing the Standard Errors

Suppose a hypothetical regression predicts completion time in minutes from a numerical score. For a sample of \(n=25\), let \(\bar{x}=50\), \(S_{xx}=2500\), and consider \(x_0=68.3193\). The fitted value is \(\hat{y}_0=32\) minutes, and the residual standard deviation is \(s=5\) minutes. Compare the standard errors at this predictor value.

First calculate the relative uncertainty in the estimated mean response:

$$ h_0=\frac{1}{25}+\frac{(68.3193-50)^2}{2500} \approx 0.17424. $$

Thus the standard error for the mean response is:

$$ SE_{\mathrm{mean}}=5\sqrt{0.17424} \approx 5(0.41742)=2.087\text{ minutes}. $$

For an individual response, add the residual-variation term before taking the square root:

$$ SE_{\mathrm{individual}}=5\sqrt{1+0.17424} =5\sqrt{1.17424} \approx 5(1.08362)=5.418\text{ minutes}. $$

The individual prediction standard error is larger because a person’s response can vary around the population mean, even if that mean were known exactly. The ratio of the two standard errors is

$$ \sqrt{\frac{1.17424}{0.17424}}\approx 2.596. $$

For this particular predictor value and sample, the individual prediction standard error is about 2.60 times the mean-response standard error. This ratio depends on \(h_0\); it is not a universal multiplier.

Why the Extra Term Appears

The difference follows from the prediction error for one new response. Let \(\mu(x_0)=\beta_0+\beta_1x_0\) be the true mean response at \(x_0\), and write the new response as \(Y_{\mathrm{new}}=\mu(x_0)+\varepsilon_{\mathrm{new}}\). Subtract the fitted value:

$$ Y_{\mathrm{new}}-\hat{y}_0 = \bigl(\mu(x_0)-\hat{y}_0\bigr)+\varepsilon_{\mathrm{new}}. $$

The first term is the error from estimating the regression mean. Under the model, its variance is \(\sigma^2h_0\). The second term is the new individual’s deviation from that mean, with variance \(\sigma^2\). The new error is independent of the data used to fit the line, so it is independent of the first term. Therefore, the variance of the sum is the sum of the variances:

$$ \begin{aligned} \operatorname{Var}(Y_{\mathrm{new}}-\hat{y}_0\mid x_0) &=\operatorname{Var}(\mu(x_0)-\hat{y}_0) +\operatorname{Var}(\varepsilon_{\mathrm{new}})\\ &=\sigma^2h_0+\sigma^2\\ &=\sigma^2(1+h_0). \end{aligned} $$

Replacing \(\sigma\) with its estimate \(s\) gives the prediction standard error \(s\sqrt{1+h_0}\). By contrast, estimating the mean response involves only the first source of uncertainty, with variance \(\sigma^2h_0\), and standard error \(s\sqrt{h_0}\). The “1” is not an arbitrary adjustment: it represents the variance of a new individual error, measured in units of \(\sigma^2\).

Result: Since \(h_0>0\), \(1+h_0>h_0\), so \(s\sqrt{1+h_0}>s\sqrt{h_0}\) when \(s>0\). At the same \(x_0\) and confidence level, the two intervals have the same center and critical value, so the prediction interval is strictly wider.

Worked Example: Comparing the Intervals in Context

Suppose a hypothetical regression predicts a measurement in units from a laboratory score. There are \(n=10\) observations, the fitted value at \(x_0\) is 32 units, \(s=4\) units, and \(h_0=0.133373\). Construct approximate 95% intervals for the mean response and an individual response at \(x_0\).

The residual degrees of freedom are \(n-2=8\), so the 95% critical value is approximately \(t^*=2.306\). The standard error for the mean response is:

$$ SE_{\mathrm{mean}}=4\sqrt{0.133373} \approx 1.4608\text{ units}. $$

Its margin of error is approximately \(2.306(1.4608)=3.369\) units. The mean response confidence interval is therefore:

$$ 32\pm 3.369=(28.63,\,35.37)\text{ units}. $$

For an individual response, the standard error includes residual variation:

$$ SE_{\mathrm{individual}} =4\sqrt{1+0.133373} \approx 4.2584\text{ units}. $$

Using the unrounded critical value and standard error, the margin of error is about \(9.820\) units. The 95% prediction interval is approximately:

$$ 32\pm 9.820=(22.18,\,41.82)\text{ units}. $$

If the rounded values \(2.306\) and \(4.258\) are multiplied directly, the product is \(9.818948\), which rounds to \(9.819\). Using the unrounded inputs gives a margin of about \(9.820\); either calculation gives the same interval endpoints to two decimal places. The mean-response interval concerns the average at this score, while the wider prediction interval allows for the additional variation of one individual measurement.

How Much Wider?

Both interval widths are twice their margin of error. Because their centers and critical values match, the ratio of their widths is the same as the ratio of their standard errors:

$$ \frac{\text{prediction-interval width}} {\text{mean-response CI width}} = \sqrt{\frac{1+h_0}{h_0}}. $$

This expression also clarifies the meaning of “wider.” The individual interval includes the mean-estimation uncertainty already present in the confidence interval, plus the variation among individual outcomes. The additional variation can be substantial when the fitted mean is estimated precisely. The ratio is not fixed, because \(h_0\) varies with the sample size and predictor design.

The addition of \(\sigma^2\) can also be seen directly in the interval variances: mean-response estimation has variance \(\sigma^2h_0\), whereas prediction error for an individual has variance \(\sigma^2h_0+\sigma^2\). It would be incorrect to add \(s\) directly to the mean standard error; independent variances add, and standard errors are square roots of variances.

Worked Example: Two 95% Intervals at One Predictor Value

In a hypothetical study, a regression predicts a concentration in milligrams per liter. Suppose \(n=16\), \(h_0=0.25\), \(s=3\) milligrams per liter, and \(\hat{y}_0=70\) milligrams per liter. Find both 95% intervals and explain their different purposes.

Here \(df=16-2=14\), with \(t^*\approx2.145\). The standard error for the mean response is:

$$ SE_{\mathrm{mean}}=3\sqrt{0.25}=1.5\text{ milligrams per liter}. $$

The margin is \(2.145(1.5)=3.2175\) milligrams per liter, giving the mean-response confidence interval:

$$ 70\pm3.2175=(66.78,\,73.22)\text{ milligrams per liter}. $$

For one individual response:

$$ SE_{\mathrm{individual}}=3\sqrt{1.25} \approx3.3541\text{ milligrams per liter}. $$

The prediction margin is approximately \(2.145(3.3541)=7.195\) milligrams per liter, so the 95% prediction interval is:

$$ 70\pm7.195=(62.81,\,77.19)\text{ milligrams per liter}. $$

The first interval estimates the population mean concentration for this predictor value. The second gives a range for one individual concentration under the model. Neither interval says that an individual measurement must lie inside the interval; it is a procedure that, under the assumptions, covers the target at the stated long-run rate.

Common Mistake: Treating an Individual Like a Mean

A common mistake is to use the mean-response confidence interval when asked to predict an individual outcome, on the grounds that both questions use the same fitted line. The fitted center is indeed the same, but the uncertainty is not. The confidence interval accounts for uncertainty in estimating the mean; the individual prediction interval accounts for that uncertainty and the individual’s residual variation.

A second mistake is to interpret the prediction interval as a range for the mean response. That interpretation reverses the target: the mean is a population average, while an individual response is one outcome around that average. Always identify whether the question asks about an average or one future observation before selecting the interval.

Key takeaway: At a fixed predictor value, use \(s\sqrt{h_0}\) for uncertainty in the estimated mean response and \(s\sqrt{1+h_0}\) for prediction of one individual response. The added \(1\) represents residual variation, making the individual prediction interval wider.

Check Your Understanding

Use the formulas and reasoning in this tutorial to answer the questions.

  1. In words, distinguish a confidence interval for the mean response at \(x_0\) from a prediction interval for one individual response at \(x_0\).
  2. Suppose \(h_0=0.20\). Write the standard-error multipliers for a mean response and an individual response, both in terms of \(s\).
  3. Why can the variance of the individual prediction error be written as \(\sigma^2h_0+\sigma^2\) under the regression model?
  4. If the two intervals are calculated at the same \(x_0\) and confidence level, what do they share, and why is one wider?
  5. A colleague says, “The prediction interval is wider because its center is farther from the fitted line.” Identify the error in this reasoning and give the correct explanation.