heteroZygotTestCpp function

Function to check whether a window is (loosely) heterozygous or not

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)

Arguments

  • x: vector of 0/1 genotypes (from genoConvert())
  • gaps: vector of differences between consecutive positions (gaps) in bps
  • maxHom: max n. of homozygous SNP in a heterozygous window
  • maxMiss: max n. of missing in a window
  • maxGap: max distance between consecutive SNP in a window to be still considered a potential run

Returns

TRUE/FALSE (whether a window is heterozygous or NOT)