Learn R Programming
GENLIB (version 1.0.3)

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)

Arguments

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

Run this code
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)
Number of alleles received by probants
table(simu_1000)
Number of alleles received by simulation
table(colSums(simu_1000))
Number of alleles received by each probant
table(simu_1000[1,],simu_1000[2,])

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)
Number of alleles received by probants
table(simu_5000)
Number of alleles received by simulation
table(colSums(simu_5000))
Number of alleles received by each probant
table(simu_5000[1,],simu_5000[2,])

Run the code above in your browser using DataLab