Lesson
The theory — The Winner's Curse
The winner's curse in statistics is the bias of estimates kept because they are statistically significant. A study's estimate d̂ of a true effect δ is significant at α = 0.05 when |d̂| > z·SE with z = 1.96, and given that, its expected size E(|d̂| | significant) exceeds δ. Andrew Gelman and John Carlin measure it with the exaggeration ratio E(|d̂| | significant)/δ (the Type M error) and with the probability that a significant estimate has the wrong sign (the Type S error). For a normal estimate with λ = δ/SE, the power is Φ(λ − z) + Φ(−λ − z) and the Type S rate is Φ(−λ − z) divided by the power.
What each symbol means
δ- the true effect, in standard deviations of the outcome
d̂- one study's estimate of δ
n- the number of people in each of the two groups
SE- the standard error of the estimate, √(2/n)
λ- the true effect measured in standard errors, δ/SE
z- the bar in standard errors, 1.96 for p < 0.05 two-sided
- Assumes
- Two groups of n with an outcome whose values are normal with a known standard deviation of 1, so the estimate is exactly normal, N(δ, 2/n), and the test is a two-sided z-test at α = 0.05. The true effect is positive and the same in every study, the studies are independent, and a study is kept exactly when it is significant.
Problem solved in full
-
The study that found 0.45 4 steps
A study with 50 people in each group reports an effect of 0.45 standard deviations, p = 0.024. Similar interventions are known to move this outcome by about 0.10. If 0.10 is the truth, how much did this design have to overstate it, how likely was it to get the sign wrong, and how many people per group would give 80% power? Press The 0.45 study: the page sets δ = 0.10 and n = 50.
-
First the reported result. The standard error and the z statistic:
-
Take δ = 0.10 from outside the study. In standard errors that is λ = 0.5, and the power is:
-
Of all studies, 7.91% reach significance and 0.69% are significant on the wrong side, which is 8.8% of the significant ones. The significant ones average 0.479 in size:
-
For 80% power, λ must reach z + Φ⁻¹(0.80), the far tail being negligible, and λ = δ√(n/2):
Answer
1,570 people per group, 31.4 times as many. With 50 per group the design had 7.9% power, a significant estimate had to be at least 0.392, and the significant ones average 0.479, 4.79 times the truth, with 8.8% of them the wrong sign. The reported 0.45 is what this design produces when the true effect is 0.10, so it gives no reason to reject 0.10. Type 1570 into the page: power 80.0%, and the significant estimates average 0.112, 1.12 times the truth.
-
Learning path
What p < 0.05 does not promise
References (1)
- Type S and Type M errors and the exaggeration ratio, the formulas for power and the Type S rate under a normal estimate (note 2), the R function retrodesign() whose figures check:compute reproduces, the beauty-and-sex-ratio example the third insight quotes, and the circle the second insight describes (a significant estimate from a small study fed back into a power analysis makes the study look reasonable): A. Gelman and J. Carlin, "Beyond Power Calculations: Assessing Type S (Sign) and Type M (Magnitude) Errors," Perspectives on Psychological Science 9(6), 641–651, 2014.