Function to check whether a window is (loosely) heterozygous or not
This is a core function. Parameters on how to consider a window heterozygous are here (maxHom, maxMiss)
heteroZygotTestCpp(x, gaps, maxHom, maxMiss, maxGap)
x
: vector of 0/1 genotypes (from genoConvert())gaps
: vector of differences between consecutive positions (gaps) in bpsmaxHom
: max n. of homozygous SNP in a heterozygous windowmaxMiss
: max n. of missing in a windowmaxGap
: max distance between consecutive SNP in a window to be still considered a potential runTRUE/FALSE (whether a window is heterozygous or NOT)
Useful links