oneLocusJacquard: Estimating Jacquard's condensed identity coefficients

In IBDsim: Simulation of Chromosomal Regions Shared by Family Members

Description

Estimates by simulation Jacquard's 9 condensed identity coefficients for a pairwise relationship. This function is rarely needed, as exact values can be obtained by using jacquard.

Usage

oneLocusJacquard(x, ind1, ind2, Nsim, 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.
verbose	A logical.
...	Further arguments to be passed on to IBDsim.

Details

For the definition and further details about these coefficients, see Jacquard (1970).

Value

A numeric of length 9, estimating the condensed Jacquard identity coefficients Δ.

Author(s)

Magnus Dehli Vigeland

See Also

jacquard, oneLocusIBD, twoLocusIBD, twoLocusJacquard

Examples

### Siblings whose parents are full siblings.
x = fullSibMating(generations=2)
Nsim = 100 # (increase to improve accuracy) 

# Estimating the 9 identity coefficients
j_est = oneLocusJacquard(x, ind1=5, ind2=6, Nsim=Nsim)

# Exact: c(2,1,4,1,4,1,7,10,2)/32

# With the "identity" package:
## Not run: 
  j_exact = jacquard(x, 5:6)

## End(Not run)

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