odelaneau / shapeit4

Segmented HAPlotype Estimation and Imputation Tool
MIT License
90 stars 18 forks source link

shapeit4: src/phaser/phaser_algorithm.cpp:50: void phaser::phaseWindow(int, int): Assertion `threadData[id_worker].Kvec[w].size()>0' failed. #4

Closed StuartKim closed 5 years ago

StuartKim commented 5 years ago

Hi,

I am getting the same error message as mdkeehan. In my case, 89.chr20.vcf.gz has a single human sample with ~14,400 SNPs on chr20. I get the error message with the default window, and 1e6, 1e7 and 20e7 used in the --window command. I get the error message for chr2 as well.

shapeit4 \ --input 89.chr20.vcf.gz \ --map /shapeit4/maps/chr20.b37.gmap.gz \ --region 20 \ --window 20e7 \ --output phased_89.20.vcf.gz

Can you help please? Stuart Kim