Science Journaling Club Founded 2024

EXPLAINER · EXPLAINER · STATISTICS

What a P Value Actually Means, and the Four Things People Think It Means Instead

Written jointly by the Science Journaling Club

Explainer · Peer-edited by the club review board · LaTeX source · Our calculation · Interactive model

Abstract A p value measures your data against a pretend world, the one where the effect you are hunting does not exist. The number describes data, given a hypothesis. Swap those two around and every misreading in this article follows. We define the thing properly, then take apart four beliefs the club has had to argue with at its own table: that p is the probability the hypothesis is wrong, that a smaller p means a bigger effect, that 0.049 and 0.051 are different animals, and that a non-significant result means nothing happened. Our own simulation of 200,000 experiments under a true null backs the definition, with the p values coming out flat to a Kolmogorov-Smirnov distance of 0.0014 from perfect uniformity. Another 400,000 simulate a whole literature. Now the headline. In a field where 10% of tested hypotheses are true and studies have 50% power, 47.4% of everything published as significant at \(p<0.05\) is false, with no fraud, no p-hacking, no hidden analyses and no arithmetic errors anywhere. Two honest studies of one real effect disagree half the time, and one line of arithmetic gives you both numbers. No results section reads the same way afterwards.

Say the Sentence Out Loud

Read the next line twice. Out loud, if the room allows.

A p value is the probability of observing data at least as extreme as the data you observed, if the null hypothesis were exactly true.

The null hypothesis is the boring version of the world: the drug does nothing, the two groups came from the same population, the coefficient is zero, the coin is fair. You do not believe it. Nobody really does. You assume it anyway, temporarily and in bad faith, because the calculation needs somewhere to stand. Now your data are a random draw from a known distribution, and you can ask where in that distribution your particular result landed. Far out in the tail? Then the boring world finds your data surprising. The p value says how far out, as an area under a curve.

This picture is the whole article. A curve showing what your result would look like if nothing whatever were going on, plus a shaded piece of it. The shaded piece is the p value. Four smaller copies follow, one for each thing people believe instead.

Now look at what the sentence leaves out. No probability that the null hypothesis is true. No size of any effect. No verdict on whether the result matters, and no promise that the lab down the corridor, running your study again next week, will land anywhere near where you did. The thing runs one way only: given the hypothesis, how weird are the data. Read it backwards and you have invented every mistake in this article.

We wrote this down because the club kept explaining the same four things at the table, month after month, usually to somebody who had already sat through it twice. Once, properly, on paper. Then we can point.

The Coffee Machine, and Why the Definition Is Backwards From What You Want

Suppose the club decides the coffee machine in the science corridor makes people faster at the Wednesday problem set. We grab 32 members. Half get coffee, half get hot water that looks like coffee. Everybody works the same twenty questions. Coffee group mean: 14.2 minutes. Water group: 15.9 minutes. Difference of 1.7 minutes, coffee wins, \(p = 0.031\).

What did we just learn?

Most readers hear this: coffee has a 3.1% chance of doing nothing. That would be a genuinely useful number to own. Nobody computed it. Wanting a number has never once summoned it. Here is what we actually computed. Assume coffee does nothing at all. Assume our sampling and our test are what we said they were. A gap of 1.7 minutes or larger then turns up in 3.1% of studies like ours, purely from which members landed in which group.

The two sentences sound alike in English, which is why people who would never confuse anything else trade one for the other. They relate the way "the probability of clouds given rain" relates to "the probability of rain given clouds." Rain all but guarantees clouds. Clouds guarantee you almost nothing. The swap has a name, the transposed conditional, and it fathers three of the four misreadings below [6, 7]. Learn the name. You will start seeing it everywhere.

Claim card 01the coffee machine
Observed
Coffee group 1.7 minutes faster, \(n = 16\) per group, \(p = 0.031\)
Means
If coffee did nothing, data this lopsided or worse would appear in 3.1% of runs
Does not mean
Coffee has a 3.1% chance of doing nothing. Coffee works, with 96.9% confidence. The effect measures 1.7 minutes. Anybody replicating us will also get \(p<0.05\). All four are wrong, and we take them in order below.

Look at the picture before you argue with the sentence. Figure 1 draws the definition for this exact experiment, and §1 stops sounding like a technicality the moment you have seen the area it names.

