A method for the construction of Stein-type covariance identities for a nonnegative continuous random variable is proposed, using a probabilistic analogue of the mean value theorem and weighted distributions. A generalized covariance identity is obtained, and applications focused on actuarial and financial science are provided. Some characterization results for gamma and Pareto distributions are also given. Identities for risk measures which have a covariance representation are obtained; these measures are connected with the Bonferroni, De Vergottini, Gini, and Wang indices. Moreover, under some assumptions, an identity for the variance of a function of a random variable is derived, and its performance is discussed with respect to well-known upper and lower bounds.

This paper investigates the ordering properties of largest claim amounts in heterogeneous insurance portfolios in the sense of some transform orders, including the convex transform order and the star order. It is shown that the largest claim amount from a set of independent and heterogeneous exponential claims is more skewed than that from a set of independent and homogeneous exponential claims in the sense of the convex transform order. As a result, a lower bound for the coefficient of variation of the largest claim amount is established without any restrictions on the parameters of the distributions of claim severities. Furthermore, sufficient conditions are presented to compare the skewness of the largest claim amounts from two sets of independent multiple-outlier scaled claims according to the star order. Some comparison results are also developed for the multiple-outlier proportional hazard rates claims. Numerical examples are presented to illustrate these theoretical results.

In this paper we study the allocation problem of relevations in coherent systems. The optimal allocation strategies are obtained by implementing stochastic comparisons of different policies according to the usual stochastic order and the hazard rate order. As special cases of relevations, the load-sharing and minimal repair policies are further investigated. Sufficient (and necessary) conditions are established for various stochastic orderings. Numerical examples are also presented as illustrations.

We consider a class of phase-type distributions (PH-distributions), to be called the MMPP class of PH-distributions, and find bounds of their mean and squared coefficient of variation (SCV). As an application, we have shown that the SCV of the event-stationary inter-event time for Markov modulated Poisson processes (MMPPs) is greater than or equal to unity, which answers an open problem for MMPPs. The results are useful for selecting proper PH-distributions and counting processes in stochastic modeling.

Using the calculus of variations, we prove the following structure theorem for noise-stable partitions: a partition of n-dimensional Euclidean space into m disjoint sets of fixed Gaussian volumes that maximise their noise stability must be $(m-1)$-dimensional, if $m-1\leq n$. In particular, the maximum noise stability of a partition of m sets in $\mathbb {R}^{n}$ of fixed Gaussian volumes is constant for all n satisfying $n\geq m-1$. From this result, we obtain:

1. (i) A proof of the plurality is stablest conjecture for three candidate elections, for all correlation parameters $\rho $ satisfying $0<\rho <\rho _{0}$, where $\rho _{0}>0$ is a fixed constant (that does not depend on the dimension n), when each candidate has an equal chance of winning.

2. (ii) A variational proof of Borell’s inequality (corresponding to the case $m=2$).

The structure theorem answers a question of De–Mossel–Neeman and of Ghazi–Kamath–Raghavendra. Item (i) is the first proof of any case of the plurality is stablest conjecture of Khot-Kindler-Mossel-O’Donnell for fixed $\rho $, with the case $\rho \to L1^{-}$ being solved recently. Item (i) is also the first evidence for the optimality of the Frieze–Jerrum semidefinite program for solving MAX-3-CUT, assuming the unique games conjecture. Without the assumption that each candidate has an equal chance of winning in (i), the plurality is stablest conjecture is known to be false.

Copula is one of the widely used techniques to describe the dependency structure between components of a system. Among all existing copulas, the family of Archimedean copulas is the popular one due to its wide range of capturing the dependency structures. In this paper, we consider the systems that are formed by dependent and identically distributed components, where the dependency structures are described by Archimedean copulas. We study some stochastic comparisons results for series, parallel, and general $r$-out-of-$n$ systems. Furthermore, we investigate whether a system of used components performs better than a used system with respect to different stochastic orders. Furthermore, some aging properties of these systems have been studied. Finally, some numerical examples are given to illustrate the proposed results.

This paper considers logarithmic asymptotics of tails of randomly stopped sums. The stopping is assumed to be independent of the underlying random walk. First, finiteness of ordinary moments is revisited. Then the study is expanded to more general asymptotic analysis. Results are applicable to a large class of heavy-tailed random variables. The main result enables one to identify if the asymptotic behaviour of a stopped sum is dominated by its increments or the stopping variable. As a consequence, new sufficient conditions for the moment determinacy of compounded sums are obtained.

We consider a birth–death process with killing where transitions from state i may go to either state $i-1$ or state $i+1$ or an absorbing state (killing). Stochastic ordering results on the killing time are derived. In particular, if the killing rate in state i is monotone in i, then the distribution of the killing time with initial state i is stochastically monotone in i. This result is a consequence of the following one for a non-negative tri-diagonal matrix M: if the row sums of M are monotone, so are the row sums of $M^n$ for all $n\ge 2$ .

In this paper we consider a new generalized finite mixture model formed by dependent and identically distributed (d.i.d.) components. We then establish results for the comparisons of lifetimes of two such generalized finite mixture models in two different cases: (i) when the two mixture models are formed from two random vectors $\textbf{X}$ and $\textbf{Y}$ but with the same weights, and (ii) when the two mixture models are formed with the same random vectors but with different weights. Because the lifetimes of k-out-of-n systems and coherent systems are special cases of the mixture model considered, we used the established results to compare the lifetimes of k-out-of-n systems and coherent systems with respect to the reversed hazard rate and hazard rate orderings.

