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]