−4 −3 −2 −1 0 1 2 3 4 t: the gap between the two group means, in standard errors the world where coffee does nothing every result this null world produces, and how often t = 2.26 what we got t = −2.26 and as far the other way 0.0155 0.0155 shaded area, both tails = 0.031 that is the p value. It is an area, and nothing else. t distribution, 30 degrees of freedom (16 per group)
Figure 1. The whole definition in one picture, drawn for the coffee experiment. The curve is the distribution of results you would get if coffee did nothing at all: most of the time the two groups land close together, occasionally chance pushes them apart. Our observed gap sits at \(t = 2.26\), which with 16 members per group and a spread of about 2.1 minutes corresponds to the 1.7-minute difference we measured. The two shaded pieces are every outcome at least that lopsided, in either direction. Together they cover 0.031 of the total area. That area is the p value, and this shape, a distribution with a tail filled in, is the motif that returns at smaller sizes beside each misreading below.

The Property That Actually Defines It

Definitions given in words slide around under pressure. So here is the version you can put on a computer and check in ten seconds, which is the version the club trusts and goes back to whenever somebody starts hedging.

Under a true null hypothesis, p values are uniformly distributed between 0 and 1, meaning every value turns up as often as every other. A p of 0.02 is exactly as probable as a p of 0.72. And 5% of the time the number drops below 0.05, purely because random numbers behave like that. Stop calling it a side effect. Somebody specified it. The p value gets built by feeding your test statistic through the null distribution's own tail area, and that transformation flattens any distribution into a uniform one. Murdoch, Tsai and Adcock argue that simulation is how the thing should be taught in the first place, because p is a random variable, and students who never see its distribution never quite believe it [2].

We ran it. Two hundred thousand simulated two-group experiments, 32 observations per group, drawn from populations identical by construction, each one analysed with an ordinary two-sided t test. The code sits linked at the top of this page, seed included.

0 0.02 0.04 0.06 0 0.25 0.50 0.75 1.00 p value fraction of experiments per bin flat line at 0.05 this bar is "significant" 0.0493 of runs, and nothing is happening 200,000 simulated experiments, n = 32 per group true difference between groups: exactly zero chi-square vs flat: 21.99 on 19 df, p = 0.285
Figure 2. What a p value does when nothing is going on. Each bar is the fraction of 200,000 simulated experiments whose p value fell in that slice of width 0.05. The bars are flat, because that is the defining property. The dashed line marks the 0.05 each bar should sit at; the largest bar misses it by 2.1% of its own height. The dark bar on the left is every result that would be reported as statistically significant. There is no effect anywhere in this simulation. A Kolmogorov-Smirnov test against a perfect uniform distribution gives D = 0.00142, which on 200,000 draws is as close to flat as you have any right to expect.

Go back and look at the far-left bar of Figure 2. In a world holding no effects whatsoever, one experiment in twenty comes out with a publishable p value. The method has not failed you. That rate was agreed in advance, in writing, by everybody who ever wrote a methods section, as the price of having any rule at all. The rest of this article exists because of it.

Figure 3 runs the same code with a real effect underneath. Same test, same sample size. Only the true difference between the two populations changes.

0 0.2 0.4 0.6 0.8 d = 0 d = 0.2 d = 0.5 d = 0.8 no effect small medium large 5.1% 12.3% 50.4% 88.2% fraction of experiments Within each panel the horizontal axis runs p = 0 on the left to p = 1 on the right, in 20 bins. 100,000 experiments per panel, n = 32 per group. Dark bar = p below 0.05. Percentages are the club's simulated power; the exact noncentral-t values are 5.00, 12.35, 50.36 and 88.30.
Figure 3. The same experiment with a real effect underneath it. \(d\) is the true difference between the two populations in units of their standard deviation. At \(d = 0\) the distribution is flat, as Figure 2 showed. As the effect grows, mass slides leftward and piles against zero, and the dark leftmost bar, which is the fraction of studies that will be called significant, grows from 5.1% to 88.2%. That fraction has a name: power. Two things to take away. A real effect does not guarantee a small p value, and at \(d = 0.5\) it fails to deliver one about half the time. And the whole shape depends on sample size as much as on effect size, which is the subject of §5.

Wrong Thing One: The Probability the Hypothesis Is True

Start with the big one. This misreading will not die, because it answers the question every person running the study actually wanted answered.

