Tutorials › Statistical Genetics › Causal Estimands For Mendelian Randomization

Mendelian Randomization Foundations · Tutorial 834 of 989

Causal Estimands For Mendelian Randomization

Specify an MR causal estimand by stating the exposure contrast, outcome scale, and population to which the effect refers.

Advanced 11 min read

What You'll Learn

  • Define a causal estimand using potential outcomes and a stated target population
  • Distinguish effects of an exposure from associations of a genetic variant with exposure or outcome
  • Specify the exposure contrast and units for a continuous exposure
  • Compare risk difference, risk ratio, and odds ratio estimands for a binary outcome
  • Explain how changing the target population can change an average causal effect
  • Recognize when a reported effect scale does not match the scientific question

Why the Estimand Must Come First

Mendelian randomization (MR) uses genetic variation to investigate whether an exposure has a causal effect on an outcome. Before asking how an estimate is calculated, specify which causal effect is being targeted. “The effect of the exposure on the outcome” is incomplete unless the exposure contrast, outcome scale, and target population are clear.

A genetic variant can be associated with both an exposure and an outcome. Those are associations involving the variant; neither association, on its own, defines the causal effect of changing the exposure. The causal estimand is a property of a specified exposure intervention in a specified population. Whether genetic data can identify that estimand is a separate question.

Definition: A causal estimand is a precisely specified quantity that describes the effect of an exposure intervention on an outcome in a target population. Its specification includes the exposure contrast or change, the outcome summary and scale, and the population over which the effect is defined.

Potential Outcomes and a Target Population

Let \(X\) denote an exposure and \(Y\) an outcome. For a person who could receive exposure level \(x\), write \(Y(x)\) for the outcome that would occur under that exposure level. These are potential outcomes: for any one person, only the outcome under the exposure actually received is observed. A causal contrast compares potential outcomes under alternative exposure interventions, rather than simply comparing people who happen to have different observed exposures.

The subscript \(T\) in an expectation indicates the target population. For example, \(E_T[Y(x)]\) is the mean potential outcome under exposure level \(x\) among members of that population. The target might be a defined population of adults, people with a particular clinical condition, or another group to which the scientific question refers. Specifying it matters: an average effect in one population need not be the average effect in another.

Potential outcomes must refer to sufficiently well-defined interventions. If “set exposure to 10” could mean several materially different interventions with different consequences, the notation \(Y(10)\) may not identify a single outcome. The interpretation also presumes that the outcome for one person is not changed by another person’s exposure assignment (no interference). These considerations concern the definition of the causal question, not whether an MR analysis can answer it.

Definition: For two exposure levels \(x_1\) and \(x_0\), the average causal effect on the mean-outcome scale in target population \(T\) is
$$ E_T[Y(x_1)]-E_T[Y(x_0)]. $$
It is the difference between the target-population mean outcomes under the two exposure interventions. Its units are the units of \(Y\).

Specify the Exposure Contrast

For a continuous exposure, an estimand should state the two exposure levels being compared or the size of the change. A contrast from 8 to 10 units is not automatically the same scientific quantity as a contrast from 10 to 12 units. The average change per exposure unit over a specified interval can be written as

$$ \frac{E_T[Y(x+\delta)]-E_T[Y(x)]}{\delta}, $$

where \(\delta\) is the stated exposure change. This is an average contrast per unit over that interval; it is not necessarily a constant effect at every exposure level. Calling it a “per-unit effect” without saying which interval or assuming a suitable constant-effect model can conceal this distinction.

The exposure’s units also matter. An effect per 1 mmHg is numerically different from an effect per 10 mmHg, although both may describe the same underlying contrast when correctly converted. For a transformed exposure, such as log concentration, the unit is on that transformed scale; it should not be described as a one-unit change in the original concentration.

Worked Example: A Two-Unit Exposure Contrast

Suppose, in a hypothetical target population, the mean potential outcome under exposure level 8 is 120 units and the mean potential outcome under level 10 is 126 units. Define the average change per exposure unit over this interval.

The exposure contrast is \(\delta=10-8=2\) exposure units. The difference in target-population mean outcomes is \(126-120=6\) outcome units. Therefore,

$$ \frac{E_T[Y(10)]-E_T[Y(8)]}{10-8} =\frac{126-120}{2} =3 $$

The estimand is an average increase of 3 outcome units per exposure unit from level 8 to level 10 in this hypothetical population. It does not, by itself, establish that a one-unit increase has the same effect at other exposure levels.

Specify the Outcome Scale

The same exposure intervention can have distinct causal estimands on different outcome scales. For a binary outcome, let \(Y(x)\) equal 1 if the outcome occurs under exposure level \(x\), and 0 otherwise. Then \(E_T[Y(x)]\) is the risk of the outcome under that intervention. Comparing two risks can produce a risk difference, a risk ratio, or an odds ratio; these quantities answer different questions.

Definition: If \(r_1=P_T(Y(x_1)=1)\) and \(r_0=P_T(Y(x_0)=1)\), then the risk difference is \(r_1-r_0\), the risk ratio is \(r_1/r_0\), and the odds ratio is \(\{r_1/(1-r_1)\}/\{r_0/(1-r_0)\}\), when the denominators are defined. The risk difference is an absolute change in probability; the ratio measures are relative changes.

Worked Example: Three Scales for One Binary Outcome

In a hypothetical target population, suppose the outcome risk under exposure level 1 is 0.15 and the risk under exposure level 0 is 0.10. Calculate the risk difference, risk ratio, and odds ratio for level 1 versus level 0.

Risk difference. Subtract the risks: \(0.15-0.10=0.05\). This is an increase of 0.05 in risk, or 5 percentage points.

