how to calculate q_vaule of GWAS Heritability Enrichment as your paper mentioned #7
Description
Hello,I'm glad to see your article,but feel some confused about how to evaluate significant heritability for each trait/cluster pairs , which descibes as below in your paper.
"P-values were calculated from z-scores assigned to coefficients reported by ldsc.py and coefficients were divided by the average per-SNP heritability for trait associated with a given test (as calculated from the number of SNPs and overall heritability reported in the .log files from ldsc.py), as recommended by the LDSC authors, producing scaled coefficients. These scaled coefficients are provided in supplemental data files but are not reported in figures here. Tests were corrected for multiple hypothesis testing using the Benjamini-Hochberg method and only tests with a q-value of 0.05 or lower were deemed to be significant."
I don't konw how to calculate 'coefficient_pval', 'coefficient_qval' and 'scaled_coefficients' . Could you please show the formula or some simple example? I would feel grateful if you can give me your advice at your convenience.