Copulas provide a powerful and flexible tool for modeling the dependence structure of random vectors, and they have many applications in finance, insurance, engineering, hydrology, and other fields. One well-known class of copulas in two dimensions is the Farlie–Gumbel–Morgenstern (FGM) copula, since its simple analytic shape enables closed-form solutions to many problems in applied probability. However, the classical definition of the high-dimensional FGM copula does not enable a straightforward understanding of the effect of the copula parameters on the dependence, nor a geometric understanding of their admissible range. We circumvent this issue by analyzing the FGM copula from a probabilistic approach based on multivariate Bernoulli distributions. This paper examines high-dimensional exchangeable FGM copulas, a subclass of FGM copulas. We show that the dependence parameters of exchangeable FGM copulas can be expressed as a convex hull of a finite number of extreme points. We also leverage the probabilistic interpretation to develop efficient sampling and estimating procedures and provide a simulation study. Throughout, we discover geometric interpretations of the copula parameters that assist one in decoding the dependence of high-dimensional exchangeable FGM copulas.

Suppose that a system is affected by a sequence of random shocks that occur over certain time periods. In this paper we study the discrete censored $\delta$ -shock model, $\delta \ge 1$ , for which the system fails whenever no shock occurs within a $\delta$ -length time period from the last shock, by supposing that the interarrival times between consecutive shocks are described by a first-order Markov chain (as well as under the binomial shock process, i.e., when the interarrival times between successive shocks have a geometric distribution). Using the Markov chain embedding technique introduced by Chadjiconstantinidis et al. (Adv. Appl. Prob. 32, 2000), we study the joint and marginal distributions of the system’s lifetime, the number of shocks, and the number of periods in which no shocks occur, up to the failure of the system. The joint and marginal probability generating functions of these random variables are obtained, and several recursions and exact formulae are given for the evaluation of their probability mass functions and moments. It is shown that the system’s lifetime follows a Markov geometric distribution of order $\delta$ (a geometric distribution of order $\delta$ under the binomial setup) and also that it follows a matrix-geometric distribution. Some reliability properties are also given under the binomial shock process, by showing that a shift of the system’s lifetime random variable follows a compound geometric distribution. Finally, we introduce a new mixed discrete censored $\delta$ -shock model, for which the system fails when no shock occurs within a $\delta$ -length time period from the last shock, or the magnitude of the shock is larger than a given critical threshold $\gamma >0$ . Similarly, for this mixed model, we study the joint and marginal distributions of the system’s lifetime, the number of shocks, and the number of periods in which no shocks occur, up to the failure of the system, under the binomial shock process.

In this paper, we consider a financial or insurance system with a finite number of individual risks described by real-valued random variables. We focus on two kinds of risk measures, referred to as the tail moment (TM) and the tail central moment (TCM), which are defined as the conditional moment and conditional central moment of some individual risk in the event of system crisis. The first-order TM and the second-order TCM coincide with the popular risk measures called the marginal expected shortfall and the tail variance, respectively. We derive asymptotic expressions for the TM and TCM with any positive integer orders, when the individual risks are pairwise asymptotically independent and have distributions from certain classes that contain both light-tailed and heavy-tailed distributions. The formulas obtained possess concise forms unrelated to dependence structures, and hence enable us to estimate the TM and TCM efficiently. To demonstrate the wide application of our results, we revisit some issues related to premium principles and optimal capital allocation from the asymptotic point of view. We also give a numerical study on the relative errors of the asymptotic results obtained, under some specific scenarios when there are two individual risks in the system. The corresponding asymptotic properties of the degenerate univariate versions of the TM and TCM are discussed separately in an appendix at the end of the paper.

In this paper, we establish some stochastic comparison results for largest claim amounts of two sets of independent and also for interdependent portfolios under the setup of the proportional odds model. We also establish stochastic comparison results for aggregate claim amounts of two sets of independent portfolios. Further, stochastic comparisons for largest claim amounts from two sets of independent multiple-outlier claims have also been studied. The results we obtained apply to the whole family of extended distributions, also known as the Marshall–Olkin family of distributions. We have given many numerical examples to illustrate the results obtained.

The principle of maximum entropy is a well-known approach to produce a model for data-generating distributions. In this approach, if partial knowledge about the distribution is available in terms of a set of information constraints, then the model that maximizes entropy under these constraints is used for the inference. In this paper, we propose a new three-parameter lifetime distribution using the maximum entropy principle under the constraints on the mean and a general index. We then present some statistical properties of the new distribution, including hazard rate function, quantile function, moments, characterization, and stochastic ordering. We use the maximum likelihood estimation technique to estimate the model parameters. A Monte Carlo study is carried out to evaluate the performance of the estimation method. In order to illustrate the usefulness of the proposed model, we fit the model to three real data sets and compare its relative performance with respect to the beta generalized Weibull family.

We examine urn models under random replacement schemes, and the related distributions, by using generating functions. A fundamental isomorphism between urn models and a certain system of differential equations has previously been established. We study the joint distribution of the numbers of balls in the urn, and determined recurrence relations for the probability generating functions. The associated partial differential equation satisfied by the generating function is derived. We develop analytical methods for the study of urn models that can lead to perspectives on urn-related problems from analytic combinatorics. The results presented here provide a broader framework for the study of exactly solvable urns than the existing framework. Finally, we examine several applications and their numerical results in order to demonstrate how our theoretical results can be employed in the study of urn models.

