In MRCIEU/TwoSampleMR: Two Sample MR Functions and Interface to MR Base Database

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.align = 'center'
)

Background

Two sample Mendelian randomisation (2SMR) is a method to estimate the causal effect of an exposure on an outcome using only summary statistics from genome wide association studies (GWAS). Though conceptually straightforward, there are a number of steps that are required to perform the analysis properly, and they can be cumbersome. The TwoSampleMR package aims to make this easy by combining three important components

• data management and harmonisation
• the statistical routines to estimate the causal effects
• connection to a large repository of the actual GWAS summary statistics needed to perform the analyses.

The general principles [@DaveySmith2003; @DaveySmithHemani2014], and statistical methods [@Pierce2013; @Bowden2015] can be found elsewhere, here we will just outline how to use the R package.

This package uses the ieugwasr package to connect to the database of thousands of complete GWAS summary data.

---

Installation

To install directly from the GitHub repository do the following:

library(remotes)
install_github("MRCIEU/TwoSampleMR")

If you don't have the remotes package install it from CRAN using install.packages("remotes").

---

Overview

The workflow for performing MR is as follows:

1. Select instruments for the exposure (perform LD clumping if necessary)
2. Extract the instruments from the IEU GWAS database for the outcomes of interest
3. Harmonise the effect sizes for the instruments on the exposures and the outcomes to be each for the same reference allele
4. Perform MR analysis, sensitivity analyses, create plots, compile reports

A diagrammatic overview is shown here:

knitr::include_graphics("img/twosamplemr_schematic_long-01.png")

A basic analysis, e.g. the causal effect of body mass index on coronary heart disease, looks like this:

library(TwoSampleMR)

# List available GWASs
ao <- available_outcomes()

# Get instruments
exposure_dat <- extract_instruments("ieu-a-2")

# Get effects of instruments on outcome
outcome_dat <- extract_outcome_data(snps=exposure_dat$SNP, outcomes = "ieu-a-7")

# Harmonise the exposure and outcome data
dat <- harmonise_data(exposure_dat, outcome_dat)

# Perform MR
res <- mr(dat)

Each step is documented on other pages in the documentation.

Authentication

The statistical methods in TwoSampleMR can be used on any data, but there are a number of functions that connect to the OpenGWAS database for data extraction. These OpenGWAS data access functions require authentication.

Authentication is changing The main differences are that:

1. Authentication is required for most queries to OpenGWAS for everyone (i.e. no more anonymous usage)
2. We are no longer using Google Oauth2. This has been replaced by a simple API key system.

Detailed information is given here: https://mrcieu.github.io/ieugwasr/articles/guide.html#authentication.

References

MRCIEU/TwoSampleMR documentation built on May 2, 2024, 10:22 p.m.