estimating covariance matrix for saturated multivariate ACE