Adjusted maximum likelihood method for multivariate Fay–Herriot model