high_logl_pairs: return the high-sib-pair-logl indivs for each indiv

Description Usage Arguments

Description

I wrote this to see how much faster I could implement the functionality of find_high_logl_sib_pairs using Rcpp. find_high_logl_sib_pairs is painfully slow, and I think this should be much faster

Usage

1
high_logl_pairs(FSP, UPF, G, loglV)

Arguments

FSP

a 9 x L matrix (3 x 3 x L) of expected genotype frequencies of a full sibling pair, but we pass it in is a vector.

UPF

a 9 x L matrix of expected genotype frequencies of an unrelated pair, which we also pass in as a vector.

G

an L x N matrix of 0, 1, or 2 or NA giving the genotypes of the individuals

loglV

the cutoff point above which you will accept the pairs.


eriqande/fullsniplings documentation built on May 16, 2019, 8:45 a.m.