gemma2 is an implementation in R of the GEMMAv 0.97 EM algorithm that is part of the GEMMA algorithm for REML estimation of multivariate linear mixed effects models of the form:
where E is a n by 2 matrix of random effects that follows the matrix-variate normal distribution
G ∼ M**N(0, K, Vg)
where K is a relatedness matrix and Vg is a 2 by 2 covariance matrix for the two traits of interest.
Additionally, the random errors matrix E follows the distribution:
E ∼ M**N(0, In, Ve)
and G and E are independent.
gemma2, use the
devtools R package from CRAN. If you haven't installed
devtools, please run this line of code:
Then, run this line of code to install