Package phase
Class MarkerCluster
java.lang.Object
phase.MarkerCluster
Class MarkerCluster
represents a partition of markers into
contiguous marker clusters.
Instances of class MarkerCluster
are immutable.
-
Constructor Summary
ConstructorsConstructorDescriptionMarkerCluster
(PhaseData phaseData, int sample) Constructs a newMarkerCluster
instance from the specified data. -
Method Summary
Modifier and TypeMethodDescriptionint
clusterEnd
(int index) Returns the exclusive end marker for the cluster.int
clusterStart
(int index) Returns the inclusive start marker for the cluster.boolean
isHet
(int cluster) Returnstrue
if the cluster is a heterozygous genotype, and returnsfalse
otherwise.boolean
isMaskedHet
(int cluster) Returnstrue
if the cluster is a masked heterozygote genotype, and returnsfalse
otherwise.boolean
isMissingGT
(int cluster) Returnstrue
if the cluster has a missing genotype, and returnsfalse
otherwise.boolean
isMissingGtOrMaskedHet
(int cluster) Returnstrue
if the cluster has a missing genotype or a masked heterozygote genotype, and returnsfalse
otherwise.boolean
isPhasedHet
(int cluster) Returnstrue
if the cluster has a phased heterozygous genotype, and returnsfalse
otherwise.boolean
isUnphasedHet
(int cluster) Returnstrue
if the cluster has an unphased heterozygous genotype, and returnsfalse
otherwise.int
Returns the number of clusterspRecomb()
Return aFloatArray
of sizethis.nClusters()
whosek
-th element is the probability of transitioning to a random HMM state between thek
-th cluster and the previous cluster.Return the estimated haplotypes.Returns a sorted list of cluster indices in increasing order for which the cluster contains an unphased heterozygote.
-
Constructor Details
-
MarkerCluster
Constructs a newMarkerCluster
instance from the specified data.- Parameters:
phaseData
- the input data for the next genotype phasing iterationsample
- a sample index- Throws:
IndexOutOfBoundsException
- ifsample < 0 || sample >= phaseData.targGT().nSamples()
NullPointerException
- ifphaseData == null
-
-
Method Details
-
samplePhase
Return the estimated haplotypes.- Returns:
- the estimated haplotypes
-
nClusters
public int nClusters()Returns the number of clusters- Returns:
- the number of clusters
-
clusterStart
public int clusterStart(int index) Returns the inclusive start marker for the cluster.- Parameters:
index
- a cluster index- Returns:
- the inclusive start marker for the cluster
- Throws:
IndexOutOfBoundsException
- ifindex < 0 || index >= this.nClusteres()
-
clusterEnd
public int clusterEnd(int index) Returns the exclusive end marker for the cluster.- Parameters:
index
- a cluster index- Returns:
- the exclusive marker for the cluster
- Throws:
IndexOutOfBoundsException
- ifindex < 0 || index >= this.nClusteres()
-
pRecomb
Return aFloatArray
of sizethis.nClusters()
whosek
-th element is the probability of transitioning to a random HMM state between thek
-th cluster and the previous cluster.- Returns:
- a
FloatArray
of sizethis.nClusters()
whosek
-th element is the probability of transitioning to a random HMM state between thek
-th cluster and the previous cluster
-
unphasedHetClusters
Returns a sorted list of cluster indices in increasing order for which the cluster contains an unphased heterozygote.- Returns:
- a sorted list of cluster indices in increasing order for which the cluster contains an unphased heterozygote
-
isUnphasedHet
public boolean isUnphasedHet(int cluster) Returnstrue
if the cluster has an unphased heterozygous genotype, and returnsfalse
otherwise.- Parameters:
cluster
- a cluster index- Returns:
true
if the cluster has an unphased heterozygous genotype, and returnsfalse
otherwise- Throws:
IndexOutOfBoundsException
- ifcluster < 0 || cluster >= this.nClusters()
-
isPhasedHet
public boolean isPhasedHet(int cluster) Returnstrue
if the cluster has a phased heterozygous genotype, and returnsfalse
otherwise.- Parameters:
cluster
- a cluster index- Returns:
true
if the cluster has a phased heterozygous genotype, and returnsfalse
otherwise- Throws:
IndexOutOfBoundsException
- ifcluster < 0 || cluster >= this.nClusters()
-
isHet
public boolean isHet(int cluster) Returnstrue
if the cluster is a heterozygous genotype, and returnsfalse
otherwise.- Parameters:
cluster
- a cluster index- Returns:
true
if the cluster is a heterozygous genotype, andfalse
otherwise- Throws:
IndexOutOfBoundsException
- ifcluster < 0 || cluster >= this.nClusters()
-
isMissingGT
public boolean isMissingGT(int cluster) Returnstrue
if the cluster has a missing genotype, and returnsfalse
otherwise.- Parameters:
cluster
- a cluster index- Returns:
true
if the cluster has a missing genotype- Throws:
IndexOutOfBoundsException
- ifcluster < 0 || cluster >= this.nClusters()
-
isMaskedHet
public boolean isMaskedHet(int cluster) Returnstrue
if the cluster is a masked heterozygote genotype, and returnsfalse
otherwise.- Parameters:
cluster
- a cluster index- Returns:
true
if the cluster is a masked heterozygote genotype- Throws:
IndexOutOfBoundsException
- iftcluster < 0 || cluster >= this.nClusters()
-
isMissingGtOrMaskedHet
public boolean isMissingGtOrMaskedHet(int cluster) Returnstrue
if the cluster has a missing genotype or a masked heterozygote genotype, and returnsfalse
otherwise.- Parameters:
cluster
- a cluster index- Returns:
true
if the cluster has a missing genotype or a masked heterozygote genotype, and returnsfalse
otherwise- Throws:
IndexOutOfBoundsException
- ifcluster < 0 || cluster >= this.nClusters()
-