We define the generalized equilibrium distribution, that is the equilibrium distribution of a random variable with support in $\mathbb{R}$. This concept allows us to prove a new probabilistic generalization of Taylor’s theorem. Then, the generalized equilibrium distribution of two ordered random variables is considered and a probabilistic analog of the mean value theorem is proved. Results regarding distortion-based models and mean-median-mode relations are illustrated as well. Conditions for the unimodality of such distributions are obtained. We show that various stochastic orders and aging classes are preserved through the proposed equilibrium transformations. Further applications are provided in actuarial science, aiming to employ the new unimodal equilibrium distributions for some risk measures, such as Value-at-Risk and Conditional Tail Expectation.

Exchangeable partitions of the integers and their corresponding mass partitions on $\mathcal{P}_{\infty} = \{\mathbf{s} = (s_{1},s_{2},\ldots)\colon s_{1} \geq s_{2} \geq \cdots \geq 0$ and $\sum_{k=1}^{\infty}s_{k} = 1\}$ play a vital role in combinatorial stochastic processes and their applications. In this work, we continue our focus on the class of Gibbs partitions of the integers and the corresponding stable Poisson–Kingman-distributed mass partitions generated by the normalized jumps of a stable subordinator with an index $\alpha \in (0,1)$, subject to further mixing. This remarkable class of infinitely exchangeable random partitions is characterized by probabilities that have Gibbs (product) form. These partitions have practical applications in combinatorial stochastic processes, random tree/graph growth models, and Bayesian statistics. The most notable class consists of random partitions generated from the two-parameter Poisson–Dirichlet distribution $\mathrm{PD}(\alpha,\theta)$. While the utility of Gibbs partitions is recognized, there is limited understanding of the broader class. Here, as a continuation of our work, we address this gap by extending the dual coagulation/fragmentation results of Pitman (1999), developed for the Poisson–Dirichlet family, to all Gibbs models and their corresponding Poisson–Kingman mass partitions, creating nested families of Gibbs partitions and mass partitions. We focus primarily on fragmentation operations, identifying which classes correspond to these operations and providing significant calculations for the resulting Gibbs partitions. Furthermore, for completion, we provide definitive results for dual coagulation operations using dependent processes. We demonstrate the applicability of our results by establishing new findings for Brownian excursion partitions, Mittag-Leffler, and size-biased generalized gamma models.

Sharp, nonasymptotic bounds are obtained for the relative entropy between the distributions of sampling with and without replacement from an urn with balls of $c\geq 2$ colors. Our bounds are asymptotically tight in certain regimes and, unlike previous results, they depend on the number of balls of each color in the urn. The connection of these results with finite de Finetti-style theorems is explored, and it is observed that a sampling bound due to Stam (1978) combined with the convexity of relative entropy yield a new finite de Finetti bound in relative entropy, which achieves the optimal asymptotic convergence rate.

In this paper we adopt the probabilistic mean value theorem in order to study differences of the variances of transformed and stochastically ordered random variables, based on a suitable extension of the equilibrium operator. We also develop a rigorous approach aimed at expressing the variance of transformed random variables. This is based on a joint distribution which, in turn, involves the variance of the original random variable, as well as its mean residual lifetime and mean inactivity time. Then we provide applications to the additive hazards model and to some well-known random variables of interest in actuarial science. These deal with a new notion, called the ‘centred mean residual lifetime’, and a suitably related stochastic order. Finally, we also address the analysis of the differences of the variances of transformed discrete random variables thanks to the use of a discrete version of the equilibrium operator.

The present paper develops a unified approach when dealing with short- or long-range dependent processes with finite or infinite variance. We are concerned with the convergence rate in the strong law of large numbers (SLLN). Our main result is a Marcinkiewicz–Zygmund law of large numbers for $S_{n}(f)= \sum_{i=1}^{n}f(X_{i})$, where $\{X_i\}_{i\geq 1}$ is a real stationary Gaussian sequence and $f(\!\cdot\!)$ is a measurable function. Key technical tools in the proofs are new maximal inequalities for partial sums, which may be useful in other problems. Our results are obtained by employing truncation alongside new maximal inequalities. The result can help to differentiate the effects of long memory and heavy tails on the convergence rate for limit theorems.

We consider a stochastic model, called the replicator coalescent, describing a system of blocks of k different types that undergo pairwise mergers at rates depending on the block types: with rate $C_{ij}\geq 0$ blocks of type i and j merge, resulting in a single block of type i. The replicator coalescent can be seen as a generalisation of Kingman’s coalescent death chain in a multi-type setting, although without an underpinning exchangeable partition structure. The name is derived from a remarkable connection between the instantaneous dynamics of this multi-type coalescent when issued from an arbitrarily large number of blocks, and the so-called replicator equations from evolutionary game theory. By dilating time arbitrarily close to zero, we see that initially, on coming down from infinity, the replicator coalescent behaves like the solution to a certain replicator equation. Thereafter, stochastic effects are felt and the process evolves more in the spirit of a multi-type death chain.

