This could be done entirely in R, but it is hard to think it through and it seemed it would be hard to maintain. Super easy in C.
1 | C_full_sibling_pair_gfreqs(L, G, PP, TP, GE)
|
L |
the number of loci |
G |
the number of genotypic states |
PP |
parent pair probs as returned by unrelated_pair_gfreqs() for example |
TP |
transmision probs as returned by trans_probs() |
GE |
prob of observed genotypes given true genotypes, as returned by lik_array_from_simple_geno_err() for example |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.