oneLocusIBD: Estimating pairwise IBD coefficients

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Estimates by simulation the IBD coefficients of a non-inbred pairwise relationship.

Usage

1
oneLocusIBD(x, ind1, ind2, Nsim, Xchrom=FALSE, verbose=TRUE, ...)

Arguments

x

A pedigree in the form of a linkdat object.

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 IBDsim.

Details

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.

Value

A numeric of length 3 (autosomal) or 2 (X-linked), estimating κ.

Author(s)

Magnus Dehli Vigeland

See Also

twoLocusIBD, oneLocusJacquard, twoLocusJacquard

Examples

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)

IBDsim documentation built on May 2, 2019, 6:01 a.m.

Related to oneLocusIBD in IBDsim...