Risk ratio. Divide the risks: \(0.15/0.10=1.50\). The risk under level 1 is 1.50 times the risk under level 0 in this target population.

Odds ratio. The odds under level 1 are \(0.15/0.85=0.1765\); the odds under level 0 are \(0.10/0.90=0.1111\). Thus,

$$ \frac{0.15/0.85}{0.10/0.90}=1.588. $$

The three answers are not interchangeable. The risk difference is an absolute probability contrast, whereas the risk ratio and odds ratio are relative contrasts on different scales. A report must name the scale rather than describing all three simply as “the effect.”

Genetic Associations Are Not the Causal Estimand

Let \(G\) denote a genetic variant. A genetic association with the exposure describes how the distribution or mean of \(X\) differs across values of \(G\). A genetic association with the outcome describes how the distribution or mean of \(Y\) differs across values of \(G\). These are statements about genotype groups in observed data. The causal estimand, by contrast, compares outcomes under interventions on \(X\).

For example, a difference in mean outcome between two genotype groups is not, just by being a genetic association, the mean outcome difference that would result from intervening to change the exposure. Genetic data can provide a basis for learning about an exposure’s causal effect only under additional identification conditions. Defining the target estimand does not establish those conditions, and this tutorial does not treat an association coefficient as a causal effect.

Worked Example: Reading a Genetic Association Correctly

Imagine hypothetical data in which the mean exposure is 4.2 units among people with genotype \(G=1\) and 3.7 units among people with \(G=0\). The corresponding mean outcomes are 12.4 and 11.9 units. Calculate and interpret these group differences.

The observed genotype-group difference in mean exposure is \(4.2-3.7=0.5\) exposure units. The observed genotype-group difference in mean outcome is \(12.4-11.9=0.5\) outcome units. Each is a genetic association: one concerns \(G\) and \(X\), and the other concerns \(G\) and \(Y\).

These calculations do not give the effect on \(Y\) of intervening to increase \(X\) by one unit. The first difference describes how the exposure varies across genotype groups; the second describes how the outcome varies across those groups. A causal effect of \(X\) requires its own defined intervention contrast and a justification that genetic information identifies that contrast.

The Target Population Changes the Average

An average causal effect combines the effects for people in the target population. If effects differ across subgroups, changing the population’s subgroup proportions can change the average, even when each subgroup-specific effect stays the same. Thus, an estimand should say whose average is intended, not merely name the exposure and outcome.

For two subgroups, suppose the same exposure contrast and outcome scale are used in both. If subgroup 1 has proportion \(w\) in target population \(T\), and subgroup 2 has proportion \(1-w\), the population average effect is the weighted average of their subgroup-specific effects:

$$ \tau_T=w\tau_1+(1-w)\tau_2. $$

This expression follows by partitioning the target-population mean under each exposure intervention by subgroup, then subtracting the two means. The weights are population proportions, not sample sizes unless the sample represents the target population in those proportions.

Worked Example: Changing the Target Population

Suppose, hypothetically, that the causal effect is 2 outcome units per exposure unit in subgroup A and 8 outcome units per exposure unit in subgroup B, for the same specified exposure contrast. A target population consists of 75% subgroup A and 25% subgroup B. Calculate its average effect, then calculate the average for a population with equal subgroup proportions.

For the first target population, \(w=0.75\), so

$$ \tau_T=(0.75)(2)+(0.25)(8)=1.5+2=3.5 $$

The average effect is 3.5 outcome units per exposure unit. With equal subgroup proportions, the calculation is

$$ \tau_T=(0.50)(2)+(0.50)(8)=1+4=5. $$

The subgroup effects did not change, but the average changed from 3.5 to 5 because the target population changed. Reporting only “the average effect” would leave the relevant population unclear.

A Practical Specification Checklist

Before interpreting an MR result as addressing a causal question, write down the intended estimand in plain language and, where useful, as a formula. Check that the words, units, and mathematical contrast agree.

1
Name the intervention contrast.
State the exposure levels or change, including the exposure scale and units.
2
Name the outcome scale.
Specify a mean difference, risk difference, risk ratio, odds ratio, or another clearly defined summary.
3
Name the target population.
State whose average effect is intended, and consider whether subgroup composition could matter.
4
Separate the target from the evidence.
Ask whether the genetic associations and the assumptions of the analysis support identification of that estimand; do not relabel a genetic association as the causal effect.

Common Mistake

Wrong reasoning: “The variant is associated with the outcome, so this association is the causal effect of the exposure on the outcome.” This confuses the variable defining the genetic groups with the exposure being intervened on. The association is about genotype and outcome; the causal estimand is about outcomes under alternative exposure interventions.

Correction: State the exposure contrast, outcome scale, and target population first. Then assess whether the genetic analysis can identify that causal quantity under its assumptions. A carefully defined estimand is necessary for interpretation, but it is not evidence that identification assumptions hold.

Key takeaway: An MR causal estimand is not just a coefficient or a genetic association. It is a specified causal contrast: an exposure change on a stated scale, an outcome summary on a stated scale, and an average for a stated target population.

Check Your Understanding

For each question, identify the estimand feature that needs to be made explicit.

  1. Why is “the effect of a continuous exposure” incomplete if the exposure contrast and units are not given?
  2. For a binary outcome, what does a risk difference express that a risk ratio does not?
  3. A variant is associated with both an exposure and an outcome. Why does that alone not define the causal effect of the exposure?
  4. Two populations have different proportions of subgroups with different causal effects. Why can their average effects differ?
  5. What three core features should be stated when describing an MR causal estimand?
Clinical Biostats