Odds ratio

The odds ratio is a measure of effect size, describing the strength of association or non-independence between two binary data values. It is used as a descriptive statistic, and plays an important role in logistic regression. Unlike other measures of association for paired binary data such as the relative risk, the odds ratio treats the two variables being compared symmetrically, and can be estimated using some types of non-random samples.

Definition

Definition in terms of group-wise odds

The odds ratio is the ratio of the odds of an event occurring in one group to the odds of it occurring in another group. The term is also used to refer to sample-based estimates of this ratio. These groups might be men and women, an experimental group and a control group, or any other dichotomous classification. If the probabilities of the event in each of the groups are p₁ (first group) and p₂ (second group), then the odds ratio is:

:$\{\; p\_1/(1-p\_1)\; \backslash over\; p\_2/(1-p\_2)\}=\{\; p\_1/q\_1\; \backslash over\; p\_2/q\_2\}=\backslash frac\{\backslash ;p\_1q\_2\backslash ;\}\{\backslash ;p\_2q\_1\backslash ;\},$

where q_x = 1 − p_x. An odds ratio of 1 indicates that the condition or event under study is equally likely to occur in both groups. An odds ratio greater than 1 indicates that the condition or event is more likely to occur in the first group. And an odds ratio less than 1 indicates that the condition or event is less likely to occur in the first group. The odds ratio must be greater than or equal to zero if it is defined. It is undefined if p₂q₁ equals zero.

Definition in terms of joint and conditional probabilities

The odds ratio can also be defined in terms of the joint probability distribution of two binary random variables. The joint distribution of binary random variables X and Y can be written

{| class="wikitable" style="text-align: center; background: #FFFFFF;" |----- | || Y = 1 || Y = 0 |----- | X = 1 || $p\_\{11\}$ || $p\_\{10\}$ |----- | X = 0 || $p\_\{01\}$ || $p\_\{00\}$ |}

where p₁₁, p₁₀, p₀₁ and p₀₀ are non-negative "cell probabilities" that sum to one. The odds for Y within the two subpopulations defined by X = 1 and X = 0 are defined in terms of the conditional probabilities given X:

{| class="wikitable" style="text-align: center; background: #FFFFFF;" |----- | || Y = 1 || Y = 0 |----- | X = 1 || $p\_\{11\}/(p\_\{11\}+p\_\{10\})$ || $p\_\{10\}/(p\_\{11\}+p\_\{10\})$ |----- | X = 0 || $p\_\{01\}/(p\_\{01\}+p\_\{00\})$ || $p\_\{00\}/(p\_\{01\}+p\_\{00\})$ |}

Thus the odds ratio is

:$\{\; \backslash dfrac\{p\_\{11\}/(p\_\{11\}+p\_\{10\})\}\{p\_\{10\}/(p\_\{11\}+p\_\{10\})\}\; \backslash bigg\; /\; \backslash dfrac\{p\_\{01\}/(p\_\{01\}+p\_\{00\})\}\{p\_\{00\}/(p\_\{01\}+p\_\{00\})\}\}\; =\; \backslash dfrac\{p\_\{11\}p\_\{00\}\}\{p\_\{10\}p\_\{01\}\}.$

The simple expression on the right, above, is easy to remember as the product of the probabilities of the "concordant cells" (X = Y) divided by the product of the probabilities of the "discordant cells" (X ≠ Y). However note that in some applications the labeling of categories as zero and one is arbitrary, so there is nothing special about concordant versus discordant values in these applications.

Symmetry

If we had calculated the odds ratio based on the conditional probabilities given Y,

{| class="wikitable" style="text-align: center; background: #FFFFFF;" |----- | || Y = 1 || Y = 0 |----- | X = 1 || $p\_\{11\}/(p\_\{11\}+p\_\{01\})$ || $p\_\{10\}/(p\_\{10\}+p\_\{00\})$ |----- | X = 0 || $p\_\{01\}/(p\_\{11\}+p\_\{01\})$ || $p\_\{00\}/(p\_\{10\}+p\_\{00\})$ |}

we would have gotten the same result

:$\{\; \backslash dfrac\{p\_\{11\}/(p\_\{11\}+p\_\{01\})\}\{p\_\{01\}/(p\_\{11\}+p\_\{01\})\}\; \backslash bigg\; /\; \backslash dfrac\{p\_\{10\}/(p\_\{10\}+p\_\{00\})\}\{p\_\{00\}/(p\_\{10\}+p\_\{00\})\}\}\; =\; \backslash dfrac\{p\_\{11\}p\_\{00\}\}\{p\_\{10\}p\_\{01\}\}.$

Other measures of effect size for binary data such as the relative risk do not have this symmetry property.

Relation to statistical independence

If X and Y are independent, their joint probabilities can be expressed in terms of their marginal probabilities p_x =  P(X = 1) and p_y =  P(Y = 1), as follows