You want \(P(\text{hypothesis} \mid \text{data})\). The p value hands you \(P(\text{data or worse} \mid \text{hypothesis})\). Different quantities. To get from one to the other you need something the p value has never carried and cannot carry, namely how plausible the hypothesis looked before you ran the study, because Bayes' rule takes a prior as an input. Go looking for a prior inside a p value. You will not find one. So no posterior is on offer.

How wrong does it get? Berger and Sellke worked out the best case, the case most flattering to whoever is doing the transposing [10]. Pick the prior that flatters the evidence against the null as far as mathematics allows. Do your worst. A p value of exactly 0.05 still supports odds against the null of no better than about 2.5 to 1. Not 19 to 1, which is what "95% confident" sounds like it should buy you. Sellke, Bayarri and Berger later reduced the whole thing to a formula you can do on a phone: the strongest possible Bayes factor against the null is \(-e\,p\ln p\) [11]. At \(p = 0.05\) that comes to 0.407. Odds of 2.5 to 1, about the evidence in a coin landing heads twice.

p = 0.05 → best-case odds against the null 2.5 to 1
p = 0.01 → 8.0 to 1
p = 0.005 → 13.9 to 1
p = 0.001 → 53.3 to 1

Those are ceilings, not estimates. Whatever your real evidence is, it sits below every line above. A p value just under 0.05, the number that publishes papers and funds laboratories, amounts to a shrug.

Wrong Thing Two: A Small P Means a Big Effect

Invent two studies. Study A tests a new tutoring scheme on 24 students and finds an improvement of 11 marks, \(p = 0.04\); study B tests a phone app on 40,000 students and finds an improvement of 0.3 marks, \(p = 0.0000002\). Which effect is bigger?

Obviously the tutoring. Eleven marks against a third of a mark. And yet study B's p value is smaller by five orders of magnitude, so skim the two abstracts at the speed everybody reads abstracts and you come away certain that B found something enormous and A something marginal.

A p value jams effect size, sample size and noise together. You cannot unjam it. Roughly, the test statistic scales as the effect divided by the standard error, and the standard error shrinks as \(1/\sqrt{n}\), so:

$$t \;\approx\; \frac{\text{effect}}{\text{noise}} \times \sqrt{\frac{n}{2}}$$

So you can manufacture any p value you like out of any effect you like, by choosing how many people to recruit. Our simulation makes this concrete. An effect of \(d = 0.5\) at 32 per group reaches significance 50.4% of the time; the same \(d = 0.5\) at 400 per group reaches it 100.0% of the time, rounded to four decimal places, across 100,000 runs. Identical effect. Completely different p values. Nothing about the world changed. Only the number of people we asked.

So when somebody reports \(p < 0.001\) and stops there, they have told you exactly one thing, and it is this: the combination of effect size and sample size came out large. Which of the two? Not saying. Ask for the effect size and its interval, every time, including of the papers whose conclusions you like. A paper that reports neither has still told you something. The omission is the finding [5].

Wrong Thing Three, as Plain Arithmetic: 0.049 Against 0.051

Take the two-group comparison from §3. Thirty-two observations per group, sixty-two degrees of freedom, a two-sided t test. Both studies are honest. One reports \(p = 0.049\), the other \(p = 0.051\). Here is what separates them.

p = 0.049 requires t = 2.0081
p = 0.051 requires t = 1.9900
difference in t: 0.0181

observed effect at p = 0.049: 0.5020 standard deviations
observed effect at p = 0.051: 0.4975 standard deviations
difference: 0.0045 standard deviations

standard error of that estimate: 0.2500 standard deviations
the gap, in standard errors: 0.018

Eighteen thousandths of one standard error separate them. Nothing else does. One gets written up as a discovery, submitted, reviewed, cited, taught. The other gets a folder nobody opens again.

Nothing in the mathematics changes as you cross the threshold. Somebody drew a line on a continuous scale. Rosenthal and Gaito asked nineteen psychologists how much confidence they placed in results at various p levels, and the drop came out steepest between 0.05 and 0.10. A cliff in belief, standing where the evidence has none [12]. Gelman and Stern put the general case plainly: the difference between significant and not significant is not itself statistically significant [13].

Two more consequences follow from the same arithmetic. A study at \(p = 0.06\) and a study at \(p = 0.04\) are not in conflict. And test each of two groups against zero, announce that one cleared the line and the other did not, and you have compared nothing at all, an error Nieuwenhuis and colleagues found in half the neuroscience papers where it could occur [14].

