Tutorials › High-Dimensional Genomic Data Analysis › Interpreting Mendelian Randomization With Binary Outcomes

Mendelian Randomization · Tutorial 693 of 795

Interpreting Mendelian Randomization With Binary Outcomes

Interpret binary-outcome MR estimates on the log-odds scale while keeping their assumptions, non-collapsibility, and causal meaning distinct.

Advanced 10 min read

What You'll Learn

  • Calculate an MR ratio estimate from SNP–exposure and SNP–outcome associations.
  • Convert a log-odds estimate to an odds ratio and, with a baseline risk, to an illustrative risk.
  • State the assumptions needed to interpret the ratio as a causal log-odds effect.
  • Explain why conditional and marginal odds ratios can differ even without confounding.
  • Distinguish a genetically induced liability contrast from a direct intervention effect.
  • Identify common interpretation errors and assess when an odds-ratio interpretation is not justified.

What a Binary-Outcome MR Ratio Estimates

When the outcome is binary, its association with a genetic variant is commonly reported as a log odds ratio. A Mendelian randomization (MR) ratio divides that association by the variant’s association with the exposure. The result is on a log-odds scale per unit of exposure, provided that the exposure association and outcome association have compatible definitions and scales. Exponentiating the ratio gives an odds ratio, not a risk ratio.

That calculation alone does not establish a causal interpretation. As in the earlier material on pleiotropy and interpretation in multivariable MR, causal interpretation depends on the instrumental-variable assumptions and on what causal contrast the model represents. With a binary outcome, it is especially important to distinguish the odds ratio’s mathematical scale from a risk change, and a genetic liability effect from the effect of assigning an intervention.

Definition: For variant \(j\), let \(\hat{\beta}_{Xj}\) be its estimated association with the exposure, and let \(\hat{\beta}_{Yj}\) be its estimated association with the binary outcome on the log-odds scale. The single-variant ratio estimate is \(\hat{\theta}_j=\hat{\beta}_{Yj}/\hat{\beta}_{Xj}\). If the exposure is measured in standard-deviation units, \(\hat{\theta}_j\) is a log-odds contrast per one standard-deviation increase in exposure under the model’s assumptions. Its exponent, \(\exp(\hat{\theta}_j)\), is the corresponding odds ratio.

A common working relation for the genetic associations is \(\beta_{Yj}=\theta\beta_{Xj}+\alpha_j\), where \(\alpha_j\) represents a variant’s outcome association not operating through the exposure. The ratio has the intended causal interpretation when the instrument is relevant, independent of causes of the outcome (subject to the study design and adjustment), and has no pathway to the outcome outside the exposure, so that \(\alpha_j=0\). In addition, the exposure and outcome associations must be suitably aligned, and the causal effect must be representable by the chosen model and scale. For multiple variants, the same scale and interpretation issues apply to the combined estimate; this tutorial focuses on interpreting the ratio rather than choosing an estimator.

From Log Odds to Odds and Risk

For an odds ratio \(R\), the log-odds effect is \(\theta=\log(R)\), so \(R=\exp(\theta)\). An odds ratio compares odds, \(p/(1-p)\), rather than probabilities \(p\). To translate an odds ratio into a probability at a specified baseline risk \(p_0\), multiply the baseline odds by the odds ratio and convert the resulting odds back to risk:

$$ p_1 = \frac{R\left[p_0/(1-p_0)\right]} {1+R\left[p_0/(1-p_0)\right]}. $$

The resulting probability is an illustration for the specified baseline risk and a contrast that supports this calculation. It is not a universal risk change: the same odds ratio yields different absolute risk differences at different baseline risks. This conversion also does not, by itself, establish that a genetic contrast corresponds to an intervention that could be assigned.

Worked Example: Converting a Ratio Estimate to an Odds Ratio

Suppose a hypothetical analysis estimates a variant’s association with an exposure at \(\hat{\beta}_{X}=0.20\) standard deviations per effect allele, and its association with a binary outcome at \(\hat{\beta}_{Y}=0.08\) log odds per allele. Assume the estimates are aligned to the same allele and the ratio-model assumptions are appropriate.

The ratio estimate is \(\hat{\theta}=0.08/0.20=0.40\) log odds per standard-deviation increase in exposure. Exponentiating gives \(\exp(0.40)=1.492\). Under the assumed causal model, this is an estimated odds ratio of 1.492 per standard-deviation increase in exposure: the odds are estimated to be about 49.2% higher, not the risk.

