Problems concerning estimation of parameters and determination the statistic, when it is known a priori that some of these parameters are subject to certain order restrictions, are of considerable interest. In the present paper, we consider the estimators of the monotonic mean vectors for two dimensional normal distributions and compare those with the unrestricted maximum likelihood estimators under two different cases. One case is that covariance matrices are known, the other one is that covariance matrices are completely unknown and unequal. We show that when the covariance matrices are known, under the squared error loss function which is similar to the mahalanobis distance, the obtained multivariate isotonic regression estimators, motivated by estimators given in Robertson et al. (1988), which are the estimators given by Sasabuchi et al. (1983) and Sasabuchi et al. (1992), have the smaller risk than the unrestricted maximum likelihood estimators uniformly, but when the covariance matrices are unknown and unequal, the estimators have the smaller risk than the unrestricted maximum likelihood estimators only over some special sets which are defined on the covariance matrices. To illustrate the results two numerical examples are presented.