Let $\{Y_{1},\ldots ,Y_{n}\}$ be a collection of interdependent nonnegative random variables, with $Y_{i}$ having an exponentiated location-scale model with location parameter $\mu _i$, scale parameter $\delta _i$ and shape (skewness) parameter $\beta _i$, for $i\in \mathbb {I}_{n}=\{1,\ldots ,n\}$. Furthermore, let $\{L_1^{*},\ldots ,L_n^{*}\}$ be a set of independent Bernoulli random variables, independently of $Y_{i}$'s, with $E(L_{i}^{*})=p_{i}^{*}$, for $i\in \mathbb {I}_{n}.$ Under this setup, the portfolio of risks is the collection $\{T_{1}^{*}=L_{1}^{*}Y_{1},\ldots ,T_{n}^{*}=L_{n}^{*}Y_{n}\}$, wherein $T_{i}^{*}=L_{i}^{*}Y_{i}$ represents the $i$th claim amount. This article then presents several sufficient conditions, under which the smallest claim amounts are compared in terms of the usual stochastic and hazard rate orders. The comparison results are obtained when the dependence structure among the claim severities are modeled by (i) an Archimedean survival copula and (ii) a general survival copula. Several examples are also presented to illustrate the established results.

There are two types of tempered stable (TS) based Ornstein–Uhlenbeck (OU) processes: (i) the OU-TS process, the OU process driven by a TS subordinator, and (ii) the TS-OU process, the OU process with TS marginal law. They have various applications in financial engineering and econometrics. In the literature, only the second type under the stationary assumption has an exact simulation algorithm. In this paper we develop a unified approach to exactly simulate both types without the stationary assumption. It is mainly based on the distributional decomposition of stochastic processes with the aid of an acceptance–rejection scheme. As the inverse Gaussian distribution is an important special case of TS distribution, we also provide tailored algorithms for the corresponding OU processes. Numerical experiments and tests are reported to demonstrate the accuracy and effectiveness of our algorithms, and some further extensions are also discussed.

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.

Two-sided bounds are explored for concentration functions and Rényi entropies in the class of discrete log-concave probability distributions. They are used to derive certain variants of the entropy power inequalities.

Pareto distribution is an important distribution in extreme value theory. In this paper, we consider parallel systems with Pareto components and study the effect of heterogeneity on skewness of such systems. It is shown that, when the lifetimes of components have different shape parameters, the parallel system with heterogeneous Pareto component lifetimes is more skewed than the system with independent and identically distributed Pareto components. However, for the case when the lifetimes of components have different scale parameters, the result gets reversed in the sense of star ordering. We also establish the relation between star ordering and dispersive ordering by extending the result of Deshpande and Kochar [(1983). Dispersive ordering is the same as tail ordering. Advances in Applied Probability 15(3): 686–687] from support $(0, \infty )$ to general supports $(a, \infty )$, $a > 0$. As a consequence, we obtain some new results on dispersion of order statistics from heterogeneous Pareto samples with respect to dispersive ordering.

We show that the sequence of moments of order less than 1 of averages of i.i.d. positive random variables is log-concave. For moments of order at least 1, we conjecture that the sequence is log-convex and show that this holds eventually for integer moments (after neglecting the first $p^2$ terms of the sequence).

We prove the sharp bound for the probability that two experts who have access to different information, represented by different $\sigma$-fields, will give radically different estimates of the probability of an event. This is relevant when one combines predictions from various experts in a common probability space to obtain an aggregated forecast. The optimizer for the bound is explicitly described. This paper was originally titled ‘Contradictory predictions’.

We consider a space-time random field on ${{\mathbb{R}^d} \times {\mathbb{R}}}$ given as an integral of a kernel function with respect to a Lévy basis with a convolution equivalent Lévy measure. The field obeys causality in time and is thereby not continuous along the time axis. For a large class of such random fields we study the tail behaviour of certain functionals of the field. It turns out that the tail is asymptotically equivalent to the right tail of the underlying Lévy measure. Particular examples are the asymptotic probability that there is a time point and a rotation of a spatial object with fixed radius, in which the field exceeds the level x, and that there is a time interval and a rotation of a spatial object with fixed radius, in which the average of the field exceeds the level x.

Extended gamma processes have been seen as a flexible extension of standard gamma processes in the recent reliability literature, for the purpose of cumulative deterioration modeling. The probabilistic properties of the standard gamma process have been well explored since the 1970s, whereas those of its extension remain largely unexplored. In particular, stochastic comparisons between degradation levels modeled by standard gamma processes and ageing properties for the corresponding level-crossing times are now well understood. The aim of this paper is to explore similar properties for extended gamma processes and see which ones can be broadened to this new context. As a by-product, new stochastic comparisons for convolutions of gamma random variables are also obtained.

Comparison results for Markov processes with respect to function-class-induced (integral) stochastic orders have a long history. The most general results so far for this problem have been obtained based on the theory of evolution systems on Banach spaces. In this paper we transfer the martingale comparison method, known for the comparison of semimartingales to Markovian semimartingales, to general Markov processes. The basic step of this martingale approach is the derivation of the supermartingale property of the linking process, giving a link between the processes to be compared. This property is achieved using the characterization of Markov processes by the associated martingale problem in an essential way. As a result, the martingale comparison method gives a comparison result for Markov processes under a general alternative but related set of regularity conditions compared to the evolution system approach.

We give a new method of proof for a result of D. Pierre-Loti-Viaud and P. Boulongne which can be seen as a generalization of a characterization of Poisson law due to Rényi and Srivastava. We also provide explicit formulas, in terms of Bell polynomials, for the moments of the compound distributions occurring in the extended collective model in non-life insurance.