p-value

From HandWiki
Short description: Function of the observed sample results

In null-hypothesis significance testing, the p-value[note 1] is the probability of obtaining test results at least as extreme as the result actually observed, under the assumption that the null hypothesis is correct.[2][3] A very small p-value means that such an extreme observed outcome would be very unlikely under the null hypothesis. Even though reporting p-values of statistical tests is common practice in academic publications of many quantitative fields, misinterpretation and misuse of p-values is widespread and has been a major topic in mathematics and metascience.[4][5] In 2016, the American Statistical Association (ASA) made a formal statement that "p-values do not measure the probability that the studied hypothesis is true, or the probability that the data were produced by random chance alone" and that "a p-value, or statistical significance, does not measure the size of an effect or the importance of a result" or "evidence regarding a model or hypothesis."[6] That said, a 2019 task force by ASA has issued a statement on statistical significance and replicability, concluding with: "p-values and significance tests, when properly applied and interpreted, increase the rigor of the conclusions drawn from data."[7]

Basic concepts

In statistics, every conjecture concerning the unknown probability distribution of a collection of random variables representing the observed data [math]\displaystyle{ X }[/math] in some study is called a statistical hypothesis. If we state one hypothesis only and the aim of the statistical test is to see whether this hypothesis is tenable, but not to investigate other specific hypotheses, then such a test is called a null hypothesis test.

As our statistical hypothesis will, by definition, state some property of the distribution, the null hypothesis is the default hypothesis under which that property does not exist. The null hypothesis is typically that some parameter (such as a correlation or a difference between means) in the populations of interest is zero. Our hypothesis might specify the probability distribution of [math]\displaystyle{ X }[/math] precisely, or it might only specify that it belongs to some class of distributions. Often, we reduce the data to a single numerical statistic, e.g., [math]\displaystyle{ T }[/math], whose marginal probability distribution is closely connected to a main question of interest in the study.

The p-value is used in the context of null hypothesis testing in order to quantify the statistical significance of a result, the result being the observed value of the chosen statistic [math]\displaystyle{ T }[/math].[note 2] The lower the p-value is, the lower the probability of getting that result if the null hypothesis were true. A result is said to be statistically significant if it allows us to reject the null hypothesis. All other things being equal, smaller p-values are taken as stronger evidence against the null hypothesis.

Loosely speaking, rejection of the null hypothesis implies that there is sufficient evidence against it.

As a particular example, if a null hypothesis states that a certain summary statistic [math]\displaystyle{ T }[/math] follows the standard normal distribution N(0,1), then the rejection of this null hypothesis could mean that (i) the mean of [math]\displaystyle{ T }[/math] is not 0, or (ii) the variance of [math]\displaystyle{ T }[/math] is not 1, or (iii) [math]\displaystyle{ T }[/math] is not normally distributed. Different tests of the same null hypothesis would be more or less sensitive to different alternatives. However, even if we do manage to reject the null hypothesis for all 3 alternatives, and even if we know the distribution is normal and variance is 1, the null hypothesis test does not tell us which non-zero values of the mean are now most plausible. The more independent observations from the same probability distribution one has, the more accurate the test will be, and the higher the precision with which one will be able to determine the mean value and show that it is not equal to zero; but this will also increase the importance of evaluating the real-world or scientific relevance of this deviation.

Definition and interpretation

Definition

The p-value is the probability under the null hypothesis of obtaining a real-valued test statistic at least as extreme as the one obtained. Consider an observed test-statistic [math]\displaystyle{ t }[/math] from unknown distribution [math]\displaystyle{ T }[/math]. Then the p-value [math]\displaystyle{ p }[/math] is what the prior probability would be of observing a test-statistic value at least as "extreme" as [math]\displaystyle{ t }[/math] if null hypothesis [math]\displaystyle{ H_0 }[/math] were true. That is:

  • [math]\displaystyle{ p = \Pr(T \geq t \mid H_0) }[/math] for a one-sided right-tail test-statistic distribution,
  • [math]\displaystyle{ p = \Pr(T \leq t \mid H_0) }[/math] for a one-sided left-tail test-statistic distribution,
  • [math]\displaystyle{ p = 2\min\{\Pr(T \geq t \mid H_0),\Pr(T \leq t \mid H_0)\} }[/math] for a two-sided test-statistic distribution. If the distribution of [math]\displaystyle{ T }[/math] is symmetric about zero, then [math]\displaystyle{ p =\Pr(|T| \geq |t| \mid H_0) }[/math]

