Residual maximum likelihood estimation of (co)variance components in multivariate mixed linear models using average information