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.
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/).
The output will be available from its corresponding .tmp subdirectory.