The 180-Day Descent
Download this day:EPUBPDF

Block I · Foundations of Knowledge & Reasoning · Day 006 / 180

Statistics & the Art of Not Fooling Yourself

Given enough defensible choices, many datasets contain a path to p < .05. The whole discipline is the struggle not to walk it by accident.

THE GARDEN OF FORKING PATHSsame data, many defensible specifications, one tempting headlinesame datasetnoise onlyreasonable choicesoutlierskeep / trimcovariatenone / Wscaleraw / ranksampleall / subsetreported resultsp = .03p = .61p = .44p = .78reported branch
The garden of forking paths. The data are pure noise — yet enough branches guarantee one ends in a "discovery." The danger is choosing after the data points the way.

In 2011 three psychologists set out to prove something impossible. They sat twenty undergraduates down, played half of them the Beatles’ “When I’m Sixty-Four” and the other half a control tune, then asked a battery of questions — including each person’s date of birth. After a perfectly ordinary statistical analysis, they announced their finding: listening to “When I’m Sixty-Four” made people a year and a half younger. Not feel younger. Be younger — their birth dates said so, at .

This is, of course, impossible. A song cannot reach back and edit the year you were born. And that was exactly the point. Joseph Simmons, Leif Nelson, and Uri Simonsohn had used only standard, respectable tools, made only choices that any working scientist makes every week — and produced a statistically significant impossibility. Their paper was a deliberately absurd demonstration designed to expose a vulnerability. Today we learn what it exposed, and how a science learns to stop fooling itself.

Where we are

On Day 2 we met the replication crisis as a set of bruising numbers (97% of original psychology studies reported significant effects; only 36% replicated) and named p-hacking without yet opening it up. Today we open it up. This is the statistical engine room beneath that crisis — and beneath Day 1’s Gettier worry (a belief that’s true by luck, not connection) scaled up to entire literatures. We’ll lean hard on Day 4’s lesson that a -value is not the probability the null is true (that confusion is the base-rate fallacy in disguise), and on Day 5’s warning that which variables you control for is itself a fork in the road. Threads today: information (signal vs noise) and computation (the lab as a fallible inference machine), with a quiet preview of emergence — science as an error-correcting system bigger than any one analyst.

The hook

The easiest person to fool

“The first principle is that you must not fool yourself — and you are the easiest person to fool.” — Richard Feynman, Caltech commencement, 1974

Richard Feynman standing at a blackboard in 1959.

Richard Feynman at a blackboard in 1959. The quote above comes from his 1974 Caltech commencement address on “cargo cult science.” Image: The Big T 1959, public domain in the United States via Wikimedia Commons.

Feynman was talking about a kind of integrity that no formula can supply. The cruel twist of modern statistics is that the formulas can still be correct, but after flexible searching they no longer protect the nominal error rate. They lend the warm authority of mathematics to whatever you were already hoping to find.

The villain has a clinical name: researcher degrees of freedom. Every real analysis hides dozens of small, defensible decisions. Which outliers to exclude? Should you control for age? Sex? Both? Do you stop collecting data now, or run twenty more participants? You measured three outcomes — which do you report? Each choice, taken alone, is innocent. Taken together, made after you’ve glanced at how the numbers are leaning, they become a machine for minting false discoveries.

Simmons, Nelson & Simonsohn proved this with arithmetic before they proved it with the Beatles. By simulation, they showed that bundling just four ordinary degrees of freedom — peeking at the data and adding more if needed; choosing between two related outcome measures; splitting the sample by a variable like gender; dropping one of three experimental conditions — can inflate your chance of a false positive from the advertised 5% to as high as 61%. Flip enough innocent switches and a coin-flip’s worth of noise becomes a near-certainty of “success.”

Run the false-positive factory: every simulated study has no true effect, so any apparent discovery comes from the choices you permit. Watch analytic freedom turn noise into a headline.

Interactive · live Monte-Carlo simulation

The False-Positive Factory

Each run simulates ~1,500 experiments in which there is no real effect at all — two groups drawn from the same distribution. If the null is true, the model assumptions hold, and you run one pre-specified test, about 5% will still cross by chance. Now hand yourself some researcher degrees of freedom and re-run. The factory simply reports the prettiest result it can find among the analyses you've allowed.

5%
false-positive rate · "findings" from pure noise
0%↑ nominal 5%~65%
No tricks enabled. If the null is true and the assumptions hold, one pre-specified test gives the error rate the system promises.

The core model

What a -value actually says (and the six lies people tell about it)

To not be fooled, you have to know precisely what the instrument reports. Many readers, including trained scientists and a worrying number of the statistics instructors teaching them, get this wrong. So, carefully:

Definition · the p-value

The p-value is the probability — computed under a specified null model, including the null hypothesis and its sampling and modeling assumptions — that the chosen test statistic would be at least as extreme as its observed value.

Concrete version: suppose the null model says two groups have the same mean. If a difference this large or larger would occur 3% of the time under that model, then . That is not a 3% chance the null is true.

Course-styled diagram of a null-model bell curve with a shaded tail area labeled p-value, emphasizing that a tail probability under the null is not the probability that the null is true.

The shaded tail is the p-value: under the null model, the probability of test-statistic values this extreme or more extreme. Adapted from Repapetilto and Chen-Pan Liao, CC BY-SA 3.0 via Wikimedia Commons.

Read it twice. It is a tail-area statement about a chosen test statistic under the null model — schematically, . It is emphatically not a statement about the hypothesis given the data, . Confusing those two is the exact inversion Day 4 warned about with the disease-test problem: is not , and no amount of staring turns one into the other without a prior. Here are the six misreadings to burn out of your reflexes — drawn from Greenland and colleagues’ field guide to twenty-five of them:

Three probabilities people confuse

Alpha is the long-run Type I error rate of a testing procedure if the null model is true: with , about 5% of true-null tests will falsely cross the line. The false discovery rate asks a different question: over repeated use of a specified discovery procedure, what is the expected proportion of false claims among those it declares significant? That depends on base rates, power, and selection, not just . A posterior probability, such as , is different again: it needs a prior. A 5% is neither a 5% false discovery rate nor a 5% chance the null is true.

  • means there’s a 3% chance the null is true.”No. That’s — a Bayesian quantity needing a prior. The -value already assumes the null.
  • means a 5% chance the result is just chance.”No. Chance (the null) is the very assumption the calculation is built on; it can’t also be the thing in doubt.
  • means there is no effect.”No. Absence of evidence is not evidence of absence — an underpowered study fails to detect real effects all the time.
  • is the probability the alternative is true.”No. Neither nor its complement is a probability about any hypothesis.
  • p tells you the result will replicate.”No. A low from one noisy study says little about the next study.
  • ”A significant result is a large or important one.”No. With a big enough sample, a trivially small effect is “significant.” Significance ≠ size ≠ importance.

That last one is the quiet killer, so it gets its own model.

Significance is not size: meet the effect size

A -value entangles effect size, noise, sample size, and model assumptions. Pour in enough participants and even a microscopic, meaningless difference crosses the line. So grown-up statistics insists on reporting the effect size separately: Cohen’s d (a difference in means measured in standard deviations) or a correlation r. Jacob Cohen offered rough labels — d ≈ 0.2 small, 0.5 medium, 0.8 large — while cheerfully admitting they were arbitrary: “Although arbitrary, the proposed conventions will be found to be reasonable by reasonable people.” The point is not the labels. The point is that a number with no effect size attached has told you almost nothing worth knowing.

Interactive · significance versus size

The Effect-Size Dial

Set a true mean difference, a per-group sample size, and a noise level. A group means one side of the two-group comparison, such as treatment versus control; n = 50 per group means 50 observations in each side, 100 total. The plot uses fixed x and y scales: mean difference from -3 to +3, and bucket frequency from repeated estimates. Finer bars show the estimate distribution; the curve is the matching sampling distribution; the horizontal bar is the 95% confidence interval on the same x scale.

95% CI widthbucket frequency
p-value.14
95% CI[-0.10, 0.70]
Cohen's d0.30
Practical sizesmall

The interval you’ve been misreading too

Confidence intervals are sold as the humane alternative to -values, and they are better — but they hide their own trap.

Definition · the 95% confidence interval

A 95% CI is produced by a procedure that, across many hypothetical repeat experiments, brackets the true value 95% of the time. The “95%” describes the long-run reliability of the method, not the odds for your one particular interval.

So the natural sentence — “there’s a 95% chance the true value is in this interval” — is, in the frequentist world, simply false. Your interval either contains the truth or it doesn’t; the dice were thrown when you ran the study. (A Day 4 Bayesian credible interval does license that sentence — but only after you’ve committed to a prior. Different tool, different promise.) Nor does a CI mean 95% of your data live inside it, nor that values just outside are impossible. What it gives you, refreshingly, is a range of effect sizes compatible with your data — and that reframing, from “yes/no” to “how much, how precisely,” is the heart of what reformers call the new statistics.

Interactive · long-run coverage

The 100-interval test

Generate 100 repeated samples from a world where the true mean is known to be 0. The x-axis is fixed from -4 to +4 so width changes stay visually comparable. The confidence-level dial changes the long-run target: higher confidence covers more often but makes wider intervals. Sample size and noise mostly change precision. The highlighted interval is one observed study: it either covers or it does not.

covered truth95 / 100
missed truth5 / 100
highlighted intervalcovers
typical width0.74

Why small studies can exaggerate what they detect

Two ways to be wrong: a Type I error (false alarm — declaring an effect that isn’t there) and a Type II error (a miss — failing to see one that is). Statistical power is your chance of catching a real effect; more data buys more power. You’d assume an underpowered study just yields shrugs and “no result.” Worse: it actively poisons the literature. When power is low, the only estimates lucky enough to clear the significance bar are the substantially inflated ones — the winner’s curse. Pair that with publication practices that preferentially reward significant findings, and you manufacture a published record of effects that are both exaggerated and fragile. A 2013 audit (“Power failure,” Button et al.) put the median power across neuroscience at roughly 21%; a 2017 reanalysis by Nord and colleagues usefully qualifies the headline, arguing that power is not uniformly low across every neuroscience subfield. The lesson survives the caveat: where studies are underpowered, real effects are missed and the effects that surface are often inflated.

