Tutorials › Regression Analysis (Undergraduate) › Why Selected Coefficients Are Biased

Subset and Stepwise Selection · Tutorial 657 of 991

Why Selected Coefficients Are Biased

Learn how selecting the strongest association exaggerates its apparent size, and why an independent sample can reveal the difference.

Advanced 9 min read

What You'll Learn

  • Define winner’s-curse bias in a predictor search.
  • Explain why the largest selected association can look impressive when every true slope is zero.
  • Distinguish exaggeration in absolute magnitude from bias in the signed coefficient.
  • Calculate the chance of finding an extreme estimate across many null predictors.
  • Compare a selected coefficient with an estimate from independent data.
  • Identify when an independent comparison does and does not remove selection bias.

Why the Strongest Result Can Be Misleading

A search across many candidate predictors is designed to find a strong association. But even if none of the predictors has a real relationship with the outcome, random sampling variation will make some estimated slopes larger than others. Choosing the strongest estimate tends to select one that is unusually far from its true value. This selection effect is called the winner’s curse: the reported association is likely to overstate the size of the underlying association.

The issue is not that least squares produces a biased slope for a predictor chosen in advance. Rather, the search and the estimate use the same data. The selected estimate has been chosen partly because its random error happened to be favorable. The effect can be particularly striking when many predictors are searched and the largest coefficient is reported without the context of that search.

Definition: Winner’s-curse bias is the tendency for an estimate selected because it is unusually large (or unusually far from zero) to exaggerate the magnitude of the corresponding underlying effect. In predictor selection, it arises when the same data are used both to choose a predictor and to estimate its coefficient. With symmetric null predictors, the selected coefficient’s absolute magnitude is biased upward, even though its signed average over repeated searches can still be zero.

A Simple Model for Seeing the Effect

Suppose an outcome \(Y\) is modeled using a set of candidate predictors \(X_1,\ldots,X_p\):

$$ Y_i=\beta_0+\beta_1X_{i1}+\cdots+\beta_pX_{ip}+\varepsilon_i. $$

Here \(\beta_j\) is the population slope for predictor \(X_j\), holding the other included predictors fixed. For a clean illustration, assume that the predictors are fixed and mutually orthogonal in the sample, the errors are independent with mean zero and common variance \(\sigma^2\), and the errors are normally distributed. Under these assumptions, the least-squares estimate for predictor \(j\) has a normal sampling distribution centered at its true slope. In particular, if all the candidate predictors are null, then \(\beta_j=0\) for every \(j\), and each estimate fluctuates around zero.

Orthogonality is a simplifying condition, not a requirement for selection bias. With predictor columns that are mutually orthogonal after removing the intercept (and any other regressors), and normal errors, the slope estimates are independent, which makes the selection effect especially easy to calculate. In real searches, predictors can be correlated and selection rules can be more complicated, but choosing an estimate partly because it is extreme still creates a risk of exaggeration.

A useful simplified description is \(\hat{\beta}_j\sim N(0,\tau^2)\) for every null predictor, where \(\tau\) is the common standard error. The search may select the largest positive estimate or the estimate with the largest absolute value. Those rules have an important difference: selecting the largest positive estimate produces a positive selected slope on average, whereas selecting the largest absolute estimate and keeping its original sign produces a signed average of zero under complete symmetry. In the latter case, the selected estimate is still unusually large in absolute value.

Worked Example: A Simulated Search Among Null Predictors

Consider a hypothetical simulation with 12 candidate predictors. The data are generated so that every true slope is zero. In one screening sample, a search produces the following estimated slopes, in outcome units per predictor unit:

PredictorEstimated slope
\(X_1\)0.06
\(X_2\)-0.14
\(X_3\)0.09
\(X_4\)0.31
\(X_5\)-0.18
\(X_6\)0.11
\(X_7\)-0.07
\(X_8\)0.22
\(X_9\)-0.12
\(X_{10}\)0.15
\(X_{11}\)-0.27
\(X_{12}\)0.04

