Description Usage Arguments Details Value Author(s) See Also Examples
Estimates by simulation the IBD coefficients of a non-inbred pairwise relationship.
1 | oneLocusIBD(x, ind1, ind2, Nsim, Xchrom=FALSE, verbose=TRUE, ...)
|
x |
A pedigree in the form of a |
ind1, ind2 |
Numeric ID labels of the two individuals. |
Nsim |
The number of simulations to be performed. |
Xchrom |
A logical indicating if the locus is X-linked (if TRUE) or autosomal (FALSE). |
verbose |
A logical. |
... |
Further arguments to be passed on to |
For any pair of non-inbred individuals, the IBD coefficients κ=(κ[0], κ[1], κ[2]) associated with the relationship, are defined as the probabilities
κ[i] = Pr(i alleles shared identically by descent).
For an X-chromosomal locus, and if at least one of the individuals is male, κ[i] is defined only for i=0,1.
A numeric of length 3 (autosomal) or 2 (X-linked), estimating κ.
Magnus Dehli Vigeland
twoLocusIBD
, oneLocusJacquard
, twoLocusJacquard
1 2 3 4 5 6 7 8 9 | ### Example 1: Full siblings
x <- nuclearPed(2)
Nsim <- 100 # Should be increased substantially
# Autosomal kappa estimate (exact = c(0.25, 0.5, 0.25))
oneLocusIBD(x, ind1=3, ind2=4, Nsim=Nsim)
# X-chromosomal kappa estimate (exact = c(0.5, 0.5))
oneLocusIBD(x, ind1=3, ind2=4, Nsim=Nsim, Xchrom=TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.