IBDsim (version 0.9-8)

oneLocusJacquard: Estimating Jacquard's condensed identity coefficients

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.

Value

A numeric of length 9, estimating the condensed Jacquard identity coefficients \(\Delta\).

Details

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

See Also

jacquard, oneLocusIBD, twoLocusIBD, twoLocusJacquard

Examples

# NOT RUN {
### 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)
# }