polyploid-genotyping icon indicating copy to clipboard operation
polyploid-genotyping copied to clipboard

Improper log-sum of exponentials in log likelihood

Open pblischak opened this issue 7 years ago • 0 comments

The log likelihoods of the different models are computed using log sums of exponentials, which can have nasty overflow/underflow problems (cf. this link). I've run into this problem when read counts are high (on the order of ~1000x). Since sequencing coverage isn't typically this high, it hopefully should be an issue for most data sets. Nevertheless, I'll need to fix it.

pblischak avatar Jan 20 '17 15:01 pblischak