Learn R Programming

cpgen (version 0.1)

cgrm.A: Additive Genomic Relationship Matrix

Description

Based on a marker matrix $\mathbf{X}$ with {-1,0,1} - coding that will be centered column-wise and a shrinkage parameter $\lambda$, cgrm.A returns the following additive genomic relationship matrix according to VanRaden (2008): $$\mathbf{G} = (1-\lambda) \frac{\mathbf{X X}^{'}}{\sum\limits_{i=1}^n 2 p_i q_i} + \mathbf{I}\lambda$$

Usage

cgrm.A(X, lambda=0, yang=FALSE)

Arguments

X
marker matrix
lambda
numeric scalar, shrinkage parameter
yang
boolean, diagonal elements of A according to Yang et al. (2010)

Value

  • Additive genomic relationship matrix with dimension nrow(X)

Details

...

References

VanRaden, P.M. "Efficient Methods to Compute Genomic Predictions". Journal of Dairy Science 91, no. 11 (November 2008): 4414-23. doi:10.3168/jds.2007-0980.

Yang, Jian, Beben Benyamin, Brian P McEvoy, Scott Gordon, Anjali K Henders, Dale R Nyholt, Pamela A Madden, et al. "Common SNPs Explain a Large Proportion of the Heritability for Human Height". Nature Genetics 42, no. 7 (July 2010): 565-69. doi:10.1038/ng.608.

See Also

cgrm, cgrm.D

Examples

Run this code
# generate random data
rand_data(500,5000)

### compute the additive genomic relationship matrix
A <- cgrm.A(M,lambda=0.01)

Run the code above in your browser using DataLab