Wrong Thing Four: Not Significant Means No Effect

Absence of evidence is not evidence of absence. Altman and Bland put that sentence into the BMJ in 1995 [15]. People have quoted it ever since, then gone off and done the thing anyway.

A p value above your threshold means one of two things and cannot tell you which: either no effect exists, or one exists and your study was too small to see it. Look back at Figure 3. At \(d = 0.5\), a real and respectable effect, our simulation returns \(p > 0.05\) in 49.6% of runs; cut the sample to 16 per group, leave the effect sitting exactly where it was at 0.5, and now 72.3% of studies miss it altogether. Every one of those misses gets written up as "no significant difference was found."

The word "found" is doing heavy lifting there. Nothing was found. The study looked with a torch too dim for the job and reported darkness.

Button and colleagues put median power across neuroscience near 21% [16]. Read that number next to the paragraph above, slowly. A field running at that power learns almost nothing from a non-significant result, because even real effects get missed four times out of five, so ask what effect sizes the study could actually have detected. If the confidence interval stretches from "harmful" through "nothing" to "miraculous", the honest summary is that nobody knows yet.

Claim card 02the null result
Reported
"No significant effect of the intervention (\(p = 0.31\))"
Compatible with
No effect at all. A small effect. A large effect in a small study. All of these, simultaneously.
Ask for
The confidence interval. If it runs from −8% to +14%, the study has ruled out almost nothing and should say so.

Base Rates, and How an Honest Field Can Still Be Mostly Wrong

Now the part that changes how you read journals. It took the club longest to swallow. One of us went home and re-derived it alone before she would accept it.

Suppose a field runs 1,000 studies. Some fraction of the hypotheses tested are actually true; call that fraction \(R\), the prior, and note that a true hypothesis gets detected with probability equal to the power. A false one gets called significant with probability \(\alpha\), the threshold, since that is what \(\alpha\) means, so out of 1,000 studies:

$$\text{true positives} = 1000\,R\,\text{power}, \qquad \text{false positives} = 1000\,(1-R)\,\alpha$$

and the share of significant findings that are real, the positive predictive value, is the first over the sum of both. Ioannidis built his 2005 paper on exactly this arithmetic and gave it the title everybody remembers [3], and here is the club's version, with numbers picked to be unflattering rather than absurd: a field where 10% of tested hypotheses are true, running studies at 50% power.

1,000 studies
100 hypotheses true, 900 false

true and detected: 100 × 0.50 = 50
false and called anyway: 900 × 0.05 = 45

significant results: 50 + 45 = 95
of which real: 50 / 95 = 52.6%
of which false: 47.4%

Forty-seven percent of the published significant findings in that field are wrong. Nobody cheated. Every study used the right test, reported what it found, hit the alpha it promised and would survive an audit. The errors are structural. They come out of the base rate, and no single study can see its own base rate.

Algebra alone did not satisfy us, so we simulated it: 400,000 studies, each with a hypothesis drawn true with probability 0.10 and an effect size tuned to give exactly 50% power. Of the 37,996 that came out significant, 20,058 were real: simulated PPV 0.5279 against the analytic 0.5263, a difference of 0.0016, and Figure 4 turns the same arithmetic into a curve. Find your own field on it.

0 25% 50% 75% 100% 0.5% 1% 2% 5% 10% 20% 50% 100% prior: share of tested hypotheses that are actually true (log scale) share of significant findings that are real power 0.80 power 0.50 power 0.20 club headline: 52.6% prior 10%, power 50% screening a long list: 3.9% true coin flip α = 0.05 throughout
Figure 4. From the club's own calculation. The horizontal axis is the prior probability that a hypothesis pulled out of a field's research programme is true; note the log scale, which is where the action is. The vertical axis is the fraction of significant findings that are real. Three curves for three levels of power. The filled circle is our headline case: a field with a 10% prior running 50%-powered studies gets 52.6% of its significant claims right, which is barely better than a coin. The pale circle at the left is a screening programme with a 1% prior and 20% power, where 96.1% of significant findings are false. The dashed horizontal line is 50%. Everything below it is a field publishing more wrong claims than right ones.
Share of significant findings that are false, at α = 0.05. Rows are the prior; columns are power.
priorpower 0.20power 0.35power 0.50power 0.80power 0.95
0.0196.1%93.4%90.8%86.1%83.9%
0.0582.6%73.1%65.5%54.3%50.0%
0.1069.2%56.2%47.4%36.0%32.1%
0.2542.9%30.0%23.1%15.8%13.6%
0.5020.0%12.5%9.1%5.9%5.0%
0.805.9%3.4%2.4%1.5%1.3%