Interpretations

The error that a practising statistician would consider the more important to avoid (which is a subjective judgment) is called the error of the first kind. The first demand of the mathematical theory is to deduce such test criteria as would ensure that the probability of committing an error of the first kind would equal (or approximately equal, or not exceed) a preassigned number α, such as α = 0.05 or 0.01, etc. This number is called the level of significance.

In a significance test, the null hypothesis [math]\displaystyle{ H_0 }[/math] is rejected if the p-value is less than or equal to a predefined threshold value [math]\displaystyle{ \alpha }[/math], which is referred to as the alpha level or significance level. [math]\displaystyle{ \alpha }[/math] is not derived from the data, but rather is set by the researcher before examining the data. [math]\displaystyle{ \alpha }[/math] is commonly set to 0.05, though lower alpha levels are sometimes used. In 2018, a group of statisticians led by Daniel Benjamin proposed the adoption of the 0.005 value as standard value for statistical significance worldwide.[8]

Different p-values based on independent sets of data can be combined, for instance using Fisher's combined probability test.

Distribution

The p-value is a function of the chosen test statistic [math]\displaystyle{ T }[/math] and is therefore a random variable. If the null hypothesis fixes the probability distribution of [math]\displaystyle{ T }[/math] precisely (e.g. [math]\displaystyle{ H_0: \theta = \theta_0 }[/math] where [math]\displaystyle{ \theta }[/math] is the only parameter), and if that distribution is continuous, then when the null-hypothesis is true the p-value is uniformly distributed between 0 and 1. Regardless of the truth of the [math]\displaystyle{ H_0 }[/math], the p-value is not fixed; if the same test is repeated independently with fresh data one will typically obtain a different p-value in each iteration.

Usually only a single p-value relating to a hypothesis is observed, so the p-value is interpreted by a significance test and no effort is made to estimate the distribution it was drawn from. When a collection of p-values are available (e.g. when considering a group of studies on the same subject), the distribution of p-values is sometimes called a p-curve.[9] A p-curve can be used to assess the reliability of scientific literature, such as by detecting publication bias or p-hacking. [9][10]

Distribution for composite hypothesis

In parametric hypothesis testing problems, a simple or point hypothesis refers to a hypothesis where the parameter's value is assumed to be a single number. In contrast, in a composite hypothesis the parameter's value is given by a set of numbers. When the null-hypothesis is composite (or the distribution of the statistic is discrete), then when the null-hypothesis is true the probability of obtaining a p-value less than or equal to any number between 0 and 1 is still less than or equal to that number. In other words, it remains the case that very small values are relatively unlikely if the null-hypothesis is true, and that a significance test at level [math]\displaystyle{ \alpha }[/math] is obtained by rejecting the null-hypothesis if the p-value is less than or equal to [math]\displaystyle{ \alpha }[/math].[11][12]

For example, when testing the null hypothesis that a distribution is normal with a mean less than or equal to zero against the alternative that the mean is greater than zero ([math]\displaystyle{ H_0: \mu \leq 0 }[/math], variance known), the null hypothesis does not specify the exact probability distribution of the appropriate test statistic. In this example that would be the Z-statistic belonging to the one-sided one-sample Z-test. For each possible value of the theoretical mean, the Z-test statistic has a different probability distribution. In these circumstances the p-value is defined by taking the least favorable null-hypothesis case, which is typically on the border between null and alternative. This definition ensures the complementarity of p-values and alpha-levels: [math]\displaystyle{ \alpha = 0.05 }[/math] means one only rejects the null hypothesis if the p-value is less than or equal to [math]\displaystyle{ 0.05 }[/math], and the hypothesis test will indeed have a maximum type-1 error rate of [math]\displaystyle{ 0.05 }[/math].

Usage

