sim_pheno_lmm: Simulate bivariate phenotype with incorporation of random...
In fboehm/pleio: Likelihood ratio test to distinguish pleiotropy from linkage

Simulate bivariate phenotype with incorporation of random effects

1	sim_pheno_lmm(genomat, bmat, Sigma)

`genomat`	a single locus's genotypes probability matrix
`bmat`	a matrix of effect sizes
`Sigma`	a 2n by 2n covariance matrix

Kmat <- diag(1, 8)
Vg <- matrix(c(1:4 / 10), nrow = 2, ncol = 2)
Ve <- matrix(c(2:5) / 10, nrow = 2, ncol = 2)
In <- diag(1, 8)
sim_pheno_lmm(genomat = matrix(c(1, 0, 0, 1, 2, 1, 0, 1), ncol = 1), 
bmat = matrix(c(2.5, 3), nrow = 1), 
Sigma = Vg %x% Kmat + Ve %x% In)