We study the tail asymptotics of two functionals (the maximum and the sum of the marks) of a generic cluster in two sub-models of the marked Poisson cluster process, namely the renewal Poisson cluster process and the Hawkes process. Under the hypothesis that the governing components of the processes are regularly varying, we extend results due to [6, 19], notably relying on Karamata’s Tauberian Theorem to do so. We use these asymptotics to derive precise large-deviation results in the fashion of [32] for the just-mentioned processes.

Motivated by applications to COVID dynamics, we describe a model of a branching process in a random environment $\{Z_n\}$ whose characteristics change when crossing upper and lower thresholds. This introduces a cyclical path behavior involving periods of increase and decrease leading to supercritical and subcritical regimes. Even though the process is not Markov, we identify subsequences at random time points $\{(\tau_j, \nu_j)\}$—specifically the values of the process at crossing times, viz. $\{(Z_{\tau_j}, Z_{\nu_j})\}$—along which the process retains the Markov structure. Under mild moment and regularity conditions, we establish that the subsequences possess a regenerative structure and prove that the limiting normal distributions of the growth rates of the process in supercritical and subcritical regimes decouple. For this reason, we establish limit theorems concerning the length of supercritical and subcritical regimes and the proportion of time the process spends in these regimes. As a byproduct of our analysis, we explicitly identify the limiting variances in terms of the functionals of the offspring distribution, threshold distribution, and environmental sequences.

Corrections are made to formulations and proofs of some theorems about convolution equivalence closure for random sum distributions. These arise because of the falsity of a much used asymptotic equivalence lemma, and they impinge on the convolution equivalence closure theorem for general infinitely divisible laws.

Let {X k , k ≥ 1} be a sequence of independent, identically distributed nonnegative random variables with common distribution function F and finite expectation μ > 0. Under the assumption that the tail probability is consistently varying as x tends to infinity, this paper investigates precise large deviations for both the partial sums S n and the random sums S N(t), where N(·) is a counting process independent of the sequence {X k , k ≥ 1}. The obtained results improve some related classical ones. Applications to a risk model with negatively associated claim occurrences and to a risk model with a doubly stochastic arrival process (extended Cox process) are proposed.

In this paper, we investigate asymptotic properties of the tail probabilities of the maxima of partial sums of independent random variables. For some large classes of heavy-tailed distributions, we show that the tail probabilities of the maxima of the partial sums asymptotically equal to the sum of the tail probabilities of the individual random variables. Then we partially extend the result to the case of random sums. Applications to some commonly used risk processes are proposed. All heavy-tailed distributions involved in this paper are supposed on the whole real line.

In this paper we study random variables related to a shock reliability model. Our models can be used to study systems that fail when k consecutive shocks with critical magnitude (e.g. above or below a certain critical level) occur. We obtain properties of the distribution function of the random variables involved and we obtain their limit behaviour when k tends to infinity or when the probability of entering a critical set tends to zero. This model generalises the Poisson shock model.

Numerical evaluation of compound distributions is one of the central numerical tasks in insurance mathematics. Two widely used techniques are Panjer recursion and transform methods. Many authors have pointed out that aliasing errors imply the need to consider the whole distribution if transform methods are used, a potential drawback especially for heavy-tailed distributions. We investigate the magnitude of aliasing errors and show that this problem can be solved by a suitable change of measure.

We derive upper bounds for the total variation distance, d, between the distributions of two random sums of non-negative integer-valued random variables. The main results are then applied to some important random sums, including cluster binomial and cluster multinomial distributions, to obtain bounds on approximating them to suitable Poisson or compound Poisson distributions. These bounds are generally better than the known results on Poisson and compound Poisson approximations. We also obtain a lower bound for d and illustrate it with an example.

In this paper we investigate operators that map one or more probability distributions on the positive real line into another via their Laplace–Stieltjes transforms. Our goal is to make it easier to construct new transforms by manipulating known transforms. We envision the results here assisting modelling in conjunction with numerical transform inversion software. We primarily focus on operators related to infinitely divisible distributions and Lévy processes, drawing upon Feller (1971). We give many concrete examples of infinitely divisible distributions. We consider a cumulant-moment-transfer operator that allows us to relate the cumulants of one distribution to the moments of another. We consider a power-mixture operator corresponding to an independently stopped Lévy process. The special case of exponential power mixtures is a continuous analog of geometric random sums. We introduce a further special case which is remarkably tractable, exponential mixtures of inverse Gaussian distributions (EMIGs). EMIGs arise naturally as approximations for busy periods in queues. We show that the steady-state waiting time in an M/G/1 queue is the difference of two EMIGs when the service-time distribution is an EMIG. We consider several transforms related to first-passage times, e.g. for the M/M/1 queue, reflected Brownian motion and Lévy processes. Some of the associated probability density functions involve Bessel functions and theta functions. We describe properties of the operators, including how they transform moments.