Berkson's paradox

An example of Berkson's paradox:

• Top: a graph where talent and attractiveness are uncorrelated in the population.
• Bottom: The same graph truncated to only include celebrities (where a person must be both talented and attractive, in some combination, to have become a celebrity). Someone sampling this population may wrongly infer that talent is negatively correlated with attractiveness.

Berkson's paradox, also known as Berkson's bias, collider bias, or Berkson's fallacy, is a result in conditional probability and statistics which is often found to be counterintuitive, and hence a veridical paradox. It is a complicating factor arising in statistical tests of proportions. Specifically, it arises when there is an ascertainment bias inherent in a study design. The effect is related to the explaining away phenomenon in Bayesian networks, and conditioning on a collider in graphical models.

It is often described in the fields of medical statistics or biostatistics, as in the original description of the problem by Joseph Berkson.

Examples

Overview

An illustration of Berkson's Paradox. The top graph represents the actual distribution, in which a positive correlation between quality of burgers and fries is observed. However, an individual who does not eat at any location where both are bad observes only the distribution on the bottom graph, which appears to show a negative correlation.

The most common example of Berkson's paradox is a false observation of a negative correlation between two desirable traits, i.e., that members of a population which have some desirable trait tend to lack a second. Berkson's paradox occurs when this observation appears true when in reality the two properties are unrelated—or even positively correlated—because members of the population where both are absent are not equally observed. For example, a person may observe from their experience that fast food restaurants in their area which serve good hamburgers tend to serve bad fries and vice versa; but because they would likely not eat anywhere where both were bad, they fail to allow for the large number of restaurants in this category which would weaken or even flip the correlation.

Original illustration

Berkson's original illustration involves a retrospective study examining a risk factor for a disease in a statistical sample from a hospital in-patient population. Because samples are taken from a hospital in-patient population, rather than from the general public, this can result in a spurious negative association between the disease and the risk factor. For example, if the risk factor is diabetes and the disease is cholecystitis, a hospital patient without diabetes is more likely to have cholecystitis than a member of the general population, since the patient must have had some non-diabetes (possibly cholecystitis-causing) reason to enter the hospital in the first place. That result will be obtained regardless of whether there is any association between diabetes and cholecystitis in the general population.

Ellenberg example

An example presented by Jordan Ellenberg: Suppose Alex will only date a man if his niceness plus his handsomeness exceeds some threshold. Then nicer men do not have to be as handsome to qualify for Alex's dating pool. So, among the men that Alex dates, Alex may observe that the nicer ones are less handsome on average (and vice versa), even if these traits are uncorrelated in the general population. Note that this does not mean that men in the dating pool compare unfavorably with men in the population. On the contrary, Alex's selection criterion means that Alex has high standards. The average nice man that Alex dates is actually more handsome than the average man in the population (since even among nice men, the ugliest portion of the population is skipped). Berkson's negative correlation is an effect that arises within the dating pool: the rude men that Alex dates must have been even more handsome to qualify.

Quantitative example

As a quantitative example, suppose a collector has 1000 postage stamps, of which 300 are pretty and 100 are rare, with 30 being both pretty and rare. 30% of all his stamps are pretty and 10% of his pretty stamps are rare, so prettiness tells nothing about rarity. He puts the 370 stamps which are pretty or rare on display. Just over 27% of the stamps on display are rare (100/370), but still only 10%(30/300) of the pretty stamps are rare (and 100% of the 70 not-pretty stamps on display are rare). If an observer only considers stamps on display, they will observe a spurious negative relationship between prettiness and rarity as a result of the selection bias (that is, not-prettiness strongly indicates rarity in the display, but not in the total collection).

Statement

Two independent events become conditionally dependent given that at least one of them occurs. Symbolically:

If ${\displaystyle P(A\cap B)=P(A)P(B)}$ and ${\displaystyle P(A\cup B)<1}$ then

${\displaystyle P(A\cap B|A\cup B)<P(A|A\cup B)P(B|A\cup B)}$

Proof: Note that ${\displaystyle P(A|A\cup B)=P(A)/P(A\cup B)}$ and ${\displaystyle P(B|A\cup B)=P(B)/P(A\cup B)}$ which, together with ${\displaystyle P(A\cap B)=P(A)P(B)}$ and ${\displaystyle P(A\cup B)<1}$ (so ${\displaystyle {\frac {1}{P(A\cup B)}}<{\frac {1}{[P(A\cup B)]^{2}}}\ }$) implies that

