Bayes' theorem

From HandWiki
Short description: Probability based on prior knowledge

In probability theory and statistics, Bayes' theorem (alternatively Bayes' law or Bayes' rule), named after Thomas Bayes, describes the probability of an event, based on prior knowledge of conditions that might be related to the event.[1] For example, if the risk of developing health problems is known to increase with age, Bayes' theorem allows the risk to an individual of a known age to be assessed more accurately by conditioning it relative to their age, rather than simply assuming that the individual is typical of the population as a whole.

One of the many applications of Bayes' theorem is Bayesian inference, a particular approach to statistical inference. When applied, the probabilities involved in the theorem may have different probability interpretations. With Bayesian probability interpretation, the theorem expresses how a degree of belief, expressed as a probability, should rationally change to account for the availability of related evidence. Bayesian inference is fundamental to Bayesian statistics, being considered by one authority as "to the theory of probability what Pythagoras's theorem is to geometry."[2]

History

Bayes' theorem is named after the Reverend Thomas Bayes (/bz/), also a statistician and philosopher. Bayes used conditional probability to provide an algorithm (his Proposition 9) that uses evidence to calculate limits on an unknown parameter. His work was published in 1763 as An Essay towards solving a Problem in the Doctrine of Chances. Bayes studied how to compute a distribution for the probability parameter of a binomial distribution (in modern terminology). On Bayes's death his family transferred his papers to a friend, the minister, philosopher, and mathematician Richard Price.

Over two years, Richard Price significantly edited the unpublished manuscript, before sending it to a friend who read it aloud at the Royal Society on 23 December 1763.[3] Price edited[4] Bayes's major work "An Essay towards solving a Problem in the Doctrine of Chances" (1763), which appeared in Philosophical Transactions,[5] and contains Bayes' theorem. Price wrote an introduction to the paper which provides some of the philosophical basis of Bayesian statistics and chose one of the two solutions offered by Bayes. In 1765, Price was elected a Fellow of the Royal Society in recognition of his work on the legacy of Bayes.[6][7] On 27 April a letter sent to his friend Benjamin Franklin was read out at the Royal Society, and later published, where Price applies this work to population and computing 'life-annuities'.[8]

Independently of Bayes, Pierre-Simon Laplace in 1774, and later in his 1812 Théorie analytique des probabilités, used conditional probability to formulate the relation of an updated posterior probability from a prior probability, given evidence. He reproduced and extended Bayes's results in 1774, apparently unaware of Bayes's work.[note 1][9] The Bayesian interpretation of probability was developed mainly by Laplace.[10]

About 200 years later, Sir Harold Jeffreys put Bayes's algorithm and Laplace's formulation on an axiomatic basis, writing in a 1973 book that Bayes' theorem "is to the theory of probability what the Pythagorean theorem is to geometry".[2]

Stephen Stigler used a Bayesian argument to conclude that Bayes' theorem was discovered by Nicholas Saunderson, a blind English mathematician, some time before Bayes;[11][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's work, recognized its importance, corrected it, contributed to the article, and found a use for it. The modern convention of employing Bayes's name alone is unfair but so entrenched that anything else makes little sense.[15]

Statement of theorem

Bayes' theorem is stated mathematically as the following equation:[16]

[math]\displaystyle{ P(A\vert B) = \frac{P(B \vert A) P(A)}{P(B)} }[/math]

where [math]\displaystyle{ A }[/math] and [math]\displaystyle{ B }[/math] are events and [math]\displaystyle{ P(B) \neq 0 }[/math].

  • [math]\displaystyle{ P(A\vert B) }[/math] is a conditional probability: the probability of event [math]\displaystyle{ A }[/math] occurring given that [math]\displaystyle{ B }[/math] is true. It is also called the posterior probability of [math]\displaystyle{ A }[/math] given [math]\displaystyle{ B }[/math].
  • [math]\displaystyle{ P(B\vert A) }[/math] is also a conditional probability: the probability of event [math]\displaystyle{ B }[/math] occurring given that [math]\displaystyle{ A }[/math] is true. It can also be interpreted as the likelihood of [math]\displaystyle{ A }[/math] given a fixed [math]\displaystyle{ B }[/math] because [math]\displaystyle{ P(B\vert A)=L(A\vert B) }[/math].
  • [math]\displaystyle{ P(A) }[/math] and [math]\displaystyle{ P(B) }[/math] are the probabilities of observing [math]\displaystyle{ A }[/math] and [math]\displaystyle{ B }[/math] respectively without any given conditions; they are known as the prior probability and marginal probability.

Proof

Visual proof of Bayes' theorem

For events

Bayes' theorem may be derived from the definition of conditional probability:

[math]\displaystyle{ P(A\vert B)=\frac{P(A \cap B)}{P(B)}, \text{ if } P(B) \neq 0, }[/math]

where [math]\displaystyle{ P(A \cap B) }[/math] is the probability of both A and B being true. Similarly,

[math]\displaystyle{ P(B\vert A)=\frac{P(A \cap B)}{P(A)}, \text{ if } P(A) \neq 0. }[/math]

Solving for [math]\displaystyle{ P(A \cap B) }[/math] and substituting into the above expression for [math]\displaystyle{ P(A\vert B) }[/math] yields Bayes' theorem:

