Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

MuTect2 LOD calculations

$
0
0

I am interested in how NLOD and TLOD are calculated in MuTect2. The MuTect2.java code (https://github.com/broadgsa/gatk-protected/blob/master/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/cancer/m2/MuTect2.java#L856) leads me to a function calcGenotypeLikelihoodsOfRefVsAny(), in which when a non-ref is observed:
genotypeLikelihoods[AB] += Math.log10(fpobs + (1-f)pobs/3.0d);
However, my understanding is that the real likelihood should be
genotypeLikelihoods[AB] += Math.log10(fpobs + (1-f)(1-pobs)/3.0d);

Is this an error, or could you clarify how these likelihood values, NLOD, TLOD are calculated?

I have read the algorithm guide here https://software.broadinstitute.org/gatk/guide/article?id=4442, but couldn't find my answer.
I have also read the MuTect1 paper https://www.ncbi.nlm.nih.gov/pubmed/23396013, but I don't know if they are the same as in MuTect2.


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>