From 0a4cf02a523f747ea076b59329f7369f3e1f338a Mon Sep 17 00:00:00 2001 From: rpoplin Date: Sun, 10 Oct 2010 17:35:15 +0000 Subject: [PATCH] Fix for index out of bounds exception in VR. git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4473 348d0f76-0448-11de-a6fe-93d51630548a --- .../VariantGaussianMixtureModel.java | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantGaussianMixtureModel.java b/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantGaussianMixtureModel.java index a8f282bc9..b436fba3b 100755 --- a/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantGaussianMixtureModel.java +++ b/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantGaussianMixtureModel.java @@ -472,11 +472,11 @@ public final class VariantGaussianMixtureModel extends VariantOptimizationModel double QUAL_STEP2 = (MAX_QUAL * 1.0 / 10.0) / ((double) NUM_BINS / 2.0); if( QUAL_STEP2 < 0.01 ) { QUAL_STEP2 = 0.01; } // QUAL field in VCF file is rounded to two decimal places - final int numKnownAtCut[] = new int[NUM_BINS+1]; - final int numNovelAtCut[] = new int[NUM_BINS+1]; - final double knownTiTvAtCut[] = new double[NUM_BINS+1]; - final double novelTiTvAtCut[] = new double[NUM_BINS+1]; - final double theCut[] = new double[NUM_BINS+1]; + final int numKnownAtCut[] = new int[NUM_BINS+2]; + final int numNovelAtCut[] = new int[NUM_BINS+2]; + final double knownTiTvAtCut[] = new double[NUM_BINS+2]; + final double novelTiTvAtCut[] = new double[NUM_BINS+2]; + final double theCut[] = new double[NUM_BINS+2]; final double fdrCutAsTiTv[] = new double[FDRtranches.length]; for( int iii = 0; iii < FDRtranches.length; iii++ ) {