gen.simuSampleFreq {GENLIB}

Gene dropping simulations - Frequencies

Description

Returns the number of alleles transmitted to specified probands from ancestors

Usage

gen.simuSampleFreq( gen, pro, ancestors, stateAncestors, simulNo = "5000")

Arguments

gen

An object of class GLgen obtained with gen.genealogy, gen.lineages or gen.branching. Required.

pro

Vector of proband id numbers to be included. Required.

ancestors

Vector of ancestors id numbers to be included. Required.

stateAncestors

Required vector indicating, for each ancestor in ancestors, the genotype state:

0 = no disease allele present

1 = 1 disease allele present

2 = 2 disease alleles present

simulNo

Number of simulations to perform. Default is 5000.

Value

A data frame containing for each row (one row per proband): Alleles.transmitted.0 Alleles.transmitted.1 Alleles.transmitted.2

See Also

gen.genealogy gen.simuProb gen.simuSet gen.simuSample

Examples

data(geneaJi) 
genJi<-gen.genealogy(geneaJi) 
gen.simuSampleFreq(genJi, pro=c(1,29), ancestors=c(20,25),stateAncestors=c(2,1), simulNo = 1000)