gen.simuProb: Gene dropping simulations - Probabilities

Description Usage Arguments Value See Also Examples

Description

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

Usage

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

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.

statePro

Required vector indicating, for each proband in pro, the probability to be calculated:

0 = no disease allele is transmitted

1 = 1 disease allele is transmitted

2 = 2 disease alleles are transmitted

3 = 1 or 2 disease alleles are transmitted

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.

probRecomb

Recombination probabilities for males and females. Default is no recombination.

probSurvival

Survival probability for homozygotes. Default is 1.

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
data(geneaJi) 
genJi<-gen.genealogy(geneaJi) 
## Not run: Probability that subjects 1 and 29 get 1 and 2 alleles from ancestors 20 and 25,
## Not run:  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)
## Not run: Probability that subjects 1 and 29 get 1 segment from ancestors 25, 
## Not run:  knowing the segment has a male 
## Not run:  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))
## Not run: Probability that subjects 1 and 29 get 1 and 2 alleles from ancestors 20 and 25, 
## Not run:  that have themselves 2 and 1 and knowing that homozygous people have a survival rate 
## Not run:  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)


Search within the GENLIB package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.