hess-pipeline

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

The scripts 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 is done to make it generic including examples for height and QUS meta-analyses. These only involve SNP ids the as sorting key. The scripts are hess.sh, hess.subs and hess.awk.

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

The syntax is as follows,

hess.sh GWAS-summary-statistics

where GWAS-summary-statistics has these columns: rsID, pos, A1, A2, Z, N.

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