- Bayes' theorem
In
probability theory , Bayes' theorem (often called Bayes' law afterThomas Bayes ) relates the conditional and marginal probabilities of tworandom event s. It is often used to computeposterior probabilities given observations. For example, a patient may be observed to have certain symptoms. Bayes' theorem can be used to compute the probability that a proposed diagnosis is correct, given that observation. (See example 2)As a formal
theorem , Bayes' theorem is valid in all common interpretations of probability. However, it plays a central role in the debate around thefoundations of statistics :frequentist and Bayesian interpretations disagree about the ways in which probabilities should be assigned in applications. Frequentists assign probabilities to random events according to their frequencies of occurrence or to subsets of populations as proportions of the whole, while Bayesians describe probabilities in terms of beliefs and degrees of uncertainty. The articles onBayesian probability and frequentist probability discuss these debates at greater length.tatement of Bayes' theorem
Bayes' theorem relates the conditional and marginal probabilities of events "A" and "B", where "B" has a non-vanishing probability:
:
Each term in Bayes' theorem has a conventional name:
* P("A") is theprior probability ormarginal probability of "A". It is "prior" in the sense that it does not take into account any information about "B".
* P("A"|"B") is theconditional probability of "A", given "B". It is also called theposterior probability because it is derived from or depends upon the specified value of "B".
* P("B"|"A") is the conditional probability of "B" given "A".
* P("B") is the prior or marginal probability of "B", and acts as anormalizing constant .Intuitively, Bayes' theorem in this form describes the way in which one's beliefs about observing 'A' are updated by having observed 'B'.
An example
Suppose there is a co-ed school having 60% boys and 40% girls as students. The girl students wear trousers or skirts in equal numbers; the boys all wear trousers. An observer sees a (random) student from a distance; all they can see is that this student is wearing trousers. What is the probability this student is a girl?
It is clear that the probability is less than 40%, but by how much? Is it half that, since only half the girls are wearing trousers? The correct answer can be computed using Bayes' theorem.
The event "A" is that the student observed is a girl, and the event "B" is that the student observed is wearing trousers. To compute P("A"|"B"), we first need to know:
* P("A"), or the probability that the student is a girl regardless of any other information. Since the observers sees a random student, meaning that all students have the same probability of being observed, and the fraction of girls among the students is 40%, this probability equals 0.4.
* P("A"' ), or the probability that the student is a boy regardless of any other information ("A"' is the complementary event to "A"). This is 60%, or 0.6.
* P("B"|"A"), or the probability of the student wearing trousers given that the student is a girl. As they are as likely to wear skirts as trousers, this is 0.5.
* P("B"|"A"' ), or the probability of the student wearing trousers given that the student is a boy. This is given as 1.
* P("B"), or the probability of a (randomly selected) student wearing trousers regardless of any other information. Since P("B") = P("B"|"A")P("A") + P("B"|"A"' )P("A"' ), this is nowrap|1= 0.5×0.4 + 1×0.6 = 0.8.Given all this information, the probability of the observer having spotted a girl given that the observed student is wearing trousers can be computed by substituting these values in the formula:
:
As expected, it is less than 40%, but more than half that.
Another, essentially equivalent way of obtaining the same result is as follows. Assume, for concreteness, that there are 100 students, 60 boys and 40 girls. Among these, 60 boys and 20 girls wear trousers. All together there are 80 trouser-wearers, of which 20 are girls. Therefore the chance that a random trouser-wearer is a girl equals nowrap|1=20/80 = 0.25.
It is often helpful when calculating conditional probabilities to create a simple table containing the number of occurrences of each outcome, or the relative frequencies of each outcome, for each of the independent variables. The table below illustrates the use of this method for the above girl-or-boy example.
Bayes' theorem in terms of likelihood
Bayes' theorem can also be interpreted in terms of likelihood:
:
Here "L"("A"|"b") is the likelihood of "A" given fixed "b". The rule is then an immediate consequence of the relationship .
With this terminology, the theorem may be paraphrased as
:(where is a normalising constant equal to "P"("B")).
In words: the posterior probability is proportional to the product of the prior probability and the likelihood.
Derivation from conditional probabilities
To derive the theorem, we start from the definition of
conditional probability . The probability of event "A" given event "B" is:
Equivalently, the probability of event "B" given event "A" is
:
Rearranging and combining these two equations, we find
:
This lemma is sometimes called the product rule for probabilities. Dividing both sides by P("B"), providing that it is non-zero, we obtain Bayes' theorem:
:
Alternative forms of Bayes' theorem
Bayes' theorem is often embellished by noting that
:
where "A""C" is the complementary event of "A" (often called "not A"). So the theorem can be restated as the following formula
:
More generally, where {"A""i"} forms a partition of the event space,
:
for any "A""i" in the partition.
See also the
law of total probability .Bayes' theorem in terms of odds and likelihood ratio
Bayes' theorem can also be written neatly in terms of a likelihood ratio Λ and
odds "O" as:
where are the "odds" of "A" given "B",
and are the odds of "A" by itself,
while is the likelihood ratio.
Bayes' theorem for probability densities
There is also a version of Bayes' theorem for
continuous distribution s.It is somewhat harder to derive, since probability densities, strictly speaking, are not probabilities,so Bayes' theorem has to be established by a limit process;see Papoulis (citation below), Section 7.3 for an elementary derivation.Bayes' theorem for probability densities is formally similar to the theorem for probabilities::
There is an analogous statement of the
law of total probability , which is used in the denominator::As in the discrete case, the terms have standard names. : is the joint distribution of "X" and "Y",: is the posterior distribution of "X" given "Y"="y",: is (as a function of "x") the likelihood function of "X" given "Y"="y",and :and :
are the marginal distributions of "X" and "Y" respectively, with being the prior distribution of "X".
Abstract Bayes' theorem
Given two
absolutely continuous probability measures on theprobability space and a sigma-algebra , the abstract Bayes theorem for a -measurable random variable becomes:.
Proof :
by definition of conditional probability,
and we have also
This formulation is used in
Kalman filtering to findZakai equation s. It is also used infinancial mathematics for change ofnumeraire techniques.Extensions of Bayes' theorem
Theorems analogous to Bayes' theorem hold in problems with more than two variables. For example:
:
This can be derived in a few steps from Bayes' theorem and the definition of conditional probability:
:
Similarly,
:
which can be regarded as a conditional Bayes' Theorem, and can be derived by as follows:
:
A general strategy is to work with a decomposition of the
joint probability , and to marginalize (integrate) over the variables that are not of interest.Depending on the form of the decomposition,it may be possible to prove that some integrals must be 1,and thus they fall out of the decomposition;exploiting this property can reduce the computations very substantially.ABayesian network , for example, specifies a factorization of ajoint distribution of several variables in which the conditional probability of any one variable given the remaining ones takes a particularly simple form (seeMarkov blanket ).Further examples
Example 1: Drug testing
Bayes' theorem is useful in evaluating the result of
drug test s. Suppose a certain drug test is 99% sensitive and 99% specific, that is, the test will correctly identify a drug user as testing positive 99% of the time, and will correctly identify a non-user as testing negative 99% of the time. This would seem to be a relatively accurate test, but Bayes' theorem will reveal a potential flaw. Let's assume a corporation decides to test its employees foropium use, and 0.5% of the employees use the drug. We want to know theprobability that, given a positive drug test, an employee is actually a drug user. Let "D" be the event of being a drug user and "N" indicate being a non-user. Let "+" be the event of a positive drug test. We need to know the following:
* "P"("D"), or the probability that the employee is a drug user, regardless of any other information. This is 0.005, since 0.5% of the employees are drug users. This is the "prior probability" of D.
* "P"("N"), or the probability that the employee is not a drug user. This is nowrap|1 − "P"("D"), or 0.995.
* "P"(+|"D"), or the probability that the test is positive, given that the employee is a drug user. This is 0.99, since the test is 99% accurate.
* "P"(+|"N"), or the probability that the test is positive, given that the employee is not a drug user. This is 0.01, since the test will produce afalse positive for 1% of non-users.
* "P"(+), or the probability of a positive test event, regardless of other information. This is 0.0149 or 1.49%, which is found by adding the probability that a true positive result will appear (= 99% x 0.5% = 0.495%) plus the probability that a false positive will appear (= 1% x 99.5% = 0.995%). This is the prior probability of +.Given this information, we can compute the posterior probability "P"("D"|+) of an employee who tested positive actually being a drug user::
Despite the apparently high accuracy of the test, the probability that an employee who tested positive actually did use drugs is only about 33%, so it is "actually more likely" that the employee is not a drug user. The rarer the condition for which we are testing, the greater the percentage of positive tests that will be false positives.
Example 2: Bayesian inference
Applications of Bayes' theorem often assume the philosophy underlying
Bayesian probability that uncertainty and degrees of belief can be measured as probabilities. One such example follows. For additional worked out examples, including simpler examples, please see the article on the examples of Bayesian inference.We describe the marginal probability distribution of a variable "A" as the
prior probability distribution or simply the 'prior'. The conditional distribution of "A" given the "data" "B" is theposterior probability distribution or just the 'posterior'.Suppose we wish to know about the proportion r of voters in a large population who will vote "yes" in a referendum. Let n be the number of voters in a random sample (chosen with replacement, so that we have
statistical independence ) and let m be the number of voters in that random sample who will vote "yes". Suppose that we observe "n" = 10 voters and "m" = 7 say they will vote yes. From Bayes' theorem we can calculate the probability distribution function for "r" using:
From this we see that from the prior probability density function "f"("r") and the likelihood function "L"("r") = "f"("m" = 7|"r", "n" = 10), we can compute the posterior probability density function "f"("r"|"n" = 10, "m" = 7).
The prior probability density function "f"("r") summarizes what we know about the distribution of "r" in the absence of any observation. We provisionally assume in this case that the prior distribution of "r" is uniform over the interval [0, 1] . That is, "f"("r") = 1. If some additional background information is found, we should modify the prior accordingly. However before we have any observations, all outcomes are equally likely.
Under the assumption of random sampling, choosing voters is just like choosing balls from an urn. The likelihood function "L"("r") = "P"("m" = 7|"r", "n" = 10,) for such a problem is just the probability of 7 successes in 10 trials for a
binomial distribution .:
As with the prior, the likelihood is open to revision -- more complex assumptions will yield more complex likelihood functions. Maintaining the current assumptions, we compute the normalizing factor,
:
and the posterior distribution for "r" is then
:
for "r" between 0 and 1, inclusive.
One may be interested in the probability that more than half the voters will vote "yes". The prior probability that more than half the voters will vote "yes" is 1/2, by the symmetry of the
uniform distribution . In comparison, the posterior probability that more than half the voters will vote "yes", i.e., the conditional probability given the outcome of the opinion poll – that seven of the 10 voters questioned will vote "yes" – is:
which is about an "89% chance".
Example 3: The Monty Hall problem
We are presented with three doors - red, green, and blue - one of which has a prize. We choose the red door, which is not opened until the presenter performs an action. The presenter "who knows what door the prize is behind, and who must open a door, but is not permitted to open the door we have picked or the door with the prize", opens the "blue" door and reveals that there is no prize behind it and subsequently asks if we wish to change our mind about our initial selection of red. What is the probability that the prize is behind each of the green and red doors?
Let us call the situation that the prize is behind a given door Ar, Ag, and Ab.
To start with, , and to make things simpler we shall assume that we have already picked the red door.
Let us call B "the presenter opens the blue door". Without any prior knowledge, we would assign this a probability of 50%.
* In the situation where the prize is behind the red door, the host is free to pick between the green or the blue door at random. Thus,
* In the situation where the prize is behind the green door, the host must pick the blue door. Thus,
* In the situation where the prize is behind the blue door, the host must pick the green door. Thus,Thus,
So, we should always choose the green door.
Note how this depends on the value of P(B).
Historical remarks
An investigation by a statistics professor (Stigler 1983) suggests that Bayes' theorem was discovered by
Nicholas Saunderson some time before Bayes.Bayes' theorem is named after the Reverend
Thomas Bayes (1702–1761), who studied how to compute a distribution for the parameter of abinomial distribution (to use modern terminology). His friend,Richard Price , edited and presented the work in 1763, after Bayes' death, as "An Essay towards solving a Problem in the Doctrine of Chances".Pierre-Simon Laplace replicated and extended these results in an essay of 1774, apparently unaware of Bayes' work.One of Bayes' results (Proposition 5) gives a simple description of
conditional probability , and shows that it can be expressed independently of the order in which things occur::If there be two subsequent events, the probability of the second b/N and the probability of both together P/N, and it being first discovered that the second event has also happened, from hence I guess that the first event has also happened, the probability I am right [i.e., the conditional probability of the first event being true given that the second has also happened] is P/b.
Note that the expression says nothing about the "order" in which the events occurred; it measures correlation, not causation. His preliminary results, in particular Propositions 3, 4, and 5, imply the result now called Bayes' Theorem (as described above), but it does not appear that Bayes himself emphasized or focused on that result.
Bayes' main result (Proposition 9 in the essay) is the following: assuming a
uniform distribution for theprior distribution of thebinomial parameter "p", the probability that "p" is between two values "a" and "b" is:
where "m" is the number of observed successes and "n" the number of observed failures.
What is "Bayesian" about Proposition 9 is that Bayes presented it as a probability for the parameter "p". So, one can compute probability for an experimental outcome, but also for the parameter which governs it, and the same algebra is used to make inferences of either kind.
Bayes states his question in a way that might make the idea of assigning a probability distribution to a parameter palatable to a frequentist. He supposes that a billiard ball is thrown at random onto a billiard table, and that the probabilities "p" and "q" are the probabilities that subsequent billiard balls will fall above or below the first ball.
Stephen Fienberg [http://ba.stat.cmu.edu/journal/2006/vol01/issue01/fienberg.pdf] describes the evolution of the field from "inverse probability" at the time of Bayes and Laplace, and even ofHarold Jeffreys (1939) to "Bayesian" in the 1950's. The irony is that this label was introduced byR.A. Fisher in a derogatory sense. So, historically, Bayes was not a "Bayesian". It is actually unclear whether or not he was a Bayesian in the modern sense of the term, i.e. whether or not he was interested in inference or merely in probability: the 1763 essay is more of a probability paper.See also
*Bayesian inference
*Bayesian network
*Bayesian probability
*Bayesian spam filtering
*Thomas Bayes
*Bogofilter
*Conjugate prior
*Empirical Bayes method
*Monty Hall problem
*Occam's razor
*Prosecutor's fallacy
*Raven paradox
*Recursive Bayesian estimation
*Revising opinions in statistics
*Sequential bayesian filtering
*Borel's paradox
*Naive Bayes classifier References
Versions of the essay
* Thomas Bayes (1763), "An Essay towards solving a Problem in the Doctrine of Chances. By the late Rev. Mr. Bayes, F. R. S. communicated by Mr. Price, in a letter to John Canton, A. M. F. R. S.", "
Philosophical Transactions , Giving Some Account of the Present Undertakings, Studies and Labours of the Ingenious in Many Considerable Parts of the World" 53:370–418.
* Thomas Bayes (1763/1958) "Studies in the History of Probability and Statistics: IX. Thomas Bayes' Essay Towards Solving a Problem in the Doctrine of Chances", "Biometrika " 45:296–315. "(Bayes' essay in modernized notation)"
* Thomas Bayes [http://www.stat.ucla.edu/history/essay.pdf "An essay towards solving a Problem in the Doctrine of Chances"] . "(Bayes' essay in the original notation)"Commentaries
* G. A. Barnard (1958) "Studies in the History of Probability and Statistics: IX. Thomas Bayes' Essay Towards Solving a Problem in the Doctrine of Chances", "Biometrika" 45:293–295. "(biographical remarks)"
* Daniel Covarrubias. [http://www.stat.rice.edu/~blairc/seminar/Files/danTalk.pdf "An Essay Towards Solving a Problem in the Doctrine of Chances"] . "(an outline and exposition of Bayes' essay)"
* Stephen M. Stigler (1982). "Thomas Bayes' Bayesian Inference," "Journal of the Royal Statistical Society", Series A, 145:250–258. "(Stigler argues for a revised interpretation of the essay; recommended)"
*Isaac Todhunter (1865). "A History of the Mathematical Theory of Probability from the time of Pascal to that of Laplace", Macmillan. Reprinted 1949, 1956 by Chelsea and 2001 by Thoemmes.
* [http://yudkowsky.net/bayes/bayes.html An Intuitive Explanation of Bayesian Reasoning] (includes biography)Additional material
* Pierre-Simon Laplace (1774/1986), "Memoir on the Probability of the Causes of Events", "Statistical Science" 1(3):364–378.
* Stephen M. Stigler (1986), "Laplace's 1774 memoir on inverse probability", "Statistical Science" 1(3):359–378.
* Stephen M. Stigler (1983), "Who Discovered Bayes' Theorem?" "The American Statistician" 37(4):290–296.
* Jeff Miller, "et al.", [http://members.aol.com/jeff570/mathword.html Earliest Known Uses of Some of the Words of Mathematics (B)] . ("very informative; recommended")
*Athanasios Papoulis (1984), "Probability, Random Variables, and Stochastic Processes", second edition. New York: McGraw-Hill.
* [http://www.inference.phy.cam.ac.uk/mackay/itila/ The on-line textbook: Information Theory, Inference, and Learning Algorithms] , byDavid J. C. MacKay provides an up to date overview of the use of Bayes' theorem in information theory and machine learning.
* Provides a comprehensive introduction to Bayes' theorem.
* [http://plato.stanford.edu/entries/logic-inductive/ Stanford Encyclopedia of Philosophy: Inductive Logic] provides a comprehensive Bayesian treatment of Inductive Logic and Confirmation Theory.
*
*
*Eliezer S. Yudkowsky (2003), " [http://yudkowsky.net/bayes/bayes.html An Intuitive Explanation of Bayesian Reasoning] "
* [http://www.celiagreen.com/charlesmccreery/statistics/bayestutorial.pdf A tutorial on probability and Bayes’ theorem devised for Oxford University psychology students]
* [http://faculty-staff.ou.edu/H/James.A.Hawthorne-1/Hawthorne%20--%20Confirmation%20Theory.pdf Confirmation Theory] An extensive presentation of Bayesian Confirmation Theory
Wikimedia Foundation. 2010.