Frequent Links
Bayes' theorem
In probability theory and statistics, Bayes' theorem (alternatively Bayes' law or Bayes' rule) describes the probability of an event, based on conditions that might be related to the event. For example, suppose one is interested in whether Addison has cancer. Furthermore, suppose that Addison is age 65. If cancer is related to age, information about Addison's age can be used to more accurately assess his or her chance of having cancer using Bayes' Theorem.
When applied, the probabilities involved in Bayes' theorem may have different interpretations. In one of these interpretations, the theorem is used directly as part of a particular approach to statistical inference. In particular, with the Bayesian interpretation of probability, the theorem expresses how a subjective degree of belief should rationally change to account for evidence: this is Bayesian inference, which is fundamental to Bayesian statistics. However, Bayes' theorem has applications in a wide range of calculations involving probabilities, not just in Bayesian inference.
Bayes' theorem is named after Rev. Thomas Bayes (/ˈbeɪz/; 1701–1761), who first^{[citation needed]} showed how to use new evidence to update beliefs. It was further developed by Pierre-Simon Laplace, who first published the modern formulation in his 1812 Théorie analytique des probabilités. Sir Harold Jeffreys put Bayes' algorithm and Laplace's formulation on an axiomatic basis. Jeffreys wrote that Bayes' theorem "is to the theory of probability what Pythagoras's theorem is to geometry".^{[1]}
Contents
Statement of theorem
Bayes' theorem is stated mathematically as the following equation:^{[2]}
- <math>P(A|B) = \frac{P(B | A)\, P(A)}{P(B)},</math>
where A and B are events.
- P(A) and P(B) are the probabilities of A and B without regard to each other.
- P(A | B), a conditional probability, is the probability of A given that B is true.
- P(B | A), is the probability of B given that A is true.
Examples
Cancer at age 65
Suppose we want to know Addison's probability of having cancer, but we don't know anything about him or her (apart from name, which we assume to be unrelated to having cancer). Despite not knowing anything about Addison, a probability can be assigned based on the general prevalence of cancer. For the sake of this example, suppose it is 1%. This is known as the base rate or prior probability of having cancer. "Prior" refers to the time before being informed about the particular case at hand.
Next, suppose we find out that Addison is 65 years old. If we assume that cancer and age are related, this new piece of information can be used to better assess Addison's chance of having cancer. More precisely, we'd like to know the probability that a person has cancer when it is known that he or she is 65 years old. This quantity is known as the current probability, where "current" refers to upon finding out information about the particular case at hand.
In order to apply knowledge of Addison's age in conjunction with Bayes' Theorem, two additional pieces of information are needed. Note, however, that the additional information is not specific to Addison. The needed information is as follows:
- The probability of being 65 years old. Suppose it is 0.2%
- The probability that a person with cancer is 65 years old. Suppose it is 0.5%. Note that this is greater than the previous value. This reflects that people with cancer are disproportionately 65 years old.
Knowing this, along with the base rate, we can calculate that a person (such as Addison) who is age 65 has a probability of having cancer equal to
<math>0.5% \div 0.2% * 1% = 2.5%</math>
It may come as a surprise that even though being 65 years old increases the risk of having cancer, Addison's probability of having cancer is still fairly low. This is because the base rate of cancer (regardless of age) is low. This illustrates both the importance of base rate, as well as that it is commonly neglected.^{[3]} Base rate neglect leads to serious misinterpretation of statistics; therefore, special care should be taken to avoid such mistakes. Becoming familiar with Bayes' Theorem is one way to combat the natural tendency to neglect base rates.
Drug testing
Suppose a drug test is 99% sensitive and 99% specific. That is, the test will produce 99% true positive results for drug users and 99% true negative results for non-drug users. Suppose that 0.5% of people are users of the drug. If a randomly selected individual tests positive, what is the probability he or she is a user?
- <math>
\begin{align} P(\text{User}\mid\text{+}) &= \frac{P(\text{+}\mid\text{User}) P(\text{User})}{P(\text{+}\mid\text{User}) P(\text{User}) + P(\text{+}\mid\text{Non-user}) P(\text{Non-user})} \\[8pt] &= \frac{0.99 \times 0.005}{0.99 \times 0.005 + 0.01 \times 0.995} \\[8pt] &\approx 33.2\% \end{align}</math>
Despite the apparent accuracy of the test, if an individual tests positive, it is more likely that they do not use the drug than that they do. This again illustrates the importance of base rates, and how the formation of policy can be eggregiously misguided if base rates are neglected.
This surprising result arises because the number of non-users is very large compared to the number of users; thus the number of false positives (0.995%) outweighs the number of true positives (0.495%). To use concrete numbers, if 1000 individuals are tested, there are expected to be 995 non-users and 5 users. From the 995 non-users, 0.01 × 995 ≃ 10 false positives are expected. From the 5 users, 0.99 × 5 ≃ 5 true positives are expected. Out of 15 positive results, only 5, about 33%, are genuine.
Note: The importance of specificity can be illustrated by showing that even if sensitivity is 100% and specificity is at 99% the probability of the person being a drug user is ≈33% but if the specificity is changed to 99.5% and the sensitivity is dropped down to 99% the probability of the person being a drug user rises to 49.8%.
More Complicated Example
The entire output of a factory is produced on three machines. The three machines account for 20%, 30%, and 50% of the output, respectively. The fraction of defective items produced is this: for the first machine, 5%; for the second machine, 3%; for the third machine, 1%. If an item is chosen at random from the total output and is found to be defective, what is the probability that it was produced by the third machine?
A solution is as follows. Let A_{i} denote the event that a randomly chosen item was made by the ith machine (for i = 1,2,3). Let B denote the event that a randomly chosen item is defective. Then, we are given the following information:
- P(A_{1}) = 0.2, P(A_{2}) = 0.3, P(A_{3}) = 0.5.
If the item was made by machine A_{1}, then the probability that it is defective is 0.05; that is, P(B | A_{1}) = 0.05. Overall, we have
- P(B | A_{1}) = 0.05, P(B | A_{2}) = 0.03, P(B | A_{3}) = 0.01.
To answer the original question, we first find P(B). That can be done in the following way:
- P(B) = Σ_{i} P(B | A_{i}) P(A_{i}) = (0.05)(0.2) + (0.03)(0.3) + (0.01)(0.5) = 0.024.
Hence 2.4% of the total output of the factory is defective.
We are given that B has occurred, and we want to calculate the conditional probability of A_{3}. By Bayes' theorem,
- P(A_{3} | B) = P(B | A_{3}) P(A_{3})/P(B) = (0.01)(0.50)/(0.024) = 5/24.
Given that the item is defective, the probability that it was made by the third machine is only 5/24. Although machine 3 produces half of the total output, it produces a much smaller fraction of the defective items. Hence the knowledge that the item selected was defective enables us to replace the prior probability P(A_{3}) = 1/2 by the smaller posterior probability P(A_{3} | B) = 5/24.
Interpretations
The interpretation of Bayes' theorem depends on the interpretation of probability ascribed to the terms. The two main interpretations are described below.
Bayesian interpretation
In the Bayesian (or epistemological) interpretation, probability measures a degree of belief. Bayes' theorem then links the degree of belief in a proposition before and after accounting for evidence. For example, suppose it is believed with 50% certainty that a coin is twice as likely to land heads than tails. If the coin is flipped a number of times and the outcomes observed, that degree of belief may rise, fall or remain the same depending on the results.
For proposition A and evidence B,
- P(A), the prior, is the initial degree of belief in A.
- P(A | B), the posterior, is the degree of belief having accounted for B.
- the quotient P(B | A)/P(B) represents the support B provides for A.
For more on the application of Bayes' theorem under the Bayesian interpretation of probability, see Bayesian inference.
Frequentist interpretation
In the frequentist interpretation, probability measures a proportion of outcomes. For example, suppose an experiment is performed many times. P(A) is the proportion of outcomes with property A, and P(B) that with property B. P(B | A) is the proportion of outcomes with property B out of outcomes with property A, and P(A | B) the proportion of those with A out of those with B.
The role of Bayes' theorem is best visualized with tree diagrams, as shown to the right. The two diagrams partition the same outcomes by A and B in opposite orders, to obtain the inverse probabilities. Bayes' theorem serves as the link between these different partitionings.
Example
An entomologist spots what might be a rare subspecies of beetle, due to the pattern on its back. In the rare subspecies, 98% have the pattern, or P(Pattern | Rare) = 98%. In the common subspecies, 5% have the pattern. The rare subspecies accounts for only 0.1% of the population. How likely is the beetle having the pattern to be rare, or what is P(Rare | Pattern)?
From the extended form of Bayes' theorem (since any beetle can be only rare or common),
<math>\begin{align}P(\text{Rare}\mid \text{Pattern}) &= \frac{P(\text{Pattern}\mid \text{Rare})P(\text{Rare})} {P(\text{Pattern}\mid \text{Rare})P(\text{Rare}) \, + \, P(\text{Pattern}\mid \text{Common})P(\text{Common})} \\[8pt] &= \frac{0.98 \times 0.001} {0.98 \times 0.001 + 0.05 \times 0.999} \\[8pt] &\approx 1.9\%. \end{align}</math>
Forms
Events
Simple form
For events A and B, provided that P(B) ≠ 0,
- <math>P(A\mid B) = \frac{P(B \mid A)\, P(A)}{P(B)}\cdot \,</math>
In many applications, for instance in Bayesian inference, the event B is fixed in the discussion, and we wish to consider the impact of its having been observed on our belief in various possible events A. In such a situation the denominator of the last expression, the probability of the given evidence B, is fixed; what we want to vary is A. Bayes' theorem then shows that the posterior probabilities are proportional to the numerator:
- <math>P(A\mid B) \propto P(A) \cdot P(B\mid A) \ </math> (proportionality over A for given B).
In words: posterior is proportional to prior times likelihood.^{[4]}
If events A_{1}, A_{2}, ..., are mutually exclusive and exhaustive, i.e., one of them is certain to occur but no two can occur together, and we know their probabilities up to proportionality, then we can determine the proportionality constant by using the fact that their probabilities must add up to one. For instance, for a given event A, the event A itself and its complement ¬A are exclusive and exhaustive. Denoting the constant of proportionality by c we have
- <math>P(A\mid B) = c \cdot P(A) \cdot P(B\mid A) \ </math> and <math>P(\neg A\mid B) = c \cdot P(\neg A) \cdot P(B\mid \neg A)\cdot </math>
Adding these two formulas we deduce that
- <math> c = \frac{1}{P(A) \cdot P(B\mid A) + P(\neg A) \cdot P(B\mid \neg A) } .</math>
Alternative form
Another form of Bayes' Theorem that is generally encountered when looking at two competing statements or hypotheses is:
- <math>P(A\mid B) = \frac{P(B\mid A)\,P(A)}{ P(B\mid A) P(A) + P(B\mid \neg A) P(\neg A)}\cdot</math>
For an epistemological interpretation:
For proposition A and evidence or background B,^{[5]}
- P(A),the prior probability, is the initial degree of belief in A.
- P(−A), is the corresponding probability of the initial degree of belief against A: 1 − P(A) = P(−A)
- P(B | A), the conditional probability or likelihood, is the degree of belief in B, given that the proposition A is true.
- P(B | -A), the conditional probability or likelihood, is the degree of belief in B, given that the proposition A is false.
- P(A | B), the posterior probability, is the probability for A after taking into account B for and against A.
Extended form
Often, for some partition {A_{j}} of the event space, the event space is given or conceptualized in terms of P(A_{j}) and P(B | A_{j}). It is then useful to compute P(B) using the law of total probability:
- <math>P(B) = {\sum_j P(B\mid A_j) P(A_j)},</math>
- <math>\Rightarrow P(A_i\mid B) = \frac{P(B\mid A_i)\,P(A_i)}{\sum\limits_j P(B\mid A_j)\,P(A_j)}\cdot</math>
In the special case where A is a binary variable:
- <math>P(A\mid B) = \frac{P(B\mid A)\,P(A)}{ P(B\mid A) P(A) + P(B\mid \neg A) P(\neg A)}\cdot</math>
Random variables
Consider a sample space Ω generated by two random variables X and Y. In principle, Bayes' theorem applies to the events A = {X = x} and B = {Y = y}. However, terms become 0 at points where either variable has finite probability density. To remain useful, Bayes' theorem may be formulated in terms of the relevant densities (see Derivation).
Simple form
If X is continuous and Y is discrete,
- <math>f_X(x\mid Y=y) = \frac{P(Y=y\mid X=x)\,f_X(x)}{P(Y=y)}.</math>
If X is discrete and Y is continuous,
- <math> P(X=x\mid Y=y) = \frac{f_Y(y\mid X=x)\,P(X=x)}{f_Y(y)}.</math>
If both X and Y are continuous,
- <math> f_X(x\mid Y=y) = \frac{f_Y(y\mid X=x)\,f_X(x)}{f_Y(y)}.</math>
Extended form
A continuous event space is often conceptualized in terms of the numerator terms. It is then useful to eliminate the denominator using the law of total probability. For f_{Y}(y), this becomes an integral:
- <math> f_Y(y) = \int_{-\infty}^\infty f_Y(y\mid X=\xi )\,f_X(\xi)\,d\xi .</math>
Bayes' rule
Bayes' rule is Bayes' theorem in odds form.
- <math>O(A_1:A_2\mid B) = O(A_1:A_2) \cdot \Lambda(A_1:A_2\mid B) </math>
where
- <math>\Lambda(A_1:A_2\mid B) = \frac{P(B\mid A_1)}{P(B\mid A_2)}</math>
is called the Bayes factor or likelihood ratio and the odds between two events is simply the ratio of the probabilities of the two events. Thus
- <math>O(A_1:A_2) = \frac{P(A_1)}{P(A_2)},</math>
- <math>O(A_1:A_2\mid B) = \frac{P(A_1\mid B)}{P(A_2\mid B)},</math>
So the rule says that the posterior odds are the prior odds times the Bayes factor, or in other words, posterior is proportional to prior times likelihood.
Derivation
For events
Bayes' theorem may be derived from the definition of conditional probability:
- <math>P(A\mid B)=\frac{P(A \cap B)}{P(B)}, \text{ if } P(B) \neq 0, \!</math>
- <math>P(B\mid A) = \frac{P(A \cap B)}{P(A)}, \text{ if } P(A) \neq 0, \!</math>
- <math>\Rightarrow P(A \cap B) = P(A\mid B)\, P(B) = P(B\mid A)\, P(A), \!</math>
- <math>\Rightarrow P(A\mid B) = \frac{P(B\mid A)\,P(A)}{P(B)}, \text{ if } P(B) \neq 0.</math>
For random variables
For two continuous random variables X and Y, Bayes' theorem may be analogously derived from the definition of conditional density:
- <math>f_X(x\mid Y=y) = \frac{f_{X,Y}(x,y)}{f_Y(y)} </math>
- <math>f_Y(y\mid X=x) = \frac{f_{X,Y}(x,y)}{f_X(x)} </math>
- <math>\Rightarrow f_X(x\mid Y=y) = \frac{f_Y(y\mid X=x)\,f_X(x)}{f_Y(y)}.</math>
History
Bayes' theorem was named after the Reverend Thomas Bayes (1701–61), who studied how to compute a distribution for the probability parameter of a binomial distribution (in modern terminology). Bayes' unpublished manuscript was significantly edited by Richard Price before it was posthumously read at the Royal Society. Price edited^{[6]} Bayes' major work An Essay towards solving a Problem in the Doctrine of Chances (1763), which appeared in Philosophical Transactions,^{[7]} and contains Bayes' Theorem. Price wrote an introduction to the paper which provides some of the philosophical basis of Bayesian statistics. In 1765 he was elected a Fellow of the Royal Society in recognition of his work on the legacy of Bayes.^{[8]}^{[9]}
The French mathematician Pierre-Simon Laplace reproduced and extended Bayes' results in 1774, apparently quite unaware of Bayes' work.^{[10]}^{[11]} Stephen Stigler suggested in 1983 that Bayes' theorem was discovered by Nicholas Saunderson some time before Bayes;^{[12]} that interpretation, however, has been disputed.^{[13]}
Martyn Hooper^{[14]} and Sharon McGrayne^{[15]} have argued that Richard Price's contribution was substantial:
By modern standards, we should refer to the Bayes–Price rule. Price discovered Bayes' work, recognized its importance, corrected it, contributed to the article, and found a use for it. The modern convention of employing Bayes' name alone is unfair but so entrenched that anything else makes little sense.—^{[15]}
See also
Notes
- ^ Jeffreys, Harold (1973). Scientific Inference (3rd ed.). Cambridge University Press. p. 31. ISBN 978-0-521-18078-8.
- ^ Stuart, A.; Ord, K. (1994), Kendall's Advanced Theory of Statistics: Volume I—Distribution Theory, Edward Arnold, §8.7.
- ^ Daniel Kahneman (25 October 2011). Thinking, Fast and Slow. Macmillan. ISBN 978-1-4299-6935-2. Retrieved 8 April 2012.
- ^ Lee, Peter M. (2012). "Chapter 1". Bayesian Statistics. Wiley. ISBN 978-1-1183-3257-3.
- ^ "Bayes Theorem: Introduction". Trinity University.
- ^ Richard Allen (1999). David Hartley on Human Nature. SUNY Press. pp. 243–4. ISBN 978-0-7914-9451-6. Retrieved 16 June 2013.
- ^ Bayes, Thomas, and Price, Richard (1763). "An Essay towards solving a Problem in the Doctrine of Chance. By the late Rev. Mr. Bayes, communicated by Mr. Price, in a letter to John Canton, A. M. F. R. S." (PDF). Philosophical Transactions of the Royal Society of London 53 (0): 370–418. doi:10.1098/rstl.1763.0053.
- ^ Holland, pp. 46–7.
- ^ Richard Price (1991). Price: Political Writings. Cambridge University Press. p. xxiii. ISBN 978-0-521-40969-8. Retrieved 16 June 2013.
- ^ Laplace refined Bayes' theorem over a period of decades:
- Laplace announced his independent discovery of Bayes' theorem in: Laplace (1774) "Mémoire sur la probabilité des causes par les événements," Mémoires de l'Académie royale des Sciences de MI (Savants étrangers), 4: 621–656. Reprinted in: Laplace, Oeuvres complètes (Paris, France: Gauthier-Villars et fils, 1841), vol. 8, pp. 27–65. Available on-line at: Gallica. Bayes' theorem appears on p. 29.
- Laplace presented a refinement of Bayes' theorem in: Laplace (read: 1783 / published: 1785) "Mémoire sur les approximations des formules qui sont fonctions de très grands nombres," Mémoires de l'Académie royale des Sciences de Paris, 423–467. Reprinted in: Laplace, Oeuvres complètes (Paris, France: Gauthier-Villars et fils, 1844), vol. 10, pp. 295–338. Available on-line at: Gallica. Bayes' theorem is stated on page 301.
- See also: Laplace, Essai philosophique sur les probabilités (Paris, France: Mme. Ve. Courcier [Madame veuve (i.e., widow) Courcier], 1814), page 10. English translation: Pierre Simon, Marquis de Laplace with F. W. Truscott and F. L. Emory, trans., A Philosophical Essay on Probabilities (New York, New York: John Wiley & Sons, 1902), page 15.
- ^ Daston, Lorraine (1988). Classical Probability in the Enlightenment. Princeton Univ Press. p. 268. ISBN 0-691-08497-1.
- ^ Stigler, Stephen M (1983). "Who Discovered Bayes' Theorem?". The American Statistician 37 (4): 290–296. doi:10.1080/00031305.1983.10483122.
- ^ Edwards, A. W. F. (1986). "Is the Reference in Hartley (1749) to Bayesian Inference?". The American Statistician 40 (2): 109–110. doi:10.1080/00031305.1986.10475370.
- ^ Hooper, Martyn (2013). "Richard Price, Bayes' theorem, and God". Significance 10 (1): 36–39. doi:10.1111/j.1740-9713.2013.00638.x.
- ^ ^{a} ^{b} McGrayne, S. B. (2011). The Theory That Would Not Die: How Bayes' Rule Cracked the Enigma Code, Hunted Down Russian Submarines & Emerged Triumphant from Two Centuries of Controversy. Yale University Press. ISBN 978-0-300-18822-6.
Further reading
- Bruss, F. Thomas (2013), "250 years of 'An Essay towards solving a Problem in the Doctrine of Chance. By the late Rev. Mr. Bayes, communicated by Mr. Price, in a letter to John Canton, A. M. F. R. S.' ", DOI 10.1365/s13291-013-0077-z, Jahresbericht der Deutschen Mathematiker-Vereinigung, Springer Verlag, Vol. 115, Issue 3-4 (2013), 129-133.
- Gelman, A, Carlin, JB, Stern, HS, and Rubin, DB (2003), "Bayesian Data Analysis", Second Edition, CRC Press.
- Grinstead, CM and Snell, JL (1997), "Introduction to Probability (2nd edition)", American Mathematical Society (free pdf available) [1].
- Hazewinkel, Michiel, ed. (2001), "Bayes formula", Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
- McGrayne, SB (2011). The Theory That Would Not Die: How Bayes' Rule Cracked the Enigma Code, Hunted Down Russian Submarines & Emerged Triumphant from Two Centuries of Controversy. Yale University Press. ISBN 978-0-300-18822-6.
- Laplace, P (1774/1986), "Memoir on the Probability of the Causes of Events", Statistical Science 1(3):364–378.
- Lee, PM (2012), "Bayesian Statistics: An Introduction", Wiley.
- Rosenthal, JS (2005), "Struck by Lightning: the Curious World of Probabilities". Harper Collings.
- Stigler, SM (1986). "Laplace's 1774 Memoir on Inverse Probability". Statistical Science 1 (3): 359–363. doi:10.1214/ss/1177013620.
- Stone, JV (2013), download chapter 1 of "Bayes' Rule: A Tutorial Introduction to Bayesian Analysis", Sebtel Press, England.
External links
- Bayes' theorem at Encyclopædia Britannica
- The Theory That Would Not Die by Sharon Bertsch McGrayne New York Times Book Review by John Allen Paulos on 5 August 2011
- Visual explanation of Bayes using trees (video)
- Bayes' frequentist interpretation explained visually (video)
- Earliest Known Uses of Some of the Words of Mathematics (B). Contains origins of "Bayesian", "Bayes' Theorem", "Bayes Estimate/Risk/Solution", "Empirical Bayes", and "Bayes Factor".
- Weisstein, Eric W., "Bayes' Theorem", MathWorld.
- Bayes' theorem at PlanetMath.org.
- Bayes Theorem and the Folly of Prediction
- A tutorial on probability and Bayes' theorem devised for Oxford University psychology students
- An Intuitive Explanation of Bayes' Theorem by Eliezer S. Yudkowsky