demuxlet
demuxlet copied to clipboard
Demuxlet Output Questions
I am trying to understand the output of demuxlet, in particular the .best file.
LLK1is forDBL.1STwhileLLK2is forDBL.2ND, but forLLK12, what does theIDmean in this case? ID of what?- I noticed that for almost all samples with the singlet as the best assignment, that
LLK12is exactly equal toSNG.LLK1. The code to pick the "best" assignment will pick singlet in this case, due to thepairLLK12 > singLLK1 + 2. But if the likelihoods are equal (or close), shouldn't the output be ambiguous (modulo some prior)? I don't see a way to control this case, as I'd like to tune my output to only have unambiguous assignments that have a minimum confidence. - Is there an easy way of getting the posterior probability of a singlet like we have for
PRB.DBL?
I am using --alpha 0.0 --alpha 0.5 --alpha 1.0 for two sample mixture.
Perhaps the --alpha 1.0 is the issue, as it's allowing a "doublet" with 100% one of the samples?