A multivariate multilevel Gaussian model with a mixed effects structure in the mean and covariance part