The deeper trap

You don’t have to cheat to be fooled

Here is the part that should worry careful researchers. You might read all of the above, swear off p-hacking, run exactly one pre-planned analysis, and report it faithfully — and still have an inflated false-positive rate. This is Andrew Gelman and Eric Loken’s garden of forking paths.

Their insight is subtle and a little haunting. The damage doesn’t require running many analyses. It only requires that the single analysis you ran was contingent on the data you happened to see. Suppose your data had come out slightly differently — a touch noisier here, a cluster there. You’d have made a different, equally reasonable choice: controlled for age instead of income, compared women instead of the whole sample, used a median instead of a mean. The other analyses are invisible because you never ran them. But the universe of paths you would have taken, in nearby worlds, is exactly the multiplicity that inflates your error rate. As Gelman and Loken put it, you can get a false positive with no “fishing expedition” and no p-hacking at all, even with the hypothesis fixed in advance.

Notice how this rhymes with Day 5. There, the lurking decision was which variables to condition on — and conditioning on a collider could manufacture a correlation out of nothing. The forking paths are the same hazard generalized: every defensible analytic choice is a junction, and the data quietly nudge you down the branch that flatters your hypothesis. A good-faith analyst and a p-hacker can produce the identical false result. The difference is only that one of them knows they did it.

The fix that names the disease

If the problem is that choices get made after seeing the data, the cleanest cure is to make them before: preregistration (write down your analysis plan, publicly, in advance) and Registered Reports (journals peer-review and accept the plan, before any results exist). We met these on Day 2; here’s why they bite. They freeze the confirmatory path; deviations become exploratory rather than hidden, converting a garden of “reasonable choices I made while peeking” into a single committed road.

The debate

Burn the threshold, or just lower it?

If is so abusable, what should replace it? Here the field splits — not into crank versus expert, but into camps of serious statisticians who genuinely disagree. It helps to lay them on a single line, from “patch it” to “torch it.”

Diagram · the reform spectrum

Four ways to live with p

Everyone on this line agrees the status quo is broken. They disagree about how radical the cure must be.

keep .05status quoredefine: .005Benjamin+ 2018(72 authors)justify your αLakens+ 2018(88 authors)retire itAmrhein+ 2019(800+ signatories)patch the rule……abolish the bright line
Redefine — Benjamin et al. (2018) keep the threshold idea but move it: call a "discovery," and .005–.05 merely "suggestive." Simple, blunt, and (critics say) treats a symptom.
Justify your alpha — Lakens et al. (2018) reject any universal number. Pick your threshold deliberately, case by case, weighing the real costs of false alarms versus misses, and show your reasoning.
Retire significance — Amrhein, Greenland & McShane (2019), with 800+ co-signers in Nature, want to scrap the "significant / not significant" dichotomy itself — the habit of treating and p = .06 as different worlds. Not a ban on -values; a ban on the bright line.
The field's statistical society — the American Statistical Association weighed in twice: a careful 2016 statement of six principles, then a more radical 2019 editorial urging us to stop saying "statistically significant" at all. (That 2019 piece was the editors' view, not formal ASA policy — a distinction that itself caused a row.)

What unites the whole spectrum, from cautious to incendiary, is a single migration: away from a yes/no verdict at a magic number, toward reporting how big an effect is, how uncertain you are, and how fragile the answer is to the choices you made. That last idea — fragility — is where the live 2020s frontier lives.

The frontier · 2024–2026

Making fragility visible — and the hype filter

If a single analysis can mislead, the modern move is brutally simple: run them all, and show the spread. Two named methods do this, and a third line of work tests the whole idea by turning loose armies of real scientists on the same data. As always, each claim gets a label.

Edge 01TransparencyFormal test

Multiverse & specification-curve analysis

Instead of defending one analysis, you enumerate every reasonable combination of choices — keep outliers or trim them, log-transform or not, control for this covariate or that — and compute the result under each. Multiverse analysis (Steegen, Tuerlinckx, Gelman & Vanpaemel, Perspectives on Psychological Science, 2016) displays the whole cloud of outcomes. Specification-curve analysis (Simonsohn, Simmons & Nelson, Nature Human Behaviour, 2020) sorts hundreds of specifications into a single curve and asks: across all the defensible ways to slice this, does the effect hold up, or does it evaporate the moment you wiggle a choice? It’s transparent, vivid, and increasingly expected by reviewers. Try the curve now, before the caveats continue.

One discipline matters more than the picture suggests: some “reasonable specifications” change the estimand. Controlling for W is not just another estimate of the same question unless W is a causally justified control; rank-transforming, trimming outliers, or dropping a subgroup can shift the target from a mean difference to an ordinal association, a robust effect, or a population-specific effect. Specification curves expose that fork, but they do not decide which question you meant to ask.

The hype filter matters here, though. These are superb transparency devices and weak inferential ones. The authors of specification-curve analysis concede the catch themselves: deciding which specifications are “reasonable” is a judgment call no algorithm can make, and “the goal to eliminate subjectivity is unattainable.” A determined arguer can still curate the multiverse. And there is no settled way to compute a single valid conclusion from a curve — which is why 2024 brought new machinery (PIMA, below) trying to supply one. Robustness display Objectivity machine

Interactive · run your own multiverse

The Specification Curve

A fixed synthetic dataset gives X and Y a weak, confounded association: the baseline model starts below the usual line, but some defensible choices can push it across. Here r is the Pearson correlation coefficient, from -1 to +1; positive r means larger X tends to come with larger Y, and values near 0 are weak. We draw the curve because one reported model can be arbitrary: the question is whether the conclusion survives reasonable analytic choices or depends on one path. The panel opens with every choice included and cherry-pick highlighted; turning off a chip narrows the multiverse. Teal dots are , grey dots are not. The check rows show which choices each specification uses, including whether it controls for the lurking Day 5 confounder W.

The full reading

Edge 02Analyst spread

Many analysts, one dataset — the clearest demonstration in the room

Here is the experiment that should change how you read every headline. Take one dataset, one clear question, and hand identical copies to dozens of expert teams. Do they converge on the same verdict? Often, no.

In Silberzahn et al. (Advances in Methods and Practices in Psychological Science, 2018), 29 teams (61 analysts) were asked a single question — do soccer referees give more red cards to dark-skin-toned players? Their estimated odds ratios ranged from 0.89 to 2.93; twenty teams found a significant positive effect, nine did not. Tellingly, the analysts’ prior beliefs and even their statistical expertise did not explain who found what. In Botvinik-Nezer et al. (Nature, 2020), 70 teams (180 researchers) analyzed the same brain-imaging dataset against nine pre-set hypotheses — and no two teams used the same pipeline; their yes/no conclusions diverged sharply even when their underlying statistical maps were highly correlated. Breznau et al. (PNAS, 2022) put the same data and hypothesis (does immigration erode support for social policy?) to 73 teams and watched estimates scatter from clearly negative to clearly positive. Even finance has its version: Menkveld et al.’s “Nonstandard Errors” (Journal of Finance, 2024) had 164 teams test the same hypotheses on the same market data, and found the cross-team variability rivaled ordinary statistical error — shrinking, encouragingly, when extra peer-review stages were added.

This is Day 1’s Gettier problem made flesh, inside a single dataset. Each team’s belief is “justified” by a competent analysis; whether it’s true — connected to reality rather than to their particular path through the garden — varies team to team. One analysis is one branch. Treat it accordingly.

Edge 03New methodsReform infra

Plugging the inference gap — and Day 4’s quiet alternative

The open problem in all of the above is valid inference across a multiverse: how do you draw one valid conclusion from a thousand entangled analyses? A 2024 entry, PIMA (Post-selection Inference in Multiverse Analysis; Girardi et al., Psychometrika, 2024), offers a sign-flipping test that aims to give the multiverse a proper joint error guarantee, reaching beyond the linear-model limits of the 2020 specification curve. It’s peer-reviewed and genuinely interesting, but new and not yet standard practice — a promising hint, not a settled tool.

Meanwhile the boring, durable fixes keep spreading: the Registered Reports format — where journals accept your plan before your results exist — is now offered at several hundred journals, and preregistration is becoming a default rather than a virtue. And recall Day 4’s e-values and “testing by betting” (Grünwald, Ramdas, Shafer): an alternative to the -value that stays valid even if you peek at your data and stop whenever you like — aimed squarely at the optional-stopping degree of freedom you toggled in the factory above. Still a research frontier rather than mainstream practice — promising, watch this space.

Open questions

What’s genuinely unsettled

  • Should “statistical significance” survive at all? Reformers are split between disciplining the threshold and abolishing it. Decisions still have to be made somewhere — a drug is approved or it isn’t — and abolitionists owe an account of how.
  • Can a multiverse ever yield a single valid verdict? Or is the choice of “reasonable specifications” an irreducibly subjective act that just relocates the forking paths one level up? PIMA and friends are early attempts; the jury is out.
  • Is variation in conclusions the same as variation in reality? A subtle 2023 reanalysis noted that in some many-analyst studies the headline-grabbing disagreement was about significance, while the underlying effect sizes were quietly consistent and small. Disagreement about a verdict can exceed disagreement about the number. Don’t over-learn “everything is hopeless.”
  • Do -values or e-values deserve the future? The betting-based tools elegantly solve optional stopping but can demand more data and import a new modeling burden (which bet?). Coexistence looks likelier than conquest.
  • And the question waiting in the AI block: when a model trained on millions of papers reports a “robust” result, has it learned to reason about evidence — or to imitate the very forking-path habits that got us here? (Days 138–145.)

The day in three sentences

Big idea
Statistical tools don’t protect you from fooling yourself — used with ordinary human flexibility they help you do it, because the dozens of defensible choices in any analysis (the “garden of forking paths”) let noise masquerade as discovery, so the cure is transparency, preregistration, effect sizes with uncertainty, and robustness across reasonable analyses rather than a verdict at a magic threshold.
Best analogy
The False-Positive Factory: hand yourself enough innocent freedoms and a coin-flip’s worth of nothing becomes a 61% chance of “success” — and the garden of forking paths, where the analyses you would have run in nearby worlds inflate your error even if you only ran one.
Live controversy
Whether to lower the threshold (), justify it case by case, or retire “statistical significance” altogether — with multiverse and many-analyst studies revealing how alarmingly a single dataset’s verdict depends on who analyzes it.