[math]\displaystyle{ P(A\vert B) = \frac{P(B\vert A) P(A)}{P(B)}, \text{ if } P(B) \neq 0. }[/math]

For continuous random variables

For two continuous random variables X and Y, Bayes' theorem may be analogously derived from the definition of conditional density:

[math]\displaystyle{ f_{X \vert Y=y} (x) = \frac{f_{X,Y}(x,y)}{f_Y(y)} }[/math]
[math]\displaystyle{ f_{Y \vert X=x}(y) = \frac{f_{X,Y}(x,y)}{f_X(x)} }[/math]

Therefore,

[math]\displaystyle{ f_{X \vert Y=y}(x) = \frac{f_{Y \vert X=x}(y) f_X(x)}{f_Y(y)}. }[/math]

General case

Let [math]\displaystyle{ P_Y^x }[/math] be the conditional distribution of [math]\displaystyle{ Y }[/math] given [math]\displaystyle{ X = x }[/math] and let [math]\displaystyle{ P_X }[/math] be the distribution of [math]\displaystyle{ X }[/math]. The joint distribution is then [math]\displaystyle{ P_{X,Y} (dx,dy) = P_Y^x (dy) P_X (dx) }[/math]. The conditional distribution [math]\displaystyle{ P_X^y }[/math] of [math]\displaystyle{ X }[/math] given [math]\displaystyle{ Y=y }[/math] is then determined by

[math]\displaystyle{ P_X^y (A) = E (1_A (X) | Y = y) }[/math]

Existence and uniqueness of the needed conditional expectation is a consequence of the Radon–Nikodym theorem. This was formulated by Kolmogorov in his famous book from 1933. Kolmogorov underlines the importance of conditional probability by writing "I wish to call attention to ... and especially the theory of conditional probabilities and conditional expectations ..." in the Preface.[17] The Bayes theorem determines the posterior distribution from the prior distribution. Bayes' theorem can be generalized to include improper prior distributions such as the uniform distribution on the real line.[18] Modern Markov chain Monte Carlo methods have boosted the importance of Bayes' theorem including cases with improper priors.[19]

Examples

Recreational mathematics

Bayes' rule and computing conditional probabilities provide a solution method for a number of popular puzzles, such as the Three Prisoners problem, the Monty Hall problem, the Two Child problem and the Two Envelopes problem.

Drug testing

Figure 1: Using a frequency box to show [math]\displaystyle{ P(\text{User}\vert \text{Positive}) }[/math] visually by comparison of shaded areas

Suppose, a particular test for whether someone has been using cannabis is 90% sensitive, meaning the true positive rate (TPR) = 0.90. Therefore, it leads to 90% true positive results (correct identification of drug use) for cannabis users.

The test is also 80% specific, meaning true negative rate (TNR) = 0.80. Therefore, the test correctly identifies 80% of non-use for non-users, but also generates 20% false positives, or false positive rate (FPR) = 0.20, for non-users.

Assuming 0.05 prevalence, meaning 5% of people use cannabis, what is the probability that a random person who tests positive is really a cannabis user?

The Positive predictive value (PPV) of a test is the proportion of persons who are actually positive out of all those testing positive, and can be calculated from a sample as:

PPV = True positive / Tested positive

If sensitivity, specificity, and prevalence are known, PPV can be calculated using Bayes theorem. Let [math]\displaystyle{ P(\text{User}\vert \text{Positive}) }[/math] mean "the probability that someone is a cannabis user given that they test positive," which is what is meant by PPV. We can write:

[math]\displaystyle{ \begin{align} P(\text{User}\vert \text{Positive}) &= \frac{P(\text{Positive}\vert \text{User}) P(\text{User})}{P(\text{Positive})} \\ &= \frac{P(\text{Positive}\vert\text{User}) P(\text{User})}{P(\text{Positive}\vert\text{User}) P(\text{User}) + P(\text{Positive}\vert\text{Non-user}) P(\text{Non-user})} \\[8pt] &= \frac{0.90 \times 0.05}{0.90 \times 0.05 + 0.20 \times 0.95} = \frac{0.045}{0.045 + 0.19} \approx 19\% \end{align} }[/math]

The fact that [math]\displaystyle{ P(\text{Positive}) = P(\text{Positive}\vert\text{User}) P(\text{User}) + P(\text{Positive}\vert\text{Non-user}) P(\text{Non-user}) }[/math] is a direct application of the Law of Total Probability. In this case, it says that the probability that someone tests positive is the probability that a user tests positive, times the probability of being a user, plus the probability that a non-user tests positive, times the probability of being a non-user. This is true because the classifications user and non-user form a partition of a set, namely the set of people who take the drug test. This combined with the definition of conditional probability results in the above statement.

In other words, even if someone tests positive, the probability that they are a cannabis user is only 19%—this is because in this group, only 5% of people are users, and most positives are false positives coming from the remaining 95%.

If 1,000 people were tested:

  • 950 are non-users and 190 of them give false positive (0.20 × 950)
  • 50 of them are users and 45 of them give true positive (0.90 × 50)

The 1,000 people thus yields 235 positive tests, of which only 45 are genuine drug users, about 19%. See Figure 1 for an illustration using a frequency box, and note how small the pink area of true positives is compared to the blue area of false positives.

Sensitivity or specificity

