Multivariate regression with covariance estimation