In this work, we consider extensions of the dual risk model with proportional gains by introducing dependence structures among gain sizes and gain interarrival times. Among others, we further consider the case where the proportionality parameter is randomly chosen, the case where it is a uniformly random variable, as well as the case where we may have upward as well as downward jumps. Moreover, we consider the case with causal dependence structure, as well as the case where the dependence is based on the generalized Farlie–Gumbel–Morgenstern copula. The ruin probability and the distribution of the time to ruin are investigated.

We consider the following problem: the drift of the wealth process of two companies, modelled by a two-dimensional Brownian motion, is controllable such that the total drift adds up to a constant. The aim is to maximize the probability that both companies survive. We assume that the volatility of one company is small with respect to the other, and use methods from singular perturbation theory to construct a formal approximation of the value function. Moreover, we validate this formal result by explicitly constructing a strategy that provides a target functional, approximating the value function uniformly on the whole state space.

Multivariate regular variation is a key concept that has been applied in finance, insurance, and risk management. This paper proposes a new dependence assumption via a framework of multivariate regular variation. Under the condition that financial and insurance risks satisfy our assumption, we conduct asymptotic analyses for multidimensional ruin probabilities in the discrete-time and continuous-time cases. Also, we present a two-dimensional numerical example satisfying our assumption, through which we show the accuracy of the asymptotic result for the discrete-time multidimensional insurance risk model.

We consider the problem of controlling the drift and diffusion rate of the endowment processes of two firms such that the joint survival probability is maximized. We assume that the endowment processes are continuous diffusions, driven by independent Brownian motions, and that the aggregate endowment is a Brownian motion with constant drift and diffusion rate. Our results reveal that the maximal joint survival probability depends only on the aggregate risk-adjusted return and on the maximal risk-adjusted return that can be implemented in each firm. Here the risk-adjusted return is understood as the drift rate divided by the squared diffusion rate.

We consider a dual risk model with constant expense rate and i.i.d. exponentially distributed gains $C_i$ ( $i=1,2,\dots$ ) that arrive according to a renewal process with general interarrival times. We add to this classical dual risk model the proportional gain feature; that is, if the surplus process just before the ith arrival is at level u, then for $a>0$ the capital jumps up to the level $(1+a)u+C_i$ . The ruin probability and the distribution of the time to ruin are determined. We furthermore identify the value of discounted cumulative dividend payments, for the case of a Poisson arrival process of proportional gains. In the dividend calculations, we also consider a random perturbation of our basic risk process modeled by an independent Brownian motion with drift.

This paper studies the joint tail asymptotics of extrema of the multi-dimensional Gaussian process over random intervals defined as $P(u)\;:\!=\; \mathbb{P}\{\cap_{i=1}^n (\sup_{t\in[0,\mathcal{T}_i]} ( X_{i}(t) +c_i t )>a_i u )\}$ , $u\rightarrow\infty$ , where $X_i(t)$ , $t\ge0$ , $i=1,2,\ldots,n$ , are independent centered Gaussian processes with stationary increments, $\boldsymbol{\mathcal{T}}=(\mathcal{T}_1, \ldots, \mathcal{T}_n)$ is a regularly varying random vector with positive components, which is independent of the Gaussian processes, and $c_i\in \mathbb{R}$ , $a_i>0$ , $i=1,2,\ldots,n$ . Our result shows that the structure of the asymptotics of P(u) is determined by the signs of the drifts $c_i$ . We also discuss a relevant multi-dimensional regenerative model and derive the corresponding ruin probability.

This paper considers a variant of the classical Cramér–Lundberg model that is particularly appropriate in the credit context, with the distinguishing feature that it corresponds to a finite number of obligors. The focus is on computing the ruin probability, i.e. the probability that the initial reserve, increased by the interest received from the obligors and decreased by the losses due to defaults, drops below zero. As well as an exact analysis (in terms of transforms) of this ruin probability, an asymptotic analysis is performed, including an efficient importance-sampling-based simulation approach.

The base model is extended in multiple dimensions: (i) we consider a model in which there may, in addition, be losses that do not correspond to defaults, (ii) then we analyze a model in which the individual obligors are coupled via a regime switching mechanism, (iii) then we extend the model so that between the losses the reserve process behaves as a Brownian motion rather than a deterministic drift, and (iv) we finally consider a set-up with multiple groups of statistically identical obligors.

