public class LearnReadOrientationModelEngine
extends java.lang.Object
Constructor and Description |
---|
LearnReadOrientationModelEngine(htsjdk.samtools.util.Histogram<java.lang.Integer> refHistogram,
java.util.List<htsjdk.samtools.util.Histogram<java.lang.Integer>> altDepthOneHistograms,
java.util.List<AltSiteRecord> altDesignMatrixForContext,
double convergenceThreshold,
int maxEMIterations,
int maxDepth,
org.apache.logging.log4j.Logger logger)
Contract: the reference contexts must be combined with its reverse complements prior to instantiating this class
|
Modifier and Type | Method and Description |
---|---|
static double[] |
computeResponsibilities(Nucleotide refAllele,
Nucleotide altAllele,
int altDepth,
int f1r2AltCount,
int depth,
double[] artifactPrior,
boolean givenNotHomRef)
Return normalized probabilities
|
double[] |
getAltDepth1Resonsibilities(int rowNum) |
double[] |
getAltResonsibilities(int rowNum) |
org.apache.commons.math3.linear.RealVector |
getEffectiveCounts() |
double |
getEffectiveCounts(ArtifactState state) |
static double[] |
getFlatPrior(Nucleotide refAllele) |
double[] |
getRefResonsibilities(int rowNum) |
ArtifactPrior |
learnPriorForArtifactStates() |
public LearnReadOrientationModelEngine(htsjdk.samtools.util.Histogram<java.lang.Integer> refHistogram, java.util.List<htsjdk.samtools.util.Histogram<java.lang.Integer>> altDepthOneHistograms, java.util.List<AltSiteRecord> altDesignMatrixForContext, double convergenceThreshold, int maxEMIterations, int maxDepth, org.apache.logging.log4j.Logger logger)
public ArtifactPrior learnPriorForArtifactStates()
public static double[] computeResponsibilities(Nucleotide refAllele, Nucleotide altAllele, int altDepth, int f1r2AltCount, int depth, double[] artifactPrior, boolean givenNotHomRef)
public double[] getRefResonsibilities(int rowNum)
public double[] getAltResonsibilities(int rowNum)
public double[] getAltDepth1Resonsibilities(int rowNum)
public org.apache.commons.math3.linear.RealVector getEffectiveCounts()
public double getEffectiveCounts(ArtifactState state)
public static double[] getFlatPrior(Nucleotide refAllele)