Calculates the variance covariance matrix for a multivariate normal vector when there are random effects.
Useful links