Consider the coupon collector problem where each box of a brand of cereal contains a coupon and there are n different types of coupons. Suppose that the probability of a box containing a coupon of a specific type is $1/n$, and that we keep buying boxes until we collect at least m coupons of each type. For $k\geq m$ call a certain coupon a k-ton if we see it k times by the time we have seen m copies of all of the coupons. Here we determine the asymptotic distribution of the number of k-tons after we have collected m copies of each coupon for any k in a restricted range, given any fixed m. We also determine the asymptotic joint probability distribution over such values of k, and the total number of coupons collected.

We propose a stochastic model for the failure times of items subject to two external random shocks occurring as events in an underlying bivariate counting process. This is a special formulation of the competing risks model, which is of interest in reliability theory and survival analysis. Specifically, we assume that a system, or an item, fails when the sum of the two types of shock reaches a critical random threshold. In detail, the two kinds of shock occur according to a bivariate space-fractional Poisson process, which is a two-dimensional vector of independent homogeneous Poisson processes time-changed by an independent stable subordinator. Various results are given, such as analytic hazard rates, failure densities, the probability that the failure occurs due to a specific type of shock, and the survival function. Some special cases and ageing notions related to the NBU characterization are also considered. In this way we generalize certain results in the literature, which can be recovered when the underlying process reduces to the homogeneous Poisson process.

Let f be the density function associated to a matrix-exponential distribution of parameters $(\boldsymbol{\alpha}, T,\boldsymbol{{s}})$ . By exponentially tilting f, we find a probabilistic interpretation which generalizes the one associated to phase-type distributions. More specifically, we show that for any sufficiently large $\lambda\ge 0$ , the function $x\mapsto \left(\int_0^\infty e^{-\lambda s}f(s)\textrm{d} s\right)^{-1}e^{-\lambda x}f(x)$ can be described in terms of a finite-state Markov jump process whose generator is tied to T. Finally, we show how to revert the exponential tilting in order to assign a probabilistic interpretation to f itself.

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 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 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.

We introduce a multivariate class of distributions with support I, a k-orthotope in $[0,\infty)^{k}$ , which is dense in the set of all k-dimensional distributions with support I. We call this new class ‘multivariate finite-support phase-type distributions’ (MFSPH). Though we generally define MFSPH distributions on any finite k-orthotope in $[0,\infty)^{k}$ , here we mainly deal with MFSPH distributions with support $[0,1)^{k}$ . The distribution function of an MFSPH variate is computed by using that of a variate in the MPH $^{*} $ class, the multivariate class of distributions introduced by Kulkarni (1989). The marginal distributions of MFSPH variates are found as FSPH distributions, the class studied by Ramaswami and Viswanath (2014). Some properties, including the mixture property, of MFSPH distributions are established. Estimates of the parameters of a particular class of bivariate finite-support phase-type distributions are found by using the expectation-maximization algorithm. Simulated samples are used to demonstrate how this class could be used as approximations for bivariate finite-support distributions.

In the literature of stochastic orders, one rarely finds results characterizing non-comparability of random variables. We prove simple tools implying the non-comparability with respect to the convex transform order. The criteria are used, among other applications, to provide a negative answer for a conjecture about comparability in a much broader scope than its initial statement.

For independent exponentially distributed random variables $X_i$ , $i\in {\mathcal{N}}$ , with distinct rates ${\lambda}_i$ we consider sums $\sum_{i\in\mathcal{A}} X_i$ for $\mathcal{A}\subseteq {\mathcal{N}}$ which follow generalized exponential mixture distributions. We provide novel explicit results on the conditional distribution of the total sum $\sum_{i\in {\mathcal{N}}}X_i$ given that a subset sum $\sum_{j\in \mathcal{A}}X_j$ exceeds a certain threshold value $t>0$ , and vice versa. Moreover, we investigate the characteristic tail behavior of these conditional distributions for $t\to\infty$ . Finally, we illustrate how our probabilistic results can be applied in practice by providing examples from both reliability theory and risk management.

We investigate stochastic comparisons of parallel systems (corresponding to the largest-order statistics) with respect to the reversed hazard rate and likelihood ratio orders for the proportional reversed hazard rate (PRHR) model. As applications of the main results, we obtain the equivalent characterizations of stochastic comparisons with respect to the reversed hazard rate and likelihood rate orders for the exponentiated generalized gamma and exponentiated Pareto distributions. Our results recover and strengthen some recent results in the literature.

We study failure rate monotonicity and generalised convex transform stochastic ordering properties of random variables, with an emphasis on applications. We are especially interested in the effect of a tail-weight iteration procedure to define distributions, which is equivalent to the characterisation of moments of the residual lifetime at a given instant. For the monotonicity properties, we are mainly concerned with hereditary properties with respect to the iteration procedure providing counterexamples showing either that the hereditary property does not hold or that inverse implications are not true. For the stochastic ordering, we introduce a new criterion, based on the analysis of the sign variation of a suitable function. This criterion is then applied to prove ageing properties of parallel systems formed with components that have exponentially distributed lifetimes.

Large samples from a light-tailed distribution often have a well-defined shape. This paper examines the implications of the assumption that there is a limit shape. We show that the limit shape determines the upper quantiles for a large class of random variables. These variables may be described loosely as continuous homogeneous functionals of the underlying random vector. They play an important role in evaluating risk in a multivariate setting. The paper also looks at various coefficients of tail dependence and at the distribution of the scaled sample points for large samples. The paper assumes convergence in probability rather than almost sure convergence. This results in an elegant theory. In particular, there is a simple characterization of domains of attraction.

The main result of this note implies that any function from the product of several vector spaces to a vector space can be uniquely decomposed into the sum of mutually orthogonal functions that are odd in some of the arguments and even in the other arguments. Probabilistic notions and facts are employed to simplify statements and proofs.