oneLocusJacquard {IBDsim}

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.

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 \Delta.

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)