diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/ApplyVariantClustersWalker.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/ApplyVariantClustersWalker.java index 67e11bbe8..37e05ff9a 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/ApplyVariantClustersWalker.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/ApplyVariantClustersWalker.java @@ -161,14 +161,8 @@ public class ApplyVariantClustersWalker extends RodWalker 0.1 ) { // only use the known prior if the value is specified (meaning not equal to zero) variantDatum.qual = 0.5 * recalQual + 0.5 * KNOWN_VAR_QUAL_PRIOR; } else { diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/VariantGaussianMixtureModel.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/VariantGaussianMixtureModel.java index 2b11b62d1..0a1d301ec 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/VariantGaussianMixtureModel.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantoptimizer/VariantGaussianMixtureModel.java @@ -721,6 +721,7 @@ public final class VariantGaussianMixtureModel extends VariantOptimizationModel final double denom = Math.pow(2.0 * 3.14159, ((double)numAnnotations) / 2.0) * Math.pow(determinant[kkk], 0.5); pVarInCluster[kkk] = (1.0 / ((double) numGaussians)) * (Math.exp( -0.5 * sum )) / denom; + /* if( isUsingTiTvModel ) { //pVarInCluster[kkk] = Math.exp( -0.5 * sum ); if( pVarInCluster[kkk] < MIN_PROB) { // Very small numbers are a very big problem @@ -733,13 +734,14 @@ public final class VariantGaussianMixtureModel extends VariantOptimizationModel //pVarInCluster[kkk] = Math.exp( -0.5 * sum ); // BUGBUG: should pCluster be the distribution from the GMM or a uniform distribution here? } + */ } - if( isUsingTiTvModel ) { - for( int kkk = 0; kkk < numGaussians; kkk++ ) { - pVarInCluster[kkk] /= sumProb; - } - } + //if( isUsingTiTvModel ) { + // for( int kkk = 0; kkk < numGaussians; kkk++ ) { + // pVarInCluster[kkk] /= sumProb; + // } + //} }