Stochastic comparisons for multivariate shock models