For illustration, take a baseline risk of 10%. The baseline odds are \(0.10/0.90=0.1111\). Multiplying by 1.492 gives odds \(0.1658\), corresponding to risk \(0.1658/(1+0.1658)=0.1422\), or about 14.2%. Thus the illustrative risk difference is about 4.2 percentage points at this baseline risk. It would not generally be 4.2 percentage points at another baseline risk.

Uncertainty on the Odds-Ratio Scale

A ratio estimate is uncertain because both association estimates are uncertain. One approximate standard error can be obtained by the delta method when the estimates are approximately normal, the denominator is sufficiently far from zero, and the two estimates are independent. If their standard errors are \(SE_Y\) and \(SE_X\), the approximation is \(SE(\hat{\theta})\approx\sqrt{SE_Y^2/\hat{\beta}_X^2+\hat{\beta}_Y^2SE_X^2/\hat{\beta}_X^4}\). With non-independent estimates, the covariance term must also be included; the independence assumption should not be used automatically.

Worked Example: A Confidence Interval for the Odds Ratio

Continue with the hypothetical associations \(\hat{\beta}_{Y}=0.08\) and \(\hat{\beta}_{X}=0.20\). Suppose \(SE_Y=0.02\), \(SE_X=0.01\), and the estimates are independent. The ratio standard error is approximately \[ \sqrt{\frac{0.02^2}{0.20^2}+ \frac{0.08^2(0.01^2)}{0.20^4}} = \sqrt{0.0100+0.0004} =0.1020. \]

A normal-approximation 95% confidence interval for the log-odds effect is \(0.40\pm1.96(0.1020)=(0.2001,0.5999)\). Exponentiating the endpoints gives an odds-ratio interval of approximately \((1.222,1.822)\). Under the model and assumptions, this interval is a range of values compatible with the data and method for the odds ratio per standard-deviation increase in exposure. It is not a 95% probability statement about a fixed parameter.

This approximation can be unreliable when the SNP–exposure association is weak: dividing by an estimate near zero can produce a skewed or unbounded ratio distribution, which a symmetric normal interval does not represent well. A reassuring-looking interval from this approximation does not repair a weak denominator.

Why Odds Ratios Are Non-Collapsible

An odds ratio is non-collapsible: a conditional odds ratio and a marginal odds ratio can differ even when the conditioning variable is not a confounder. Conditional odds compare outcome odds between exposure groups within levels of another variable. Marginal odds compare odds after averaging risks across that variable’s distribution. Averaging probabilities and then taking odds is not equivalent to taking odds within each group and averaging those odds.

The following hypothetical calculation demonstrates the distinction. Let \(Z\) define two equally common strata. In each stratum, suppose exposure doubles the outcome odds. In stratum A, the unexposed risk is 0.10, so the exposed risk is \(2/11\). In stratum B, the unexposed risk is 0.50, so the exposed risk is \(2/3\). Each stratum-specific odds ratio is exactly 2, but the risks averaged across the equally sized strata give a different marginal odds ratio.

Worked Example: Conditional and Marginal Odds Ratios

For stratum A, unexposed odds are \(0.10/0.90=1/9\). Doubling them gives \(2/9\), which corresponds to exposed risk \((2/9)/(1+2/9)=2/11\). For stratum B, unexposed odds are \(0.50/0.50=1\); doubling gives odds 2 and exposed risk \(2/(1+2)=2/3\). Thus the conditional odds ratio is 2 in both strata.

With equal stratum proportions, the marginal unexposed risk is \((0.10+0.50)/2=0.30\). The marginal exposed risk, retaining exact values, is \[ \frac{(2/11)+(2/3)}{2}=\frac{14}{33}. \] The marginal odds ratio is therefore \[ \frac{(14/33)/(19/33)}{0.30/0.70} = \frac{14/19}{3/7} = \frac{98}{57} \approx 1.719. \]

The difference between 2 and 1.719 does not demonstrate confounding in this example: the stratum-specific effect was identical, and the strata were balanced between exposure groups by construction. It demonstrates that conditional and marginal odds ratios are different summaries. Therefore, a marginal genetic association on the log-odds scale cannot automatically be substituted for a conditional causal log-odds effect. The model and target contrast must make the scales compatible.

