Sample_Haplotypes: Choose random block of haplotypes for theoretical power

In lin-lab/CEPSKAT: Finding rare variant associations for extreme phenotype samples.

Usage

Sample_Haplotypes(SNP.location, Region.length, n)

Arguments

SNP.location
Region.length
n

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (SNP.location, Region.length, n) 
{
    MAX.loc = max(SNP.location)
    MIN.loc = min(SNP.location)
    Start.loc = floor(runif(1) * (MAX.loc - MIN.loc - Region.length) + 
        MIN.loc)
    End.loc = Start.loc + Region.length
    Hap.ind = intersect(which(SNP.location > Start.loc), which(SNP.location < 
        End.loc))
    return(Hap.ind)
  }

lin-lab/CEPSKAT documentation built on May 29, 2019, 3:41 a.m.