{| class="wikitable" style="text-align: center; background: #FFFFFF;" |----- | || Y = 1 || Y = 0 |----- | X = 1 || $p\_xp\_y$ || $p\_x(1-p\_y)$ |----- | X = 0 || $(1-p\_x)p\_y$ || $(1-p\_x)(1-p\_y)$ |}

In this case, the odds ratio equals one, and conversely the odds ratio can only equal one if the joint probabilities can be factored in this way. Thus the odds ratio equals one if and only if X and Y are independent.

Recovering the cell probabilities from the odds ratio and marginal probabilities

The odds ratio is a function of the cell probabilities, and conversely, the cell probabilities can be recovered given knowledge of the odds ratio and the marginal probabilities P(X = 1) = p₁₁ + p₁₀ and P(Y = 1) = p₁₁ + p₀₁. If the odds ratio R differs from 1, then

:$p\_\{11\}\; =\; \backslash frac\{1\; +\; (p\_\{1\backslash cdot\}+p\_\{\backslash cdot\; 1\})(R-1)\; -\; S\}\{2(R-1)\}$

where p_1• = p₁₁ + p₁₀,  p_•1 = p₁₁ + p₀₁, and

:$S\; =\; \backslash sqrt\{(1+(p\_\{1\backslash cdot\}+p\_\{\backslash cdot\; 1\})(R-1))^2\; +\; 4R(1-R)p\_\{1\backslash cdot\}p\_\{\backslash cdot\; 1\}\}.$

In the case where R = 1, we have independence, so p₁₁ = p_1•p_•1.

Once we have p₁₁, the other three cell probabilities can easily be recovered from the marginal probabilities.

Example

Suppose that in a sample of 100 men, 90 have drunk wine in the previous week, while in a sample of 100 women only 20 have drunk wine in the same period. The odds of a man drinking wine are 90 to 10, or 9:1, while the odds of a woman drinking wine are only 20 to 80, or 1:4 = 0.25:1. The odds ratio is thus 9/0.25, or 36, showing that men are much more likely to drink wine than women. The detailed calculation is:

:$\{\; 0.9/0.1\; \backslash over\; 0.2/0.8\}=\backslash frac\{\backslash ;0.9\backslash times\; 0.8\backslash ;\}\{\backslash ;0.1\backslash times\; 0.2\backslash ;\}\; =\{0.72\; \backslash over\; 0.02\}\; =\; 36.$

This example also shows how odds ratios are sometimes sensitive in stating relative positions: in this sample men are 90/20 = 4.5 times more likely to have drunk wine than women, but have 36 times the odds. The logarithm of the odds ratio, the difference of the logits of the probabilities, tempers this effect, and also makes the measure symmetric with respect to the ordering of groups. For example, using natural logarithms, an odds ratio of 36/1 maps to 3.584, and an odds ratio of 1/36 maps to −3.584.

Statistical inference

Several approaches to statistical inference for odds ratios have been developed.

One approach to inference uses large sample approximations to the sampling distribution of the log odds ratio (the natural logarithm of the odds ratio). If we use the joint probability notation defined above, the population log odds ratio is

:$\{\backslash log\backslash left(\backslash frac\{p\_\{11\}p\_\{00\}\}\{p\_\{01\}p\_\{10\}\}\backslash right)\; =\; \backslash log(p\_\{11\})\; +\; \backslash log(p\_\{00\}\backslash big)\; -\; \backslash log(p\_\{10\})\; -\; \backslash log(p\_\{01\})\}.\backslash ,$

If we observe data in the form of a contingency table

{| class="wikitable" style="text-align: center; background: #FFFFFF;" |----- | || Y = 1 || Y = 0 |----- | X = 1 || $n\_\{11\}$ || $n\_\{10\}$ |----- | X = 0 || $n\_\{01\}$ || $n\_\{00\}$ |}

then the probabilities in the joint distribution can be estimated as

{| class="wikitable" style="text-align: center; background: #FFFFFF;" |----- | || Y = 1 || Y = 0 |----- | X = 1 || $\backslash hat\{p\}\_\{11\}$ || $\backslash hat\{p\}\_\{10\}$ |----- | X = 0 || $\backslash hat\{p\}\_\{01\}$ || $\backslash hat\{p\}\_\{00\}$ |}

where p̂ = n_ij / n, with n = n₁₁ + n₁₀ + n₀₁ + n₀₀ being the sum of all four cell counts. The sample log odds ratio is

:$\{L=\backslash log\backslash left(\backslash dfrac\{\backslash hat\{p\}\_\{11\}\backslash hat\{p\}\_\{00\}\}\{\backslash hat\{p\}\_\{10\}\backslash hat\{p\}\_\{01\}\}\backslash right)\; =\; \backslash log\backslash left(\backslash dfrac\{n\_\{11\}n\_\{00\}\}\{n\_\{10\}n\_\{01\}\}\backslash right)\}$.