Genetic Liability Is Not Automatically an Intervention

A genetic variant generally shifts exposure over a lifetime, often by a small amount. A ratio estimate consequently describes a contrast associated with genetically induced differences in exposure under the MR assumptions. It does not automatically describe the effect of starting a drug, changing a behaviour at a particular age, or assigning an exposure value in a clinical trial. Such an intervention may differ in timing, duration, magnitude, and affected biological pathways.

The distinction is especially important when the exposure itself is binary or when its genetic association represents a change in liability to being exposed. A ratio involving a genetic association with an exposure does not turn a liability contrast into the effect of assigning everyone from unexposed to exposed. The estimand depends on how exposure is defined and modeled, the size and nature of the genetic shift, and the assumptions linking that shift to the outcome. Interpret the estimate using the actual exposure scale, not as a generic treatment effect.

Key interpretation: A ratio of a SNP–outcome log-odds association to a SNP–exposure association is a candidate causal log-odds effect per exposure unit under stated MR and scale assumptions. Its exponential is an odds ratio. Neither quantity alone is a risk difference, a risk ratio, or proof of the effect of a direct intervention.

Diagnostics, Alternatives, and Common Mistakes

Before interpreting a binary-outcome ratio, check that the reported outcome association is a log odds ratio, the exposure association has a clear unit, and both estimates refer to the same allele and compatible population. Verify that the exposure association is not close to zero and that standard errors and covariance assumptions suit the data sources. Examine whether the causal model supports a constant effect on the chosen scale; if effects vary, one ratio may not summarize all relevant contrasts.

The odds ratio is only one effect measure. A risk ratio compares probabilities, and a risk difference subtracts them. These measures answer different questions and do not generally equal a transformed odds ratio without a baseline risk and additional assumptions. For common outcomes, the odds ratio can be farther from 1 than the risk ratio for the same comparison. When absolute clinical impact is the target, reporting an odds ratio alone may be insufficient; a risk translation requires a relevant baseline risk and a defensible target population.

  • Mistake: “An odds ratio of 1.5 means 50% more people will have the outcome.” It means 50% higher odds. The probability change depends on baseline risk.
  • Mistake: “A causal MR ratio is automatically the effect of treatment.” The estimate pertains to a genetically induced exposure contrast under its assumptions; an intervention effect may be a different estimand.
  • Mistake: “Different marginal and conditional odds ratios prove confounding.” They can differ through non-collapsibility alone. Confounding and non-collapsibility are distinct issues.
  • Mistake: “Exponentiating a ratio makes its causal interpretation valid.” Exponentiation changes the scale, not the assumptions. Instrument validity, compatible association scales, and an appropriate causal model remain necessary.

A practical interpretation should name the outcome, exposure unit, scale, and causal caveat together. For example: “Under the MR assumptions and a compatible log-odds model, a one-standard-deviation genetically induced increase in the exposure is associated with an estimated odds ratio of 1.49 for the binary outcome.” Avoid replacing “genetically induced” with a direct intervention claim unless the design and estimand justify that step.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A ratio estimate is \(-0.20\) log odds per exposure unit. What is the corresponding odds ratio, and how should its direction be described?
  2. At baseline risk 20%, calculate the risk implied by an odds ratio of 2. Use odds, not a direct multiplication of risk.
  3. In the two-stratum example, why does a marginal odds ratio of about 1.719 not contradict the stratum-specific odds ratios of 2?
  4. Give one reason a genetically induced exposure effect may differ from the effect of assigning an intervention.
  5. What makes a delta-method ratio confidence interval questionable when the exposure association is weak?

Solutions. (1) \(\exp(-0.20)=0.8187\). Under the model, the odds are about 18.1% lower per exposure unit; this is not an 18.1% lower risk. (2) Baseline odds are \(0.20/0.80=0.25\). Doubled odds are 0.50, giving risk \(0.50/(1+0.50)=0.3333\), or 33.3%. (3) Odds ratios are non-collapsible: averaging stratum-specific risks and then forming odds need not preserve the conditional odds ratio, even without confounding. (4) Genetic effects may operate over a different time span, magnitude, or pathway than a treatment assigned at a particular time. (5) Division by a denominator near zero makes the ratio distribution unstable and potentially asymmetric, so a symmetric normal approximation may misrepresent its uncertainty.

Clinical Biostats