search

liusihan / seGMM

A new tool to infer sex from massively parallel sequencing data.
MIT License
12 stars 2 forks source link

paste XH.sorted.txt with unsorted Xmap.txt results in error results #3

Open ShannonDaddy opened 1 year ago

ShannonDaddy commented 1 year ago

Hi, thanks for the wonderful tool seGMM. But I may have found a serious bug when I use it in the WES mode. In my data, I have a sample 1-XJ-XY20220727002 with Xmap value 0.0599118 in the Xmap.txt file, but in the feature.txt and seGMM_result.txt, the Xmap value become 0.031172, which is totally different. And the final gender is wrong. And I read the source code, I found that the XH.sorted.txt，Xmap.txt, Ymap.txt are directly pasted together. But the XH.sorted.txt has already sorted with a different sample order, the sample order may be totally different with Xmap.txt and Ymap.txt etc. Thus, the feature is wrong, which results in error gender calculation.

my data snapshots: image image image image image

please check that, thanks!

liusihan commented 1 year ago

Thanks for your reminder and we apologize for this bug. Please update seGMM to the newest version which we have fixed this bug.

ShannonDaddy commented 1 year ago

Thanks for your reminder and we apologize for this bug. Please update seGMM to the newest version which we have fixed this bug.

Thanks！