For two n-dimensional elliptical random vectors X and Y, we establish an identity for $\mathbb{E}[f({\bf Y})]- \mathbb{E}[f({\bf X})]$, where $f\,{:}\, \mathbb{R}^n \rightarrow \mathbb{R}$ satisfies some regularity conditions. Using this identity we provide a unified method to derive sufficient and necessary conditions for classifying multivariate elliptical random vectors according to several main integral stochastic orders. As a consequence we obtain new inequalities by applying the method to multivariate elliptical distributions. The results generalize the corresponding ones for multivariate normal random vectors in the literature.

We investigate several aspects of a self-similar evolutionary process that builds a random bipolar network from building blocks that are themselves small bipolar networks. We characterize admissible outdegrees in the history of the evolution. We obtain the limit distribution of the polar degrees (when suitably scaled) characterized by its sequence of moments. We also obtain the asymptotic joint multivariate normal distribution of the number of nodes of small admissible outdegrees. Five possible substructures arise, and each has its own parameters (mean vector and covariance matrix) in the multivariate distribution. Several results are obtained by mapping bipolar networks into Pólya urns.

In Reinert and Röllin (2009) a new approach - called the ‘embedding method’ - was introduced, which allows us to make use of exchangeable pairs for normal and multivariate normal approximations with Stein's method in cases where the corresponding couplings do not satisfy a certain linearity condition. The key idea is to embed the problem into a higher-dimensional space in such a way that the linearity condition is then satisfied. Here we apply the embedding to U-statistics as well as to subgraph counts in random graphs.

The problem of finding bounds for P(A1 ∪ ⋯ ∪ An) based on P(Ak1 ∩ ⋯ ∩ Aki) (1 ≤ k1 < ⋯ < ki ≤ n, i = 1,…,d) goes back to Boole (1854), (1868) and Bonferroni (1937). In this paper upper bounds are presented using methods in graph theory. The main theorem is a common generalization of the earlier results of Hunter, Worsley and recent results of Prékopa and the author. Algorithms are given to compute bounds. Examples for bounding values of multivariate normal distribution functions are presented.

We consider the convex ordering for random vectors and some weaker versions of it, like the convex ordering for linear combinations of random variables. First we establish conditions of stochastic equality for random vectors that are ordered by one of the convex orderings. Then we establish necessary and sufficient conditions for the convex ordering to hold in the case of multivariate normal distributions and sufficient conditions for the positive linear convex ordering (without the restriction to multi-normality).

The problem of bounding P(∪ Ai) given P(Ai) and P(AiAj) for i ≠ j = 1, …, k goes back to Boole (1854) and Bonferroni (1936). In this paper a new family of upper bounds is derived using results in graph theory. This family contains the bound of Kounias (1968), and the smallest upper bound in the family for a given application is easily derivable via the minimal spanning tree algorithm of Kruskal (1956). The properties of the algorithm and of the multivariate normal and t distributions are shown to provide considerable simplifications when approximating tail probabilities of maxima from these distributions.