Threads today › information (separating signal from noise; the -value as a misread evidence measure) · computation (the lab as a fallible inference engine; the multiverse as brute enumeration) · emergence (science as an error-correcting system larger than any analyst — preregistration, many-analyst crowds, meta-analysis), leading directly into Information Theory.

Tomorrow → Day 7

Information Theory

Statistics separates signal from noise. Tomorrow we ask what signal is: Shannon’s bit, entropy as expected surprise, channel capacity, Maxwell’s demon — and Landauer’s discovery that forgetting has a physical price.

Sources

Sources & further reading

  1. Feynman, R. P. (1974). “Cargo Cult Science.” Engineering and Science 37(7): 10–13. Caltech’s 1974 commencement address; source of the opening quote. calteches.library.caltech.edu
  2. Open Science Collaboration. (2015). “Estimating the reproducibility of psychological science.” Science 349(6251): aac4716. doi:10.1126/science.aac4716; source for the 97% original significant / 36% significant replication statistic. doi.org/10.1126/science.aac4716
  3. Simmons, J. P., Nelson, L. D. & Simonsohn, U. (2011). “False-Positive Psychology: Undisclosed Flexibility in Data Collection and Analysis Allows Presenting Anything as Significant.” Psychological Science 22(11): 1359–1366. doi:10.1177/0956797611417632. doi.org/10.1177/0956797611417632 — the “When I’m Sixty-Four” demonstration; the 5%→61% simulation.
  4. Gelman, A. & Loken, E. (2014). “The Statistical Crisis in Science.” American Scientist 102(6): 460–465. — the “garden of forking paths”; earlier as a 2013 Columbia working paper. americanscientist.org
  5. Greenland, S., Senn, S. J., Rothman, K. J., Carlin, J. B., Poole, C., Goodman, S. N. & Altman, D. G. (2016). “Statistical tests, P values, confidence intervals, and power: a guide to misinterpretations.” European Journal of Epidemiology 31(4): 337–350. doi:10.1007/s10654-016-0149-3. doi.org/10.1007/s10654-016-0149-3 — the canonical list of 25 misreadings.
  6. Haller, H. & Krauss, S. (2002). “Misinterpretations of Significance: A Problem Students Share with Their Teachers?” Methods of Psychological Research Online 7(1): 1–20. — survey evidence that p-value misunderstandings persist even among instructors teaching statistics. epub.uni-regensburg.de/34338
  7. Badenes-Ribera, L., Frias-Navarro, D., Iotti, B., Bonilla-Campos, A. & Longobardi, C. (2016). “Misconceptions of the p-value among Chilean and Italian Academic Psychologists.” Frontiers in Psychology 7: 1247. doi:10.3389/fpsyg.2016.01247. doi.org/10.3389/fpsyg.2016.01247
  8. Wasserstein, R. L. & Lazar, N. A. (2016). “The ASA Statement on p-Values: Context, Process, and Purpose.” The American Statistician 70(2): 129–133. doi:10.1080/00031305.2016.1154108. doi.org/10.1080/00031305.2016.1154108 — the six principles.
  9. Wasserstein, R. L., Schirm, A. L. & Lazar, N. A. (2019). “Moving to a World Beyond ‘p < 0.05’.” The American Statistician 73(sup1): 1–19. doi:10.1080/00031305.2019.1583913. doi.org/10.1080/00031305.2019.1583913 — editorial to a 43-article special issue; “stop saying statistically significant” (editors’ view, not ASA policy).
  10. Benjamin, D. J., Berger, J. O., Johannesson, M., et al. (2018). “Redefine statistical significance.” Nature Human Behaviour 2(1): 6–10. doi:10.1038/s41562-017-0189-z. doi.org/10.1038/s41562-017-0189-z
  11. Lakens, D., Adolfi, F. G., Albers, C. J., et al. (2018). “Justify your alpha.” Nature Human Behaviour 2(3): 168–171. doi:10.1038/s41562-018-0311-x. doi.org/10.1038/s41562-018-0311-x
  12. Amrhein, V., Greenland, S. & McShane, B. (2019). “Scientists rise up against statistical significance.” Nature 567(7748): 305–307. doi:10.1038/d41586-019-00857-9; 800+ signatories. doi.org/10.1038/d41586-019-00857-9
  13. Button, K. S., Ioannidis, J. P. A., Mokrysz, C., et al. (2013). “Power failure: why small sample size undermines the reliability of neuroscience.” Nature Reviews Neuroscience 14(5): 365–376. doi:10.1038/nrn3475; median power ≈ 21%. doi.org/10.1038/nrn3475
  14. Nord, C. L., Valton, V., Wood, J. & Roiser, J. P. (2017). “Power-up: A Reanalysis of ‘Power Failure’ in Neuroscience Using Mixture Modeling.” The Journal of Neuroscience 37(34): 8051–8061. doi:10.1523/JNEUROSCI.3592-16.2017; caveat on the uniformity of low power across neuroscience. doi.org/10.1523/JNEUROSCI.3592-16.2017
  15. Center for Open Science. “Registered Reports.” Official COS overview and participating-journal list; source for “over 300 journals.” cos.io/initiatives/registered-reports
  16. Steegen, S., Tuerlinckx, F., Gelman, A. & Vanpaemel, W. (2016). “Increasing Transparency Through a Multiverse Analysis.” Perspectives on Psychological Science 11(5): 702–712. doi:10.1177/1745691616658637. doi.org/10.1177/1745691616658637
  17. Simonsohn, U., Simmons, J. P. & Nelson, L. D. (2020). “Specification curve analysis.” Nature Human Behaviour 4(11): 1208–1214. doi:10.1038/s41562-020-0912-z. doi.org/10.1038/s41562-020-0912-z
  18. Silberzahn, R., Uhlmann, E. L., Martin, D. P., et al. (2018). “Many Analysts, One Data Set.” Advances in Methods and Practices in Psychological Science 1(3): 337–356. doi:10.1177/2515245917747646; 29 teams, ORs 0.89–2.93. doi.org/10.1177/2515245917747646
  19. Botvinik-Nezer, R., Holzmeister, F., Camerer, C. F., et al. (2020). “Variability in the analysis of a single neuroimaging dataset by many teams.” Nature 582(7810): 84–88. doi:10.1038/s41586-020-2314-9; 70 teams. doi.org/10.1038/s41586-020-2314-9
  20. Breznau, N., Rinke, E. M., Wuttke, A., et al. (2022). “Observing many researchers using the same data and hypothesis reveals a hidden universe of uncertainty.” PNAS 119(44): e2203150119. doi:10.1073/pnas.2203150119; 73 teams. doi.org/10.1073/pnas.2203150119
  21. Mathur, M. B., Covington, C. & VanderWeele, T. J. (2023). “Variation across analysts in statistical significance, yet consistently small effect sizes.” PNAS 120(3): e2218957120. doi:10.1073/pnas.2218957120; caveat on reading many-analyst disagreement through effect sizes and intervals, not significance labels alone. doi.org/10.1073/pnas.2218957120
  22. Menkveld, A. J., Dreber, A., Holzmeister, F., et al. (2024). “Nonstandard Errors.” The Journal of Finance 79(3): 2339–2390. doi:10.1111/jofi.13337; 164 teams. doi.org/10.1111/jofi.13337
  23. Girardi, P., Vesely, A., Lakens, D., et al. (2024). “Post-selection Inference in Multiverse Analysis (PIMA): An Inferential Framework Based on the Sign Flipping Score Test.” Psychometrika 89(2): 542–568. doi:10.1007/s11336-024-09973-6. doi.org/10.1007/s11336-024-09973-6
  24. Ramdas, A., Grünwald, P., Vovk, V. & Shafer, G. (2023). “Game-theoretic statistics and safe anytime-valid inference.” Statistical Science 38(4): 576–601. doi:10.1214/23-STS894; e-processes, testing by betting, and optional-stopping validity. doi.org/10.1214/23-STS894
  25. Cohen, J. (1988). Statistical Power Analysis for the Behavioral Sciences (2nd ed.). Lawrence Erlbaum. — effect-size conventions, offered as deliberately arbitrary.
Deep dive appendixThe Deeper MachineryOptional extension.

The main lesson left a debt unpaid. We said, almost in passing, that a -value is one of the most misunderstood numbers in science — and then we hurried on to the famous misreadings, as if the confusion were everyone’s fault but the number’s. It isn’t. The deeper truth is that the tool itself is a chimera: stitched together from two rival philosophies whose founders despised each other, taught as if it were one coherent idea, and quietly incapable of meaning what almost everyone takes it to mean. This appendix is the autopsy. We’ll find out why repeated uncorrected testing makes false positives increasingly likely — approaching certainty in the idealized case of indefinite optional stopping; how wrong a “significant” result really is once you do the full arithmetic; how a small army of statistical detectives now hunts for the fingerprints; and what the people who actually know this stuff reach for instead. We close on a single, lovely idea that turns a -value into a number of coin flips — and hands us straight to tomorrow.

Continuing from the main lesson

You’ve met the False-Positive Factory, the garden of forking paths, the reform spectrum, and the many-analysts duels. Hold onto three threads in particular: Day 4’s Bayes (a -value is not P(hypothesis | data), and here we’ll finally compute what it costs to forget that), Day 2’s Popper and Ioannidis (severity and “most findings false”), and Day 1’s calibration instinct. Everything below is the substructure those ideas were standing on.

§1 The original sin

A tool built from two incompatible blueprints

Here is a fact that explains an astonishing amount of misery: the standard procedure every science student is taught as “hypothesis testing” is a hybrid. It welds Fisherian evidence to Neyman-Pearson long-run decision rules, and that hybrid often causes confusion.

