qmarcou / IGoR

IGoR is a C++ software designed to infer V(D)J recombination related processes from sequencing data. Find full documentation at:
https://qmarcou.github.io/IGoR/
GNU General Public License v3.0
47 stars 25 forks source link

All 0 output in infer and question about L_thresh #71

Open zhuyixin43 opened 1 year ago

zhuyixin43 commented 1 year ago

Hi, I noticed that when I am running -infer I always get all zero parameters in the _finalmarginals.txt. This leads to no difference in the generated sequence for all my samples. The documentation suggested me to lower my --L_thresh. If I'm not mistaken, --L_thresh is set to default 1e-60 right? I tried to change it to different values such as 1e-50 or 1e-70, but the problem persisted and some samples' running processes gets killed. I also check the "mean number of errors" in the _inferencelogs.txt, which is all zeros, so I don't think it's the alignment issue? Could you give me some guidance on how to select a reasonable --L_thresh value? Should I be looking at the _inferencelogs.txt "seq_best_senario" value and use it to get a general idea of my sequence's probability?

Thanks a lot!