hess-pipeline

Introduction

The scripts here are developed to facilitate heritability estimation from summary statistics (HESS) as available from https://github.com/huwenboshi/hess.

They were initially tested using EPIC-Norfolk GWAS, for which test.sh, test.subs and test.awk are created to work on /genetics/bin/hess/tests directory and later on changed to /genetics/data/gwas/8-7-15/QUS.

Further work was done to make it generic with GIANT data on height meta-analysis.

Recently there has been considerable changes in HESS, documented in https://huwenboshi.github.io/hess/; ealier version is available from the release as well as the oldrelease directory.

Setup

To set it up, modify hess.sh as appropriate for your system.

The syntax is as follows,

hess.sh GWAS-summary-statistics

where GWAS-summary-statistics is a file in the current directory with these columns: rsID, pos, A1, A2, Z, N, e.g., code as in height.sh, since it only involve rsids chromosomal positions were added on.

NOTE that HESS describes bed file in line with UCSC (https://genome.ucsc.edu/).

Output

The output will be available from its corresponding .tmp subdirectory.