On one side stood Ronald Fisher. In Statistical Methods for Research Workers (1925), he offered the -value as an informal, continuous gauge of evidence against a single null hypothesis — a measure of how strange your data would look if nothing were going on. He needed no alternative hypothesis, no rejection region, no long run. He suggested only as a homely rule of thumb, the rough point past which a coincidence starts to feel implausible, and explicitly said you could pick a different line if you preferred. For Fisher, a small p was a prompt to think, not a verdict to obey. The null, he wrote, is never proved — only, occasionally, embarrassed.

On the other side stood Jerzy Neyman and Egon Pearson. In a series of papers culminating in 1933, they rebuilt the whole enterprise as a decision rule, not a measure of evidence. You specify two hypotheses up front, fix in advance your tolerable rate of false alarms (α, the Type I error rate) and misses (β, the Type II rate; power = ), draw a line, and thereafter behave: reject or accept, accepting that over the long run you’ll be wrong at the rates you chose. There is, pointedly, no measure of evidence for the case in front of you — only a procedure with known long-run properties. Neyman thought Fisher’s inductive “evidence” talk was mush; Fisher thought Neyman had turned science into factory quality-control. The feud was lifelong and genuinely vicious.

Diagram · two philosophies, forced into one box

Fisher vs Neyman–Pearson

They answer different questions. The modern ” test” pretends to answer both at once — which is exactly why it answers neither cleanly.

Fisher’s significance test

1925 · evidential / inductive

  • One hypothesis (the null) — no alternative required.
  • p is a continuous measure of evidence against it.
  • .05 is a flexible convention, not a law.
  • A prompt to investigate, judged in context.
  • No “accept the null,” no fixed long-run error rate.

Neyman–Pearson test

1928–1933 · behavioural / decision

  • Two hypotheses fixed in advance ( and ).
  • , , and power set before seeing data.
  • .05 is a committed threshold, chosen for its costs.
  • A rule to act on — reject or accept, and move on.
  • No evidence claim about the single case; only long-run rates.
— welded together —

The null ritual (what’s actually taught)

Compute . Compare it to . Treat crossing it as both a Fisherian “the evidence is strong” and a Neyman–Pearson “the decision is made” — while reporting the exact as if it graded the evidence, yet thresholding it as if it were a fixed- decision. The psychologist Gerd Gigerenzer calls this incoherent mash-up the “null ritual,” performed by rote and understood by few. The confusion you met in the main lesson isn’t a failure of student attention. It’s baked into a tool that was never designed.

Once you see the seam, the whole replication crisis reads differently. We took a flexible Fisherian prompt-to-think, froze it into a Neyman–Pearson gate, stripped out the part of Neyman–Pearson that made the gate valid (pre-committed , adequate power, a real alternative), and then published only what got through. The wonder is not that so much fails to replicate. It’s that anything survives.

§2 Why repeated testing raises false alarms

The two facts that let repeated analysis manufacture false positives

The False-Positive Factory wasn’t a trick of clever programming. It rests on two mathematical facts so simple they’re almost insulting — and once you hold both, you understand why enough uncorrected looks can turn noise into apparently significant evidence.

Fact one: under the null, small p-values still arrive by chance

Run an exact test with continuous data where nothing is happening, under the assumptions of the null model, and the -value you get is a uniform random number between 0 and 1. For discrete, composite, or conservative tests, valid -values are often super-uniform: , not necessarily exactly equal. The practical warning survives that qualification. In the calibrated continuous case, exactly 5% of the interval [0, 1] lies below .05; in conservative cases the rate is no larger than that. But flip it around and the menace appears: even by-the-book, effect-free experiments sometimes hand you a “significant” result by raw luck. Run the experiment twenty ways — twenty outcomes, twenty covariates, twenty subgroups — and you’ve effectively bought many tickets on a small-probability prize. The expected number of “discoveries” in a drawer full of nothing is driven by how many times you looked.

Use the live panel to watch an idealized continuous exact test. With the true effect at zero, -values lie flat, with 5% falling below the line. Turn up a real effect and the distribution leans toward smaller values: genuine signal makes very small significant -values more common than barely significant ones. That surplus is a literature’s pulse; its absence is a warning sign.

Interactive · live simulation

The shape of a p-value

