A Multivariate Fay-Herriot model is used to aid the prediction of small area parameters of dependent variables with sample data aggregated to area level. The empirical best linear unbiased predictor of the parameter vector is used, and an approximation of the elements of the mean cross product...