In this paper, a new multivariate counting process model (called Multivariate Poisson Generalized Gamma Process) is developed and its main properties are studied. Some basic stochastic properties of the number of events in the new multivariate counting process are initially derived. It is shown that this new multivariate counting process model includes the multivariate generalized Pólya process as a special case. The dependence structure of the multivariate counting process model is discussed. Some results on multivariate stochastic comparisons are also obtained.

Many wear processes used for modeling accumulative deterioration in a reliability context are nonhomogeneous Lévy processes and, hence, have independent increments, which may not be suitable in an application context. In this work we consider Lévy processes transformed by monotonous functions to overcome this restriction, and provide a new state-dependent wear model. These transformed Lévy processes are first observed to remain tractable Markov processes. Some distributional properties are derived. We investigate the impact of the current state on the future increment level and on the overall accumulated level from a stochastic monotonicity point of view. We also study positive dependence properties and stochastic monotonicity of increments.

In this paper we extend some recent results on the comparison of multivariate risk vectors with respect to supermodular and related orderings. We introduce a dependence notion called the ‘weakly conditional increasing in sequence order’ that allows us to conclude that ‘more dependent’ vectors in this ordering are also comparable with respect to the supermodular ordering. At the same time, this ordering allows us to compare two risks with respect to the directionally convex order if the marginals increase convexly. We further state comparison criteria with respect to the directionally convex order for some classes of risk vectors which are modelled by functional influence factors. Finally, we discuss Fréchet bounds with respect to Δ-monotone functions when multivariate marginals are given. It turns out that, in the case of multivariate marginals, comonotone vectors no longer yield necessarily the largest risks but, in some cases, may even be vectors which minimize risk.

Simple approximation techniques are developed exploiting relationships between generalized convex orders and appropriate probability metrics. In particular, the distance between s-convex ordered random variables is investigated. Results connecting positive or negative dependence concepts and convex ordering are also presented. These results lead to approximations and bounds for the distributions of sums of positively or negatively dependent random variables. Applications and extensions of the main results pertaining to compound Poisson, normal and exponential approximation are provided as well.

Univariate probability inequalities have received extensive attention. It has been shown that under certain conditions, product-type bounds are valid and sharper than summation-type bounds. Although results concerning multivariate inequalities have appeared in the literature, product-type bounds in a multivariate setting have not yet been studied. This note explores an approach using graph theory and linear programming techniques to construct product-type lower bounds for the probability of the intersection among unions of k sets of events.

We will state a general version of Simpson's paradox, which corresponds to the loss of some dependence properties under marginalization. We will then provide conditions under which the paradox is avoided. Finally we will relate these Simpson-type paradoxes to some well-known paradoxes concerning the loss of ageing properties when the level of information changes.

Let P(c) = P(X1 ≦ c1, · ··, Xp ≦ cp) for a random vector (X1, · ··, Xp). Bounds are considered of the form where T is a spanning tree corresponding to the bivariate probability structure and di is the degree of the vertex i in T. An optimized version of this inequality is obtained. The main result is that alwayṡ dominates certain second-order Bonferroni bounds. Conditions on the covariance matrix of a N(0,Σ) distribution are given so that this bound applies, and various applications are given.