The distribution of the log odds ratio is approximately normal with:

$X\backslash \; \backslash sim\backslash \; \backslash mathcal\{N\}(\backslash log\; (OR),\backslash ,\backslash sigma^2).\; \backslash ,$

The standard error for the log odds ratio is approximately

:$+\; \backslash dfrac\{1\}\{n\_\{10\}\}\; +\; \backslash dfrac\{1\}\{n\_\{01\}\}\; +\; \backslash dfrac\{1\}\{n\_\{00\}\}\}\}$.

This is an asymptotic approximation, and will not give a meaningful result if any of the cell counts are very small. If L is the sample log odds ratio, an approximate 95% confidence interval for the population log odds ratio is L ± 1.96SE. This can be mapped to exp(L − 1.96SE), exp(L + 1.96SE) to obtain a 95% confidence interval for the odds ratio. If we wish to test the hypothesis that the population odds ratio equals one, the two-sided p-value is 2P(Z< −|L|/SE), where P denotes a probability, and Z denotes a standard normal random variable.

An alternative approach to inference for odds ratios looks at the distribution of the data conditionally on the marginal frequencies of X and Y. An advantage of this approach is that the sampling distribution of the odds ratio can be expressed exactly.

Role in logistic regression

Logistic regression is one way to generalize the odds ratio beyond two binary variables. Suppose we have a binary response variable Y and a binary predictor variable X, and in addition we have other predictor variables Z₁, ..., Z_p that may or may not be binary. If we use multiple logistic regression to regress Y on X, Z₁, ..., Z_p, then the estimated coefficient $\backslash hat\{\backslash beta\}\_x$ for X is related to a conditional odds ratio. Specifically, at the population level

:$\backslash exp(\backslash beta\_x)\; =\; \backslash frac\{P(Y=1|X=1,\; Z\_1,\; \backslash ldots,\; Z\_p)/P(Y=0|X=1,\; Z\_1,\; \backslash ldots,\; Z\_p)\}\{P(Y=1|X=0,\; Z\_1,\; \backslash ldots,\; Z\_p)/P(Y=0|X=0,\; Z\_1,\; \backslash ldots,\; Z\_p)\},$

so $\backslash exp(\backslash hat\{\backslash beta\}\_x)$ is an estimate of this conditional odds ratio. The interpretation of $\backslash exp(\backslash hat\{\backslash beta\}\_x)$ is as an estimate of the odds ratio between Y and X when the values of Z₁, ..., Z_p are held fixed.

Insensitivity to the type of sampling

If the data form a "population sample", then the cell probabilities p̂_ij are interpreted as the frequencies of each of the four groups in the population as defined by their X and Y values. In many settings it is impractical to obtain a population sample, so a selected sample is used. For example, we may choose to sample units with X = 1 with a given probability f, regardless of their frequency in the population (which would necessitate sampling units with X = 0 with probability 1 − f). In this situation, our data would follow the following joint probabilities:

{| class="wikitable" style="text-align: center; background: #FFFFFF;" |----- | || Y = 1 || Y = 0 |----- | X = 1 || $fp\_\{11\}/(p\_\{11\}+p\_\{10\})$ || $fp\_\{10\}(p\_\{11\}+p\_\{10\})$ |----- | X = 0 || $(1-f)p\_\{01\}/(p\_\{01\}+p\_\{00\})$ || $(1-f)p\_\{00\}/(p\_\{01\}+p\_\{00\})$ |}

The odds ratio p₁₁p₀₀ / p₀₁p₁₀ for this distribution does not depend on the value of f. This shows that the odds ratio (and consequently the log odds ratio) is invariant to non-random sampling based on one of the variables being studied. Note however that the standard error of the log odds ratio does depend on the value of f. This fact is exploited in two important situations:

* Suppose it is inconvenient or impractical to obtain a population sample, but it is practical to obtain a convenience sample of units with different X values, such that within the X = 0 and X = 1 subsamples the Y values are representative of the population (i.e. they follow the correct conditional probabilities).

* Suppose the marginal distribution of one variable, say X, is very skewed. For example, if we are studying the relationship between high alcohol consumption and pancreatic cancer in the general population, the incidence of pancreatic cancer would be very low, so it would require a very large population sample to get a modest number of pancreatic cancer cases. However we could use data from hospitals to contact most or all of their pancreatic cancer patients, and then randomly sample an equal number of subjects without pancreatic cancer (this is called a "case-control study").

In both these settings, the odds ratio can be calculated from the selected sample, without biasing the results relative to what would have been obtained for a population sample.

Use in quantitative research

