An Identity for Expectations and Characteristic Function of Matrix Variate Skew-normal Distribution with Applications to Associated Stochastic Orderings

We establish an identity for EfY-EfX$$Ef\left( \varvec{Y}\right) - Ef\left( \varvec{X}\right) $$, when X$$\varvec{X}$$ and Y$$\varvec{Y}$$ both have matrix variate skew-normal distributions and the function f satisfies some weak conditions. The characteristic function of matrix variate skew normal distribution is then derived. We then make use of it to derive some necessary and sufficient conditions for the comparison of matrix variate skew-normal distributions under six different orders, such as usual stochastic order, convex order, increasing convex order, upper orthant order, directionally convex order and supermodular order.