maskTransitionSites: Masking of double recombinations at homozygous-heterozygous...

Description Usage Arguments Examples

Description

Sets short double recombination calls to missing, e.g. in the sequence HHAABB, AA is replaced by missing data, NN.

Usage

1
maskTransitionSites(inputGenos, maxHapLength = 1)

Arguments

inputGenos

ABH genotype file to mask.

maxHapLength

Max length (in number of markers) of double recombination haplotype to set to missing. Default is 1.

Examples

1
2
3
4
# replaces double recombination transitions up to 4 marker lengths
data(abhgeno)
g <- maskTransitionSites(abhgeno, 4)
g

johncarlosignacio/gtmask documentation built on Dec. 21, 2021, 2:12 a.m.