gemma2 uses restricted maximum likelihood methods to infer variance components in a multivariate linear mixed effects model for QTL mapping.qtl2pleio performs a d-variate, d-QTL scan over a genomic region of interest before calculating a likelihood ratio test statistic for the competing hypotheses of pleiotropy and separate QTL.