GENLIB (version 1.0.1)

gen.simuSample: Gene dropping simulations - Sample

Description

Returns the number of alleles transmitted to specified probands from ancestors

Usage

gen.simuSample(gen, pro, ancestors, stateAncestors, simulNo = 5000)

Value

A matrix with number of columns equal to the number of simulations and number of rows equal to the number of probands.

See Also

gen.genealogy gen.simuProb gen.simuSet gen.simuSampleFreq

Examples

data(geneaJi) 
genJi<-gen.genealogy(geneaJi) 
simu_1000<-gen.simuSample(genJi, pro=c(1,29), ancestors=c(20,25), stateAncestors=c(2,1),
                          simulNo=1000)
table(simu_1000) #Number of alleles received by proposants
table(colSums(simu_1000)) #Number of alleles received by simulation
table(simu_1000[1,],simu_1000[2,]) #Number of alleles received by each proposant


data(geneaJi) 
genJi<-gen.genealogy(geneaJi) 
simu_5000<-gen.simuSample(genJi, pro=c(1,29), ancestors=c(20,25), stateAncestors=c(2,1),
                          simulNo=5000)
table(simu_5000) #Number of alleles received by proposants
table(colSums(simu_5000)) #Number of alleles received by simulation
table(simu_5000[1,],simu_5000[2,]) #Number of alleles received by each proposant