Hi everyone,
I have a question about the difference between --h2 and --h2-cts flags. I know that h2-cts is used to partition heritability to cell-type-specific annotations and LDSC performs a one-sided test when this flag is used, so all heritability estimates are positive. However, I run part. heritability with custom annotations, so I'm using the h2 flag, which runs a two-sided test. As a result, sometimes I get negative heritability estimates. So, I was wondering if there is a way to run a one-sided test using the h2 flag and custom annotations.
I also use a custom baseline model in addition to the generic baseline model. Here is the command:
./ldsc.py \
--h2 /data/munged/ancreg_munged.txt.sumstats.gz \
--ref-ld-chr /data/active_merged/active_merged.,/data/active_marks/active_marks.,/data/baselineLD/baselineLD. \
--out /data/results/active_merged/ancreg_munged.txt.sumstats.gz \
--overlap-annot \
--frqfile-chr /data/LDscores/Phase3/1000G_Phase3_frq/1000G.EUR.QC. \
--w-ld-chr /data/LDscores/Phase3/1000G_Phase3_weights_hm3_no_MHC/weights.hm3_noMHC. \
--print-coefficients
Here is the log file:
Beginning analysis at Thu Dec 3 12:06:19 2020 Reading summary statistics from /data/ancreg_munged.txt.sumstats.gz ... Read summary statistics for 1098387 SNPs. Reading reference panel LD Score from /data/active_merged/active_merged.,/data/active_marks/active_marks.,/data/baselineLD/baselineLD.[1-22] ... (ldscore_fromlist) Read reference panel LD Scores for 1190321 SNPs. Removing partitioned LD Scores with zero variance. Reading regression weight LD Score from /data/LDscores/Phase3/1000G_Phase3_weights_hm3_no_MHC/weights.hm3_noMHC.[1-22] ... (ldscore_fromlist) Read regression weight LD Scores for 1187349 SNPs. After merging with reference panel LD, 1094449 SNPs remain. After merging with regression SNP LD, 1094449 SNPs remain. Removed 29 SNPs with chi^2 > 80 (1094420 SNPs remain) Total Observed scale h2: 0.3553 (0.0509) Categories: L2_0 L2_1 baseL2_2 Coding_UCSCL2_2 Coding_UCSC.flanking.500L2_2 Conserved_LindbladTohL2_2 Conserved_LindbladToh.flanking.500L2_2 CTCF_HoffmanL2_2 CTCF_Hoffman.flanking.500L2_2 DGF_ENCODEL2_2 DGF_ENCODE.flanking.500L2_2 DHS_peaks_TrynkaL2_2 DHS_TrynkaL2_2 DHS_Trynka.flanking.500L2_2 Enhancer_AnderssonL2_2 Enhancer_Andersson.flanking.500L2_2 Enhancer_HoffmanL2_2 Enhancer_Hoffman.flanking.500L2_2 FetalDHS_TrynkaL2_2 FetalDHS_Trynka.flanking.500L2_2 H3K27ac_HniszL2_2 H3K27ac_Hnisz.flanking.500L2_2 H3K27ac_PGC2L2_2 H3K27ac_PGC2.flanking.500L2_2 H3K4me1_peaks_TrynkaL2_2 H3K4me1_TrynkaL2_2 H3K4me1_Trynka.flanking.500L2_2 H3K4me3_peaks_TrynkaL2_2 H3K4me3_TrynkaL2_2 H3K4me3_Trynka.flanking.500L2_2 H3K9ac_peaks_TrynkaL2_2 H3K9ac_TrynkaL2_2 H3K9ac_Trynka.flanking.500L2_2 Intron_UCSCL2_2 Intron_UCSC.flanking.500L2_2 PromoterFlanking_HoffmanL2_2 PromoterFlanking_Hoffman.flanking.500L2_2 Promoter_UCSCL2_2 Promoter_UCSC.flanking.500L2_2 Repressed_HoffmanL2_2 Repressed_Hoffman.flanking.500L2_2 SuperEnhancer_HniszL2_2 SuperEnhancer_Hnisz.flanking.500L2_2 TFBS_ENCODEL2_2 TFBS_ENCODE.flanking.500L2_2 Transcr_HoffmanL2_2 Transcr_Hoffman.flanking.500L2_2 TSS_HoffmanL2_2 TSS_Hoffman.flanking.500L2_2 UTR_3_UCSCL2_2 UTR_3_UCSC.flanking.500L2_2 UTR_5_UCSCL2_2 UTR_5_UCSC.flanking.500L2_2 WeakEnhancer_HoffmanL2_2 WeakEnhancer_Hoffman.flanking.500L2_2 GERP.NSL2_2 GERP.RSsup4L2_2 MAFbin1L2_2 MAFbin2L2_2 MAFbin3L2_2 MAFbin4L2_2 MAFbin5L2_2 MAFbin6L2_2 MAFbin7L2_2 MAFbin8L2_2 MAFbin9L2_2 MAFbin10L2_2 MAF_Adj_Predicted_Allele_AgeL2_2 MAF_Adj_LLD_AFRL2_2 Recomb_Rate_10kbL2_2 Nucleotide_Diversity_10kbL2_2 Backgrd_Selection_StatL2_2 CpG_Content_50kbL2_2 MAF_Adj_ASMCL2_2 GTEx_eQTL_MaxCPPL2_2 BLUEPRINT_H3K27acQTL_MaxCPPL2_2 BLUEPRINT_H3K4me1QTL_MaxCPPL2_2 BLUEPRINT_DNA_methylation_MaxCPPL2_2 synonymousL2_2 non_synonymousL2_2 Conserved_Vertebrate_phastCons46wayL2_2 Conserved_Vertebrate_phastCons46way.flanking.500L2_2 Conserved_Mammal_phastCons46wayL2_2 Conserved_Mammal_phastCons46way.flanking.500L2_2 Conserved_Primate_phastCons46wayL2_2 Conserved_Primate_phastCons46way.flanking.500L2_2 BivFlnkL2_2 BivFlnk.flanking.500L2_2 Human_Promoter_VillarL2_2 Human_Promoter_Villar.flanking.500L2_2 Human_Enhancer_VillarL2_2 Human_Enhancer_Villar.flanking.500L2_2 Ancient_Sequence_Age_Human_PromoterL2_2 Ancient_Sequence_Age_Human_Promoter.flanking.500L2_2 Ancient_Sequence_Age_Human_EnhancerL2_2 Ancient_Sequence_Age_Human_Enhancer.flanking.500L2_2 Human_Enhancer_Villar_Species_Enhancer_CountL2_2 Human_Promoter_Villar_ExACL2_2 Human_Promoter_Villar_ExAC.flanking.500L2_2 Lambda GC: 1.1207 Mean Chi^2: 1.1573 Intercept: 1.0206 (0.0104) Ratio: 0.1308 (0.0662) Reading annot matrix from /data/active_merged/active_merged.,/data/active_marks/active_marks.,/data/baselineLD/baselineLD.[1-22] ... (annot) Results printed to /data/active_merged/ancreg_munged.txt.sumstats.gz.results Analysis finished at Thu Dec 3 12:15:19 2020 Total time elapsed: 8.0m:59.99s