Consider a multidimensional risk model, in which an insurer simultaneously confronts m (m ≥ 2) types of claims sharing a common non-stationary and non-renewal arrival process. Assuming that the claims arrival process satisfies a large deviation principle and the claim-size distributions are heavy-tailed, asymptotic estimates for two common types of ruin probabilities for this multidimensional risk model are obtained. As applications, we give two examples of the non-stationary point process: a Hawkes process and a Cox process with shot noise intensity, and asymptotic ruin probabilities are obtained for these two examples.

We find explicit estimates for the exponential rate of long-term convergence for the ruin probability in a level-dependent Lévy-driven risk model, as time goes to infinity. Siegmund duality allows us to reduce the problem to long-term convergence of a reflected jump-diffusion to its stationary distribution, which is handled via Lyapunov functions.

The De Vylder and Goovaerts conjecture is an open problem in risk theory, stating that the finite-time ruin probability in a standard risk model is greater than or equal to the corresponding ruin probability evaluated in an associated model with equalized claim amounts. Equalized means here that the jump sizes of the associated model are equal to the average jump in the initial model between 0 and a terminal time T.In this paper, we consider the diffusion approximations of both the standard risk model and its associated risk model. We prove that the associated model, when conveniently renormalized, converges in distribution to a Gaussian process satisfying a simple SDE. We then compute the probability that this diffusion hits the level 0 before time T and compare it with the same probability for the diffusion approximation for the standard risk model. We conclude that the De Vylder and Goovaerts conjecture holds for the diffusion limits.

Consider a particular bidimensional risk model, in which two insurance companies divide between them in different proportions both the premium income and the aggregate claims. In practice, it can be interpreted as an insurer–reinsurer scenario, where the reinsurer takes over a proportion of the insurer's losses. Under the assumption that the claim sizes and inter-arrival times form a sequence of independent and identically distributed random pairs, with each pair obeying a dependence structure, an asymptotic expression for the ruin probability of this bidimensional risk model with constant interest rates is established.

Transform inversions, in which density and survival functions are computed from their associated moment generating function $\mathcal{M}$, have largely been based on methods which use values of $\mathcal{M}$ in its convergence region. Prominent among such methods are saddlepoint approximations and Fourier-series inversion methods, including the fast Fourier transform. In this paper we propose inversion methods which make use of values for $\mathcal{M}$ which lie outside of its convergence region and in its analytic continuation. We focus on the simplest and perhaps richest setting for applications in which $\mathcal{M}$ is either a meromorphic function in its analytic continuation, so that all of its singularities are poles, or else the singularities are isolated essential. Asymptotic expansions of finite- and infinite-orders are developed for density and survival functions using the poles of $\mathcal{M}$ in its analytic continuation. For finite-order expansions, the expansion error is a contour integral in the analytic continuation, which we approximate using the saddlepoint method based on following the path of steepest descent. Such saddlepoint error approximations accurately determine expansion errors and, thus, provide the means for determining the order of the expansion needed to achieve some preset accuracy. They also provide an additive correction term which increases accuracy of the expansion. Further accuracy is achieved by computing the expansion errors numerically using a contour path which ultimately tracks the steepest descent direction. Important applications include Wilks’ likelihood ratio test in MANOVA, compound distributions, and the Sparre Andersen and Cramér–Lundberg ruin models.

In this study, we show how expressions for the probability of ultimate ruin can be obtained from the probability function of the time of ruin in a particular compound binomial risk model, and from the density of the time of ruin in a particular Sparre Andersen risk model. In each case evaluation of generalised binomial series is required, and the argument of each series has a common form. We evaluate these series by creating an identity based on the generalised negative binomial distribution. We also show how the same ideas apply to the probability function of the number of claims in a particular Sparre Andersen model.

We analyse ruin probabilities for an insurance risk process with a more generalised dependence structure compared to the one introduced in Constantinescu et al. (2016). In this paper, we assume that a random threshold window is generated every time after a claim occurs. By comparing the previous inter-claim time with the threshold window, the distributions of the current threshold window and the inter-arrival time are determined. Furthermore, the statuses for the previous and current inter-arrival times give rise to the current claim size distribution as well. Like Constantinescu et al. (2016), we first identify the embedded Markov additive process where all the randomness takes a general form. Inspired by the Erlangisation technique, the key message of this paper is to analyse such risk process using a Markov fluid flow model where the underlying random variables follow phase-type distributions. This would further allow us to approximate the fixed observation windows by Erlang random variables. Then ruin probabilities under the process with Erlang(n) observation windows are proved to be Erlangian approximations for those related to the process with fixed threshold windows at the limit. An exact form of the limit can be obtained whose application will be illustrated further by a numerical example.

