Instrumental variable analysis with a nonlinear exposure-outcome relationship

Instrumental variable methods can estimate the causal effect of an exposure on an outcome using observational data. Many instrumental variable methods assume that the exposure-outcome relation is linear, but in practice this assumption is often in doubt, or perhaps the shape of the relation is a target for investigation. We investigate this issue in the context of Mendelian randomization, the use of genetic variants as instrumental variables.

Using simulations, we demonstrate the performance of a simple linear instrumental variable method when the true shape of the exposure-outcome relation is not linear. We also present a novel method for estimating the effect of the exposure on the outcome within strata of the exposure distribution. This enables the estimation of localized average causal effects within quantile groups of the exposure or as a continuous function of the exposure using a sliding window approach.

Our simulations suggest that linear instrumental variable estimates approximate a population-averaged causal effect. This is the average difference in the outcome if the exposure for every individual in the population is increased by a fixed amount. Estimates of localized average causal effects reveal the shape of the exposure-outcome relation for a variety of models. These methods are used to investigate the relations between body mass index and a range of cardiovascular risk factors.

Nonlinear exposure-outcome relations should not be a barrier to instrumental variable analyses. When the exposure-outcome relation is not linear, either a population-averaged causal effect or the shape of the exposure-outcome relation can be estimated.