Imputation of missing covariates under a multivariate linear mixed model