The p-value is widely used in statistical hypothesis testing, specifically in null hypothesis significance testing. In this method, before conducting the study, one first chooses a model (the null hypothesis) and the alpha level α (most commonly 0.05). After analyzing the data, if the p-value is less than α, that is taken to mean that the observed data is sufficiently inconsistent with the null hypothesis for the null hypothesis to be rejected. However, that does not prove that the null hypothesis is false. The p-value does not, in itself, establish probabilities of hypotheses. Rather, it is a tool for deciding whether to reject the null hypothesis.[13]

Misuse

Main page: Misuse of p-values

According to the ASA, there is widespread agreement that p-values are often misused and misinterpreted.[3] One practice that has been particularly criticized is accepting the alternative hypothesis for any p-value nominally less than .05 without other supporting evidence. Although p-values are helpful in assessing how incompatible the data are with a specified statistical model, contextual factors must also be considered, such as "the design of a study, the quality of the measurements, the external evidence for the phenomenon under study, and the validity of assumptions that underlie the data analysis".[3] Another concern is that the p-value is often misunderstood as being the probability that the null hypothesis is true.[3][14]

Some statisticians have proposed abandoning p-values and focusing more on other inferential statistics,[3] such as confidence intervals,[15][16] likelihood ratios,[17][18] or Bayes factors,[19][20][21] but there is heated debate on the feasibility of these alternatives.[22]Cite error: Closing </ref> missing for <ref> tag John Arbuthnot studied this question in 1710,[23][24][25][26] and examined birth records in London for each of the 82 years from 1629 to 1710. In every year, the number of males born in London exceeded the number of females. Considering more male or more female births as equally likely, the probability of the observed outcome is 1/282, or about 1 in 4,836,000,000,000,000,000,000,000; in modern terms, the p-value. This is vanishingly small, leading Arbuthnot that this was not due to chance, but to divine providence: "From whence it follows, that it is Art, not Chance, that governs." In modern terms, he rejected the null hypothesis of equally likely male and female births at the p = 1/282 significance level. This and other work by Arbuthnot is credited as "… the first use of significance tests …"[27] the first example of reasoning about statistical significance,[28] and "… perhaps the first published report of a nonparametric test …",[24] specifically the sign test; see details at Sign test § History.

The same question was later addressed by Pierre-Simon Laplace, who instead used a parametric test, modeling the number of male births with a binomial distribution:[29]

In the 1770s Laplace considered the statistics of almost half a million births. The statistics showed an excess of boys compared to girls. He concluded by calculation of a p-value that the excess was a real, but unexplained, effect.

The p-value was first formally introduced by Karl Pearson, in his Pearson's chi-squared test,[30] using the chi-squared distribution and notated as capital P.[30] The p-values for the chi-squared distribution (for various values of χ2 and degrees of freedom), now notated as P, were calculated in (Elderton 1902), collected in (Pearson 1914).

The use of the p-value in statistics was popularized by Ronald Fisher,[31][full citation needed] and it plays a central role in his approach to the subject.[32] In his influential book Statistical Methods for Research Workers (1925), Fisher proposed the level p = 0.05, or a 1 in 20 chance of being exceeded by chance, as a limit for statistical significance, and applied this to a normal distribution (as a two-tailed test), thus yielding the rule of two standard deviations (on a normal distribution) for statistical significance (see 68–95–99.7 rule).[33][note 3][34]

He then computed a table of values, similar to Elderton but, importantly, reversed the roles of χ2 and p. That is, rather than computing p for different values of χ2 (and degrees of freedom n), he computed values of χ2 that yield specified p-values, specifically 0.99, 0.98, 0.95, 0,90, 0.80, 0.70, 0.50, 0.30, 0.20, 0.10, 0.05, 0.02, and 0.01.[35] That allowed computed values of χ2 to be compared against cutoffs and encouraged the use of p-values (especially 0.05, 0.02, and 0.01) as cutoffs, instead of computing and reporting p-values themselves. The same type of tables were then compiled in (Fisher Yates), which cemented the approach.[34]

As an illustration of the application of p-values to the design and interpretation of experiments, in his following book The Design of Experiments (1935), Fisher presented the lady tasting tea experiment,[36] which is the archetypal example of the p-value.