We consider phase-type scale mixture distributions which correspond to distributions of a product of two independent random variables: a phase-type random variable Y and a non-negative but otherwise arbitrary random variable S called the scaling random variable. We investigate conditions for such a class of distributions to be either light- or heavy-tailed, we explore subexponentiality and determine their maximum domains of attraction. Particular focus is given to phase-type scale mixture distributions where the scaling random variable S has discrete support – such a class of distributions has been recently used in risk applications to approximate heavy-tailed distributions. Our results are complemented with several examples.

Motor insurance is a very competitive business where insurers operate with quite large portfolios, often decisions must be taken under short horizons and therefore ruin probabilities should be calculated in finite time. The probability of ruin, in continuous and finite time, is numerically evaluated under the classical Cramér–Lundberg risk process framework for a large motor insurance portfolio, where we allow for a posteriori premium adjustments, according to the claim record of each individual policyholder. Focusing on the classical model for bonus-malus systems, we propose that the probability of ruin can be interpreted as a measure to decide between different bonus-malus scales or even between different bonus-malus rules. In our work, the required initial surplus can also be evaluated. We consider an application of a bonus-malus system for motor insurance to study the impact of experience rating in ruin probabilities. For that, we used a real commercial scale of an insurer operating in the Portuguese market, and we also work on various well-known optimal bonus-malus scales estimated with real data from that insurer. Results involving these scales are discussed.

In this paper we identify three questions concerning the management of risk networks with a central branch, which may be solved using the extensive machinery available for one-dimensional risk models. First, we propose a criterion for judging whether a subsidiary is viable by its readiness to pay dividends to the central branch, as reflected by the optimality of the zero-level dividend barrier. Next, for a deterministic central branch which must bailout a single subsidiary each time its surplus becomes negative, we determine the optimal bailout policy, as well as the ruin probability and other risk measures, in closed form. Moreover, we extend these results to the case of hierarchical networks. Finally, for nondeterministic central branches with one subsidiary, we compute approximate risk measures by applying rational approximations, and by using the recently developed matrix scale methodology.

We consider a Cramér–Lundberg insurance risk process with the added feature of reinsurance. If an arriving claim finds the reserve below a certain threshold γ, or if it would bring the reserve below that level, then a reinsurer pays part of the claim. Using fluctuation theory and the theory of scale functions of spectrally negative Lévy processes, we derive expressions for the Laplace transform of the time to ruin and of the joint distribution of the deficit at ruin and the surplus before ruin. We specify these results in much more detail for the threshold set-up in the case of proportional reinsurance.

This paper focuses on the evaluation of the probability that both components of a bivariate stochastic process ever simultaneously exceed some large level; a leading example is that of two Markov fluid queues driven by the same background process ever reaching the set (u, ∞)×(u, ∞), for u>0. Exact analysis being prohibitive, we resort to asymptotic techniques and efficient simulation, focusing on large values of u. The first contribution concerns various expressions for the decay rate of the probability of interest, which are valid under Gärtner–Ellis-type conditions. The second contribution is an importance-sampling-based rare-event simulation technique for the bivariate Markov modulated fluid model, which is capable of asymptotically efficiently estimating the probability of interest; the efficiency of this procedure is assessed in a series of numerical experiments.

We study an M/G/1-type queueing model with the following additional feature. The server works continuously, at fixed speed, even if there are no service requirements. In the latter case, it is building up inventory, which can be interpreted as negative workload. At random times, with an intensity ω(x) when the inventory is at level x>0, the present inventory is removed, instantaneously reducing the inventory to 0. We study the steady-state distribution of the (positive and negative) workload levels for the cases ω(x) is constant and ω(x) = a x. The key tool is the Wiener–Hopf factorization technique. When ω(x) is constant, no specific assumptions will be made on the service requirement distribution. However, in the linear case, we need some algebraic hypotheses concerning the Laplace–Stieltjes transform of the service requirement distribution. Throughout the paper, we also study a closely related model arising from insurance risk theory.