We implement an expectation-maximization algorithm for multivariate variance components after the GEMMA software's algorithm.