To evaluate a lady's claim that she (Muriel Bristol) could distinguish by taste how tea is prepared (first adding the milk to the cup, then the tea, or first tea, then milk), she was sequentially presented with 8 cups: 4 prepared one way, 4 prepared the other, and asked to determine the preparation of each cup (knowing that there were 4 of each). In that case, the null hypothesis was that she had no special ability, the test was Fisher's exact test, and the p-value was [math]\displaystyle{ 1/\binom{8}{4} = 1/70 \approx 0.014, }[/math] so Fisher was willing to reject the null hypothesis (consider the outcome highly unlikely to be due to chance) if all were classified correctly. (In the actual experiment, Bristol correctly classified all 8 cups.)

Fisher reiterated the p = 0.05 threshold and explained its rationale, stating:[37]

It is usual and convenient for experimenters to take 5 per cent as a standard level of significance, in the sense that they are prepared to ignore all results which fail to reach this standard, and, by this means, to eliminate from further discussion the greater part of the fluctuations which chance causes have introduced into their experimental results.

He also applies this threshold to the design of experiments, noting that had only 6 cups been presented (3 of each), a perfect classification would have only yielded a p-value of [math]\displaystyle{ 1/\binom{6}{3} = 1/20 = 0.05, }[/math] which would not have met this level of significance.[37] Fisher also underlined the interpretation of p, as the long-run proportion of values at least as extreme as the data, assuming the null hypothesis is true.

In later editions, Fisher explicitly contrasted the use of the p-value for statistical inference in science with the Neyman–Pearson method, which he terms "Acceptance Procedures".[38] Fisher emphasizes that while fixed levels such as 5%, 2%, and 1% are convenient, the exact p-value can be used, and the strength of evidence can and will be revised with further experimentation. In contrast, decision procedures require a clear-cut decision, yielding an irreversible action, and the procedure is based on costs of error, which, he argues, are inapplicable to scientific research.

Related indices

The E-value can refer to two concepts, both of which are related to the p-value and both of which play a role in multiple testing. First, it corresponds to a generic, more robust alternative to the p-value that can deal with optional continuation of experiments. Second, it is also used to abbreviate "expect value", which is the expected number of times that one expects to obtain a test statistic at least as extreme as the one that was actually observed if one assumes that the null hypothesis is true.[39] This expect-value is the product of the number of tests and the p-value.

The q-value is the analog of the p-value with respect to the positive false discovery rate.[40] It is used in multiple hypothesis testing to maintain statistical power while minimizing the false positive rate.[41]

The Probability of Direction (pd) is the Bayesian numerical equivalent of the p-value.[42] It corresponds to the proportion of the posterior distribution that is of the median's sign, typically varying between 50% and 100%, and representing the certainty with which an effect is positive or negative.

Second-generation p-values extend the concept of p-values by not considering extremely small, practically irrelevant effect sizes as significant.[43]

See also

Notes

  1. Italicisation, capitalisation and hyphenation of the term vary. For example, AMA style uses "P value", APA style uses "p value", and the American Statistical Association uses "p-value". In all cases, the "p" stands for probability.[1]
  2. The statistical significance of a result does not imply that the result also has real-world relevance. For instance, a medicine might have a statistically significant effect that is too small to be interesting.
  3. To be more specific, the p = 0.05 corresponds to about 1.96 standard deviations for a normal distribution (two-tailed test), and 2 standard deviations corresponds to about a 1 in 22 chance of being exceeded by chance, or p ≈ 0.045; Fisher notes these approximations.