The importance of specificity can be seen by showing that even if sensitivity is raised to 100% and specificity remains at 80%, the probability of someone testing positive really being a cannabis user only rises from 19% to 21%, but if the sensitivity is held at 90% and the specificity is increased to 95%, the probability rises to 49%.

Test
Actual
Positive Negative Total
User 45 5 50
Non-user 190 760 950
Total 235 765 1000
90% sensitive, 80% specific, PPV=45/235 ≈ 19%
Test
Actual
Positive Negative Total
User 50 0 50
Non-user 190 760 950
Total 240 760 1000
100% sensitive, 80% specific, PPV=50/240 ≈ 21%
Test
Actual
Positive Negative Total
User 45 5 50
Non-user 47 903 950
Total 92 908 1000
90% sensitive, 95% specific, PPV=45/92 ≈ 49%

Cancer rate

Even if 100% of patients with pancreatic cancer have a certain symptom, when someone has the same symptom, it does not mean that this person has a 100% chance of getting pancreatic cancer. Assuming the incidence rate of pancreatic cancer is 1/100000, while 10/99999 healthy individuals have the same symptoms worldwide, the probability of having pancreatic cancer given the symptoms is only 9.1%, and the other 90.9% could be "false positives" (that is, falsely said to have cancer; "positive" is a confusing term when, as here, the test gives bad news).

Based on incidence rate, the following table presents the corresponding numbers per 100,000 people.

Symptom
Cancer
Yes No Total
Yes 1 0 1
No 10 99989 99999
Total 11 99989 100000

Which can then be used to calculate the probability of having cancer when you have the symptoms:

[math]\displaystyle{ \begin{align} P(\text{Cancer}|\text{Symptoms}) &= \frac{P(\text{Symptoms}|\text{Cancer}) P(\text{Cancer})}{P(\text{Symptoms})} \\ &= \frac{P(\text{Symptoms}|\text{Cancer}) P(\text{Cancer})}{P(\text{Symptoms}|\text{Cancer}) P(\text{Cancer}) + P(\text{Symptoms}|\text{Non-Cancer}) P(\text{Non-Cancer})} \\[8pt] &= \frac{1 \times 0.00001}{1 \times 0.00001 + (10/99999) \times 0.99999} = \frac1{11} \approx 9.1\% \end{align} }[/math]

Defective item rate

Condition

Machine
Defective Flawless Total
A 10 190 200
B 9 291 300
C 5 495 500
Total 24 976 1000

A factory produces items using three machines—A, B, and C—which account for 20%, 30%, and 50% of its output respectively. Of the items produced by machine A, 5% are defective; similarly, 3% of machine B's items and 1% of machine C's are defective. If a randomly selected item is defective, what is the probability it was produced by machine C?

Once again, the answer can be reached without using the formula by applying the conditions to a hypothetical number of cases. For example, if the factory produces 1,000 items, 200 will be produced by Machine A, 300 by Machine B, and 500 by Machine C. Machine A will produce 5% × 200 = 10 defective items, Machine B 3% × 300 = 9, and Machine C 1% × 500 = 5, for a total of 24. Thus, the likelihood that a randomly selected defective item was produced by machine C is 5/24 (~20.83%).

This problem can also be solved using Bayes' theorem: Let Xi denote the event that a randomly chosen item was made by the i th machine (for i = A,B,C). Let Y denote the event that a randomly chosen item is defective. Then, we are given the following information:

[math]\displaystyle{ P(X_A) = 0.2, \quad P(X_B) = 0.3, \quad P(X_C) = 0.5. }[/math]

If the item was made by the first machine, then the probability that it is defective is 0.05; that is, P(Y | XA) = 0.05. Overall, we have

[math]\displaystyle{ P(Y| X_A) = 0.05, \quad P(Y |X_B) = 0.03, \quad P(Y| X_C) = 0.01. }[/math]

To answer the original question, we first find P(Y). That can be done in the following way:

[math]\displaystyle{ P(Y) = \sum_i P(Y| X_i) P(X_i) = (0.05)(0.2) + (0.03)(0.3) + (0.01)(0.5) = 0.024. }[/math]

Hence, 2.4% of the total output is defective.

We are given that Y has occurred, and we want to calculate the conditional probability of XC. By Bayes' theorem,

[math]\displaystyle{ P(X_C|Y) = \frac{P(Y | X_C) P(X_C)}{P(Y)} = \frac{0.01 \cdot 0.50}{0.024} = \frac{5}{24} }[/math]

Given that the item is defective, the probability that it was made by machine C is 5/24. Although machine C 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(XC) = 1/2 by the smaller posterior probability P(XC | Y) = 5/24.

Interpretations

Figure 2: A geometric visualisation of Bayes' theorem

The interpretation of Bayes' rule depends on the interpretation of probability ascribed to the terms. The two predominant interpretations are described below. Figure 2 shows a geometric visualization.

Bayesian interpretation

In the Bayesian (or epistemological) interpretation, probability measures a "degree of belief". Bayes' theorem 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 will probably rise or fall, but might even 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 after incorporating news that B is true.
  • 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

Figure 3: Illustration of frequentist interpretation with tree diagrams

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 (the prior) and P(B) is the proportion with property B. P(B | A) is the proportion of outcomes with property B out of outcomes with property A, and P(A | B) is the proportion of those with A out of those with B (the posterior).