Two things jump out of that table. Power helps, and it helps less than you would hope: at a 1% prior, going from 20% power to 95% power takes you from 96.1% wrong to 83.9% wrong. The prior does the heavy lifting. A field testing well-motivated hypotheses grounded in earlier work stands somewhere entirely different from one screening ten thousand candidate genes, and the same p value means something different in each.

Which is why the proposal to move the default threshold to 0.005 has teeth [8]: hold our headline case fixed, change only alpha, and the false share drops from 47.4% to 8.3%. Lakens and eighty-seven co-authors object that one new default merely relocates the problem, and that researchers should justify their alpha against the specific costs of being wrong in their own field [9]. Both papers are worth your time. They agree entirely about the arithmetic.

Working Notes From the Club Table

Session 1 · somebody asks the obvious question

Nadia: "If two labs run the same study on the same real effect, how often do they agree?"

Nobody knew. Somebody guessed 90%. Onto the board it went, with the arithmetic under it. Ninety seconds, start to finish: give each study power \(w\), keep the two independent, and they land on opposite sides of the threshold with probability \(2w(1-w)\).

Plug in \(w = 0.5\). Two times a half times a half. One half. Two competent studies of one genuine effect, both run correctly by people who know exactly what they are doing, will contradict each other half the time, and neither has done a single thing wrong.

Four seconds of silence. Then everybody talked at once.

Session 1 · we refused to believe it and simulated it

200,000 pairs of studies, \(n = 32\) per group, true \(d = 0.5\). Simulated disagreement rate 0.5018. Analytic 0.5000. Correct, and always going to be. Cumming named this the dance of the p values [4]. Follow a first study at \(p = 0.05\) with a replication, and the second p value lands anywhere across a range roughly two orders of magnitude wide. Study size does not rescue you.

Our own version of his number, taken from the pairs where the first study came out significant: the replication's p value falls below 0.0004 five percent of the time and above 0.6791 five percent of the time. Same effect. Same design. Same code.

Session 2 · the gallery, which is the thing that finally landed

Rob said the tables were not working. He wanted twenty-four individual studies, printed out, one at a time. He was right. Three separate people who did not believe the 50% number have since been shown this figure.

1 0.1 0.01 0.001 0.0001 study number, in the order the simulation produced them p value (log scale) p = 0.05 24 honest replications of one real effect (d = 0.5, n = 32 per group) 14 significant, 10 not. Filled = would be published as a finding. study 5: p = 0.87 study 18: p = 0.0002
Figure 5. The dance. Twenty-four simulated studies, all measuring the same effect of exactly \(d = 0.5\), all with 32 observations per group, all analysed identically. The only difference between them is the random numbers. Fourteen cross the line and would be written up as discoveries; ten do not and would be written up as null results or never written up at all. Study 18 lands at \(p = 0.0002\) and study 5 at \(p = 0.87\), a spread of three and a half orders of magnitude with nothing underneath it but chance. Study 4 (\(p = 0.0392\)) and study 12 (\(p = 0.0577\)) sit on opposite sides of the threshold and are separated by less than two hundredths of a standard error.
Session 2 · what we decided the figure means

A single p value is one draw from a wide distribution. Reporting it to three decimal places claims a precision the quantity does not have. We round to two now. And we say what the interval was.

Why the Wrong Versions Are So Appealing

None of these misreadings are stupid. Each is a reasonable answer to a reasonable question, attached to the wrong number. Understanding why they feel right is most of what stops you making them again.

The first one answers the question you asked. You did not run the study to learn about hypothetical data in a world you disbelieve; you ran it to find out whether the thing is real. \(P(\text{hypothesis} \mid \text{data})\) is the quantity your brain ordered. So a number turns up looking like a probability, sitting right beside a conclusion, and your brain signs for the delivery without reading the label. Goodman calls the resulting muddle the p-value fallacy, and traces it to a historical mashing-together of two incompatible systems: Fisher's evidential p and the Neyman-Pearson error-rate framework. Two systems never meant to appear in the same sentence [7, 17].