Due to the widespread use of logistic regression, the odds ratio is widely used in many fields of medical and social science research. The odds ratio is commonly used in survey research, in epidemiology, and to express the results of some clinical trials, such as in case-control studies. It is often abbreviated "OR" in reports. When data from multiple surveys is combined, it will often be expressed as "pooled OR".

Relation to relative risk

In clinical studies, as well as in some other settings, the parameter of greatest interest is often the relative risk rather than the odds ratio. The relative risk is best estimated using a population sample, but if the rare disease assumption holds, the odds ratio is a good approximation to the relative risk — the odds is p / (1 − p), so when p moves towards zero, 1 − p moves towards 1, meaning that the odds approaches the risk, and the odds ratio approaches the relative risk. When the rare disease assumption does not hold, the odds ratio can overestimate the relative risk.

If the absolute risk in the control group is available, conversion between the two is calculated by:. Another alternative estimator is the Mantel-Haenszel estimator.

Numerical examples

The following four contingency tables contain observed cell counts, along with the corresponding sample odds ratio (OR) and sample log odds ratio (LOR):

{| cellpadding="5" cellspacing="0" align="center" |- ! rowspan=2 | ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;" colspan=2 | OR = 1, LOR = 0 ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;" colspan=2 | OR = 1, LOR = 0 ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;" colspan=2 | OR = 4, LOR = 1.39 ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;border-right:1px solid black;" colspan=2 | OR = 0.25, LOR = −1.39 |- ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;" | Y = 0 ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;" | Y = 0 ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;" | Y = 0 ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;border-right:1px solid black;" | Y = 0 |- ! style="background:#ffdead;border-left:1px solid black;border-top:1px solid black;" | X = 1 ! style="border-left:1px solid black;" | 10 ! 10 ! style="border-left:1px solid black;" | 100 ! 100 ! style="border-left:1px solid black;" | 20 ! 10 ! style="border-left:1px solid black;" | 10 ! style="border-right:1px solid black;" |20 |- ! style="background:#ffdead;border-bottom:1px solid black;border-left:1px solid black;" | X = 0 ! style="border-bottom:1px solid black;border-left:1px solid black;" | 5 ! style="border-bottom:1px solid black;" | 5 ! style="border-left:1px solid black;border-bottom:1px solid black;" | 50 ! style="border-bottom:1px solid black;" | 50 ! style="border-left:1px solid black;border-bottom:1px solid black;" | 10 ! style="border-bottom:1px solid black;" | 20 ! style="border-left:1px solid black;border-bottom:1px solid black;" | 20 ! style="border-right:1px solid black;border-bottom:1px solid black;" | 10 |}

The following joint probability distributions contain the population cell probabilities, along with the corresponding population odds ratio (OR) and population log odds ratio (LOR):

{| cellpadding="5" cellspacing="0" align="center" |- ! rowspan=2 | ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;" colspan=2 | OR = 1, LOR = 0 ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;" colspan=2 | OR = 1, LOR = 0 ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;" colspan=2 | OR = 16, LOR = 2.77 ! style="background:#efefef;border-left:1px solid black;border-top:1px solid black;border-right:1px solid black;" colspan=2 | OR = 0.67, LOR = −0.41 |- ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;" | Y = 0 ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;" | Y = 0 ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;" | Y = 0 ! style="background:#ffdead;border-left:1px solid black;" | Y = 1 ! style="background:#ffdead;border-right:1px solid black;" | Y = 0 |- ! style="background:#ffdead;border-left:1px solid black;border-top:1px solid black;" | X = 1 ! style="border-left:1px solid black;" | 0.2 ! 0.2 ! style="border-left:1px solid black;" | 0.4 ! 0.4 ! style="border-left:1px solid black;" | 0.4 ! 0.1 ! style="border-left:1px solid black;" | 0.1 ! style="border-right:1px solid black;" | 0.3 |- ! style="background:#ffdead;border-bottom:1px solid black;border-left:1px solid black;" | X = 0 ! style="border-bottom:1px solid black;border-left:1px solid black;" | 0.3 ! style="border-bottom:1px solid black;" | 0.3 ! style="border-left:1px solid black;border-bottom:1px solid black;" | 0.1 ! style="border-bottom:1px solid black;" | 0.1 ! style="border-left:1px solid black;border-bottom:1px solid black;" | 0.1 ! style="border-bottom:1px solid black;" | 0.4 ! style="border-left:1px solid black;border-bottom:1px solid black;" | 0.2 ! style="border-right:1px solid black;border-bottom:1px solid black;" | 0.4 |}

Worked example

See also

Hazard ratio

External links

Odds Ratio Calculator — website

Odds ratio definition and examples

OpenEpi, a web-based program that calculates the odds ratio, both unmatched and pair-matched

References

Category:Epidemiology Category:Medical statistics Category:Statistical terminology Category:Bayesian statistics