hakyimlab/PrediXcan

predict.py error

Opened this issue · 0 comments

Hello,
I have a question related to predixcan software.
I am trying to run predict.py script. But its showing error for me.
The command use:
python $PXCN_TOOLS/predict.py --model_db_path $MODELS/en_Whole_Blood.db --model_db_snp_key rsid --vcf_mode genotyped --vcf_genotypes $VCF_FILES/*.vcf --prediction_output $OUTPUT/GVDS_PrediXcan_Test_2021.txt
the error:
[E::bcf_hdr_parse] Could not parse the header, sample line not found
Segmentation fault
I checked the vcf files but they have the sample line:
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT HG00105 HG00107 HG00115 HG00132 HG00145 HG00157 HG00181 HG00308 HG00365 HG00371 HG00379 HG00380 HG01789 HG01790 HG01791 HG02215 NA06985 NA07346 NA11832 NA11840 NA11881 NA11918 NA12005 NA12156 NA12234 NA12760 NA12762 NA12776 NA12813 NA18488 NA19092 NA19141 NA19143 NA19144 NA19153 NA19159 NA19184 NA19201 NA19206 NA19210 NA19214 NA20514 HG00096 HG00097 HG00099 HG00100 HG00101 HG00102 HG00103 HG00104 HG00106 HG00108 HG00109 HG00110 HG00111 HG00112 HG00114 HG00116 HG00117 HG00118 HG00119 HG00120 HG00121 HG00122 HG00123 HG00124 HG00125 HG00126 HG00127
I searched this error in google and it says to make sure that its tab delimited and not space. I checked the vcf files and its sample line is separated by tab.
Can you help me out with this?
Thank you.