Competing and complementary risk (CCR) problems are often modelled using a class of distributions of the maximum, or minimum, of a random number of independent and identically distributed random variables, called the CCR class of distributions. While CCR distributions generally do not have an easy-to-calculate density or probability mass function, two special cases, namely the Poisson–exponential and exponential–geometric distributions, can easily be calculated. Hence, it is of interest to approximate CCR distributions with these simpler distributions. In this paper, we develop Stein’s method for the CCR class of distributions to provide a general comparison method for bounding the distance between two CCR distributions, and we contrast this approach with bounds obtained using a Lindeberg argument. We detail the comparisons for Poisson–exponential, and exponential–geometric distributions.

We show that $\alpha $-stable Lévy motions can be simulated by any ergodic and aperiodic probability-preserving transformation. Namely we show that: for $0<\alpha <1$ and every $\alpha $-stable Lévy motion ${\mathbb {W}}$, there exists a function f whose partial sum process converges in distribution to ${\mathbb {W}}$; for $1\leq \alpha <2$ and every symmetric $\alpha $-stable Lévy motion, there exists a function f whose partial sum process converges in distribution to ${\mathbb {W}}$; for $1< \alpha <2$ and every $-1\leq \beta \leq 1$ there exists a function f whose associated time series is in the classical domain of attraction of an $S_\alpha (\ln (2), \beta ,0)$ random variable.

Motivated by the investigation of probability distributions with finite variance but heavy tails, we study infinitely divisible laws whose Lévy measure is characterized by a radial component of geometric (tempered) stable type. We closely investigate the univariate case: characteristic exponents and cumulants are calculated, as well as spectral densities; absolute continuity relations are shown, and short- and long-time scaling limits of the associated Lévy processes analyzed. Finally, we derive some properties of the involved probability density functions.

We introduce a comprehensive method for establishing stochastic orders among order statistics in the independent and identically distributed case. This approach relies on the assumption that the underlying distribution is linked to a reference distribution through a transform order. Notably, this method exhibits broad applicability, particularly since several well-known nonparametric distribution families can be defined using relevant transform orders, including the convex and the star transform orders. Moreover, for convex-ordered families, we show that an application of Jensen’s inequality gives bounds for the probability that a random variable exceeds the expected value of its corresponding order statistic.

Conditional risk measures and their associated risk contribution measures are commonly employed in finance and actuarial science for evaluating systemic risk and quantifying the effects of risk interactions. This paper introduces various types of contribution ratio measures based on the multivariate conditional value-at-risk (MCoVaR), multivariate conditional expected shortfall (MCoES), and multivariate marginal mean excess (MMME) studied in [34] (Ortega-Jiménez, P., Sordo, M., & Suárez-Llorens, A. (2021). Stochastic orders and multivariate measures of risk contagion. Insurance: Mathematics and Economics, vol. 96, 199–207) and [11] (Das, B., & Fasen-Hartmann, V. (2018). Risk contagion under regular variation and asymptotic tail independence. Journal of Multivariate Analysis 165(1), 194–215) to assess the relative effects of a single risk when other risks in a group are in distress. The properties of these contribution risk measures are examined, and sufficient conditions for comparing these measures between two sets of random vectors are established using univariate and multivariate stochastic orders and statistically dependent notions. Numerical examples are presented to validate these conditions. Finally, a real dataset from the cryptocurrency market is used to analyze the spillover effects through our proposed contribution measures.

For a continuous-time phase-type (PH) distribution, starting with its Laplace–Stieltjes transform, we obtain a necessary and sufficient condition for its minimal PH representation to have the same order as its algebraic degree. To facilitate finding this minimal representation, we transform this condition equivalently into a non-convex optimization problem, which can be effectively addressed using an alternating minimization algorithm. The algorithm convergence is also proved. Moreover, the method we develop for the continuous-time PH distributions can be used directly for the discrete-time PH distributions after establishing an equivalence between the minimal representation problems for continuous-time and discrete-time PH distributions.

We propose Rényi information generating function (RIGF) and discuss its properties. A connection between the RIGF and the diversity index is proposed for discrete-type random variables. The relation between the RIGF and Shannon entropy of order q > 0 is established and several bounds are obtained. The RIGF of escort distribution is derived. Furthermore, we introduce the Rényi divergence information generating function (RDIGF) and discuss its effect under monotone transformations. We present nonparametric and parametric estimators of the RIGF. A simulation study is carried out and a real data relating to the failure times of electronic components is analyzed. A comparison study between the nonparametric and parametric estimators is made in terms of the standard deviation, absolute bias, and mean square error. We have observed superior performance for the newly proposed estimators. Some applications of the proposed RIGF and RDIGF are provided. For three coherent systems, we calculate the values of the RIGF and other well-established uncertainty measures, and similar behavior of the RIGF is observed. Further, a study regarding the usefulness of the RDIGF and RIGF as model selection criteria is conducted. Finally, three chaotic maps are considered and then used to establish a validation of the proposed information generating function.

