Package vcf
Class RefGT
java.lang.Object
vcf.RefGT
- All Implemented Interfaces:
GT
Class RefGT stores a list of samples and a
haplotype pair for each sample.
Instances of class BasicRefGT are immutable.
-
Constructor Summary
Constructors -
Method Summary
Modifier and TypeMethodDescriptionintallele(int marker, int haplotype) Returns the allele on the specified haplotype for the specified marker or return -1 if the allele is missing.get(int marker) Returns theRefGTRecfor the specified marker.booleanisPhased()Returnstrueif the genotype for each marker and sample has non-missing alleles and is either haploid or diploid with a phased allele separator, and returnsfalseotherwise.booleanReturnstrueif the markers are ordered by decreasing chromosome base position, and returnsfalseotherwise.marker(int marker) Returns the specified marker.markers()Returns the list of markers in order of increasing chromosome position.intnHaps()Returns the number of haplotypes.intnMarkers()Returns the number of markers.intnSamples()Returns the number of samples.restrict(int start, int end) Returns a newGTinstance restricted to genotype data for the specified markers.Returns aGTinstance restricted to genotype data for the specified markers.static RefGTReturns aRefGTinstance restricted to genotype data for the specified markers.samples()Returns the list of samples.
-
Constructor Details
-
RefGT
Constructs a newRefGTinstance.- Parameters:
markers- the sequence of markerssamples- the sequence of samplesrefVcfRecs- the sequence of per-marker genotype data- Throws:
IllegalArgumentException- ifmarkers.nMarkers() != refVcfRecs.lengthIllegalArgumentException- ifrefVcfRecs[k].samples().equals(samples) == falsefor anyksatisfying0 <= k && k < refVcfRecs.lengthIllegalArgumentException- ifrefVcfRecs[k].marker().equals(markers.marker(k)) == falsefor anyksatisfying0 <= k && k < refVcfRecs.lengthIllegalArgumentException- ifrefVcfRecs[k].isPhased() == falsefor anyksatisfying0 <= k && k < refVcfRecs.lengthNullPointerException- ifmarkers == null || samples == null || refVcfRecs == null || refVcfRecs[k] == nullfor anyksatisfying0 <= k && k <= refVcfRecs.length
-
RefGT
Constructs a newRefHapPairsinstance.- Parameters:
refVcfRecs- the sequence of per-marker genotype data- Throws:
IllegalArgumentException- ifrefVcfRecs.length == 0IllegalArgumentException- ifrefVcfRecs[k].samples().equals(samples) == falsefor anyksatisfying0 <= k && k < refVcfRecs.lengthIllegalArgumentException- ifrefVcfRecs[k].isPhased() == falsefor anyksatisfying0 <= k && k < refVcfRecs.lengthNullPointerException- ifsamples == null || refVcfRecs == nullNullPointerException- if(refVcfRecs[k] == null)for anyksatisfying(0 <= k && k <= refVcfRecs.length)
-
-
Method Details
-
isReversed
public boolean isReversed()Description copied from interface:GTReturnstrueif the markers are ordered by decreasing chromosome base position, and returnsfalseotherwise.- Specified by:
isReversedin interfaceGT- Returns:
trueif the markers are ordered by decreasing chromosome base position
-
nMarkers
public int nMarkers()Description copied from interface:GTReturns the number of markers. -
marker
Description copied from interface:GTReturns the specified marker. -
markers
Description copied from interface:GTReturns the list of markers in order of increasing chromosome position. If(this.isReversed() == false)then(this.markers().marker(j).equals(this.marker(j)) == true)for all(0 <= j && j < this.nMarkers()). If(this.isReversed() == true)then(this.markers().marker(this.nMarkers() - 1 - j).equals(this.marker(j)) == true)for all(0 <= j && j < this.nMarkers()) -
nHaps
public int nHaps()Description copied from interface:GTReturns the number of haplotypes. The returned value is equal to2*this.nSamples(). -
nSamples
public int nSamples()Description copied from interface:GTReturns the number of samples. -
samples
Description copied from interface:GTReturns the list of samples. -
isPhased
public boolean isPhased()Description copied from interface:GTReturnstrueif the genotype for each marker and sample has non-missing alleles and is either haploid or diploid with a phased allele separator, and returnsfalseotherwise. -
allele
public int allele(int marker, int haplotype) Description copied from interface:GTReturns the allele on the specified haplotype for the specified marker or return -1 if the allele is missing. The order of the two alleles is unspecified ifthis.isPhased() == false. -
restrict
Returns aRefGTinstance restricted to genotype data for the specified markers.- Parameters:
refGT- theRefGTinstance to be restrictedindices- a list of distinct marker indices (fromthis.markers())in increasing order- Returns:
- a
RefGTinstance restricted to genotype data for the specified markers - Throws:
IndexOutOfBoundsException- if there existsjsuch that(0 <= j && j < indices.length)such that(indices[j] < 0 || indices[j] >= gt.nMarkers())IllegalArgumentException- if there existsjsuch that(1 <= j && j < indices.length)such that(indices[j] <= indice[j - 1])NullPointerException- ifgt == null || indices == null
-
restrict
Description copied from interface:GTReturns aGTinstance restricted to genotype data for the specified markers. -
restrict
Description copied from interface:GTReturns a newGTinstance restricted to genotype data for the specified markers. -
get
Returns theRefGTRecfor the specified marker.- Parameters:
marker- the marker index- Returns:
- the
RefGTRecfor the specified marker - Throws:
IndexOutOfBoundsException- ifmarker < 0 || marker >= this.nMarkers()
-