Each run simulates 4,000 experiments (two groups, n = 30 each) and bins their -values. Slide the true effect from 0 (nothing is real) upward. In this continuous simulation, the null distribution is flat; as signal appears, small -values become more common and the share below climbs from about 5% (that's ) to your study's actual power.

Fact two: keep looking without correction, and the false-alarm risk grows

The optional-stopping toggle in the main factory hid a result that ought to be more famous than it is. Suppose you test after every new subject and stop the instant dips below .05. If there’s a real effect, fine. But if there isn’t — if the null is exactly true — and you’re allowed to keep sampling indefinitely, then the probability that you eventually cross .05 is not 5%; in the mathematical idealization, it approaches 100%. This is “sampling to a foregone conclusion,” proved by Anscombe in 1954 and hammered home by Armitage and colleagues in 1969: the wandering -value, given enough patience, is liable to stumble below any fixed line at some point, purely by chance, even when nothing is there.

This is why peeking is not a venial sin but a structural exploit — and why Day 4’s e-values and anytime-valid tests are such a clean fix: they’re built so that the guarantee survives unlimited peeking. It’s also why real clinical trials, which must monitor accumulating data for safety, use alpha-spending and group-sequential designs — they pay a small statistical tax at each interim look so the total false-alarm budget still adds up to 5%. Optional stopping done lawlessly is cheating; done with a spending rule, it’s one of the triumphs of 20th-century statistics.

§3 Reading the fingerprints

The forensic turn

If a real effect leaves many very small significant -values, while p-hacking leaves a pile-up just below .05, then you can audit a whole literature from its published numbers alone — without rerunning a single experiment. That realization launched a small, slightly gleeful field of statistical forensics.

The flagship tool is p-curve (Simonsohn, Nelson & Simmons, 2014 — yes, the same trio from the main lesson). Collect the significant -values from a set of studies on one effect and look at their shape. A healthy p-curve has more very small significant -values than barely significant ones: many .01s, fewer .04s, because real signal mostly yields convincingly small -values. A literature with no underlying effect, propped up entirely by p-hacking, shows the suspicious pattern instead — an unnatural surplus of .04s and .049s, the telltale residue of researchers nudging results just over the line and stopping the moment they cleared it. A pile-up just below .05 is suspicious.

And then there are the sting operations — demonstrations built to make the multiple-comparisons problem impossible to wave away:

🐟 The dead salmon

Bennett et al., 2009 · Ig Nobel 2012

Neuroscientists slid a dead Atlantic salmon into an fMRI scanner and “showed” it photographs of human social situations. With no correction for the ~130,000 voxels tested, a cluster of the deceased fish’s brain lit up with “significant” activity. The point: test enough things and noise will counterfeit a discovery — in a corpse.

🪄 Voodoo correlations

Vul et al., 2009

Many celebrated brain-and-emotion studies reported correlations so high they exceeded the reliability of their own measures — impossible. The culprit was circular analysis: pick the voxels that correlate most, then report that correlation as if it were independent. The selection is the result. A forking path wearing a lab coat.

🔢 GRIM & the sleuths

Brown & Heathers, 2017

The GRIM test checks whether a reported average is even arithmetically possible for the stated sample size and integer answers. Astonishingly often, it isn’t. A loose collective of “data detectives” now uses such tools to surface impossible numbers in published papers.

🤖 statcheck

Nuijten et al., 2016

An algorithm that re-computes the -values in psychology papers from the reported test statistics. Scanning a quarter-million results, it found that roughly half of papers contained at least one inconsistent -value, and about one in eight an error large enough to potentially flip a conclusion.

Hype filter: p-curve and these forensic tools are established as diagnostics — but they detect patterns, not guilt; an innocent literature with publication bias can mimic some signatures, so they flag suspicion, not proof.

§4 The full arithmetic

How wrong is a “significant” result, really?

Now for the number that should be printed on the first page of every statistics textbook and almost never is. You ran a clean study, no hacking, and got . Intuition whispers: “about a 5% chance this is a fluke.” Intuition is catastrophically wrong, and we can prove it with nothing more than Day 4’s Bayes.

A -value lives in the world of . What you actually care about is — the chance your shiny finding is bogus. Bridging the two requires a prior: how plausible was your hypothesis before the data? And the moment you supply a realistic prior, the picture darkens. The engine is the same one Ioannidis used in his famous Day 2 essay “Why Most Published Research Findings Are False”: combine the prior odds, the study’s power, and , and read off the false positive risk — the share of your “significant” results that are false alarms.

Use the calculator to see how often significant findings are wrong—even with a coin-flip prior and decent power, and especially when the hypothesis is a long shot.

Interactive · the false-positive risk engine

Of 100 "discoveries," how many are false?

Imagine many hypotheses tested at the usual . Set how often your hypotheses are really true (the prior) and your study's power. The grid shows 100 results you'd call "significant," coloured by whether they're genuine (teal) or false alarms (red).

of your "significant" findings are false
true discovery false alarm

Two refinements make the news worse, not better — and they’re the real essence here.

It’s grimmer than the grid suggests (Colquhoun). The calculator above uses the whole .05 tail as the false-positive rate. But David Colquhoun (2014, 2017) pointed out that a result sitting right at is barely more likely under a real effect than under the null, so its calibrated false-positive risk is higher than the tail suggests. His verdict, for an even-odds prior: observe and the minimum chance you’re wrong is around 26% — not 5%. Demand a genuine discovery rate, and simply isn’t a high enough bar. This is the quantitative heart of the Benjamin et al. push toward .005 you met in the main lesson.

The null is almost never exactly true anyway (Meehl’s crud factor). Paul Meehl observed that in large social-science datasets, everything correlates a little with everything — birth order with politics, height with humor — because the world is a dense web of faint real connections he called the crud factor. So a point null of “exactly zero” is a straw man; with a big enough sample you can reject it for any pair of variables you like. “Statistically significant” then certifies only that your sample was large, not that you found anything that matters. It’s the deepest reason significance and importance are different animals.

When Bayes and frequentism openly fight: the Jeffreys–Lindley paradox

Here’s a genuinely unsettling one. With a very large sample, a result at exactly can be evidence for the null on a Bayesian analysis — the two frameworks don’t just differ, they point in opposite directions (Lindley, 1957). The frequentist rejects; the Bayesian shrugs and leans toward . The resolution isn’t that one side is dumb; it’s that they’re answering different questions, and at large n the gap between “surprising under the null” and “more probable under the alternative” yawns wide. The modern Bayesian answer is the Bayes factor — a direct ratio of how well each hypothesis predicted the data, popularized in accessible software like JASP. It can do something a -value constitutionally cannot: gather evidence for the null, not just fail to reject it.

§5 The better toolkit

What the people who know this reach for

Demolition is easy; the appendix owes you the rebuild. The estimation-era toolkit isn’t exotic — it’s a handful of moves that quietly fix the specific failures above.

Better errors: Type S and Type M

Andrew Gelman and John Carlin (2014) argue that in the low-power world we actually inhabit, the classic Type I / Type II framing asks the wrong question. Two sharper failures matter more. A Type S (“sign”) error is getting the direction wrong — concluding a treatment helps when it harms. A Type M (“magnitude”) error is the exaggeration ratio — how many times too large your “significant” estimate is. Their punchline ties straight to the main lesson’s winner’s curse: when power is low, the few estimates that reach significance aren’t just uncertain, they’re systematically inflated, often two- or threefold. The published effect isn’t a slightly noisy truth; it’s a funhouse-mirror enlargement of it.

Proving a negative: equivalence testing

The main lesson warned that "" does not mean “no effect.” But sometimes you genuinely want to conclude there’s nothing meaningful here — and you can, with equivalence testing (the TOST procedure; Lakens, 2017). You define the smallest effect you’d care about, then test whether the true effect is confidently smaller than that. It flips absence-of-evidence into evidence-of-absence, legitimately. A non-significant result plus a passed equivalence test is a real conclusion; a non-significant result alone is just a shrug.

Facing Multiplicity: Controlling the False Discovery Rate

When you really must run many tests (every gene, every brain region), the old fix was Bonferroni — divide by the number of tests — which is correct but so brutal it leaves you no power. Benjamini and Hochberg’s false discovery rate (1995), one of the most-cited papers in all of science, changed the target: instead of guaranteeing zero false positives among all tests, control the expected proportion of false positives among the ones you flag. It’s the pragmatic bargain that made genomics and neuroimaging statistically livable — accept a known, small fraction of duds in exchange for actually finding things.

Renaming the interval to mean what it is

Finally, a small move with large clarifying power: Amrhein, Greenland and colleagues urge calling a confidence interval a compatibility interval — the range of effect sizes reasonably compatible with your data and model. The word kills two errors at once. It stops you treating values inside as proven and values outside as impossible, and it nudges you to report the whole interval as a statement about magnitude and uncertainty rather than a yes/no membership test for zero.

§6 The case for the defence

Significance testing, done severely

It would be too tidy to end with “p-values bad, estimation good.” Evenhandedness — and accuracy — demand a serious counterweight, and the philosopher Deborah Mayo supplies it. In Statistical Inference as Severe Testing (2018), she defends error-statistical reasoning against both the abolitionists and the Bayesians, with a single demanding principle.

You have evidence for a claim only to the degree it has passed a test it would probably have failed, had the claim been false.

That is severity, and notice it is Day 2’s Popper rendered quantitative: a result corroborates a hypothesis only if the test genuinely risked refuting it. On this view, the rot isn’t the -value — it’s the abandonment of severity. A -value squeezed out of a forking-paths fishing trip is worthless not because it’s a -value but because the “test” was never severe: it could hardly have come out any other way. Mayo’s warning to the reformers is pointed: drop all thresholds and error-control, and you don’t get more nuanced science — you get less accountable science, with no public standard for when a claim has actually earned its keep. The synthesis the main lesson reached for lives here: not “abolish” or “obey,” but demand severity, report magnitude and uncertainty, and make the test one the hypothesis could have flunked.

§7 The bridge to tomorrow

A p-value, measured in coin flips

One last idea, and it’s a beauty — because it quietly turns the page to Day 7. Sander Greenland (2019) suggests we stop staring at the bare -value and instead take its surprisal, borrowed straight from information theory: the S-value, defined as . This converts the p-value into bits of surprisal, expressing how incompatible the data are with the tested model under its assumptions — and bits, as we’ll see tomorrow, are the universal currency of surprise.

The translation is wonderfully concrete. Ask: how many fair coins would I have to flip and get all heads to be this surprised? That count is the S-value. A result at carries about 4.3 bits — roughly the surprise of calling four or five coin tosses in a row. Stated that way, the spell breaks. Four heads in a row is mildly interesting; it is not a discovery, and you would never bet your reputation on it.

Move the slider to translate any -value into Day 7’s information units.

Interactive · surprise in bits

The S-value

Drag to set a -value. The S-value, , tells you how many coins-all-heads worth of surprise it represents — the real evidence against the tested model, in bits.

4.3bits of surprise  ≈  4–5 coins, all heads

What to actually do

The working scientist’s defence kit

Distilled from everything above and in the main lesson — the moves that convert luck into knowledge:

  • Preregister the confirmatory test. Freeze the confirmatory path before you see the data; label everything else exploratory.
  • Report the effect size with a compatibility interval, always. The magnitude and its uncertainty are the result; the -value is a footnote.
  • Power the study for the smallest effect you'd care about — and remember that a “significant” result from an underpowered study is probably inflated (Type M).
  • Translate your p into bits. If it’s only a few coin flips’ worth of surprise, say so, and don’t oversell.
  • Check robustness across reasonable analyses (a specification curve), and ask whether the result survived a severe test or merely a convenient one.
  • When you must test many things, control the false discovery rate — not the family-wise error into oblivion.
  • If you want to claim "no effect," run an equivalence test. Absence of evidence isn’t evidence of absence until you’ve actually looked for absence.

The appendix in three sentences

Deepest idea
The ” test” is an incoherent weld of Fisher’s evidence-measure and Neyman–Pearson’s decision-rule, and the misreadings everyone makes are symptoms of that original sin. In the continuous exact-test idealization, null p-values are uniform; under indefinite, uncorrected optional stopping, the chance of eventually crossing a fixed threshold can approach one. Those conditions explain how repeated analytic flexibility inflates false positives without making every analysis a guaranteed false discovery.
Best new analogy
A -value translated into coin flips: is the surprise of ~4–5 heads in a row — mildly interesting, never a discovery. Under Colquhoun’s specific even-odds prior, power, and model setup, an observed can correspond to a false-positive risk around 26%; that posterior risk is not a fixed meaning of p.
The constructive turn
The fix is a toolkit, not a ban: severity (Mayo), effect sizes with compatibility intervals, Type S/M thinking, equivalence tests, false-discovery-rate control, and Bayes factors — plus the forensic methods (p-curve, the dead salmon, statcheck) that audit a literature from its own numbers.

Sources · appendix

Sources & further reading

  1. Fisher, R. A. (1925). Statistical Methods for Research Workers. Oliver & Boyd. — the p-value as a continuous evidential measure; .05 as a flexible convention.
  2. Neyman, J. & Pearson, E. S. (1933). “On the Problem of the Most Efficient Tests of Statistical Hypotheses.” Philosophical Transactions of the Royal Society A 231: 289–337. — testing as a decision rule with α, β, and power.
  3. Gigerenzer, G. (2004). “Mindless statistics.” The Journal of Socio-Economics 33(5): 587–606. — the “null ritual” as an incoherent Fisher/Neyman–Pearson hybrid.
  4. Goodman, S. N. (1999). “Toward Evidence-Based Medical Statistics. 1: The P Value Fallacy.” Annals of Internal Medicine 130(12): 995–1004. doi:10.7326/0003-4819-130-12-199906150-00008.
  5. Anscombe, F. J. (1954). “Fixed-Sample-Size Analysis of Sequential Observations.” Biometrics 10(1): 89–100. — sampling to a foregone conclusion.
  6. Armitage, P., McPherson, C. K. & Rowe, B. C. (1969). “Repeated Significance Tests on Accumulating Data.” Journal of the Royal Statistical Society A 132(2): 235–244. — optional stopping crosses any threshold under H₀ with probability → 1.
  7. Simonsohn, U., Nelson, L. D. & Simmons, J. P. (2014). “P-curve: A Key to the File-Drawer.” Journal of Experimental Psychology: General 143(2): 534–547. doi:10.1037/a0033242.
  8. Bennett, C. M., Baird, A. A., Miller, M. B. & Wolford, G. L. (2010). “Neural Correlates of Interspecies Perspective Taking in the Post-Mortem Atlantic Salmon.” Journal of Serendipitous and Unexpected Results 1: 1–5. (HBM 2009 poster; Ig Nobel 2012) — multiple comparisons in fMRI.
  9. Vul, E., Harris, C., Winkielman, P. & Pashler, H. (2009). “Puzzlingly High Correlations in fMRI Studies of Emotion, Personality, and Social Cognition.” Perspectives on Psychological Science 4(3): 274–290. — “voodoo correlations”; non-independence error.
  10. Brown, N. J. L. & Heathers, J. A. J. (2017). “The GRIM Test: A Simple Technique Detects Numerous Anomalies in the Reporting of Results in Psychology.” Social Psychological and Personality Science 8(4): 363–369. doi:10.1177/1948550616673876.
  11. Nuijten, M. B., Hartgerink, C. H. J., van Assen, M. A. L. M., Epskamp, S. & Wicherts, J. M. (2016). “The prevalence of statistical reporting errors in psychology (1985–2013).” Behavior Research Methods 48: 1205–1226. doi:10.3758/s13428-015-0664-2 — the statcheck audit.
  12. Ioannidis, J. P. A. (2005). “Why Most Published Research Findings Are False.” PLoS Medicine 2(8): e124. — the positive-predictive-value engine (revisited from Day 2).
  13. Colquhoun, D. (2014). “An investigation of the false discovery rate and the misinterpretation of p-values.” Royal Society Open Science 1: 140216. doi:10.1098/rsos.140216.  Colquhoun, D. (2017). “The reproducibility of research and the misinterpretation of p-values.” Royal Society Open Science 4: 171085. — observed false-positive risk .
  14. Lindley, D. V. (1957). “A Statistical Paradox.” Biometrika 44(1–2): 187–192. — the Jeffreys–Lindley divergence of p-values and Bayes factors at large n.
  15. Meehl, P. E. (1990). “Why Summaries of Research on Psychological Theories Are Often Uninterpretable.” Psychological Reports 66(1): 195–244. — the “crud factor”; the point null is almost never exactly true.
  16. Gelman, A. & Carlin, J. (2014). “Beyond Power Calculations: Assessing Type S (Sign) and Type M (Magnitude) Errors.” Perspectives on Psychological Science 9(6): 641–651. doi:10.1177/1745691614551642.
  17. Lakens, D. (2017). “Equivalence Tests: A Practical Primer for t Tests, Correlations, and Meta-Analyses.” Social Psychological and Personality Science 8(4): 355–362. doi:10.1177/1948550617697177 — the TOST procedure.
  18. Benjamini, Y. & Hochberg, Y. (1995). “Controlling the False Discovery Rate.” Journal of the Royal Statistical Society B 57(1): 289–300. — FDR vs family-wise error.
  19. Mayo, D. G. (2018). Statistical Inference as Severe Testing: How to Get Beyond the Statistics Wars. Cambridge University Press. — severity; the error-statistical defence.
  20. Greenland, S. (2019). “Valid P-Values Behave Exactly as They Should: Some Misleading Criticisms of P-Values and Their Resolution With S-Values.” The American Statistician 73(sup1): 106–114. doi:10.1080/00031305.2018.1529625 — the S-value / surprisal, in bits.
Deep dive appendixThe Incoming WaveOptional extension.

The main lesson taught the old danger: a p-value plus human flexibility can turn noise into literature. This appendix looks forward. The sharpest post-2020 work is not merely “use AI carefully” or “report confidence intervals.” It is a deeper shift: validity is being moved outside the model, outside the analyst’s discretion, and sometimes outside the raw data itself. The new tools wrap black boxes in guarantees, split exploration from inference, keep tests valid while people peek, correct AI labels with human truth, quantify privacy noise, and audit the research machine itself. They are powerful. They are also easy to oversell.

Scope and ground rules

This version keeps the Day 6 promise: no method is treated as magic. Each wave gets a maturity label, a misuse warning, and a callback to the p-hacking / base-rate / causation lessons from Days 1–5. Older foundations are included only when the post-2020 wave is the part that changed practice.

Field map

Ten waves, sorted by the failure mode they attack

Filter the atlas by maturity. Green means the core theory is established; amber means serious and promising but still settling; red means potentially transformative and currently dangerous.

Wave 01

Anytime-valid inference

Tests and intervals that survive optional stopping, continuous monitoring, and sequential accumulation.

established
Wave 02

Conformal risk control

Distribution-free wrappers for prediction sets, abstention, and general risk metrics around black-box models.

core established
Wave 03

Prediction-powered inference

Use cheap AI predictions without pretending they are truth; correct them with a small gold-standard sample.

Science / PNAS
Wave 04

Orthogonal causal ML

Let flexible ML learn nuisance functions while orthogonal scores protect the causal estimate.

maturing
Wave 05

Inference after discovery

Post-selection inference, data fission, and causal-discovery corrections for exploration without double dipping.

fast frontier
Wave 06

Private-data inference

Confidence intervals that include the uncertainty induced by differential privacy and synthetic-data release.

growing
Wave 07

Bayesian workflow

Prior predictive checks, posterior predictive checks, SBC, model criticism: inference as an audited workflow.

practice shift
Wave 08

Forensics at scale

Data detectives, decision markets, ML replicability forecasts, paper-mill detection, and journal-level triage.

promising
Wave 09

Silicon samples

LLMs as simulated survey respondents or experimental participants: useful for piloting, dangerous as replacement.

contested
Wave 10

AI analysts & benchmarks

Automated analysis agents, benchmark contamination, reproducibility drift, and “AI scientist” pipelines.

unsettled

§1 The p-value that survives the refresh button

Anytime-valid inference: the end of illegal peeking

Wave 01Theory establishedAdoption still uneven

The old fixed-sample test is fragile because it assumes you decided the sample size before seeing the data. But modern science does not behave like that. A/B tests stream in hourly. Platform teams refresh dashboards. Clinical trials have interim looks. Labs add subjects when the first batch is “almost significant.” The ordinary p-value is not built for this world.

Anytime-valid inference changes the contract. Instead of a p-value that is valid only at a pre-fixed stopping time, it uses e-values, e-processes, always-valid p-values, and confidence sequences: evidence measures and intervals designed to remain valid at arbitrary stopping times. You can look early, look often, stop when the evidence is decisive, or continue when it is ambiguous, without silently inflating the false-positive rate.

This is not a decorative upgrade. It attacks one of Day 6’s central machines for fooling yourself: optional stopping. The frontier after 2020 is the consolidation of a “game-theoretic statistics / safe anytime-valid inference” language. It treats evidence as a capital process in a game against nature. If the null is true, your evidence wealth cannot systematically grow without bound; Ville-style inequalities turn that into sequential tests and confidence sequences.

The hype filter: the mathematics is now serious and mature. The practical risk is cultural. People like dashboards but dislike changing the number they report. If a team runs a continuously monitored experiment and still reports ordinary fixed-horizon p-values, the UI may look modern while the inferential contract is 1930s paper-and-pencil.

Day 6 callback

Optional stopping was one of the False-Positive Factory switches. Anytime-valid inference does not say “never peek.” It says: if peeking is part of the workflow, put peeking into the statistical procedure rather than pretending it did not happen.

Interactive · optional stopping

The peek meter

A toy approximation: treat each look as a fresh chance to cross . Real sequential dependence is subtler, but the direction is the point. Fixed-sample p-values are brittle under repeated looks; anytime-valid procedures are designed for the stopping rule.

fixed-sample p46%
anytime-valid5%

The red bar is the chance that at least one look produces a false “discovery” under the simple independent-looks approximation. The green bar stays at the designed false-alarm level.

§2 The guarantee that wraps a black box

Conformal prediction becomes conformal risk control

Wave 02Core method establishedRisk-control frontier

Conformal prediction is the simplest profound idea in the current statistics/ML interface. Train any predictive model. Hold out calibration data. Measure how wrong the model is on those held-out cases. Use the empirical quantile of those errors to form a prediction set for the next case. Under exchangeability, the set covers the true answer at the advertised rate in finite samples, regardless of whether the model is linear, neural, brilliant, or embarrassing.

A common simplification treats conformal prediction as “intervals for any model.” That is correct but incomplete. The post-2020 frontier is broader: conformal risk control asks not merely “does the true label fall in the set?” but “can we control a chosen risk function?” In vision and language tasks, the useful error metric may be false negative rate, graph distance, abstention error, retrieval miss-rate, or token-level F1. Conformal risk control extends the conformal idea to monotone losses and gives finite-sample control of the expected risk.

That matters for AI systems because the output is often not a scalar. A chatbot answer, segmentation mask, medical image region, or structured extraction is not a single yes/no label. Risk control lets you build a calibrated ship / abstain / escalate / ask for review wrapper around a black-box model. In the safest version, the model remains a prediction engine; the conformal wrapper becomes the trust policy.

The hype filter: marginal coverage has the advertised guarantee; the harder target is conditional coverage. Under distribution shift, the classical exchangeability guarantee weakens. Recent work handles weighted and non-exchangeable settings, conditional-style guarantees, and e-value variants, but the practitioner still has to ask whether the calibration data resemble deployment. Conformal methods make uncertainty explicit. They do not make a bad deployment context safe.

Interactive · conformal bargain

Coverage is bought with set size

This is a conceptual dial, not a theorem. Lower α asks for higher coverage; weaker model skill and more shift make the prediction set larger or less trustworthy. The real guarantee depends on exchangeability or a justified shift correction.

target coverage90%
stress coverage90%
set size42%

A uselessly wide set can be valid. Validity is a safety net, not a hammock.

§3 Using AI predictions without swallowing their errors

Prediction-powered inference: the rectifier is the point

Wave 03Science 2023PNAS 2024 extensionsAdoption growing

AI can label a million galaxies, satellite tiles, pathology patches, protein records, or clinical notes. But a statistical estimate from those predictions is not automatically valid, because model errors are not random clerical typos. They are structured. They follow the model’s blind spots.

Prediction-powered inference (PPI) solves the problem by refusing to trust the model alone. You combine a small labeled sample, where humans or instruments provide gold-standard truth, with a much larger unlabeled sample on which the model predicts. The labeled sample estimates a rectifier: the correction for the model’s bias on the estimand you actually care about. The unlabeled sample supplies cheap precision; the labeled sample supplies validity.

PPI++ improves the efficiency of this idea by adapting to prediction quality. If the model is useful, the interval shrinks. If the model is weak, the method falls back toward the labeled-only interval rather than letting bad predictions poison the estimate. Cross-prediction-powered inference goes further by distributing prediction and correction across folds. The architecture is exactly the incoming-wave pattern: use the black box for power, but put the guarantee in the wrapper.

The limitation is practical, not philosophical. You still need a gold-standard labeled sample drawn from the right target population. You still need to define the estimand before the model seduces you into measuring whatever it can label. And if the labeled set misses the subgroup where the model fails, the rectifier can inherit that blind spot.

Interactive · the PPI bargain

Cheap predictions buy precision only after gold-standard correction

Move the sliders to compare three intervals: classical labeled-only, naive ML-as-truth, and PPI. The naive interval can be narrow and wrong. PPI aims to be narrower than labeled-only while retaining validity through the rectifier.

labeled-only
naive ML
PPI
labels only naive prediction PPI truth

§4 Machines for nuisance, humans for assumptions

Double / debiased machine learning for causal effects

Wave 04Theory establishedTooling mature

Day 5’s warning was blunt: association is not intervention. To estimate a causal effect from observational data, you must adjust for confounders, and you must do it without imposing a fantasy linear model on a nonlinear world.

Double/debiased machine learning (DML) uses ML where ML is strongest: learning high-dimensional nuisance functions such as outcome regressions and propensity scores. Then it uses orthogonal scores and cross-fitting where statistics is strongest: protecting a low-dimensional target estimate from first-order errors in those nuisance fits. You let the random forest, gradient booster, or neural net approximate the messy parts, then construct the causal estimator so that small nuisance mistakes wash out.

The post-2020 change is not that DML was invented after 2020; it was not. The change is that software and applied practice caught up. DoubleML implementations in Python and R, causal-forest workflows, and biomedical heterogeneous-treatment-effect reviews have turned “orthogonalize the nuisance” from a seminar phrase into an applied workflow.

The hype filter is essential: DML weakens modeling assumptions; it does not remove causal identification assumptions. If the back-door set is wrong, if a confounder is unmeasured, if treatment assignment depends on hidden variables, a beautiful orthogonal score will compute a precise answer to the wrong causal question. DML solves flexible adjustment, not the metaphysics of causation.

Day 5 callback

DML helps estimate only after the design and identification argument say that this target is reachable. It is not a machine for converting into causation by force.

§5 Explore first, infer second, pay the bill

Inference after discovery: post-selection, causal discovery, and data fission

Wave 05Fast-moving frontierCore problem settled

The garden of forking paths is not going away. Modern analysis often begins with discovery: cluster the cells, find the subgroup, select the variables, learn the causal graph, choose the image region, tune the transformation, ask the AI analyst to “find what matters.” The fatal mistake is to use the same data to discover a target and then pretend the target was specified in advance.

Post-selection inference is the statistical attempt to account for that bill directly. The classical approaches condition on the selection event, split the data into discovery and confirmation sets, or build simultaneous intervals wide enough to survive selection. The post-2020 frontier adds three especially relevant ideas.

A

Algorithmic stability. If the selection rule is randomized and stable, its sensitivity to any one observation can be quantified. That stability can be turned into post-selection corrections, borrowing ideas from differential privacy and adaptive data analysis.

B

Data fission and data thinning. Instead of splitting rows, split the statistical information inside the data into pieces. One piece explores; the other tests. For some distributions, recent thinning methods produce independent pieces without throwing away whole observations.

C

Valid inference after causal discovery. Running a graph-discovery algorithm and then estimating effects on the same data can break coverage. New work randomizes and corrects discovery so the subsequent causal intervals are not fantasy intervals.

The hype filter: this is not easy, and it is not a substitute for preregistration when preregistration is possible. But it is the right frontier for exploratory science because “never explore” is not a serious rule. The serious rule is: if the data suggested the question, the inference must know that.

Interactive · AI analyst multiverse

The garden grows faster when analysis is cheap

Assume pure noise and a nominal . Each prompt, model, outcome, and cleaning rule multiplies the number of possible analyses. The independence approximation is crude; the moral is not.

240 analysis paths
false-positive risk≈100%
preregistered5%

An AI analyst that reports the full multiverse can be a transparency machine. An AI analyst that silently selects the prettiest path is a p-hacking accelerator.

§6 When the released data are not the data

Differential privacy turns uncertainty into a public design problem

Wave 06High-stakes frontierPolicy-relevant

Open science wants data sharing. Human subjects, firms, hospitals, and governments want confidentiality. Differential privacy is the strongest modern language for this tradeoff: random noise is added so that one person’s presence or absence has bounded influence on the released output.

For Day 6, the statistical sting is this: privacy noise is not a file-format detail. If you analyze differentially private statistics or synthetic data as if they were ordinary observations, your intervals can be too narrow, your p-values too confident, and your null findings too frequent. The privacy mechanism is part of the data-generating process, so valid inference has to include it.

Post-2020 work has been building noise-aware confidence intervals, private bootstraps, inference from DP synthetic data, and guidelines for evaluating whether a claimed privacy guarantee means what users think it means. This is not glamorous, but it is central. In the next decade, many researchers will never touch the raw administrative or medical data. They will touch privacy-filtered outputs. The frontier question becomes: can we make those outputs scientifically useful without making people identifiable?

The hype filter: differential privacy protects individuals by making some analyses less precise. A valid private interval may be much wider. That is not failure; that is the price tag made visible. The dangerous failure is a narrow interval that forgot the noise.

§7 The model is a workflow, not a formula

Bayesian workflow and simulation-based calibration

Wave 07Practice shiftComputation-sensitive

Bayesian statistics is often sold as a clean alternative to p-values: put a prior on the unknowns, compute a posterior, report probabilities about parameters. That is attractive, but it can become another ritual if the prior is unexamined, the model is misspecified, and the sampler is silently broken.

The post-2020 Bayesian workflow movement says: the unit of trust is not the posterior; it is the workflow that produced and criticized the posterior. Prior predictive checks ask whether the model could generate plausible data before seeing the data. Posterior predictive checks ask what the fitted model still cannot reproduce. Cross-validation asks whether predictions survive held-out cases. Simulation-based calibration checks whether the inference algorithm can recover known parameters when the data are simulated from the model.

This belongs in the cutting-edge appendix because it is a direct answer to Day 1’s stopped-clock worry. A posterior can be numerically impressive for accidental reasons. A workflow makes the path visible: what was assumed, what was simulated, what failed, what changed, and what uncertainty remains.

The hype filter: model checking is not model proof. Posterior predictive checks can miss the discrepancy that matters. SBC validates computation under the assumed generative model, not the truth of the model itself. But compared with a lone posterior interval, a documented workflow is a much harder machine to fool.

§8 The error-correction system grows teeth

Statistical forensics and replicability triage at scale

Wave 08Forensics validatedScale still early

Appendix I’s forensic toolkit was artisanal: GRIM checks, statcheck, p-curve, digit anomalies, spreadsheet oddities. The post-2020 change is that the stakes and scale have grown.

The Francesca Gino / Data Colada case showed both the power and the danger. Data Colada published detailed anomaly analyses in 2023. Harvard’s investigation concluded research misconduct; Gino denied wrongdoing and sued. In September 2024, a federal judge dismissed all claims against the Data Colada defendants; in May 2025 Harvard revoked Gino’s tenure and terminated her employment, while some legal disputes involving Harvard continued. The scientific lesson is not “bloggers should run universities.” It is that transparent forensic arguments can become part of the formal error-correction system.

The next layer is triage. We cannot replicate every claim. Decision markets have shown that scientists’ aggregated forecasts can identify high- and low-replicability findings: in one Nature Human Behaviour study, top-market studies replicated far more often than bottom-market studies. The Center for Open Science’s Predicting Replicability Challenge adds machine-learning models trained on thousands of prior replication outcomes. The early results are mixed: Round 1 was weak by some metrics; Round 2 improved. That nuance matters. Replicability prediction is a promising smoke alarm, not a truth machine.

The hype filter: forensic tools can be misused as accusation engines. A weird dataset is not automatically fraud; a low replicability forecast is not a verdict. The frontier needs transparent code, calibrated false-alarm rates, and institutional protections for both whistleblowers and accused researchers.

Timeline · maturity and caution

The wave at a glance

2021–2024

Anytime-valid inference consolidates

Always-valid A/B testing, e-processes, confidence sequences, and safe testing become a coherent post-p-value language for sequential science.

2023–2025

Conformal leaves the classroom

Gentle tutorials, beyond-exchangeability theory, conformal risk control, conditional-style guarantees, and NLP/LLM applications widen the use case.

2023–2024

Prediction-powered inference lands

Science paper, PPI++, and cross-PPI formalize valid estimation with cheap model predictions plus a small gold-standard sample.

2023–2025

Inference after discovery gets practical

Algorithmic stability, data thinning/fission, and valid inference after causal discovery directly address double dipping in exploratory workflows.

2022–2026

Privacy-aware inference becomes unavoidable

DP synthetic data and private confidence intervals force analysts to include privacy noise in uncertainty, rather than treating released data as raw data.

2023–2026

LLMs enter the research pipeline

Silicon samples, AI data analysts, benchmark contamination, and AI-scientist systems create a new replication problem: not just what did humans choose, but what did the machine choose?

§9 Synthetic humans are not humans

Silicon samples and LLM social simulations

Wave 09Useful for pilotingContested as evidence

Silicon sampling asks an irresistible question: if a language model has absorbed enough social text, can it act as a simulated survey respondent or experimental participant? The first wave of results was exciting. Conditioned on demographic backstories, models sometimes reproduced aggregate opinion patterns. Agent-based simulations produced believable local interactions. Large-scale replications of psychology and management scenarios found that LLMs often matched the direction of main effects.

The second wave is the caution. LLMs can distort correlations even when margins look right. They often compress distributions toward modal answers. They may amplify effect sizes. They can perform worse on socially sensitive topics. Their “population” is a residue of training data, alignment tuning, prompt wording, and model version. That is not a sample frame. It is a machine-shaped mirror.

The defensible use case is therefore narrow but real: pilot instruments, debug ambiguous survey wording, estimate rough sample-size needs, test whether a manipulation is even legible, and generate hypotheses before spending human-subjects budgets. The reckless use case is confirmatory evidence without human validation.

Reporting standards should become severe: model name and version, date, system prompt, persona construction, temperature, sampling scheme, repeated-run variability, validation against human data, and a preregistered rule for when the simulation is allowed to count. Otherwise silicon samples become the garden of forking paths wearing a human mask.

§10 When the analyst is silicon too

AI analysts, automated science, and benchmark contamination

Wave 10Deeply unsettledTooling advancing fast

The most dangerous frontier is not AI as a model inside an analysis. It is AI as the entity that chooses the analysis.

Data-analysis agents can clean data, write code, choose models, run tests, generate plots, and narrate findings. New benchmarks such as StatLLM evaluate statistical code; data-agent benchmarks evaluate longer workflows. Studies of LLM data analysis show variation across models, prompts, temperatures, and repeated executions; even the same task and dataset can lead to meaningfully different estimates. At the far edge, autonomous research systems now generate ideas, write code, run experiments, plot results, draft papers, and perform self-review.

This does not make AI useless. It makes AI analysis statistically dangerous in a very familiar way. The many-analysts problem has become cheap enough to run in seconds. The “researcher degrees of freedom” are now prompt degrees of freedom, model degrees of freedom, retry degrees of freedom, tool-call degrees of freedom, and hidden chain-of-reasoning degrees of freedom.

AI benchmarking has the same pathology at field scale. Public leaderboards are repeated tests. Benchmark contamination happens when static benchmarks leak into training data. Developers can tune to the benchmark. Scores can rise because models generalize better, because they memorized more, because the harness changed, or because the benchmark itself is exhausted. The statistical object is no longer “accuracy on a test set”; it is a dynamic game between model builders, benchmark maintainers, and hidden evaluation data.

The responsible frontier looks like this: prompt logs, containerized code, repeated independent agent runs, multiverse reports, hidden holdout evaluations, contamination audits, preregistered analysis instructions, and human review of the decisions the system made. The dream is an AI analyst that exposes every fork. The nightmare is an AI analyst that silently chooses the fork that gets published.

Day 1 callback

AI can give an illusion of understanding: more outputs, more polish, more confidence, less actual contact with why the claim is true. The antidote is not less computation. It is more audit.

Interactive · benchmark games

When the leaderboard is not deployment

A toy model of benchmark inflation. Contamination and repeated submissions can raise a public score without improving hidden deployment performance.

hidden holdout72%
public board80%
generalization gap8%

Healthy evaluation rotates or hides tests, audits contamination, limits leaderboard overfitting, and reports uncertainty around rank differences.

The through-line

The incoming wave is not “trust AI.” It is “wrap discretion.”

The strongest post-2020 methods share a pattern. They accept that the model may be a black box, the analyst may explore, the experiment may be monitored continuously, the data may be privatized, and the research system may be too large for manual checking. Then they ask: where can a guarantee still be attached?

Failure modeOld temptationNew wrapperWhat still can fail
Optional stoppingPeek until p < .05Anytime-valid p-values, e-processes, confidence sequencesPower, implementation, and whether people actually use the right procedure
Black-box predictionsTreat AI outputs as truthConformal sets, conformal risk control, PPI rectifiersExchangeability, calibration coverage, bad gold labels, deployment shift
Flexible causal adjustmentChoose a convenient regressionDML, orthogonal scores, cross-fittingUnmeasured confounding and wrong identification arguments
Exploratory discoverySelect then infer as if preplannedPost-selection inference, data splitting/fission/thinning, stability correctionsPower loss, narrow applicability, unreported exploration
Privacy-filtered dataAnalyze released outputs as raw observationsNoise-aware DP inference, private bootstrap, valid DP CIsWide intervals, weak utility, misunderstood privacy budgets
Automated scienceLet the agent report the prettiest resultPrompt logs, repeated runs, hidden evals, multiverse reporting, forensic triageOpaque models, benchmark games, publication incentives

The original statistics crisis was not caused by evil people. It was caused by flexible humans using brittle procedures in incentive systems that rewarded clean stories. The next crisis will rhyme: flexible machines using brittle evaluations in incentive systems that reward impressive outputs. The cure is the same at a higher level: make the forks visible, attach guarantees where they belong, and refuse to mistake polish for contact with truth.

The appendix in three sentences

The wave
Statistics is being rebuilt for black-box, streaming, privacy-filtered, AI-assisted science: conformal wrappers, PPI rectifiers, orthogonal causal scores, anytime-valid evidence, post-selection correction, and workflow-level audits.
The hard lesson
The guarantee has to move to the part of the workflow that can actually bear it: the stopping rule, the calibration set, the gold-standard labels, the selection procedure, the privacy mechanism, or the audit log.
The immediate danger
LLM participants, AI analysts, and public benchmarks can recreate p-hacking at machine speed unless prompts, model versions, reruns, hidden tests, and full analytic multiverses become ordinary scientific records.

Sources · appendix II

Sources & further reading

  1. Ramdas, A., Grünwald, P., Vovk, V. & Shafer, G. (2023). “Game-Theoretic Statistics and Safe Anytime-Valid Inference.” Statistical Science. doi. Core survey of e-processes and confidence sequences.
  2. “Game-Theoretic Statistics and Safe Anytime-Valid Inference” editorial (2024), New England Journal of Statistics in Data Science. journal page.
  3. Johari, R., Pekelis, L. & Walsh, D. J. “Always Valid Inference: Continuous Monitoring of A/B Tests.” Operations Research. doi.
  4. Angelopoulos, A. N. & Bates, S. (2023). “Conformal Prediction: A Gentle Introduction.” Foundations and Trends in Machine Learning 16(4):494–591. doi.
  5. Barber, R. F., Candès, E. J., Ramdas, A. & Tibshirani, R. J. (2023). “Conformal Prediction Beyond Exchangeability.” Annals of Statistics. doi.
  6. Angelopoulos, A. N., Bates, S., Fisch, A., Lei, L. & Schuster, T. (2024). “Conformal Risk Control.” ICLR. OpenReview.
  7. Gibbs, I., Cherian, J. J. & Candès, E. J. (2025). “Conformal Prediction with Conditional Guarantees.” JRSS-B. doi.
  8. Gauthier, E., Bach, F. & Jordan, M. I. (2025). “E-Values Expand the Scope of Conformal Prediction.” preprint. arXiv.
  9. Campos, M. et al. (2024). “Conformal Prediction for Natural Language Processing: A Survey.” TACL. journal page.
  10. Angelopoulos, A. N., Bates, S., Fannjiang, C., Jordan, M. I. & Zrnic, T. (2023). “Prediction-Powered Inference.” Science 382:669–674. doi.
  11. Angelopoulos, A. N., Duchi, J. C. & Zrnic, T. (2023). “PPI++: Efficient Prediction-Powered Inference.” arXiv.
  12. Zrnic, T. et al. (2024). “Cross-Prediction-Powered Inference.” PNAS. doi.
  13. Chernozhukov, V. et al. (2018). “Double/Debiased Machine Learning for Treatment and Structural Parameters.” The Econometrics Journal. doi.
  14. Bach, P., Chernozhukov, V., Kurz, M. S. & Spindler, M. (2022). “DoubleML — An Object-Oriented Implementation of Double Machine Learning in Python.” JMLR. journal page.
  15. Bach, P. et al. (2024). “DoubleML — An Object-Oriented Implementation of Double Machine Learning in R.” Journal of Statistical Software. journal page.
  16. Abécassis, J. et al. (2025). “From Prediction to Prescription: Machine Learning and Causal Inference for the Heterogeneous Treatment Effect.” Annual Review of Biomedical Data Science. doi.
  17. Gradu, P., Zrnic, T., Wang, Y. & Jordan, M. I. (2025). “Valid Inference After Causal Discovery.” JASA. doi.
  18. Zrnic, T. & Jordan, M. I. (2023). “Post-Selection Inference via Algorithmic Stability.” Annals of Statistics. doi.
  19. Dharamshi, A. et al. (2025). “Generalized Data Thinning Using Sufficient Statistics.” JASA. doi.
  20. Neufeld, A., Dharamshi, A., Gao, L. L. & Witten, D. (2024). “Data Thinning for Convolution-Closed Distributions.” JMLR. pdf.
  21. Räisä, O. et al. (2023). “Noise-Aware Statistical Inference with Differentially Private Synthetic Data.” AISTATS. PMLR.
  22. Drechsler, J. et al. (2022). “Nonparametric Differentially Private Confidence Intervals for the Median.” Journal of Survey Statistics and Methodology. journal page.
  23. NIST (2025). Guidelines for Evaluating Differential Privacy Guarantees, SP 800-226. pdf.
  24. Gelman, A. et al. (2026). Bayesian Workflow. maintained edition.
  25. Modrák, M. et al. “Simulation-Based Calibration Checking for Bayesian Computation.” Bayesian Analysis. doi.
  26. Säilynoja, T., Schmitt, M., Bürkner, P.-C. & Vehtari, A. (2026). “Posterior SBC: Simulation-Based Calibration Checking Conditional on Data.” Statistics and Computing. doi.
  27. Argyle, L. P. et al. (2023). “Out of One, Many: Using Language Models to Simulate Human Samples.” Political Analysis. doi.
  28. Bisbee, J., Clinton, J., Dorff, C., Kenkel, B. & Larson, J. (2024). “Synthetic Replacements for Human Survey Data? The Perils of Large Language Models.” Political Analysis. doi.
  29. Cui, Z., Li, N. & Zhou, H. (2025). “A Large-Scale Replication of Scenario-Based Experiments in Psychology and Management Using Large Language Models.” Nature Computational Science. doi.
  30. Messeri, L. & Crockett, M. J. (2024). “Artificial Intelligence and Illusions of Understanding in Scientific Research.” Nature. doi.
  31. Holzmeister, F. et al. (2025). “Examining the Replicability of Online Experiments Selected by a Decision Market.” Nature Human Behaviour. doi.
  32. Center for Open Science (2026). “Predicting Replicability Challenge: Round 2 Results.” COS blog.
  33. Science (2024). “Honesty Researcher’s Lawsuit Against Data Sleuths Dismissed.” news.
  34. The Guardian (2025). “Harvard Professor Fired Following Claims She Falsified Ethics Research Data.” news.
  35. Lu, C. et al. (2026). “Towards End-to-End Automation of AI Research.” Nature. doi.
  36. Song, X. et al. (2026). “StatLLM: A Dataset for Evaluating the Performance of Large Language Models in Statistical Analysis.” Scientific Data. doi.
  37. Cui, J. & Alexander, R. (2026). “Same Prompt, Different Outcomes: Evaluating the Reproducibility of Data Analysis by LLMs.” arXiv.
  38. Deng, C. et al. (2024). “Investigating Data Contamination in Modern Benchmarks for Large Language Models.” NAACL. ACL Anthology.
  39. Xu, C. et al. (2024). “Benchmark Data Contamination of Large Language Models: A Survey.” arXiv.
  40. Choi, H. K. et al. (2025). “How Contaminated Is Your Benchmark? Measuring Dataset Leakage in Large Language Models with Kernel Divergence.” ICML. PMLR.

End of Day 006 · 174 descents remain