View source: R/Elston-algorithm.r
Elston | R Documentation |
Compute the probability of a pedigree, for a parameter theta
and a given modele
Elston(ped, modele, theta, mem = new.env())
ped |
a pedigree, created by |
modele |
a modele |
theta |
a parameter for the modele |
mem |
an environment used for memoization |
This function runs Elston-Stewart algorithm. If mem
is provided, some intermediate results of previous runs stored in it can be re-used.
A list with two components, result
: the probability value, and mem
: an object storing intermediate results for future runs
Likelihood
, es.pedigree
## cf Elston-Stewart vignette for more coments on this example data(conrad2) # creating an es.pedigree object genotypes <- c( rep(list(0:2), 21), 2 ) X <- es.pedigree( id = conrad2$id, father = conrad2$father, mother = conrad2$mother, sex = conrad2$sex, pheno = rep(0, 22), geno = genotypes ) # running Elston-Stewart r <- Elston(X, modele.di, list(p = 0.98)) r$result # using the memoization... r <- Elston(X, modele.di, list(p = 0.99), r$mem) r$result
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.