From 6bda5696664da40bc2baef4f4cb69e4ef1f86ce5 Mon Sep 17 00:00:00 2001 From: Ryan Poplin Date: Wed, 28 Aug 2013 12:40:08 -0400 Subject: [PATCH] One of the log10sumLog10s in the VQSR was missed in a previous bug fix. Thanks to Mike McCowan for spotting this one. --- .../gatk/walkers/variantrecalibration/GaussianMixtureModel.java | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/protected/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/GaussianMixtureModel.java b/protected/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/GaussianMixtureModel.java index 864abad7d..ed26bc17a 100644 --- a/protected/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/GaussianMixtureModel.java +++ b/protected/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/GaussianMixtureModel.java @@ -283,7 +283,7 @@ public class GaussianMixtureModel { } // add this sample's probability to the pile in order to take an average in the end - sumPVarInGaussian += Math.pow(10.0, MathUtils.log10sumLog10(pVarInGaussianLog10)); // p = 10 ^ Sum(pi_k * p(v|n,k)) + sumPVarInGaussian += Math.pow(10.0, nanTolerantLog10SumLog10(pVarInGaussianLog10)); // p = 10 ^ Sum(pi_k * p(v|n,k)) numRandomDraws++; } }