GTEx SNP-gene association statistics used in mash analysis

Last updated: 2018-06-21

workflowr checks: (Click a bullet for more information)

✔ R Markdown file: up-to-date

Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

✔ Repository version: f77307a

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility. The version displayed above was the version of the Git repository at the time these results were generated.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:


Ignored files:
    Ignored:    output/bmaonlybetasd5lfsr.txt.gz
    Ignored:    output/bmaonlybetasd5posterior.means.txt.gz
    Ignored:    output/independentsim.rds
    Ignored:    output/independentsimesd05.rds
    Ignored:    output/indsim05sdlfsr.txt.gz
    Ignored:    output/indsim05sdposterior.means.txt.gz
    Ignored:    output/noashsharedwithzerobmaalllfsr.txt.gz
    Ignored:    output/noashsharedwithzerobmaallposterior.means.txt.gz
    Ignored:    output/sharedashcutoffomega2jun15lfsr.txt.gz
    Ignored:    output/sharedashcutoffomega2jun15posterior.means.txt.gz
    Ignored:    output/simdata.rds
    Ignored:    output/univariate.ash.lfsr.txt.gz
    Ignored:    output/univariate.ash.pm.txt.gz
    Ignored:    output/univariate.ash.pmindesd.txt.gz
    Ignored:    output/univariate.ashind.lfsresd.txt.gz

Unstaged changes:
    Modified:   analysis/_site.yml

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.

Expand here to see past versions:

File	Version	Author	Date	Message
Rmd	f77307a	Peter Carbonetto	2018-06-21	wflow_publish(“gtexdata.Rmd”)
Rmd	87259b7	Peter Carbonetto	2018-06-21	wflow_publish(“gtexdata.Rmd”)

To apply multivariate adaptive shrinkage (mash) to data from the GTEx study, we created an R data set (serialized R object) containing matrices of SNP-gene association statistics. These association statistics include effect estimates, Z scores and corresponding standard errors.

See here for the scripts used to generate these statistics from the SNP-gene data that were provided by the GTEx Project.

How to download the data

These are the recommended steps for retrieving the GTEx SNP-gene association statistics:

Download or clone the git repository.
The association statistics are found in file MatrixEQTLSumStats.Portable.Z.rds.

How to load the data into R

Change the working directory in R (or RStudio) to the analysis directory of the gtexresults repository, e.g.,

R setwd("gtexresults/analysis")

Read the data object into R:

R dat <- readRDS("../data/MatrixEQTLSumStats.Portable.Z.rds")

Get an overview of the data from this file:

R names(dat)

Overview of SNP-gene association statistics

This file contains SNP-gene association statistics for 16,069 genes and 44 human tissues. These 16,069 genes were selected because they all showed some indication of being expressed in all 44 tissues. Therefore, the association statistics are stored as matrices each with 16,069 rows and 44 columns.

As input to mash, we use a matrix of expression quantitative trait loci (eQTL) effect estimate, and corresponding standard errors. (We also provide Z scores.) See the manuscript for details on how these association statistics were obtained.

These association statistics were subdivided into three subsets:

Results from a subset `strong” tests corresponding to stronger effects in your study. For example, these tests might have been identified by taking the “top” eQTL in each gene based on univariate test results, or by some other approach such as a simple meta-analysis.

Results from a random subset of all tests. It is important that these be an unbiased representation of all the tests you are considering, including null and non-null tests, because mashr uses these tests to learn about the amount of signal in the data, and to “correct” estimates for the fact that many tests are null (analagous to a kind of multiple testing correction.

This reproducible R Markdown analysis was created with workflowr 1.0.1.9000

GTEx SNP-gene association statistics used in mash analysis

Sarah Urbut, Gao Wang, Peter Carbonetto and Matthew Stephens

How to download the data

How to load the data into R

Overview of SNP-gene association statistics