This is a core function. It implements the consecutive method for detection of runs in diploid genomes (see Marras et al. 2015)
consecutiveRunsCpp(indGeno, individual, mapFile, ROHet = TRUE,
minSNP = 3L, maxOppositeGenotype = 1L, maxMiss = 1L,
minLengthBps = 1000L, maxGap = 1000000L)vector of 0/1/NAs of individual genotypes (0: homozygote; 1: heterozygote)
list of group (breed, population, case/control etc.) and ID of individual sample
Plink map file (for SNP position)
shall we detect ROHet or ROHom?
minimum number of SNP in a run
max n. of homozygous/heterozygous SNP
max. n. of missing SNP
min length of a run in bps
max distance between consecutive SNP in a window to be still considered a potential run
A data frame of runs per individual sample
The consecutive method detect runs by consecutively scanning SNP loci along the genome. No sliding windows are used. Checks on minimum n. of SNP, max n. of opposite and missing genotypes, max gap between adjacent loci and minimum length of the run are implemented (as in the sliding window method). Both runs of homozygosity (RoHom) and of heterozygosity (RoHet) can be search for (option ROHet: TRUE/FALSE)