The role of Bayes' theorem is best visualized with tree diagrams such as Figure 3. The two diagrams partition the same outcomes by A and B in opposite orders, to obtain the inverse probabilities. Bayes' theorem links the different partitionings.

Example

Figure 4: Tree diagram illustrating the beetle example. R, C, P and [math]\displaystyle{ \overline{P} }[/math] are the events rare, common, pattern and no pattern. Percentages in parentheses are calculated. Three independent values are given, so it is possible to calculate the inverse tree.

An entomologist spots what might, due to the pattern on its back, be a rare subspecies of beetle. A full 98% of the members of the rare subspecies have the pattern, so P(Pattern | Rare) = 98%. Only 5% of members of the common subspecies have the pattern. The rare subspecies is 0.1% of the total population. How likely is the beetle having the pattern to be rare: what is P(Rare | Pattern)?

From the extended form of Bayes' theorem (since any beetle is either rare or common),

[math]\displaystyle{ \begin{align} P(\text{Rare} \vert \text{Pattern}) &= \frac{P(\text{Pattern} \vert \text{Rare})P(\text{Rare})} {P(\text{Pattern})}\\ [8pt] &= \frac{P(\text{Pattern}\vert \text{Rare})P(\text{Rare})} {P(\text{Pattern} \vert \text{Rare}) P(\text{Rare}) + P(\text{Pattern}\vert \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]\displaystyle{ P(A| B) = \frac{P(B | A) P(A)}{P(B)} . }[/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, so the last equation becomes:

[math]\displaystyle{ P(A| B) \propto P(A) \cdot P(B| A) . }[/math]

In words, the posterior is proportional to the prior times the likelihood.[20]

If events A1, A2, ..., are mutually exclusive and exhaustive, i.e., one of them is certain to occur but no two can occur together, 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]\displaystyle{ P(A| B) = c \cdot P(A) \cdot P(B| A) \text{ and } P(\neg A| B) = c \cdot P(\neg A) \cdot P(B| \neg A). }[/math]

Adding these two formulas we deduce that

[math]\displaystyle{ 1 = c \cdot (P(B| A)\cdot P(A) + P(B| \neg A) \cdot P(\neg A)), }[/math]

or

[math]\displaystyle{ c = \frac{1}{P(B| A)\cdot P(A) + P(B| \neg A) \cdot P(\neg A)} = \frac 1 {P(B)}. }[/math]

Alternative form

Contingency table
  Background

Proposition
B ¬B
(not B)
Total
A P(B|A)·P(A)
= P(A|B)·P(B)
P(¬B|A)·P(A)
= P(A|¬B)·P(¬B)
P(A)
¬A
(not A)
P(B|¬A)·P(¬A)
= P(¬A|B)·P(B)
P(¬B|¬A)·P(¬A)
= P(¬A|¬B)·P(¬B)
P(¬A) =
1−P(A)
Total    P(B)    P(¬B) = 1−P(B) 1

Another form of Bayes' theorem for two competing statements or hypotheses is:

[math]\displaystyle{ P(A| B) = \frac{P(B| A) P(A)}{ P(B| A) P(A) + P(B| \neg A) P(\neg A)}. }[/math]

For an epistemological interpretation:

For proposition A and evidence or background B,[21]

  • [math]\displaystyle{ P(A) }[/math] is the prior probability, the initial degree of belief in A.
  • [math]\displaystyle{ P(\neg A) }[/math] is the corresponding initial degree of belief in not-A, that A is false, where [math]\displaystyle{ P(\neg A) =1-P(A) }[/math]
  • [math]\displaystyle{ P(B| A) }[/math] is the conditional probability or likelihood, the degree of belief in B given that proposition A is true.
  • [math]\displaystyle{ P(B|\neg A) }[/math] is the conditional probability or likelihood, the degree of belief in B given that proposition A is false.
  • [math]\displaystyle{ P(A| B) }[/math] is the posterior probability, the probability of A after taking into account B.

Extended form

Often, for some partition {Aj} of the sample space, the event space is given in terms of P(Aj) and P(B | Aj). It is then useful to compute P(B) using the law of total probability:

[math]\displaystyle{ P(B)=\sum_{j}P(B \cap A_j), }[/math]

Or (using the multiplication rule for conditional probability),[22]

[math]\displaystyle{ P(B) = {\sum_j P(B| A_j) P(A_j)}, }[/math]
[math]\displaystyle{ \Rightarrow P(A_i| B) = \frac{P(B| A_i) P(A_i)}{\sum\limits_j P(B| A_j) P(A_j)}\cdot }[/math]

In the special case where A is a binary variable:

[math]\displaystyle{ P(A| B) = \frac{P(B| A) P(A)}{ P(B| A) P(A) + P(B| \neg A) P(\neg A)}\cdot }[/math]

Random variables

Figure 5: Bayes' theorem applied to an event space generated by continuous random variables X and Y with known probability distibutions. There exists an instance of Bayes' theorem for each point in the domain. In practice, these instances might be parametrized by writing the specified probability densities as a function of x and y.

Consider a sample space Ω generated by two random variables X and Y with known probability distributions. In principle, Bayes' theorem applies to the events A = {X = x} and B = {Y = y}.

[math]\displaystyle{ P( X{=}x | Y {=} y) = \frac{P(Y{=}y | X{=}x) P(X{=}x)}{P(Y{=}y)} }[/math]

However, terms become 0 at points where either variable has finite probability density. To remain useful, Bayes' theorem can be formulated in terms of the relevant densities (see Derivation).

Simple form

If X is continuous and Y is discrete,

[math]\displaystyle{ f_{X | Y{=}y}(x) = \frac{P(Y{=}y| X{=}x) f_X(x)}{P(Y{=}y)} }[/math]

where each [math]\displaystyle{ f }[/math] is a density function.

If X is discrete and Y is continuous,

[math]\displaystyle{ P(X{=}x| Y{=}y) = \frac{f_{Y | X{=}x}(y) P(X{=}x)}{f_Y(y)}. }[/math]

If both X and Y are continuous,

[math]\displaystyle{ f_{X| Y{=}y}(x) = \frac{f_{Y | X{=}x}(y) f_X(x)}{f_Y(y)}. }[/math]

Extended form

Figure 6: A way to conceptualize event spaces generated by continuous random variables X and Y

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 fY(y), this becomes an integral:

[math]\displaystyle{ f_Y(y) = \int_{-\infty}^\infty f_{Y| X = \xi}(y) f_X(\xi)\,d\xi . }[/math]

Bayes' rule in odds form

Bayes' theorem in odds form is:

[math]\displaystyle{ O(A_1:A_2\vert B) = O(A_1:A_2) \cdot \Lambda(A_1:A_2\vert B) }[/math]

where

[math]\displaystyle{ \Lambda(A_1:A_2\vert B) = \frac{P(B\vert A_1)}{P(B\vert A_2)} }[/math]

is called the Bayes factor or likelihood ratio. The odds between two events is simply the ratio of the probabilities of the two events. Thus

[math]\displaystyle{ O(A_1:A_2) = \frac{P(A_1)}{P(A_2)}, }[/math]
[math]\displaystyle{ O(A_1:A_2\vert B) = \frac{P(A_1\vert B)}{P(A_2\vert B)}, }[/math]

Thus, the rule says that the posterior odds are the prior odds times the Bayes factor, or in other words, the posterior is proportional to the prior times the likelihood.

In the special case that [math]\displaystyle{ A_1 = A }[/math] and [math]\displaystyle{ A_2 = \neg A }[/math], one writes [math]\displaystyle{ O(A)=O(A:\neg A) =P(A)/(1-P(A)) }[/math], and uses a similar abbreviation for the Bayes factor and for the conditional odds. The odds on [math]\displaystyle{ A }[/math] is by definition the odds for and against [math]\displaystyle{ A }[/math]. Bayes' rule can then be written in the abbreviated form

[math]\displaystyle{ O(A\vert B) = O(A) \cdot \Lambda(A\vert B) , }[/math]

or, in words, the posterior odds on [math]\displaystyle{ A }[/math] equals the prior odds on [math]\displaystyle{ A }[/math] times the likelihood ratio for [math]\displaystyle{ A }[/math] given information [math]\displaystyle{ B }[/math]. In short, posterior odds equals prior odds times likelihood ratio.

For example, if a medical test has a sensitivity of 90% and a specificity of 91%, then the positive Bayes factor is [math]\displaystyle{ \Lambda_+ = P(\text{True Positive})/P(\text{False Positive}) = 90\%/(100\%-91\%)=10 }[/math]. Now, if the prevalence of this disease is 9.09%, and if we take that as the prior probability, then the prior odds is about 1:10. So after receiving a positive test result, the posterior odds of actually having the disease becomes 1:1, which means that the posterior probability of having the disease is 50%. If a second test is performed in serial testing, and that also turns out to be positive, then the posterior odds of actually having the disease becomes 10:1, which means a posterior probability of about 90.91%. The negative Bayes factor can be calculated to be 91%/(100%-90%)=9.1, so if the second test turns out to be negative, then the posterior odds of actually having the disease is 1:9.1, which means a posterior probability of about 9.9%.

The example above can also be understood with more solid numbers: Assume the patient taking the test is from a group of 1000 people, where 91 of them actually have the disease (prevalence of 9.1%). If all these 1000 people take the medical test, 82 of those with the disease will get a true positive result (sensitivity of 90.1%), 9 of those with the disease will get a false negative result (false negative rate of 9.9%), 827 of those without the disease will get a true negative result (specificity of 91.0%), and 82 of those without the disease will get a false positive result (false positive rate of 9.0%). Before taking any test, the patient's odds for having the disease is 91:909. After receiving a positive result, the patient's odds for having the disease is

[math]\displaystyle{ \frac{91}{909}\times\frac{90.1\%}{9.0\%}=\frac{91\times90.1\%}{909\times9.0\%}=1:1 }[/math]

which is consistent with the fact that there are 82 true positives and 82 false positives in the group of 1000 people.

Correspondence to other mathematical frameworks

Propositional logic

Using [math]\displaystyle{ P(\neg B\vert A)=1-P(B\vert A) }[/math] twice, one may use Bayes' theorem to also express [math]\displaystyle{ P(\neg B\vert \neg A) }[/math] in terms of [math]\displaystyle{ P(A\vert B) }[/math] and without negations:

[math]\displaystyle{ P(\neg B\vert \neg A) = 1 - \left(1-P(A\vert B)\right)\frac{P(B)}{P(\neg A)} }[/math],

when [math]\displaystyle{ P(\neg A) = 1 - P(A) \neq 0 }[/math]. From this we can read off the inference

[math]\displaystyle{ P(A\vert B) = 1\implies P(\neg B\vert \neg A) = 1 }[/math].

In words: If certainly [math]\displaystyle{ B }[/math] implies [math]\displaystyle{ A }[/math], we infer that certainly [math]\displaystyle{ \neg A }[/math] implies [math]\displaystyle{ \neg B }[/math]. Where [math]\displaystyle{ P(B)\neq 0 }[/math], the two implications being certain are equivalent statements. In the probability formulas, the conditional probability [math]\displaystyle{ P(A \vert B) }[/math] generalizes the logical implication [math]\displaystyle{ B \implies A }[/math], where now beyond assigning true or false, we assign probability values to statements. The assertion of [math]\displaystyle{ B \implies A }[/math] is captured by certainty of the conditional, the assertion of [math]\displaystyle{ P(A\vert B) = 1 }[/math]. Relating the directions of implication, Bayes' theorem represents a generalization of the contraposition law, which in classical propositional logic can be expressed as:

[math]\displaystyle{ (B \implies A)\iff(\neg A \implies \neg B) }[/math].

In this relation between implications, the positions of [math]\displaystyle{ A }[/math] resp. [math]\displaystyle{ B }[/math] get flipped.

The corresponding formula in terms of probability calculus is Bayes' theorem, which in its expanded form involving the prior probability/base rate [math]\displaystyle{ a }[/math] of only [math]\displaystyle{ A }[/math], is expressed as:[23]

[math]\displaystyle{ P(A \vert B) = P(B \vert A) \frac{a(A)}{P(B\vert A)\,a(A)+P(B \vert \neg A)\,a(\neg A)} }[/math].

Subjective logic

Bayes' theorem represents a special case of deriving inverted conditional opinions in subjective logic expressed as:

[math]\displaystyle{ (\omega^S_{A\tilde{|}B},\omega^S_{A\tilde{|}\lnot B}) = (\omega^S_{B\vert A}, \omega^S_{B\vert\lnot A}) \widetilde{\phi} a_A, }[/math]

where [math]\displaystyle{ \widetilde{\phi} }[/math] denotes the operator for inverting conditional opinions. The argument [math]\displaystyle{ (\omega^S_{B\vert A},\omega^S_{B\vert\lnot A}) }[/math] denotes a pair of binomial conditional opinions given by source [math]\displaystyle{ S }[/math], and the argument [math]\displaystyle{ a_{A} }[/math] denotes the prior probability (aka. the base rate) of [math]\displaystyle{ A }[/math]. The pair of derivative inverted conditional opinions is denoted [math]\displaystyle{ (\omega^S_{A\tilde{|}B},\omega^{S}_{A\tilde{|}\lnot B}) }[/math]. The conditional opinion [math]\displaystyle{ \omega^S_{A\vert B} }[/math] generalizes the probabilistic conditional [math]\displaystyle{ P(A \vert B) }[/math], i.e. in addition to assigning a probability the source [math]\displaystyle{ S }[/math] can assign any subjective opinion to the conditional statement [math]\displaystyle{ (A\vert B) }[/math]. A binomial subjective opinion [math]\displaystyle{ \omega^{S}_{A} }[/math] is the belief in the truth of statement [math]\displaystyle{ A }[/math] with degrees of epistemic uncertainty, as expressed by source [math]\displaystyle{ S }[/math]. Every subjective opinion has a corresponding projected probability [math]\displaystyle{ P(\omega^{S}_{A}) }[/math]. The application of Bayes' theorem to projected probabilities of opinions is a homomorphism, meaning that Bayes' theorem can be expressed in terms of projected probabilities of opinions:

[math]\displaystyle{ P(\omega^S_{A \tilde{|} B}) = \frac{P(\omega^S_{B \vert A}) a(A)}{P(\omega^S_{B\vert A}) a(A) + P(\omega^S_{B \vert \lnot A}) a(\lnot A)}. }[/math]

Hence, the subjective Bayes' theorem represents a generalization of Bayes' theorem.[24]

Generalizations

Bayes theorem for 3 events

A version of Bayes' theorem for 3 events[25] results from the addition of a third event [math]\displaystyle{ C }[/math], with [math]\displaystyle{ P(C)\gt 0, }[/math] on which all probabilities are conditioned:

[math]\displaystyle{ P(A \vert B \cap C) = \frac{P(B \vert A \cap C) \, P(A \vert C)}{P(B \vert C)} }[/math]

Derivation

Using the chain rule

[math]\displaystyle{ P(A \cap B \cap C) = P(A \vert B \cap C) \, P(B \vert C) \, P(C) }[/math]

And, on the other hand

[math]\displaystyle{ P(A \cap B \cap C) = P(B \cap A \cap C) = P(B \vert A \cap C) \, P(A \vert C) \, P(C) }[/math]

The desired result is obtained by identifying both expressions and solving for [math]\displaystyle{ P(A \vert B \cap C) }[/math].

Use in genetics

In genetics, Bayes' rule can be used to estimate the probability of an individual having a specific genotype. Many people seek to approximate their chances of being affected by a genetic disease or their likelihood of being a carrier for a recessive gene of interest. A Bayesian analysis can be done based on family history or genetic testing, in order to predict whether an individual will develop a disease or pass one on to their children. Genetic testing and prediction is a common practice among couples who plan to have children but are concerned that they may both be carriers for a disease, especially within communities with low genetic variance.[26]

The first step in Bayesian analysis for genetics is to propose mutually exclusive hypotheses: for a specific allele, an individual either is or is not a carrier. Next, four probabilities are calculated: Prior Probability (the likelihood of each hypothesis considering information such as family history or predictions based on Mendelian Inheritance), Conditional Probability (of a certain outcome), Joint Probability (product of the first two), and Posterior Probability (a weighted product calculated by dividing the Joint Probability for each hypothesis by the sum of both joint probabilities). This type of analysis can be done based purely on family history of a condition or in concert with genetic testing.[citation needed]

Using pedigree to calculate probabilities

Hypothesis Hypothesis 1: Patient is a carrier Hypothesis 2: Patient is not a carrier
Prior Probability 1/2 1/2
Conditional Probability that all four offspring will be unaffected (1/2) · (1/2) · (1/2) · (1/2) = 1/16 About 1
Joint Probability (1/2) · (1/16) = 1/32 (1/2) · 1 = 1/2
Posterior Probability (1/32) / (1/32 + 1/2) = 1/17 (1/2) / (1/32 + 1/2) = 16/17

Example of a Bayesian analysis table for a female individual's risk for a disease based on the knowledge that the disease is present in her siblings but not in her parents or any of her four children. Based solely on the status of the subject's siblings and parents, she is equally likely to be a carrier as to be a non-carrier (this likelihood is denoted by the Prior Hypothesis). However, the probability that the subject's four sons would all be unaffected is 1/16 (​12·​12·​12·​12) if she is a carrier, about 1 if she is a non-carrier (this is the Conditional Probability). The Joint Probability reconciles these two predictions by multiplying them together. The last line (the Posterior Probability) is calculated by dividing the Joint Probability for each hypothesis by the sum of both joint probabilities.[27]

Using genetic test results

Parental genetic testing can detect around 90% of known disease alleles in parents that can lead to carrier or affected status in their child. Cystic fibrosis is a heritable disease caused by an autosomal recessive mutation on the CFTR gene,[28] located on the q arm of chromosome 7.[29]

Bayesian analysis of a female patient with a family history of cystic fibrosis (CF), who has tested negative for CF, demonstrating how this method was used to determine her risk of having a child born with CF:

Because the patient is unaffected, she is either homozygous for the wild-type allele, or heterozygous. To establish prior probabilities, a Punnett square is used, based on the knowledge that neither parent was affected by the disease but both could have been carriers:

Mother


Father
W

Homozygous for the wild-
type allele (a non-carrier)

M

Heterozygous
(a CF carrier)

W

Homozygous for the wild-
type allele (a non-carrier)

WW MW
M

Heterozygous (a CF carrier)

MW MM

(affected by cystic fibrosis)

Given that the patient is unaffected, there are only three possibilities. Within these three, there are two scenarios in which the patient carries the mutant allele. Thus the prior probabilities are ​23 and ​13.

Next, the patient undergoes genetic testing and tests negative for cystic fibrosis. This test has a 90% detection rate, so the conditional probabilities of a negative test are 1/10 and 1.  Finally, the joint and posterior probabilities are calculated as before.

Hypothesis Hypothesis 1: Patient is a carrier Hypothesis 2: Patient is not a carrier
Prior Probability 2/3 1/3
Conditional Probability of a negative test 1/10 1
Joint Probability 1/15 1/3
Posterior Probability 1/6 5/6

After carrying out the same analysis on the patient's male partner (with a negative test result), the chances of their child being affected is equal to the product of the parents' respective posterior probabilities for being carriers times the chances that two carriers will produce an affected offspring (​14).

Genetic testing done in parallel with other risk factor identification

Bayesian analysis can be done using phenotypic information associated with a genetic condition, and when combined with genetic testing this analysis becomes much more complicated. Cystic fibrosis, for example, can be identified in a fetus through an ultrasound looking for an echogenic bowel, meaning one that appears brighter than normal on a scan. This is not a foolproof test, as an echogenic bowel can be present in a perfectly healthy fetus. Parental genetic testing is very influential in this case, where a phenotypic facet can be overly influential in probability calculation. In the case of a fetus with an echogenic bowel, with a mother who has been tested and is known to be a CF carrier, the posterior probability that the fetus actually has the disease is very high (0.64). However, once the father has tested negative for CF, the posterior probability drops significantly (to 0.16).[27]

Risk factor calculation is a powerful tool in genetic counseling and reproductive planning, but it cannot be treated as the only important factor to consider. As above, incomplete testing can yield falsely high probability of carrier status, and testing can be financially inaccessible or unfeasible when a parent is not present.

See also

Notes

  1. Laplace refined Bayes's 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), p. 15.

References

  1. Joyce, James (2003), Zalta, Edward N., ed., Bayes' Theorem (Spring 2019 ed.), Metaphysics Research Lab, Stanford University, https://plato.stanford.edu/archives/spr2019/entries/bayes-theorem/, retrieved 2020-01-17 
  2. 2.0 2.1 Jeffreys, Harold (1973). Scientific Inference (3rd ed.). Cambridge University Press. p. 31. ISBN 978-0521180788. https://archive.org/details/scientificinfere0000jeff. 
  3. Frame, Paul (2015) (in en). Liberty's Apostle. Wales: University of Wales Press. pp. 44. ISBN 978-1783162161. https://www.uwp.co.uk/book/libertys-apostle-richard-price-his-life-and-times/. Retrieved 23 February 2021. 
  4. Allen, Richard (1999). David Hartley on Human Nature. SUNY Press. pp. 243–244. ISBN 978-0791494516. https://books.google.com/books?id=NCu6HhGlAB8C&pg=PA243. Retrieved 16 June 2013. 
  5. Bayes, Thomas; 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.". Philosophical Transactions of the Royal Society of London 53: 370–418. doi:10.1098/rstl.1763.0053. 
  6. Holland, pp. 46–7.
  7. Price, Richard (1991). Price: Political Writings. Cambridge University Press. p. xxiii. ISBN 978-0521409698. https://books.google.com/books?id=xdH-gjy2vzUC&pg=PR23. Retrieved 16 June 2013. 
  8. Mitchell 1911, p. 314.
  9. Daston, Lorraine (1988). Classical Probability in the Enlightenment. Princeton Univ Press. p. 268. ISBN 0691084971. https://books.google.com/books?id=oq8XNbKyUewC&pg=PA268. 
  10. Stigler, Stephen M. (1986). "Inverse Probability". The History of Statistics: The Measurement of Uncertainty Before 1900. Harvard University Press. pp. 99–138. ISBN 978-0674403413. https://books.google.com/books?id=M7yvkERHIIMC&pg=PA99. 
  11. Stigler, Stephen M. (1983). "Who Discovered Bayes' Theorem?". The American Statistician 37 (4): 290–296. doi:10.1080/00031305.1983.10483122. 
  12. de Vaux, Richard; Velleman, Paul; Bock, David (2016). Stats, Data and Models (4th ed.). Pearson. pp. 380–381. ISBN 978-0321986498. 
  13. 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. 
  14. Hooper, Martyn (2013). "Richard Price, Bayes' theorem, and God". Significance 10 (1): 36–39. doi:10.1111/j.1740-9713.2013.00638.x. 
  15. 15.0 15.1 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-0300188226. https://archive.org/details/theorythatwouldn0000mcgr. 
  16. Stuart, A.; Ord, K. (1994), Kendall's Advanced Theory of Statistics: Volume I – Distribution Theory, Edward Arnold, §8.7 
  17. Kolmogorov, A.N. (1933). Foundations of the Theory of Probability. Chelsea Publishing Company. 
  18. Taraldsen, Gunnar; Tufto, Jarle; Lindqvist, Bo H. (2021-07-24). "Improper priors and improper posteriors" (in en). Scandinavian Journal of Statistics 49 (3): 969–991. doi:10.1111/sjos.12550. ISSN 0303-6898. 
  19. Robert, Christian P.; Casella, George (2004). Monte Carlo Statistical Methods. Springer. ISBN 978-1475741452. OCLC 1159112760. http://worldcat.org/oclc/1159112760. 
  20. Lee, Peter M. (2012). "Chapter 1". Bayesian Statistics. John Wiley & Sons. ISBN 978-1-1183-3257-3. http://www-users.york.ac.uk/~pml1/bayes/book.htm. 
  21. "Bayes' Theorem: Introduction". http://www.trinity.edu/cbrown/bayesweb/. 
  22. "Bayes Theorem - Formula, Statement, Proof | Bayes Rule" (in en). https://www.cuemath.com/data/bayes-theorem/. 
  23. Audun Jøsang, 2016, Subjective Logic; A formalism for Reasoning Under Uncertainty. Springer, Cham, ISBN:978-3-319-42337-1
  24. Audun Jøsang, 2016, Generalising Bayes' Theorem in Subjective Logic. IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems (MFI 2016), Baden-Baden, September 2016
  25. Koller, D.; Friedman, N. (2009). Probabilistic Graphical Models. Massachusetts: MIT Press. pp. 1208. ISBN 978-0-262-01319-2. http://pgm.stanford.edu/. 
  26. Kraft, Stephanie A; Duenas, Devan; Wilfond, Benjamin S; Goddard, Katrina AB (24 September 2018). "The evolving landscape of expanded carrier screening: challenges and opportunities". Genetics in Medicine 21 (4): 790–797. doi:10.1038/s41436-018-0273-4. PMID 30245516. 
  27. 27.0 27.1 Ogino, Shuji; Wilson, Robert B; Gold, Bert; Hawley, Pamela; Grody, Wayne W (October 2004). "Bayesian analysis for cystic fibrosis risks in prenatal and carrier screening". Genetics in Medicine 6 (5): 439–449. doi:10.1097/01.GIM.0000139511.83336.8F. PMID 15371910. 
  28. "Types of CFTR Mutations". Cystic Fibrosis Foundation, www.cff.org/What-is-CF/Genetics/Types-of-CFTR-Mutations/.
  29. "CFTR Gene – Genetics Home Reference". U.S. National Library of Medicine, National Institutes of Health, ghr.nlm.nih.gov/gene/CFTR#location.

Bibliography

Further reading

External links