GENLIB (version 1.0.3)

gen.simuProb: Gene dropping simulations - Probabilities

Description

Returns the probabilities that specified probands inherit disease alleles from ancestors.

Usage

gen.simuProb(gen, pro, statePro, ancestors, stateAncestors, simulNo=5000, 
                    probRecomb=c(0,0), probSurvival=1.0)

Value

A list containing the following: the joint probability of specified statePro for all probands, the probability of specified statePro for each proband, the probability that, 0, 1, ..., and all probands inherit the specified number of disease alleles.

See Also

gen.genealogy gen.simuSample gen.simuSet gen.simuSampleFreq

Examples

data(geneaJi) 
genJi<-gen.genealogy(geneaJi) 
Probability that subjects 1 and 29 get 1 and 2 alleles from ancestors 20 and 25,
that have themselves 2 and 1.
gen.simuProb(genJi, pro=c(1,29), statePro=c(1,2), ancestors=c(20,25), stateAncestors=c(2,1),
             simulNo=10000)
Probability that subjects 1 and 29 get 1 segment from ancestors 25,
knowing the segment has a male
recombination rate of 0.02 and a female recombination rate of 0.04.
gen.simuProb(genJi, pro=c(1,29), statePro=c(1,1), ancestors=c(25), stateAncestors=c(1),
             simulNo=10000, probRecomb = c(0.02, 0.04))
Probability that subjects 1 and 29 get 1 and 2 alleles from ancestors 20 and 25,
that have themselves 2 and 1 and knowing that homozygous people have a survival rate
of 0.50.
gen.simuProb(genJi, pro=c(1,29), statePro=c(1,2), ancestors=c(20,25), stateAncestors=c(2,1),
             simulNo=10000, probSurvival=0.5)