For coherent systems with components and active redundancies having heterogeneous and dependent lifetimes, we prove that the lifetime of system with redundancy at component level is stochastically larger than that with redundancy at system level. In particular, in the setting of homogeneous components and redundancy lifetimes linked by an Archimedean survival copula, we develop sufficient conditions for the reversed hazard rate order, the hazard rate order and the likelihood ratio order between two system lifetimes, respectively. The present results substantially generalize some related results in the literature. Several numerical examples are presented to illustrate the findings as well.

In this paper, we define weighted failure rate and their means from the stand point of an application. We begin by emphasizing that the formation of n independent component series system having weighted failure rates with sum of weight functions being unity is same as a mixture of n distributions. We derive some parametric and non-parametric characterization results. We discuss on the form invariance property of baseline failure rate for a specific choice of weight function. Some bounds on means of aging functions are obtained. Here, we establish that weighted increasing failure rate average (IFRA) class is not closed under formation of coherent systems unlike the IFRA class. An interesting application of the present work is credited to the fact that the quantile version of means of failure rate is obtained as a special case of weighted means of failure rate.

The pth ($p\geq 1$) moment exponential stability, almost surely exponential stability and stability in distribution for stochastic McKean–Vlasov equation are derived based on some distribution-dependent Lyapunov function techniques.

A series of papers by Hickey (1982, 1983, 1984) presents a stochastic ordering based on randomness. This paper extends the results by introducing a novel methodology to derive models that preserve stochastic ordering based on randomness. We achieve this by presenting a new family of pseudometric spaces based on a majorization property. This class of pseudometrics provides a new methodology for deriving the randomness measure of a random variable. Using this, the paper introduces the Gini randomness measure and states its essential properties. We demonstrate that the proposed measure has certain advantages over entropy measures. The measure satisfies the value validity property, provides an adequate extension to continuous random variables, and is often more appropriate (based on sensitivity) than entropy in various scenarios.

Random bridges have gained significant attention in recent years due to their potential applications in various areas, particularly in information-based asset pricing models. This paper aims to explore the potential influence of the pinning point’s distribution on the memorylessness and stochastic dynamics of the bridge process. We introduce Lévy bridges with random length and random pinning points, and analyze their Markov property. Our study demonstrates that the Markov property of Lévy bridges depends on the nature of the distribution of their pinning points. The law of any random variables can be decomposed into singular continuous, discrete, and absolutely continuous parts with respect to the Lebesgue measure (Lebesgue’s decomposition theorem). We show that the Markov property holds when the pinning points’ law does not have an absolutely continuous part. Conversely, the Lévy bridge fails to exhibit Markovian behavior when the pinning point has an absolutely continuous part.

We propose a monotone approximation scheme for a class of fully nonlinear degenerate partial integro-differential equations which characterize nonlinear $\alpha$-stable Lévy processes under a sublinear expectation space with $\alpha\in(1,2)$. We further establish the error bounds for the monotone approximation scheme. This in turn yields an explicit Berry–Esseen bound and convergence rate for the $\alpha$-stable central limit theorem under sublinear expectation.

In this paper, we explore the applications of Tail Variance (TV) as a measure of tail riskiness and the confidence level of using Tail Conditional Expectation (TCE)-based risk capital. While TCE measures the expected loss of a risk that exceeds a certain threshold, TV measures the variability of risk along its tails. We first derive analytical formulas of TV and TCE for a large variety of probability distributions. These formulas are useful instruments for relevant research works on tail risk measures. We then propose a distribution-free approach utilizing TV to estimate the lower bounds of the confidence level of using TCE-based risk capital. In doing so, we introduce sharpened conditional probability inequalities, which halve the bounds of conventional Markov and Cantelli inequalities. Such an approach is easy to implement. We further investigate the characterization of tail risks by TV through an exploration of TV’s asymptotics. A distribution-free limit formula is derived for the asymptotics of TV. To further investigate the asymptotic properties, we consider two broad distribution families defined on tails, namely, the polynomial-tailed distributions and the exponential-tailed distributions. The two distribution families are found to exhibit an asymptotic equivalence between TV and the reciprocal square of the hazard rate. We also establish asymptotic relationships between TCE and VaR for the two families. Our asymptotic analysis contributes to the existing research by unifying the asymptotic expressions and the convergence rate of TV for Student-t distributions, exponential distributions, and normal distributions, which complements the discussion on the convergence rate of univariate cases in [28]. To show the usefulness of our results, we present two case studies based on real data from the industry. We first show how to use conditional inequalities to assess the confidence of using TCE-based risk capital for different types of insurance businesses. Then, for financial data, we provide alternative evidence for the relationship between the data frequency and the tail categorization by the asymptotics of TV.