Package picard.fingerprint
Class HaplotypeProbabilitiesFromSequence
- java.lang.Object
-
- picard.fingerprint.HaplotypeProbabilities
-
- picard.fingerprint.HaplotypeProbabilitiesFromSequence
-
- Direct Known Subclasses:
HaplotypeProbabilitiesFromContaminatorSequence
public class HaplotypeProbabilitiesFromSequence extends HaplotypeProbabilities
Represents the probability of the underlying haplotype given the data. By convention the alleles stored for each SNP are in phase.
-
-
Nested Class Summary
-
Nested classes/interfaces inherited from class picard.fingerprint.HaplotypeProbabilities
HaplotypeProbabilities.Genotype
-
-
Field Summary
Fields Modifier and Type Field Description protected int
obsAllele1
protected int
obsAllele2
protected int
obsAlleleOther
-
Constructor Summary
Constructors Constructor Description HaplotypeProbabilitiesFromSequence(HaplotypeBlock haplotypeBlock)
-
Method Summary
All Methods Instance Methods Concrete Methods Modifier and Type Method Description void
addToProbs(Snp snp, byte base, byte qual)
Adds a base observation with the observed quality to the evidence for this haplotype based on the fact that the SNP is part of the haplotype.double
getFractionUnexpectedAlleleObs()
double[]
getLikelihoods()
getter for likelihoods;double
getLodMostProbableGenotype()
Overridden to calculate the LOD from the loglikelihoods instead of the probabilities because it will allow for more accurate calculation before overflowing.double[]
getLogLikelihoods()
Since this class uses loglikelihoods natively, we override and return the native variableint
getObsAllele1()
Returns the number of bases/reads that support the first allele.int
getObsAllele2()
Returns the number of bases/reads that support the second allele.double[]
getPosteriorProbabilities()
getter for posteriorsSnp
getRepresentativeSnp()
Simple returns the SNP from the haplotype that has the lowest genome coordinate.int
getTotalObs()
Gets the total number of observations presented at this locus.boolean
hasEvidence()
Returns true if evidence has been added, false if the probabilities are just the priors.void
merge(HaplotypeProbabilities other)
Merges information from another haplotype probabilities object for the same haplotype into this object.void
setLogLikelihoods(double[] ll)
-
Methods inherited from class picard.fingerprint.HaplotypeProbabilities
getHaplotype, getMostLikelyGenotype, getMostLikelyHaplotype, getPriorProbablities, scaledEvidenceProbabilityUsingGenotypeFrequencies, shiftedLogEvidenceProbability, shiftedLogEvidenceProbabilityGivenOtherEvidence, shiftedLogEvidenceProbabilityUsingGenotypeFrequencies
-
-
-
-
Constructor Detail
-
HaplotypeProbabilitiesFromSequence
public HaplotypeProbabilitiesFromSequence(HaplotypeBlock haplotypeBlock)
-
-
Method Detail
-
hasEvidence
public boolean hasEvidence()
Description copied from class:HaplotypeProbabilities
Returns true if evidence has been added, false if the probabilities are just the priors.
-
addToProbs
public void addToProbs(Snp snp, byte base, byte qual)
Adds a base observation with the observed quality to the evidence for this haplotype based on the fact that the SNP is part of the haplotype.- Parameters:
snp
- The snp in the HaplotypeBlock to which evidence is being addedbase
- the base observedqual
- the quality of the observed base
-
merge
public void merge(HaplotypeProbabilities other)
Merges information from another haplotype probabilities object for the same haplotype into this object. Useful for when probabilities need to be merged to levels higher than the read group, e.g. the sample or individual.- Parameters:
other
- Another haplotype probabilities object to merge in
-
getObsAllele1
public int getObsAllele1()
Returns the number of bases/reads that support the first allele.- Overrides:
getObsAllele1
in classHaplotypeProbabilities
- Returns:
- int
-
getObsAllele2
public int getObsAllele2()
Returns the number of bases/reads that support the second allele.- Overrides:
getObsAllele2
in classHaplotypeProbabilities
- Returns:
- int
-
getTotalObs
public int getTotalObs()
Gets the total number of observations presented at this locus.- Overrides:
getTotalObs
in classHaplotypeProbabilities
- Returns:
- int
-
getFractionUnexpectedAlleleObs
public double getFractionUnexpectedAlleleObs()
-
getRepresentativeSnp
public Snp getRepresentativeSnp()
Simple returns the SNP from the haplotype that has the lowest genome coordinate.- Specified by:
getRepresentativeSnp
in classHaplotypeProbabilities
-
getPosteriorProbabilities
public double[] getPosteriorProbabilities()
getter for posteriors- Overrides:
getPosteriorProbabilities
in classHaplotypeProbabilities
-
getLikelihoods
public double[] getLikelihoods()
getter for likelihoods;- Specified by:
getLikelihoods
in classHaplotypeProbabilities
-
getLogLikelihoods
public double[] getLogLikelihoods()
Since this class uses loglikelihoods natively, we override and return the native variable- Overrides:
getLogLikelihoods
in classHaplotypeProbabilities
-
setLogLikelihoods
public void setLogLikelihoods(double[] ll)
-
getLodMostProbableGenotype
public double getLodMostProbableGenotype()
Overridden to calculate the LOD from the loglikelihoods instead of the probabilities because it will allow for more accurate calculation before overflowing.- Overrides:
getLodMostProbableGenotype
in classHaplotypeProbabilities
-
-