${\displaystyle {\begin{aligned}P(A\cap B|A\cup B)={\frac {P(A\cap B)}{P(A\cup B)}}={\frac {P(A)P(B)}{P(A\cup B)}}<{\frac {P(A)P(B)}{[P(A\cup B)]^{2}}}=P(A|A\cup B)P(B|A\cup B).\end{aligned}}}$

One can see this in tabular form as follows: the yellow regions are the outcomes where at least one event occurs (and ~A means "not A").

	A	~A
B	A & B	~A & B
~B	A & ~B	~A & ~B

For instance, if one has a sample of ${\displaystyle 100}$, and both ${\displaystyle A}$ and ${\displaystyle B}$ occur independently half the time ( ${\displaystyle P(A)=P(B)=1/2}$ ), one obtains:

	A	~A
B	25	25
~B	25	25

So in ${\displaystyle 75}$ outcomes, either ${\displaystyle A}$ or ${\displaystyle B}$ occurs, of which ${\displaystyle 50}$ have ${\displaystyle A}$ occurring. By comparing the conditional probability of ${\displaystyle A}$ to the unconditional probability of ${\displaystyle A}$:

${\displaystyle P(A|A\cup B)=50/75=2/3>P(A)=50/100=1/2}$

We see that the probability of ${\displaystyle A}$ is higher (${\displaystyle 2/3}$) in the subset of outcomes where (${\displaystyle A}$or${\displaystyle B}$) occurs, than in the overall population (${\displaystyle 1/2}$). On the other hand, the probability of ${\displaystyle A}$ given both ${\displaystyle B}$ and (${\displaystyle A}$ or ${\displaystyle B}$) is simply the unconditional probability of ${\displaystyle A}$, ${\displaystyle P(A)}$, since ${\displaystyle A}$ is independent of ${\displaystyle B}$. In the numerical example, we have conditioned on being in the top row:

	A	~A
B	25	25
~B	25	25

Here the probability of ${\displaystyle A}$ is ${\displaystyle 25/50=1/2}$.

Berkson's paradox arises because the conditional probability of ${\displaystyle A}$ given ${\displaystyle B}$within the three-cell subset equals the conditional probability in the overall population, but the unconditional probability within the subset is inflated relative to the unconditional probability in the overall population, hence, within the subset, the presence of ${\displaystyle B}$ decreases the conditional probability of ${\displaystyle A}$ (back to its overall unconditional probability):

${\displaystyle P(A|B,A\cup B)=P(A|B)=P(A)}$

${\displaystyle P(A|A\cup B)>P(A)}$

Because the effect of conditioning on ${\displaystyle (A\cup B)}$ derives from the relative size of ${\displaystyle P(A|A\cup B)}$ and ${\displaystyle P(A)}$ the effect is particularly large when ${\displaystyle A}$ is rare (${\displaystyle P(A)<<1}$) but very strongly correlated to ${\displaystyle B}$ (${\displaystyle P(A|B)\approx 1}$). For example, consider the case below where N is very large:

	A	~A
B	1	0
~B	0	N

For the case without conditioning on ${\displaystyle (A\cup B)}$ we have

${\displaystyle P(A)=1/(N+1)}$

${\displaystyle P(A|B)=1}$

So A occurs rarely, unless B is present, when A occurs always. Thus B is dramatically increasing the likelihood of A.

For the case with conditioning on ${\displaystyle (A\cup B)}$ we have

${\displaystyle P(A|A\cup B)=1}$

${\displaystyle P(A|B,A\cup B)=P(A|B)=1}$

Now A occurs always, whether B is present or not. So B has no impact on the likelihood of A. Thus we see that for highly correlated data a huge positive correlation of B on A can be effectively removed when one conditions on ${\displaystyle (A\cup B)}$.

See also

• Simpson's paradox
• Survivorship bias

References

• Berkson, Joseph (June 1946). "Limitations of the Application of Fourfold Table Analysis to Hospital Data". Biometrics Bulletin . 2 (3): 47–53. doi:10.2307/3002000. JSTOR   3002000. PMID   21001024. (The paper is frequently miscited as Berkson, J. (1949) Biological Bulletin 2, 47–53.)
• Jordan Ellenberg, "Why are handsome men such jerks?"