Calculate the phenotypes covariance matrix Sigma
calc_Sigma(Vg, Ve, kinship = NULL, n_mouse = nrow(kinship))
d by d genetic covariance matrix for the d phenotypes
d by d error covariance matrix for the d phenotypes
optional n by n kinship matrix. if NULL, Vg is not used.
number of subjects
dn by dn covariance matrix