Updated
In this lesson
Wide tables, equations, and code scroll sideways. Swipe, or Tab to focus them and use the left and right arrow keys.
Regression to the mean
You will learn: Explain why selecting extreme noisy measurements changes the next average.
Start with: Joint and marginal distributions
Select people with unusually high scores on a noisy test, then measure the same fixed underlying quantity with independent fresh measurement noise. Their average repeat score may be less extreme because the first selection also favored positive measurement noise. Compare standardized measurements on the same scale; a change in raw height across generations is not the model here.
Filled points meet the first-score selection rule; open points do not. The large outlined point marks the selected-group average. The dashed line is the cutoff; the solid line is the model conditional mean.
Selected 49/200 observations. The population probability of this selection is 0.250000. Selection uses X only; every selected unit receives a fresh measurement error for Y. The plot clips 0 pairs outside ±4, but their values remain in group calculations.
| Quantity | Population conditional mean | Observed selected mean |
|---|---|---|
| First measurement X | 1.27111 | 1.21072 |
| Repeat measurement Y | 0.88977 | 0.76409 |
| Repeat minus first | -0.38133 | -0.44663 |
Keep the underlying quantity fixed; change the measurement noise
Each unit has a fixed latent value A with variance r=0.7. Its two independent errors have variance 1−r=0.3000. The simulator sets X=A+e₁ and Y=A+e₂. Each measurement has variance one, and their correlation is r. The errors are independent of A and one another; X and Y share A and need not be independent.
In the selected sample, average A=0.75057, first error=0.46014, and repeat error=0.01352. Selection favors the first error in the chosen direction. It does not select the fresh repeat error, although its finite-sample average can differ from zero.
Inspect the selected measurements and their errors
| A | e₁ | e₂ | X | Y |
|---|---|---|---|---|
| 0.8082 | 0.6151 | -0.0270 | 1.4233 | 0.7812 |
| 0.5419 | 0.9253 | -0.3642 | 1.4673 | 0.1777 |
| 0.4605 | 1.1048 | -0.0802 | 1.5653 | 0.3804 |
| 0.3328 | 0.4529 | 0.6652 | 0.7857 | 0.9980 |
| 1.5087 | -0.3825 | 0.3765 | 1.1262 | 1.8852 |
| 0.1758 | 0.7492 | 1.4550 | 0.9250 | 1.6308 |
| 0.1471 | 0.7831 | -0.8973 | 0.9302 | -0.7502 |
| 1.2374 | 0.0770 | 0.7231 | 1.3144 | 1.9605 |
| 0.2959 | 0.5330 | 0.4608 | 0.8289 | 0.7566 |
| 0.5359 | 0.6839 | -0.2166 | 1.2198 | 0.3193 |
| 0.1169 | 1.1184 | 0.0720 | 1.2354 | 0.1889 |
| 0.9745 | 0.2529 | 0.5815 | 1.2274 | 1.5561 |
| 0.9643 | 0.7231 | 0.5222 | 1.6874 | 1.4865 |
| 0.3079 | 0.4222 | -0.2367 | 0.7302 | 0.0713 |
| 0.2481 | 0.8670 | -0.3593 | 1.1152 | -0.1112 |
| 1.0710 | -0.3861 | -0.7447 | 0.6849 | 0.3262 |
| 0.5987 | 0.5813 | -0.5992 | 1.1799 | -0.0005 |
| 1.4775 | -0.7337 | -0.9508 | 0.7438 | 0.5267 |
| 0.4842 | 0.1977 | -0.1411 | 0.6819 | 0.3431 |
| 0.3786 | 1.1440 | 0.3479 | 1.5226 | 0.7265 |
| 0.5749 | 1.0158 | -1.4805 | 1.5907 | -0.9056 |
| 0.5619 | 0.1229 | -0.9507 | 0.6848 | -0.3888 |
| 0.9630 | 0.5821 | -0.0578 | 1.5451 | 0.9052 |
| 0.7248 | 0.1457 | 0.3452 | 0.8705 | 1.0700 |
| 0.4227 | 0.6069 | -1.1733 | 1.0296 | -0.7506 |
| 0.8629 | 0.1018 | -0.3436 | 0.9648 | 0.5194 |
| 1.8781 | 1.1646 | 0.1272 | 3.0427 | 2.0053 |
| 1.1717 | 0.0217 | -0.4472 | 1.1934 | 0.7245 |
| 1.4990 | 0.4159 | -0.1732 | 1.9149 | 1.3258 |
| 0.6799 | 0.8661 | -0.4054 | 1.5460 | 0.2745 |
| 0.3777 | 0.4725 | 0.4849 | 0.8502 | 0.8626 |
| 0.2658 | 0.5783 | 0.1952 | 0.8441 | 0.4609 |
| 0.8695 | 0.0151 | -0.2662 | 0.8846 | 0.6033 |
| 0.5735 | 0.6832 | -0.6950 | 1.2567 | -0.1215 |
| 1.2698 | 0.4367 | 0.7513 | 1.7066 | 2.0212 |
| 1.2690 | 0.0827 | 0.2977 | 1.3517 | 1.5667 |
| 1.1595 | 1.1220 | -0.3825 | 2.2816 | 0.7770 |
| 0.7810 | 0.9045 | 0.4473 | 1.6855 | 1.2283 |
| 0.3229 | 0.4123 | 0.8174 | 0.7352 | 1.1403 |
| 0.6884 | 0.2645 | 0.5736 | 0.9529 | 1.2620 |
| 0.4526 | 0.7351 | -0.2215 | 1.1877 | 0.2311 |
| 0.3177 | 0.4618 | 0.1002 | 0.7794 | 0.4179 |
| 1.3218 | 0.6113 | -0.1233 | 1.9331 | 1.1984 |
| 0.2591 | 0.6368 | 1.0819 | 0.8960 | 1.3411 |
| 0.7841 | 0.2135 | -0.2380 | 0.9975 | 0.5461 |
| 1.3268 | 0.0505 | 0.3094 | 1.3773 | 1.6363 |
| 0.9088 | -0.2116 | 0.4240 | 0.6972 | 1.3328 |
| 0.9321 | 0.1160 | 0.3460 | 1.0481 | 1.2781 |
| 0.8932 | 0.1894 | 0.7319 | 1.0826 | 1.6251 |
The red line is E[Y|X=x]=rx under this normal model. At r=1 there is no measurement noise and Y=X; at r=0 the two measurements are independent. A selected group's observed mean need not equal its prediction in one sample. No intervention is applied here, so a before/after change alone cannot establish an intervention effect.
What to notice
- Orange dots pass the selected cutoff on the first measurement. The larger dot is their observed first/repeat group mean.
- The population conditional means are compared with the actual selected sample. A small selected group can fluctuate above or below its prediction; an empty group has no fabricated sample mean.
- r=1 removes measurement noise and makes both readings identical. r=0 makes the two readings independent, so the selected repeat population mean is zero.
- The red line is the conditional expectation under the standardized normal model, not a rule requiring every individual point to lie on it.
Why it happens
The math is simple: in a bivariate normal with correlation r, the conditional expectation is . For |r| < 1, the conditional mean of Y is less extreme in magnitude than X on these standardized scales. The linear conditional-mean formula follows from the joint normal model. Imperfect correlation alone is not a universal guarantee of this exact line.
The practical danger
Regression to the mean produces many false causal stories:
- A student scores unusually low on a test, gets tutoring, scores higher — tutoring gets the credit.
- A sports team has an exceptional season, regresses the next year — the coach gets blamed.
- A company tries an unusual intervention when sales are worst, sales recover — intervention gets the credit.
Such changes can occur without an intervention, but this pattern alone cannot establish whether a particular intervention helped. A suitable comparison group is needed. Selection, repeated-measurement noise, and an intervention are different parts of a causal question; the scatterplot alone cannot separate every possible explanation in real data.
Construct repeated measurements of the same quantity
The simulator draws a latent value A with variance r, then independent errors e₁ and e₂ with variance 1−r. It sets X=A+e₁ and Y=A+e₂. The latent value does not change between measurements. Each reading has variance one, and their covariance is Var(A)=r, so their correlation is r as well.
The errors are independent conditional on the shared A; the readings themselves are correlated whenever r>0. The data table exposes the pieces for every selected unit. Selecting high X tends to select positive e₁ as well as high A. The fresh e₂ does not inherit that selection.
Work the conditional prediction
At r=0.6, a first standardized reading x=2 gives predicted repeat reading 1.2, hence a predicted change of −0.8. That does not mean the latent quantity deteriorated. It also does not mean every unit with x near 2 will score 1.2 next time: conditional variance is 1−r²=0.64 in this standardized joint normal model.
For selection above the population upper quartile, the cutoff is approximately 0.674490 and the conditional mean first reading is approximately 1.271106. At r=0.6 the selected repeat population mean is 0.762664. The model table uses this truncated-normal calculation for the selected cutoff, rather than multiplying noisy observed group averages and labeling them exact population values.
Make a prediction
If every selected low-scoring participant improves after an intervention, does this simulation show the intervention had no effect?
Explore the answer
No. It shows that selection and fresh measurement noise can create improvement even without an intervention. The actual intervention could still help or harm. Estimating its effect requires an appropriate comparison and assumptions about how participants received it.
References
Random Services: multivariate normal distributions gives the joint-normal conditional expectation and variance. The latent-plus-error construction above verifies its covariance directly. Compare joint and conditional distributions and Simpson’s paradox for other ways aggregation and selection affect interpretation.