
运行前两条染色体是正常的,运行到第三条染色体开始出现如图报错index error。
其中,所有染色体ldscore文件(Add_all.chr)和参考的文件都存在且无缺损。每条染色体单独分开运行也是同样的问题(从第三条染色体开始报错)。
代码
for i in {1..22}do python ./ldsc.py --h2 SSc_ldscnew.sumstats.gz --ref-ld-chr regions_enh_E048_folder/Add_all.chr${i} --overlap-annot --w-ld-chr 1000G_Phase3_weights_hm3_no_MHC/weights.hm3_noMHC.${i} --frqfile-chr 1000G_Phase3_frq/1000G.EUR.QC.${i} --out SSc048done