Stochastic Ordering of Multivariate Normal Distributions

We show an interesting identity for Ef(Y) − Ef(X), where X, Yare normally distributed random vectors and f is a function fulfilling some weak regularity condition. This identity will be used for a unified derivation of sufficient conditions for stochastic ordering results of multivariate normal distributions, some well known ones as well as some new ones. Moreover, we will show that many of these conditions are also necessary. As examples we will consider the usual stochastic order, convex order, upper orthant order, supermodular order and directionally convex order.