The search selects \(X_4\), because its estimated slope has the greatest absolute value: \(|0.31|>|-0.27|\), and all the other absolute estimates are smaller. The reported screening-sample coefficient is therefore \(0.31\), even though the true slope for \(X_4\), like every other slope, is zero. This single search does not prove that the estimate is biased in every possible repeated run; it illustrates how an extreme estimate can arise from the selection rule itself.

Now suppose an independent validation sample is collected from the same population, and the already-selected predictor \(X_4\) is fitted there without searching again. In this hypothetical run, its estimated slope is \(0.04\). The comparison is \(0.31\) in the screening sample versus \(0.04\) in validation. Neither estimate has to equal the true value in a particular sample, but the independent estimate was not selected for being extreme. Across repeated validation samples, its expected value is zero.

Why the Selected Estimate Is Exaggerated

The basic mechanism can be derived with just two candidate slopes. Suppose \(U\) and \(V\) are independent estimates of two null slopes, each with distribution \(N(0,\tau^2)\), and the search reports the larger estimate, \(M=\max(U,V)\), which is negative when both estimates are negative. Since \(M=(U+V+|U-V|)/2\), and both \(U\) and \(V\) have mean zero, its expectation depends on the absolute difference:

$$ \begin{aligned} E(M) &=\frac{1}{2}E\left(U+V+|U-V|\right)\\ &=\frac{1}{2}E|U-V|. \end{aligned} $$

Because \(U\) and \(V\) are independent normal variables, \(U-V\sim N(0,2\tau^2)\). A mean-zero normal variable with standard deviation \(s\) has expected absolute value \(s\sqrt{2/\pi}\). Here \(s=\sqrt{2}\tau\), so \(E|U-V|=2\tau/\sqrt{\pi}\). Consequently:

$$ E(M)=\frac{\tau}{\sqrt{\pi}}>0. $$

The true slopes are both zero, but the larger of the two estimates has a positive expectation. This is selection bias in a direct form: the estimates themselves are centered correctly before selection, but taking the larger one changes the average of what is reported.

Worked Example: Expected Size of the Larger of Two Estimates

Suppose the two independent null slope estimates in the derivation have standard error \(\tau=0.15\) outcome units per predictor unit. The search selects the larger estimate, \(M=\max(U,V)\), which may be negative. Its unconditional expected value is:

$$ E(M)=\frac{\tau}{\sqrt{\pi}} =\frac{0.15}{\sqrt{\pi}} \approx 0.0846. $$

Thus, although both true slopes are zero, the average of the selected larger estimate over repeated searches is approximately \(0.0846\) outcome units per predictor unit. The estimate for a predictor fixed in advance would instead have expectation zero. The calculation concerns the larger estimate; it is not the expected signed coefficient after selecting the largest absolute value and retaining its sign.

For a rule that selects the largest absolute estimate, write \(J\) for the selected predictor index, so that \(|\hat{\beta}_J|=\max_j|\hat{\beta}_j|\). The search has deliberately chosen an estimate with an unusually large magnitude. Under the all-null, symmetric model, a positive extreme and a negative extreme are equally likely; therefore \(E(\hat{\beta}_J)=0\) across repeated searches. But \(E|\hat{\beta}_J|\) is greater than the expected absolute value of one estimate. Once a particular search has selected, for example, a positive coefficient, its magnitude is the part most likely to be exaggerated.

Worked Example: How Often a Null Search Finds an Extreme Estimate

Suppose 20 independent null slope estimates each have standard error \(0.20\) outcome units per predictor unit. For any one estimate, a magnitude of at least \(0.52\) corresponds to a standardized value of:

$$ z=\frac{0.52}{0.20}=2.60. $$

For a standard normal variable, the two-sided probability of a value at least this far from zero is approximately \(0.00932\). Under independence, the probability that none of the 20 estimates reaches that magnitude is \((1-0.00932)^{20}\). Therefore, the chance that at least one does is:

$$ P\left(\max_{1\leq j\leq 20}|\hat{\beta}_j|\geq 0.52\right) =1-(1-0.00932)^{20} \approx 0.1708. $$