Small numbers feel like big claims, and that keeps the second one alive. Scientific writing has trained everybody to read \(p < 0.001\) as emphasis. Three zeros look like shouting. A grain of truth sits underneath, because at a fixed sample size a smaller p really does mean a bigger observed effect, and the grain of truth is exactly what makes the habit dangerous. Sample sizes are never fixed across the papers you are comparing.

Thresholds survive because decisions genuinely work that way. Somebody has to decide whether to fund the trial, a decision rule needs a line, and Neyman and Pearson designed their framework for exactly that, a procedure with known long-run error rates, inside which a threshold is the right tool [17]. The error comes from carrying the decision rule back into the evidence. The line tells you what to do next; what is true stays exactly as unsettled as it was before anybody drew it, and 0.051 remains very nearly identical evidence to 0.049.

Our language has no word for "we looked and could not tell". The fourth misreading fills the gap. Every report needs a finding. "No significant difference" sounds like a finding, reads like a finding, and gets summarised as a finding, where the accurate version says "our study was too small to distinguish a 30% benefit from a 10% harm". Terrible headline. Almost nobody writes it.

Every one of these is the mind doing something sensible with a number that was not built for it.

A fifth misreading deserves a mention, quieter and worse. People treat the p value as a property of the effect rather than of the study. The effect does not have a p value; the study has one, and a different study of the same effect would hand you a different one, as Figure 5 shows twenty-four times over.

The Strongest Case Against Everything Above

We should take the other side seriously, because the other side is not silly and includes people who have thought about this far longer than we have.

The argument runs like this. Yes, p values are misinterpreted. So is every statistic ever invented, including the ones proposed as replacements. Bayes factors depend on a prior over effect sizes, usually chosen by convention and rarely defended, and a reader who cannot interpret a p value will do no better with a quantity carrying an extra arbitrary input. Confidence intervals get misread as containing the true value with 95% probability, which is the same transposed conditional in a different coat, so the misinterpretation plainly lives in the reader. Move the reader to a new tool and the misinterpretation moves too.

More pointedly: the p value does one job honestly, and it is a job nothing else does as cheaply. It calibrates a test against a reference distribution and gives a number whose long-run behaviour under the null is known exactly, without requiring you to specify what you believe beforehand. That property is what makes regulatory decisions possible between parties who disagree about priors. A drug regulator and a drug company do not share a prior. They can share an alpha.

Now our own weakest joint. The base-rate argument in §8 depends on a number, \(R\), that nobody can measure. We chose 10% because it makes a good demonstration. A defender can reasonably answer that a well-motivated field, working from mechanism and earlier evidence, might run at \(R = 0.5\) or higher. Look that row up in the table. Only 9.1% of significant findings are wrong, which is roughly what a reasonable person expects anyway. Our headline number is a choice about \(R\) dressed up as a result. We should say so plainly. We just did.

The last objection has teeth. Reproducibility numbers are not clean evidence for the story we are telling. The Open Science Collaboration replicated 100 psychology studies and got 36% significant. Widespread false positives would produce that [18]. So would real but smaller effects. So would replications differing from the originals in ways that mattered. So would the plain power arithmetic in §9, under which honest replications of real effects fail routinely.

Where we come out. The strongest version of the defence is correct about the tool and wrong about the practice. Nothing in this article says the p value is a bad statistic. It says the p value is a narrow one, asked to carry a decision, an effect size, a measure of evidence and a claim about replication. Nobody built it for any of that. The ASA's 2016 statement lands in the same place, and it had to be negotiated among statisticians who agreed about almost nothing else. P values are useful. They do not measure the probability that a hypothesis is true, nor the size of an effect, nor the importance of a result, nor on its own how much evidence a model deserves [5]. We will keep using them. We will keep reading the confidence interval first.

What We Do Now, Concretely

Rules the club adopted for itself. We apply them to other people's papers too, which has made us harder to sit next to.

Report the estimate and its interval before the p value, in the same sentence. "The coffee group finished 1.7 minutes faster (95% CI 0.2 to 3.2, \(p = 0.031\))." Now the reader can see it. The data sit comfortably with an effect anywhere from trivial to substantial. Say so.

Never write "no effect" for a non-significant result. Write what the study could have detected. If the interval includes effects you would have cared about, the study did not settle the question.

If you do not know the study’s power, say so. Power depends on the effect size you were hoping to find, so state that too, because a study designed to detect \(d = 0.8\) has almost nothing to say about \(d = 0.2\).

