fit1_bvlmm: Fit a 2-locus, multivariate trait linear mixed effects model
In fboehm/pleio: Likelihood ratio test to distinguish pleiotropy from linkage

Fit a 2-locus, multivariate trait linear mixed effects model

1 2	fit1_bvlmm(Y, X1, X2, Kmat, tol = 1e-07, maxcyc = 100, verbose = 10, start = c(0.1, 0.1, 0, 1, 1, 0), taper = rep(1/2, maxcyc))

`Y`	a phenotypes matrix with one row per subject and one trait per column
`X1`	genotype probabilities for first locus, one row per subject
`X2`	genotype probabilities for second locus, one row per subject
`Kmat`	a kinship matrix
`tol`	value of tol to be passed to regress()
`maxcyc`	value of maxcyc passed to regress()
`verbose`	value of verbose passed to regress()
`start`	value of start passed to regress()
`taper`	value of taper passed to regress()