That is about \(17.08\%\). The result does not mean that a particular coefficient has a \(17.08\%\) chance of being nonzero; every true slope in this example is zero. It means that an estimate with magnitude at least \(0.52\) can appear somewhere in this 20-predictor null search about 17.08% of the time under the stated assumptions. The calculation helps explain why a coefficient that looks large when viewed alone may be less surprising after accounting for the number of predictors searched.

Using Independent Data to Check the Magnitude

An independent sample can separate the act of choosing a predictor from the act of estimating its coefficient. Use the screening sample to select an index \(J\), then estimate the slope for that fixed selected predictor in a new sample whose outcomes were not used in the search. If the new-sample estimate is unbiased for the same population slope, then, conditional on the selected index, its expected value is the true slope:

$$ E(\hat{\beta}_{J,\mathrm{new}}\mid J)=\beta_J. $$

The reason is that the independent estimate has not been chosen for being large in the new sample. If all candidate slopes are zero, this conditional expectation is zero. The validation estimate may still be positive, negative, or large by chance in a particular sample, but it does not inherit the screening sample’s winner’s-curse mechanism, provided the validation outcomes were kept separate from selection.

Independence matters. If the same outcomes are inspected repeatedly, or if the validation sample is used to change which predictor is reported, the new estimate is no longer a clean independent check. Nor does independence guarantee that the selected association will replicate: the validation estimate is subject to sampling variation. It changes the expectation and the interpretation of the estimate, not the fact that finite samples can be noisy.

1
Choose using screening data.
Apply the stated search rule to the screening sample and record the selected predictor and its screening estimate.
2
Hold the choice fixed.
Do not use validation outcomes to replace the selected predictor or to search for a more impressive estimate.
3
Estimate in independent data.
Fit the selected predictor in the validation sample using the same outcome and predictor definitions.
4
Compare in context.
Report both estimates with their units and explain that the validation estimate is independent of the selection step, while still subject to sampling error.

Common Mistakes and Why the Distinction Matters

  • Calling every selected signed coefficient upward-biased. If the rule chooses the largest absolute estimate but retains its sign, positive and negative selected coefficients balance under a completely symmetric all-null model. The average signed coefficient is zero, while the absolute magnitude is exaggerated. State which selection rule and which quantity—signed estimate or magnitude—you mean.
  • Treating the selected coefficient as though the predictor had been specified in advance. Its size partly reflects the search. Interpret it in light of how many candidates and which search rules were considered.
  • Assuming a small screening estimate proves the effect is absent, or a smaller validation estimate proves the first result was false. Both samples are variable. Independent data provide a less selection-affected estimate, not certainty about the population slope.
  • Using validation data to search again. Choosing whichever validation result looks strongest recreates the same selection problem in the validation sample.

The practical consequence is that a selected coefficient is not a neutral summary of evidence when the same observations drove the search. The selection criterion may be useful for building a model, but the selected coefficient can overstate the association’s magnitude. Independent estimation helps make that limitation visible. The size of the winner’s curse depends on the search, the variability of estimates, and their dependence; it is not a universal correction that can be subtracted from every reported slope.

Key takeaway: Searching among candidate predictors tends to select estimates with unusually large magnitudes. For all-null predictors, that creates upward bias in absolute size or in a direction chosen after seeing the data, even though the signed average under a symmetric largest-absolute-value rule can be zero. An independent sample estimates the already-selected coefficient without selecting it again.

Check Your Understanding

Use the assumptions and distinctions in the tutorial to answer these questions.

  1. In a search among null predictors, why can the selected absolute slope be large even though every true slope is zero?
  2. Under a symmetric all-null model, what is the difference between the expected signed coefficient selected by largest absolute value and the expected absolute value of that coefficient?
  3. Two independent null estimates each have standard error \(\tau\), and the search reports the larger estimate. What is its expected value?
  4. For 20 independent null estimates with standard error \(0.20\), what is the probability that at least one has absolute value at least \(0.52\), and what does that probability describe?
  5. Why should a predictor selected in screening data be held fixed when estimating its slope in an independent validation sample?