zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

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
This commit is contained in:
rpoplin 2010-10-10 17:35:15 +00:00
parent 116309b3c3
commit 0a4cf02a52
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantGaussianMixtureModel.java
View File

@ -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++ ) {