DEPRECATED. I keep this around to remind me how slow R is. It is better to use high_logl_pairs(), implemented in Rcpp, in practice.
1 | find_high_logl_sib_pairs(FSP, UPF, GI, loglV)
|
FSP |
full-sibling pair geno freqs as returned by full_sibling_pair_gfreqs() |
UPF |
unrelated pair geno freqs as would be returned by unrelated_pair_gfreqs() |
GI |
genotypes of all the individuals as a snp_indics array |
loglV |
the log-likelihood value that must be exceeded to be included |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.