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

Add error handling for genotypes that are called but have no PLs

This commit is contained in:
Laura Gauthier 2014-12-18 13:28:56 -05:00
parent a4a9e73ec8
commit a9694951d2
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/variantutils/FamilyLikelihoodsUtils.java
View File

@ -107,9 +107,9 @@ public class FamilyLikelihoodsUtils {
*/ */
public void getUpdatedGenotypes(final VariantContext vc, final Genotype motherGenotype, final Genotype fatherGenotype, final Genotype childGenotype, final ArrayList<Genotype> updatedGenotypes){ public void getUpdatedGenotypes(final VariantContext vc, final Genotype motherGenotype, final Genotype fatherGenotype, final Genotype childGenotype, final ArrayList<Genotype> updatedGenotypes){
//genotypes here can be no call //genotypes here can be no call
boolean fatherIsCalled = fatherGenotype != null && hasCalledGT(fatherGenotype.getType()); boolean fatherIsCalled = fatherGenotype != null && hasCalledGT(fatherGenotype.getType()) && fatherGenotype.hasLikelihoods();
boolean motherIsCalled = motherGenotype != null && hasCalledGT(motherGenotype.getType()); boolean motherIsCalled = motherGenotype != null && hasCalledGT(motherGenotype.getType()) && motherGenotype.hasLikelihoods();
boolean childIsCalled = childGenotype != null && hasCalledGT(childGenotype.getType()); boolean childIsCalled = childGenotype != null && hasCalledGT(childGenotype.getType()) && childGenotype.hasLikelihoods();
//default to posteriors equal to likelihoods (flat priors) in case input genotypes are not called //default to posteriors equal to likelihoods (flat priors) in case input genotypes are not called
double[] uninformativeLikelihoods = {ONE_THIRD, ONE_THIRD, ONE_THIRD}; double[] uninformativeLikelihoods = {ONE_THIRD, ONE_THIRD, ONE_THIRD};