The number of tests you ran belongs in the paper. Twenty independent tests under a true null give you a 64% chance of at least one result below 0.05, and the error rate that matters runs across every test you performed. Not the one you chose to write up [6].

Treat the threshold as a filing decision, not a conclusion. Amrhein, Greenland and McShane, with 800 co-signatories, argued for retiring the phrase "statistically significant" entirely [19]. We have not gone that far. We have stopped writing the word "significant" without a number next to it.

One last thing, and keep this one if you keep nothing else: the p value answers a question about data in a world where the effect is absent. Any sentence beginning "therefore the hypothesis is probably" has left that world without telling you. Go back and check what you actually computed.

The area, and only the area. Everything in this article follows from the curve being the null world and the shading being a tail of it. Hold that picture and you catch all four mistakes as they go past.

References

  1. Fisher, R. A. (1925). Statistical Methods for Research Workers. Oliver & Boyd, Edinburgh.
  2. Murdoch, D. J., Tsai, Y.-L. & Adcock, J. (2008). P-values are random variables. The American Statistician 62, 242–245. doi:10.1198/000313008X332421
  3. Ioannidis, J. P. A. (2005). Why most published research findings are false. PLoS Medicine 2, e124. doi:10.1371/journal.pmed.0020124
  4. Cumming, G. (2008). Replication and p intervals: p values predict the future only vaguely, but confidence intervals do much better. Perspectives on Psychological Science 3, 286–300. doi:10.1111/j.1745-6924.2008.00079.x
  5. Wasserstein, R. L. & Lazar, N. A. (2016). The ASA statement on p-values: context, process, and purpose. The American Statistician 70, 129–133. doi:10.1080/00031305.2016.1154108
  6. 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, 337–350. doi:10.1007/s10654-016-0149-3
  7. Goodman, S. N. (2008). A dirty dozen: twelve p-value misconceptions. Seminars in Hematology 45, 135–140. doi:10.1053/j.seminhematol.2008.04.003
  8. Benjamin, D. J., Berger, J. O., Johannesson, M. et al. (2018). Redefine statistical significance. Nature Human Behaviour 2, 6–10. doi:10.1038/s41562-017-0189-z
  9. Lakens, D., Adolfi, F. G., Albers, C. J. et al. (2018). Justify your alpha. Nature Human Behaviour 2, 168–171. doi:10.1038/s41562-018-0311-x
  10. Berger, J. O. & Sellke, T. (1987). Testing a point null hypothesis: the irreconcilability of p values and evidence. Journal of the American Statistical Association 82, 112–122. doi:10.1080/01621459.1987.10478397
  11. Sellke, T., Bayarri, M. J. & Berger, J. O. (2001). Calibration of p values for testing precise null hypotheses. The American Statistician 55, 62–71. doi:10.1198/000313001300339950
  12. Rosenthal, R. & Gaito, J. (1963). The interpretation of levels of significance by psychological researchers. The Journal of Psychology 55, 33–38. doi:10.1080/00223980.1963.9916596
  13. Gelman, A. & Stern, H. (2006). The difference between "significant" and "not significant" is not itself statistically significant. The American Statistician 60, 328–331. doi:10.1198/000313006X152649
  14. Nieuwenhuis, S., Forstmann, B. U. & Wagenmakers, E.-J. (2011). Erroneous analyses of interactions in neuroscience: a problem of significance. Nature Neuroscience 14, 1105–1107. doi:10.1038/nn.2886
  15. Altman, D. G. & Bland, J. M. (1995). Statistics notes: absence of evidence is not evidence of absence. BMJ 311, 485. doi:10.1136/bmj.311.7003.485
  16. Button, K. S., Ioannidis, J. P. A., Mokrysz, C., Nosek, B. A., Flint, J., Robinson, E. S. J. & Munafò, M. R. (2013). Power failure: why small sample size undermines the reliability of neuroscience. Nature Reviews Neuroscience 14, 365–376. doi:10.1038/nrn3475
  17. 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. doi:10.1098/rsta.1933.0009
  18. Open Science Collaboration (2015). Estimating the reproducibility of psychological science. Science 349, aac4716. doi:10.1126/science.aac4716
  19. Amrhein, V., Greenland, S. & McShane, B. (2019). Scientists rise up against statistical significance. Nature 567, 305–307. doi:10.1038/d41586-019-00857-9