CalculateTheta {Jacquard} | R Documentation |
Calculate Relatedness Paramenters from Jacquard Coefficients tA Capitalized Title (ideally limited to 65 characters)
Description
Function CalculateTheta
calculates five identifiable relatedness parameters from a
list structure containing nine matrices of Jacquard coefficients.
Usage
CalculateTheta(J, ind.sub = 1:nrow(J[[1]]))
Arguments
J |
A list structure with nine matrices of pairwise Jacquard coefficients |
ind.sub |
Index for subsetting the individuals |
Details
CalculateTheta
produces four matrices according to the expressions in
Graffelman et al. (2024).
Value
A list object with four fields:
T1 |
The pairwise coancestry or kinship coefficients (symmetric) |
T2 |
The pairwise inbreeding coefficients (non-symmetric) |
T3 |
The pairwise Least one IBD out of three (symmetric) |
T4 |
T4 (skew-symmetric) |
Author(s)
Jan Graffelman (jan.graffelman@upc.edu)
References
Csuros, M. (2014) Non-identifiability of identity coefficients at biallelic loci. Theoretical Population Biology 92, pp. 22-29. doi:10.1016/j.tpb.2013.11.001.
Graffelman, J., Weir, B.S. and Goudet, J. (2024) Estimation of Jacquard's genetic identity coefficients with bi-allelic variants by constrained least-squares. Preprint at bioRxiv doi:10.1101/2024.03.25.586682.
See Also
Examples
data(DeltaSimulatedPedigree)
Theta <- CalculateTheta(DeltaSimulatedPedigree)