MakeM {Jacquard} | R Documentation |
Create matrix of conditional joint genotype probabilities
Description
Function MakeM
creates the matrix of conditional joint genotype probabilities for biallelic markers for a given allele frequency.
Usage
MakeM(p)
Arguments
p |
the allele frequency |
Value
a 9 by 9 matrix
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. Under review.
Examples
set.seed(123)
p <- runif(1)
M <- MakeM(p)
[Package Jacquard version 1.0.2 Index]