References

  1. "ASA House Style". Amstat News. American Statistical Association. http://magazine.amstat.org/wp-content/uploads/STATTKadmin/style%5B1%5D.pdf. 
  2. "Not Even Scientists Can Easily Explain P-values". 2015-11-24. https://fivethirtyeight.com/features/not-even-scientists-can-easily-explain-p-values/. 
  3. 3.0 3.1 3.2 3.3 3.4 "The ASA's Statement on p-Values: Context, Process, and Purpose". The American Statistician 70 (2): 129–133. 7 March 2016. doi:10.1080/00031305.2016.1154108. 
  4. "Why P Values Are Not a Useful Measure of Evidence in Statistical Significance Testing". Theory & Psychology 18 (1): 69–88. 2008. doi:10.1177/0959354307086923. 
  5. "A manifesto for reproducible science". Nature Human Behaviour 1: 0021. January 2017. doi:10.1038/s41562-016-0021. PMID 33954258. 
  6. Wasserstein, Ronald L.; Lazar, Nicole A. (2016-04-02). "The ASA Statement on p -Values: Context, Process, and Purpose" (in en). The American Statistician 70 (2): 129–133. doi:10.1080/00031305.2016.1154108. ISSN 0003-1305. 
  7. Benjamini, Yoav; De Veaux, Richard D.; Efron, Bradley; Evans, Scott; Glickman, Mark; Graubard, Barry I.; He, Xuming; Meng, Xiao-Li et al. (2021-10-02). "ASA President's Task Force Statement on Statistical Significance and Replicability". Chance (Informa UK Limited) 34 (4): 10–11. doi:10.1080/09332480.2021.2003631. ISSN 0933-2480. 
  8. Benjamin, Daniel J.; Berger, James O.; Johannesson, Magnus; Nosek, Brian A.; Wagenmakers, E.-J.; Berk, Richard; Bollen, Kenneth A.; Brembs, Björn et al. (1 September 2017). "Redefine statistical significance". Nature Human Behaviour 2 (1): 6–10. doi:10.1038/s41562-017-0189-z. PMID 30980045. http://doi.org/10.1038/s41562-017-0189-z. 
  9. 9.0 9.1 "The extent and consequences of p-hacking in science". PLOS Biology 13 (3): e1002106. March 2015. doi:10.1371/journal.pbio.1002106. PMID 25768323. 
  10. "p-Curve and Effect Size: Correcting for Publication Bias Using Only Significant Results". Perspectives on Psychological Science 9 (6): 666–681. November 2014. doi:10.1177/1745691614553988. PMID 26186117. 
  11. "Median of the p value under the alternative hypothesis". The American Statistician 56 (3): 202–6. 2002. doi:10.1198/000313002146. 
  12. "The behavior of the P-value when the alternative hypothesis is true". Biometrics 53 (1): 11–22. March 1997. doi:10.2307/2533093. PMID 9147587. https://zenodo.org/record/1235121. 
  13. "Scientific method: statistical errors". Nature 506 (7487): 150–152. February 2014. doi:10.1038/506150a. PMID 24522584. Bibcode2014Natur.506..150N. 
  14. "An investigation of the false discovery rate and the misinterpretation of p-values". Royal Society Open Science 1 (3): 140216. November 2014. doi:10.1098/rsos.140216. PMID 26064558. Bibcode2014RSOS....140216C. 
  15. "Alternatives to P value: confidence interval and effect size". Korean Journal of Anesthesiology 69 (6): 555–562. December 2016. doi:10.4097/kjae.2016.69.6.555. PMID 27924194. 
  16. "Why the P-value culture is bad and confidence intervals a better alternative". Osteoarthritis and Cartilage 20 (8): 805–808. August 2012. doi:10.1016/j.joca.2012.04.001. PMID 22503814. 
  17. "Sifting the evidence. Likelihood ratios are alternatives to P values". BMJ 322 (7295): 1184–1185. May 2001. doi:10.1136/bmj.322.7295.1184. PMID 11379590. 
  18. "The Likelihood Paradigm for Statistical Evidence" (in en). The Nature of Scientific Evidence. 2004. pp. 119–152. doi:10.7208/chicago/9780226789583.003.0005. ISBN 9780226789576. 
  19. "Replacing p-values with Bayes-Factors: A Miracle Cure for the Replicability Crisis in Psychological Science". 30 April 2015. https://replicationindex.wordpress.com/2015/04/30/replacing-p-values-with-bayes-factors-a-miracle-cure-for-the-replicability-crisis-in-psychological-science/. 
  20. "Hypothesis Testing: From p Values to Bayes Factors". Journal of the American Statistical Association 95 (452): 1316–1320. December 2000. doi:10.2307/2669779. 
  21. "A Test by Any Other Name: P Values, Bayes Factors, and Statistical Inference". Multivariate Behavioral Research 51 (1): 23–29. 16 February 2016. doi:10.1080/00273171.2015.1099032. PMID 26881954. 
  22. "In defense of P values". Ecology 95 (3): 611–617. March 2014. doi:10.1890/13-0590.1. PMID 24804441. Bibcode2014Ecol...95..611M. https://zenodo.org/record/894459. 
  23. "An argument for Divine Providence, taken from the constant regularity observed in the births of both sexes". Philosophical Transactions of the Royal Society of London 27 (325–336): 186–190. 1710. doi:10.1098/rstl.1710.0011. http://www.york.ac.uk/depts/maths/histstat/arbuthnot.pdf. 
  24. 24.0 24.1 "Chapter 3.4: The Sign Test". Practical Nonparametric Statistics (Third ed.). Wiley. 1999. pp. 157–176. ISBN 978-0-471-16068-7. 
  25. Applied Nonparametric Statistical Methods (Second ed.). Chapman & Hall. 1989. ISBN 978-0-412-44980-2. 
  26. The History of Statistics: The Measurement of Uncertainty Before 1900. Harvard University Press. 1986. pp. 225–226. ISBN 978-0-67440341-3. 
  27. "John Arbuthnot". Statisticians of the Centuries. Springer. 2001. pp. 39–42. ISBN 978-0-387-95329-8. 
  28. "Chapter 4. Chance or Design: Tests of Significance". A History of Mathematical Statistics from 1750 to 1930. Wiley. 1998. pp. 65. 
  29. The History of Statistics: The Measurement of Uncertainty Before 1900. Harvard University Press. 1986. p. 134. ISBN 978-0-67440341-3. 
  30. 30.0 30.1 "On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling". Philosophical Magazine. Series 5 50 (302): 157–175. 1900. doi:10.1080/14786440009463897. http://www.economics.soton.ac.uk/staff/aldrich/1900.pdf. 
  31. Inman 2004.
  32. "Confusion Over Measures of Evidence (p′s) Versus Errors (α′s) in Classical Statistical Testing", The American Statistician 57 (3): 171–178 [p. 171], 2003, doi:10.1198/0003130031856 
  33. Fisher 1925, p. 47, Chapter III. Distributions.
  34. 34.0 34.1 Dallal 2012, Note 31: Why P=0.05?.
  35. Fisher 1925, pp. 78–79, 98, Chapter IV. Tests of Goodness of Fit, Independence and Homogeneity; with Table of χ2, Table III. Table of χ2.
  36. Fisher 1971, II. The Principles of Experimentation, Illustrated by a Psycho-physical Experiment.
  37. 37.0 37.1 Fisher 1971, Section 7. The Test of Significance.
  38. Fisher 1971, Section 12.1 Scientific Inference and Acceptance Procedures.
  39. "Definition of E-value". National Institutes of Health. https://www.ncbi.nlm.nih.gov/blast/Blast.cgi?CMD=Web&PAGE_TYPE=BlastDocs&DOC_TYPE=FAQ#expect. 
  40. "The positive false discovery rate: a Bayesian interpretation and the q-value". The Annals of Statistics 31 (6): 2013–2035. 2003. doi:10.1214/aos/1074290335. 
  41. "Statistical significance for genomewide studies". Proceedings of the National Academy of Sciences of the United States of America 100 (16): 9440–9445. August 2003. doi:10.1073/pnas.1530509100. PMID 12883005. Bibcode2003PNAS..100.9440S. 
  42. "Indices of Effect Existence and Significance in the Bayesian Framework". Frontiers in Psychology 10: 2767. 10 December 2019. doi:10.3389/fpsyg.2019.02767. PMID 31920819. 
  43. An Introduction to Second-Generation p-Values Jeffrey D. Blume,Robert A. Greevy,Valerie F. Welty,Jeffrey R. Smith &William D. Dupont https://www.tandfonline.com/doi/full/10.1080/00031305.2018.1537893

Further reading

External links

ja:有意#p値