Given that the parents are unrelated and have expected genotype frequencies given in the 3 x L array gf, and given genotypin error rates specified for each locus mu (will recycle) this function returns the probablities of the 9 possible genotypes (at all L loci) of a pair of full siblings in the population. This can be done entirely in R but it would be hard to come back to it and understand what it is doing, so I am going to do much of it with a call to an Rcpp function
1  | full_sibling_pair_gfreqs(gf, mu = 0)
 | 
gf | 
 array of genotype frequencies  | 
mu | 
 genotyping error rates  | 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.