Hubbry Logo
Ruin theoryRuin theoryMain
Open search
Ruin theory
Community hub
Ruin theory
logo
7 pages, 0 posts
0 subscribers
Be the first to start a discussion here.
Be the first to start a discussion here.
Contribute something
Ruin theory
Ruin theory
Ruin theory
View on Wikipedia
from Wikipedia

In actuarial science and applied probability, ruin theory (sometimes risk theory[1] or collective risk theory) uses mathematical models to describe an insurer's vulnerability to insolvency/ruin. In such models key quantities of interest are the probability of ruin, distribution of surplus immediately prior to ruin and deficit at time of ruin.

Classical model

[edit]
A sample path of compound Poisson risk process

The theoretical foundation of ruin theory, known as the Cramér–Lundberg model (or classical compound-Poisson risk model, classical risk process[2] or Poisson risk process) was introduced in 1903 by the Swedish actuary Filip Lundberg.[3] Lundberg's work was republished in the 1930s by Harald Cramér.[4]

The model describes an insurance company who experiences two opposing cash flows: incoming cash premiums and outgoing claims. Premiums arrive a constant rate from customers and claims arrive according to a Poisson process with intensity and are independent and identically distributed non-negative random variables with distribution and mean (they form a compound Poisson process). So for an insurer who starts with initial surplus , the aggregate assets are given by:[5]

The central object of the model is to investigate the probability that the insurer's surplus level eventually falls below zero (making the firm bankrupt). This quantity, called the probability of ultimate ruin, is defined as

,

where the time of ruin is with the convention that . This can be computed exactly using the Pollaczek–Khinchine formula as[6] (the ruin function here is equivalent to the tail function of the stationary distribution of waiting time in an M/G/1 queue[7])

where is the transform of the tail distribution of ,

and denotes the -fold convolution. In the case where the claim sizes are exponentially distributed, this simplifies to[7]

Sparre Andersen model

[edit]

E. Sparre Andersen extended the classical model in 1957[8] by allowing claim inter-arrival times to have arbitrary distribution functions.[9]

where the claim number process is a renewal process and are independent and identically distributed random variables. The model furthermore assumes that almost surely and that and are independent. The model is also known as the renewal risk model.

Expected discounted penalty function

[edit]

Michael R. Powers[10] and Gerber and Shiu[11] analyzed the behavior of the insurer's surplus through the expected discounted penalty function, which is commonly referred to as Gerber-Shiu function in the ruin literature and named after actuarial scientists Elias S.W. Shiu and Hans-Ulrich Gerber. It is arguable whether the function should have been called Powers-Gerber-Shiu function due to the contribution of Powers.[10]

In Powers' notation, this is defined as

,

where is the discounting force of interest, is a general penalty function reflecting the economic costs to the insurer at the time of ruin, and the expectation corresponds to the probability measure . The function is called expected discounted cost of insolvency by Powers.[10]

In Gerber and Shiu's notation, it is given as

,

where is the discounting force of interest and is a penalty function capturing the economic costs to the insurer at the time of ruin (assumed to depend on the surplus prior to ruin and the deficit at ruin ), and the expectation corresponds to the probability measure . Here the indicator function emphasizes that the penalty is exercised only when ruin occurs.

It is quite intuitive to interpret the expected discounted penalty function. Since the function measures the actuarial present value of the penalty that occurs at , the penalty function is multiplied by the discounting factor , and then averaged over the probability distribution of the waiting time to . While Gerber and Shiu[11] applied this function to the classical compound-Poisson model, Powers[10] argued that an insurer's surplus is better modeled by a family of diffusion processes.

There are a great variety of ruin-related quantities that fall into the category of the expected discounted penalty function.

Special case Mathematical representation Choice of penalty function
Probability of ultimate ruin
Joint (defective) distribution of surplus and deficit
Defective distribution of claim causing ruin
Trivariate Laplace transform of time, surplus and deficit
Joint moments of surplus and deficit

Other finance-related quantities belonging to the class of the expected discounted penalty function include the perpetual American put option,[12] the contingent claim at optimal exercise time, and more.

Recent developments

[edit]
  • Compound-Poisson risk model with constant interest
  • Compound-Poisson risk model with stochastic interest
  • Brownian-motion risk model
  • General diffusion-process model
  • Markov-modulated risk model
  • Accident probability factor (APF) calculator – risk analysis model (@SBH)

See also

[edit]

References

[edit]

Further reading

[edit]
Revisions and contributorsEdit on WikipediaRead on Wikipedia

Ruin theory

View on Grokipedia
from Grokipedia
Ruin theory is a branch of and applied probability that models the risk of an insurance company's , defined as the event where its surplus—the difference between accumulated premiums and claims—falls below zero at any time. In this framework, the surplus process is typically expressed as U(t)=u+ct−S(t)U(t) = u + ct - S(t), where uu is the initial surplus, cc is the constant premium income rate per unit time, and S(t)S(t) represents the aggregate claims up to time tt. The ultimate goal is to compute the ruin probability ψ(u)=Pr⁡{U(t)<0 for some t≥0}\psi(u) = \Pr\{ U(t) < 0 \text{ for some } t \geq 0 \}, which quantifies the likelihood of technical ruin starting from initial capital uu. The foundational model in ruin theory is the Cramér-Lundberg model, introduced in the early 20th century and refined by Harald Cramér and Filip Lundberg, which assumes that claims arrive according to a Poisson process with rate λ\lambda, individual claim sizes are independent and identically distributed positive random variables with mean μ\mu, and premiums are collected at a constant rate c>λμc > \lambda \mu to ensure positive safety loading. Under this model, explicit expressions for ψ(u)\psi(u) are available for certain claim distributions, such as exponential claims, where ψ(u)=λμce−(1/μ−λ/c)u\psi(u) = \frac{\lambda \mu}{c} e^{-(1/\mu - \lambda/c) u}, highlighting the exponential decay of ruin risk with increasing initial capital. More generally, the model yields the Cramér-Lundberg approximation for large uu, approximating ψ(u)≈Ce−γu\psi(u) \approx C e^{-\gamma u}, where γ\gamma is the adjustment coefficient solving the Lundberg equation λ+cr=λm^(r)\lambda + c r = \lambda \hat{m}(r) (with m^(r)\hat{m}(r) the moment generating function of claim sizes), and CC is a constant depending on the claim distribution. Key results include Lundberg's inequality, which provides an upper bound ψ(u)≤e−ρu\psi(u) \leq e^{-\rho u} for the adjustment coefficient ρ>0\rho > 0, ensuring that ruin probability diminishes exponentially under positive loading. Extensions of the classical model, such as the Sparre Andersen model, relax the Poisson assumption by allowing general renewal interarrival times for claims, while incorporating dependence structures or perturbations to better capture real-world dynamics. theory also addresses finite-time ruin probabilities, dividend strategies to optimize surplus management, and expected discounted penalties at ruin, aiding insurers in capital allocation, premium setting, and compliance. These tools remain central to non-life risk , with ongoing integrating for claim modeling and climate-related risks.

Overview

Definition and Scope

Ruin theory is a branch of that models the of an insurer's , specifically the event known as "ruin," where the insurer's surplus—defined as reserves plus accumulated premiums minus claims paid—falls below zero at some point in time. This framework quantifies the likelihood of such financial distress under various claim arrival and severity assumptions, enabling actuaries to assess long-term and set appropriate premium rates. The scope of ruin theory extends beyond pure applications to broader areas of applied probability and risk theory, where it analyzes the ruin probability, denoted ψ(u)\psi(u), which represents the probability that ruin ever occurs given an initial surplus u≥0u \geq 0. Central to these models is the assumption of a positive safety loading, meaning the premium income rate exceeds the expected claim outflow rate, ensuring a positive probability of (i.e., ψ(u)<1\psi(u) < 1) and preventing certain ruin in the long run. This condition underpins the theory's practical relevance in maintaining insurer stability while balancing risk and profitability. At its core, ruin theory employs a surplus process to track the insurer's financial position over time, typically expressed as U(t)=u+ct−S(t)U(t) = u + ct - S(t), where uu is the initial surplus, c>0c > 0 is the constant premium collection rate, and S(t)S(t) denotes the aggregate claims up to time tt. Ruin occurs if inf⁡t≥0U(t)<0\inf_{t \geq 0} U(t) < 0. The classical Cramér-Lundberg model serves as the foundational framework for these analyses, incorporating Poisson-distributed claim arrivals and independent claim sizes.

Historical Background

Ruin theory originated in the early 20th century with the foundational work of Swedish actuary Filip Lundberg, whose 1903 dissertation, Approximerad framställning af sannolikhetsfunktionen. Återförsäkring af kollektivrisker, addressed the calculation of risk premiums and provided approximate methods for determining ruin probabilities in insurance portfolios. Lundberg's analysis focused on collective risk processes, modeling the insurer's surplus as a balance between incoming premiums and outgoing claims, thereby laying the groundwork for probabilistic assessments of insolvency risk. In the 1930s, Harald Cramér refined and popularized Lundberg's ideas through his treatise On the Mathematical Theory of Risk, introducing the now-famous Cramér-Lundberg approximation for the tail of ruin probabilities under light-tailed claim distributions. Cramér's contributions integrated stochastic process theory into risk assessment, emphasizing asymptotic behavior and exponential decay rates for ruin probabilities when the premium loading is positive, which became central to the classical model. Following World War II, the theory expanded with Erik Sparre Andersen's 1957 paper "On the collective theory of risk in the case of contagion between claims," which generalized the interclaim time distributions beyond the Poisson assumption to arbitrary renewal processes, enabling broader applicability to non-exponential waiting times between claims. This Sparre Andersen model marked a significant milestone by incorporating renewal theory, allowing for more flexible modeling of claim arrival patterns. During the 1990s, the framework evolved to include discounted penalty functions, pioneered by Hans U. Gerber, Michel Goovaerts, and Rob Kaas, with Gerber and Elias S.W. Shiu's 1998 paper "On the Time Value of Ruin" formalizing the Gerber-Shiu function to capture the expected discounted cost at ruin, incorporating the time of ruin, surplus immediately prior, and deficit at ruin. This development shifted focus toward practical actuarial applications, such as valuing ruin-related penalties under time-discounted criteria. In the 1960s and 1970s, ruin theory drew substantial influence from advancements in and problems, with researchers like Lajos Takács establishing equivalences between ruin probabilities and waiting times in single-server queues or barrier-crossing probabilities in , enhancing analytical tools through combinatorial and fluctuation methods.

Classical Model

Surplus Process

In the classical Cramér-Lundberg model of ruin theory, the surplus process models the financial position of an insurance company over time, starting from an initial capital u≥0u \geq 0. It is defined as U(t)=u+ct−∑i=1N(t)Xi,t≥0,U(t) = u + c t - \sum_{i=1}^{N(t)} X_i, \quad t \geq 0, where c>0c > 0 is the constant rate of premium income, N(t)N(t) is a homogeneous Poisson process with intensity λ>0\lambda > 0 representing the number of claims arriving by time tt, and the XiX_i (for i=1,2,…i = 1, 2, \dots) are independent and identically distributed positive random variables denoting individual claim sizes, independent of the arrival process N(t)N(t). This formulation, introduced by Harald Cramér in his foundational work on collective risk theory, captures the continuous inflow of premiums offset by discrete, random claim outflows. The premium rate cc is chosen to exceed the expected claim outflow, typically as c=λμ(1+θ)c = \lambda \mu (1 + \theta), where μ=E[X1]<∞\mu = \mathbb{E}[X_1] < \infty is the claim size and θ>0\theta > 0 is the security loading factor that ensures long-term profitability. Key assumptions include the independence of claim amounts XiX_i and their distribution from the Poisson arrival , with interclaim times exponentially distributed with 1/λ1/\lambda. The aggregate claims up to time tt, S(t)=∑i=1N(t)XiS(t) = \sum_{i=1}^{N(t)} X_i, constitute a , which has stationary and independent increments, reflecting the memoryless property of claim arrivals. The process exhibits positive drift under the net profit condition c>λμc > \lambda \mu, meaning E[U(t)]=u+(c−λμ)t\mathbb{E}[U(t)] = u + (c - \lambda \mu) t increases linearly with tt. Ruin is defined as the first entrance into negative surplus, at time τ=inf⁡{t≥0:U(t)<0}\tau = \inf\{ t \geq 0 : U(t) < 0 \}, with the paths of U(t)U(t) being piecewise linear (increasing at rate cc between jumps) and right-continuous with left limits due to the compound Poisson structure of claims. This setup forms the basis for analyzing the probability of ultimate ruin, ψ(u)=Pr⁡(τ<∞∣U(0)=u)\psi(u) = \Pr(\tau < \infty \mid U(0) = u).

Ultimate Ruin Probability

In the classical Cramér-Lundberg model, the ultimate ruin probability, denoted ψ(u), is defined as the probability that ruin ever occurs given an initial surplus of u, that is, ψ(u) = P(τ < ∞ | U(0) = u), where τ denotes the time of ruin and U(t) represents the surplus process at time t. This quantity captures the long-term risk of insolvency in a surplus process driven by constant premium income and compound Poisson claim arrivals. The function ψ(u) satisfies the following integro-differential equation for u > 0: cψ′(u)=λ∫0uψ(u−x) dF(x)+λ∫u∞dF(x),c \psi'(u) = \lambda \int_0^u \psi(u - x) \, dF(x) + \lambda \int_u^\infty dF(x), where c is the constant premium rate per unit time, λ is the Poisson arrival rate of claims, and F is the cumulative distribution function of the individual claim sizes. This equation arises from conditioning on the time and size of the first claim, balancing the infinitesimal change in surplus against the expected impact of claims. Boundary conditions are ψ(0) = \frac{\lambda \mu}{c} (ruin from zero surplus under positive loading c > λ μ) and ψ(u) → 0 as u → ∞ (ruin becomes impossible with infinite capital under the net profit condition c > λ ∫_0^∞ x dF(x)). When claim sizes are exponentially distributed with mean 1 (that is, = 1 - e^{-x} for x ≥ 0), the admits an exact closed-form solution: ψ(u) = \frac{\lambda}{c} e^{-(1 - \lambda/c) u} for u ≥ 0, provided that the safety loading condition λ < c holds to ensure ψ(u) < 1. This explicit expression highlights the exponential decay of the ruin probability with increasing initial capital, reflecting the memoryless property of the exponential distribution that simplifies the renewal structure of the surplus process. For general claim size distributions F with finite mean μ (where the net profit condition requires ρ = λμ / c < 1), the ultimate ruin probability admits the Pollaczek–Khinchine formula, which expresses ψ(u) as a compound geometric tail probability: ψ(u)=(1−ρ)∑n=0∞ρnGˉ∗n(u),\psi(u) = (1 - \rho) \sum_{n=0}^\infty \rho^n \bar{G}^{*n}(u), where ρ = λμ / c, \bar{G}(u) = 1 - G(u), G^{*n} denotes the n-fold convolution of G with itself (with G^{*0} being the unit point mass at 0), and G is the ladder height distribution given by G(u) = \frac{1}{\mu} \int_0^u (1 - F(x)) , dx for u ≥ 0. This representation interprets ψ(u) probabilistically as the tail of a geometric sum of independent ladder height increments, each distributed according to the equilibrium (or integrated tail) version of F, providing a versatile framework for numerical computation and approximation even when closed forms are unavailable.

Adjustment Coefficient and Lundberg Inequality

In the classical Cramér-Lundberg model, the adjustment coefficient plays a central role in providing asymptotic estimates for the ultimate ruin probability when claim sizes are light-tailed, meaning their moment generating function exists in a neighborhood of zero. The adjustment coefficient R>0R > 0 is defined as the unique positive solution to the Lundberg equation λ+cr=λMX(r),\lambda + c r = \lambda M_X(r), where λ>0\lambda > 0 is the Poisson arrival rate of claims, c>0c > 0 is the constant premium income rate, and MX(r)=E[erX]M_X(r) = \mathbb{E}[e^{r X}] is the moment generating function of the individual claim size X>0X > 0. The existence of such an RR requires two key conditions: the claim size distribution must be light-tailed (i.e., MX(r)<∞M_X(r) < \infty for some r>0r > 0), and the net profit condition (positive safety loading) must hold, c>λμc > \lambda \mu, where μ=E[X]\mu = \mathbb{E}[X]. These ensure that the equation has exactly one positive root, as the function λ(MX(r)−1)−cr\lambda (M_X(r) - 1) - c r is convex, zero at r=0r = 0, and changes sign appropriately for r>0r > 0. The adjustment coefficient quantifies the rate at which the surplus process drifts positively, enabling exponential bounds on ruin. A fundamental result derived using martingale techniques or change-of-measure arguments is the Lundberg inequality, which states that the ultimate ruin probability satisfies ψ(u)≤e−Ru,u≥0.\psi(u) \leq e^{-R u}, \quad u \geq 0. This exponential upper bound holds under the same conditions as the existence of RR, providing a conservative estimate for solvency assessments even when exact computation is intractable. The inequality originates from the supermartingale property of e−RU(t)e^{-R U(t)}, where U(t)U(t) is the surplus process, and is tight in the sense that ψ(u)\psi(u) decays at least as fast as e−Rue^{-R u}. For light-tailed claims, a sharper asymptotic result known as the Cramér-Lundberg approximation refines this bound, stating that ψ(u)∼Ce−Ruasu→∞,\psi(u) \sim C e^{-R u} \quad \text{as} \quad u \to \infty, where CC is a positive constant less than 1 depending on the claim distribution. This approximation captures the precise leading-order behavior of ψ(u)\psi(u), derived via renewal theory applied to the ladder-height process or Lundberg exponent analysis, and is particularly useful for large initial capital uu. The value of CC reflects the sub-exponential decay relative to the bound.

Generalized Models

Sparre Andersen Model

The Sparre Andersen model generalizes the classical Cramér-Lundberg model by permitting interclaim times to follow an arbitrary distribution, rather than restricting them to the characteristic of the Poisson arrival process. In this model, the interclaim times {Ti}i=1∞\{T_i\}_{i=1}^\infty are independent and identically distributed (i.i.d.) positive random variables with a general distribution function FT(t)=P(Ti≤t)F_T(t) = P(T_i \leq t) and E[Ti]=1/λ>0\mathbb{E}[T_i] = 1/\lambda > 0, where λ>0\lambda > 0 represents the intensity of claims. The claim sizes {Xi}i=1∞\{X_i\}_{i=1}^\infty are also i.i.d. positive random variables with distribution function FX(y)=P(Xi≤y)F_X(y) = P(X_i \leq y) and μ=E[Xi]<∞\mu = \mathbb{E}[X_i] < \infty, independent of the interclaim times. The claim arrival process is defined by the renewal counting process N(t)=sup⁡{n≥0:Sn≤t}N(t) = \sup\{n \geq 0 : S_n \leq t\}, where S0=0S_0 = 0 and Sn=T1+⋯+TnS_n = T_1 + \cdots + T_n for n≥1n \geq 1. Premiums are collected at a constant rate c>λμc > \lambda \mu, ensuring a positive security loading condition θ=(c−λμ)/(λμ)>0\theta = (c - \lambda \mu)/(\lambda \mu) > 0. The surplus process is then given by U(t)=u+ct−∑i=1N(t)Xi,U(t) = u + c t - \sum_{i=1}^{N(t)} X_i, where u≥0u \geq 0 is the initial surplus. The probability of ultimate ruin is ψ(u)=P(inf⁡t≥0U(t)<0∣U(0)=u)\psi(u) = P(\inf_{t \geq 0} U(t) < 0 \mid U(0) = u). The ruin probability ψ(u)\psi(u) satisfies a defective renewal equation derived from renewal theory, where the process is analyzed in terms of ladder heights or the net loss over renewal cycles. Specifically, conditioning on the time and size of the first claim leads to an integral equation of the form ψ(u)=∫0∞∫0u+ctψ(u+ct−y) dFX(y) dFT(t)+∫0∞∫u+ct∞dFX(y) dFT(t),\psi(u) = \int_0^\infty \int_0^{u + c t} \psi(u + c t - y) \, dF_X(y) \, dF_T(t) + \int_0^\infty \int_{u + c t}^\infty dF_X(y) \, dF_T(t), which can be expressed as a convolution ψ(u)=∫0uψ(u−y) dK(y)+Kˉ(u)\psi(u) = \int_0^u \psi(u - y) \, dK(y) + \bar{K}(u), where KK is the defective distribution function of the claim surplus process over the first interclaim period (with total mass less than 1) and Kˉ(u)=1−K(u)\bar{K}(u) = 1 - K(u). This defective renewal equation is solved using the Pollaczek-Khintchine formula, yielding ψ(u)=∑k=1∞(1−ρ)ρk−1Hˉ∗k(u)\psi(u) = \sum_{k=1}^\infty (1 - \rho) \rho^{k-1} \bar{H}^{*k}(u), where ρ=λμ/c<1\rho = \lambda \mu / c < 1 is the traffic intensity and HH is the integrated tail distribution of the equilibrium claim surplus. Under the assumption that the moment generating functions E[erXi]\mathbb{E}[e^{r X_i}] and E[erTi]\mathbb{E}[e^{r T_i}] exist in a neighborhood of the origin, the adjustment coefficient R>0R > 0 is defined as the unique positive solution to the Lundberg equation E[eR(cT1−X1)]=1\mathbb{E}[e^{R(c T_1 - X_1)}] = 1, or equivalently, m^T(Rc)m^X(−R)=1\hat{m}_T(R c) \hat{m}_X(-R) = 1, where m^T\hat{m}_T and m^X\hat{m}_X are the moment generating functions of T1T_1 and X1X_1, respectively. The asymptotic behavior of the ruin probability is then ψ(u)∼Ce−Ru\psi(u) \sim C e^{-R u} as u→∞u \to \infty, where C=λμ−cE[T1eRcT1]/E[T1]RcE[T1eRcT1]−λμC = \frac{\lambda \mu - c \mathbb{E}[T_1 e^{R c T_1}] / \mathbb{E}[T_1]}{R c \mathbb{E}[T_1 e^{R c T_1}] - \lambda \mu} is a constant depending on the model parameters; this mirrors the classical case but with the general interclaim distribution incorporated into RR. The Lundberg inequality ψ(u)≤e−Ru\psi(u) \leq e^{-R u} holds for all u≥0u \geq 0. Exact closed-form expressions for ψ(u)\psi(u) are available when the interclaim times follow an , Erlang(n,βn, \beta), which is the sum of nn i.i.d. exponential random variables with rate β>0\beta > 0 (mean n/βn/\beta). In this case, the renewal structure allows the ruin probability to be expressed as a of nn exponential terms, obtained by solving an nnth-order derived from the integro-differential form of the renewal equation. For instance, when n=2n=2 (Erlang(2) interclaims) and claims are exponentially distributed, explicit formulas involve the roots of the characteristic equation associated with the Lundberg equation. This tractability arises because the rational of the Erlang density facilitates recursive or transform-based solutions.

Renewal and Diffusion Approximations

In the Sparre Andersen model, where claim interarrival times follow a general renewal process, the ultimate ruin probability ψ(u) can be expressed using Beekman's convolution series, which facilitates numerical computation through successive convolutions. Specifically, ψ(u) = \sum_{n=1}^\infty (1 - \rho) \rho^{n-1} \bar{H}^{*n}(u), where \rho = \lambda \mu / c < 1 is the traffic intensity with \lambda the claim arrival rate, \mu = E[X] the mean claim size, and c the premium rate, and \bar{H}(u) = 1 - H(u) with H the integrated tail distribution H(u) = (1/\mu) \int_0^u \bar{F}(y) , dy, where \bar{F}(y) = 1 - F(y) is the claim size survival function. This series representation arises from the renewal structure of the surplus process and the geometric tail of the number of ladder heights exceeding the initial capital u. The renewal approximation is particularly useful for practical evaluation when exact solutions are intractable, as it leverages the defective renewal equation for ψ(u) and allows truncation after a finite number of terms for sufficient accuracy in moderate initial surplus scenarios. For instance, with exponentially distributed claims, the convolutions simplify to closed forms involving the gamma distribution, demonstrating rapid convergence under the net profit condition c > \lambda \mu. Diffusion approximations provide an alternative approach by modeling the surplus process as a with positive drift \theta = c - \lambda \mu and variance parameter \sigma^2 = \lambda E[X^2], valid under conditions of high claim (large \lambda with small individual claims) or light-tailed claim distributions where the process exhibits diffusive behavior. For large initial surplus u, the ruin probability approximates ψ(u) \approx e^{-2 \theta u / \sigma^2}, matching the exact ruin probability for the approximating and capturing the driven by the safety loading. For moderate values of u, the power refines the estimate by incorporating a z, yielding ψ(u) \approx \Phi\left( -\frac{\theta u + z}{\sqrt{\sigma^2 u}} \right), where \Phi is the standard normal and z \approx 2/3 is an empirical adjustment to improve accuracy over the plain . This method performs well when the scaled surplus process converges to , as validated in numerical studies for compound Poisson risks with light tails.

Discounted Penalty Functions

Gerber-Shiu Function

The Gerber-Shiu function, named after actuaries Hans U. Gerber and Elias S. W. Shiu who introduced it in , generalizes the ultimate ruin probability in the classical compound Poisson risk model by incorporating a discounted penalty that depends on the time of ruin, the surplus immediately prior to ruin, and the deficit at ruin. Specifically, for initial surplus u≥0u \geq 0, it is defined as m(u)=E[e−δτw(U(τ−),∣U(τ)∣)1{τ<∞}∣U(0)=u],m(u) = \mathbb{E}\left[ e^{-\delta \tau} w(U(\tau^-), |U(\tau)|) \mathbf{1}_{\{\tau < \infty\}} \mid U(0) = u \right], where δ≥0\delta \geq 0 is the force of interest (discount rate), τ=inf⁡{t≥0:U(t)<0}\tau = \inf\{ t \geq 0 : U(t) < 0 \} is the time of ruin, U(τ−)U(\tau^-) is the surplus just before ruin, ∣U(τ)∣|U(\tau)| is the deficit at ruin, w(x,y)≥0w(x, y) \geq 0 is a non-negative penalty function, and 1{τ<∞}\mathbf{1}_{\{\tau < \infty\}} is the indicator that ruin occurs. This expectation captures the expected present value of a penalty incurred upon ruin, allowing for analysis of various risk measures beyond mere occurrence probability. In the Cramér-Lundberg model, where claims arrive as a Poisson process with intensity λ>0\lambda > 0 and individual claim sizes follow a distribution with K(x)K(x) and density k(x)k(x), the Gerber-Shiu function satisfies the δm(u)+λm(u)=cm′(u)+λ∫0um(u−x) dK(x)+λ∫u∞w(u,x−u) dK(x),u>0,\delta m(u) + \lambda m(u) = c m'(u) + \lambda \int_0^u m(u - x) \, dK(x) + \lambda \int_u^\infty w(u, x - u) \, dK(x), \quad u > 0, with m(u)→0m(u) \to 0 as u→∞u \to \infty under net profit conditions, and m(0)m(0) determined separately by conditioning on the time and size of the first claim from initial surplus 0. Here, c>0c > 0 is the constant premium income rate. This equation arises from conditioning on the time and size of the first claim, balancing the discounted function value with the drift due to premiums and the terms accounting for non-ruin and ruin scenarios, respectively. Special cases of the Gerber-Shiu function recover classical ruin quantities. When δ=0\delta = 0 and w(x,y)=1w(x, y) = 1, it reduces to the ultimate ruin probability ψ(u)=P(τ<∞∣U(0)=u)\psi(u) = \mathbb{P}(\tau < \infty \mid U(0) = u). For δ=0\delta = 0 and linear penalty w(x,y)=yw(x, y) = y, it yields E[∣U(τ)∣1{τ<∞}∣U(0)=u]E[ |U(\tau)| \mathbf{1}_{\{\tau < \infty\}} \mid U(0) = u ], from which the conditional expectation E[∣U(τ)∣∣τ<∞,U(0)=u]=m(u)/ψ(u)E[ |U(\tau)| \mid \tau < \infty, U(0) = u ] = m(u) / \psi(u) is obtained. Solutions to the integro-differential equation are typically obtained via Laplace transforms, which convert the equation into an algebraic form solvable for the transform, followed by inversion, or through defective renewal equations derived by integrating over ladder height distributions. These methods exploit the renewal structure of the surplus process, expressing m(u)m(u) in terms of convolutions that can be analyzed asymptotically or numerically.

Applications to Ruin Analysis

The Gerber-Shiu function provides a versatile framework for quantifying key features of the ruin event, including the severity of the deficit, the surplus level just prior to ruin, and the timing of ruin occurrence. By selecting specific forms of the penalty function w(x,y)w(x, y), where xx represents the surplus immediately before ruin U(τ−)U(\tau^-) and yy the deficit at ruin ∣U(τ)∣|U(\tau)|, analysts can derive conditional expectations and distributions given that ruin occurs (τ<∞\tau < \infty). This approach extends classical ruin probability analysis to capture the magnitude and circumstances of insolvency. One primary application involves the distribution of the deficit at ruin. Setting the penalty function to w(x,y)=yw(x, y) = y with discount factor δ=0\delta = 0 yields the Gerber-Shiu function m(u)=E[∣U(τ)∣1{τ<∞}∣U(0)=u]m(u) = E[ |U(\tau)| \mathbf{1}_{\{\tau < \infty\}} \mid U(0) = u ], which, when divided by the ultimate ruin probability ψ(u)\psi(u), provides the conditional expectation E[∣U(τ)∣∣τ<∞,U(0)=u]E[ |U(\tau)| \mid \tau < \infty, U(0) = u ]. This metric assesses the typical overshoot below zero upon ruin, aiding in evaluating potential recovery costs or capital shortfalls for insurers. For the cumulative distribution, w(x,y)=1{y≤z}w(x, y) = \mathbf{1}_{\{y \leq z\}} derives P(∣U(τ)∣≤z∣τ<∞,U(0)=u)=m(u)/ψ(u)P( |U(\tau)| \leq z \mid \tau < \infty, U(0) = u ) = m(u) / \psi(u). Similarly, the surplus immediately before ruin can be analyzed by choosing w(x,y)=xw(x, y) = x, resulting in m(u)/ψ(u)=E[U(τ−)∣τ<∞,U(0)=u]m(u) / \psi(u) = E[ U(\tau^-) \mid \tau < \infty, U(0) = u ]. This expected value highlights the proximity to ruin thresholds, informing risk management strategies around near-miss levels and dividend policies that might precipitate insolvency. Such applications reveal how initial surplus uu influences pre-ruin accumulation, with higher uu typically leading to larger expected pre-ruin surpluses in compound Poisson models. To study the timing of ruin, higher moments of the time to ruin τ\tau are obtained via a generalized Gerber-Shiu function that incorporates powers of τ\tau. With δ=0\delta = 0 and the penalty extended to include τk\tau^k, the function computes E[τk1{τ<∞}∣U(0)=u]/ψ(u)=E[τk∣τ<∞,U(0)=u]E[ \tau^k \mathbf{1}_{\{\tau < \infty\}} \mid U(0) = u ] / \psi(u) = E[ \tau^k \mid \tau < \infty, U(0) = u ] for k=1,2,…k = 1, 2, \dots. These moments quantify the variability and expected duration until ruin, essential for assessing long-term solvency horizons; for instance, the first moment E[τ∣τ<∞]E[\tau \mid \tau < \infty] grows with initial surplus in classical risk processes. Joint distributions of the surplus before and deficit at ruin further enrich the analysis. For example, setting w(x,y)=1{x≤a,y≤b}w(x, y) = \mathbf{1}_{\{x \leq a, y \leq b\}} produces m(u)/ψ(u)=P(U(τ−)≤a,∣U(τ)∣≤b∣τ<∞,U(0)=u)m(u) / \psi(u) = P( U(\tau^-) \leq a, |U(\tau)| \leq b \mid \tau < \infty, U(0) = u ), capturing the bivariate dependence between pre- and post-ruin levels. This joint perspective is crucial for understanding correlated risk exposures, such as in models where large claims drive both high pre-ruin surpluses and substantial deficits. For practical evaluation in models with intricate claim size distributions, numerical methods are employed to compute Gerber-Shiu functions. In discrete-time settings, such as the compound binomial model, Panjer recursion efficiently calculates the function by iteratively convolving claim and interclaim probabilities, offering stability and speed for finite-time horizons. For continuous or complex cases, Monte Carlo simulations approximate expectations through importance sampling or direct path generation of the surplus process, providing flexible estimates with controllable error via sample size, though at higher computational cost.

Recent Advances

Models with Investments

In ruin theory, models incorporating investment income extend the classical surplus process by allowing the insurer's assets to generate returns, thereby influencing the dynamics of ruin probability. The surplus process with constant interest force rr is typically modeled in continuous time as dU(t)=(c+rU(t)) dt−dS(t),dU(t) = (c + r U(t)) \, dt - dS(t), where U(t)U(t) is the surplus at time tt, c>0c > 0 is the constant premium rate, and S(t)S(t) represents the aggregate claims process, often a . This formulation assumes risk-free investments, such as bonds, yielding a deterministic return proportional to the current surplus. In discrete-time variants, the surplus may be allocated between and bonds, with returns updating at each period based on market performance. Under this constant interest model, the ultimate ruin probability ψ(u)\psi(u), starting from initial surplus uu, satisfies a modified Lundberg equation that incorporates rr. The adjustment coefficient RR solves an equation adjusting the classical Lundberg exponent by the interest factor, leading to the bound ψ(u)≤e−Ru\psi(u) \leq e^{-R u}, where RR is increased compared to the no-interest case due to the growth from investments. This bound tightens the classical estimate, reflecting how interest mitigates ruin risk when the net profit condition holds. For stochastic investments, the surplus process integrates financial market volatility, often modeling asset returns via (GBM) with drift μ\mu and volatility σ\sigma: dA(t)=μA(t) dt+σA(t) dW(t)dA(t) = \mu A(t) \, dt + \sigma A(t) \, dW(t), where A(t)A(t) is the asset value and W(t)W(t) is a standard . The insurer allocates a proportion α\alpha of surplus to the risky asset, yielding dU(t)=c dt+α(μU(t) dt+σU(t) dW(t))−dS(t)dU(t) = c \, dt + \alpha (\mu U(t) \, dt + \sigma U(t) \, dW(t)) - dS(t). Adaptations to Heston's stochastic volatility model replace constant σ\sigma with a mean-reverting process dσ2(t)=κ(θ−σ2(t)) dt+ξσ(t) dZ(t)d\sigma^2(t) = \kappa (\theta - \sigma^2(t)) \, dt + \xi \sigma(t) \, dZ(t), capturing in investments. Asymptotic analysis of ψ(u)\psi(u) employs large deviation principles, showing exponential decay for small volatility (σ2/2<μ\sigma^2 / 2 < \mu) but certain ruin (ψ(u)=1\psi(u) = 1) when volatility dominates drift. Post-2010 developments have yielded explicit expressions for the Gerber-Shiu function in these models, particularly with exponential claims. For instance, in the constant interest case with exponentially distributed claims (mean 1/β1/\beta), the Gerber-Shiu expected discounted penalty function m(u)=E[e−δτw(Uτ−,∣Uτ∣)1{τ<∞}∣U0=u]m(u) = \mathbb{E}[e^{-\delta \tau} w(U_{\tau^-}, |U_\tau|) \mathbf{1}_{\{\tau < \infty\}} \mid U_0 = u] admits a closed-form solution involving modified Bessel functions, solving the associated integro-differential equation. A 2015 survey highlights extensions to risky investments, confirming these results hold under mild conditions on claim sizes. Further, ruin probabilities under dividend strategies, such as barrier dividends where surplus above a level bb is paid out, or proportional reinsurance combined with investments, show that optimal allocation reduces ψ(u)\psi(u) by balancing insurance and financial risks.

Dependent Risks and Stochastic Premiums

In ruin theory, dependent risks arise when claim occurrences or sizes across multiple lines of business or portfolios exhibit correlation, often modeled using copulas to capture joint distributions beyond independence. For bivariate claim processes, copulas such as Clayton or Frank facilitate the modeling of dependence between interwaiting times and claim amounts, enabling the computation of moments and risk measures via Monte Carlo simulation. This approach reveals that negative dependence (e.g., parameter θ = -1) can lead to higher-order moments and elevated value-at-risk estimates compared to independence, particularly for overdispersed data like catastrophe claims. Lévy copulas further extend this to model dependencies in claim arrival processes, allowing for the minimization of ultimate ruin probability ψ(u) by optimizing premium loadings that account for correlated jumps; in multi-risk settings, these loadings vary with initial surplus u and dependence strength, differing from profit-maximizing strategies. Stochastic premiums introduce variability in income streams, often modulated by Markov processes or renewal structures to reflect economic fluctuations or policy adjustments. In discrete-time models, premiums can be autocorrelated via copulas, with claim numbers following a Poisson INMA(1) process parameterized by dependence β, while premiums use a Poisson INAR(1) with thinning α. The ultimate ruin probability ψ(u) is then derived using an embedded Markov chain, yielding the Lundberg exponent R as the positive solution to 1 = M_X(r) G_N(r) / G_M(r), where M_X(r) is the moment generating function of net premium per period, and G_N, G_M are probability generating functions for claims and premiums; under positive safety loading, ψ(u) ≈ C e^{-R u} for some constant C. Such models highlight how autocorrelation amplifies ruin risk when α > 0, necessitating adjusted reserving. For heavy-tailed claims under dependence, subexponential distributions paired with copulas like the Farlie-Gumbel-Morgenstern yield asymptotic expressions for finite-time sum-ruin probabilities in bidimensional renewal models. Specifically, when claims are dependent and subexponential, the equivalence leads to ψ(u, t) ∼ \int_u^\infty \bar{F}(y) , dy / \mu as u → ∞, where \bar{F} is the integrated distribution and μ the mean net profit per unit time, extending classical results to correlated environments. This approximation holds uniformly under mild dependence conditions, providing bounds crucial for high initial capitals. Diffusion-perturbed models incorporate to add continuous volatility to renewal claim processes, with dependence introduced across dimensions for multi-line . In two-dimensional settings, the surplus follows a correlated alongside jumps, yielding diffusion approximations for both finite- and infinite-time ψ(u); for example, the infinite-time ruin satisfies a system of integro-differential equations solvable via scale functions, revealing that positive between Brownian components reduces ψ(u) compared to . Recent numerical advances leverage to simulate dependent ruin paths, enhancing efficiency over traditional for complex copula structures. Neural networks, such as recurrent variants, approximate ψ(u) by learning path dependencies in heavy-tailed or Markov-modulated settings. Recent applications in optimization demonstrate improved computational efficiency. These trends, evident in frameworks for , facilitate real-time solvency assessments in dependent portfolios.

References

  1. https://onlinelibrary.wiley.com/doi/book/10.1002/9780470391341
  2. https://www.worldscientific.com/worldscibooks/10.1142/7431
  3. https://doi.org/10.18637/jss.v025.i07
  4. https://www.actuaries.org.uk/system/files/documents/pdf/e11-paper.pdf
  5. https://cran.r-project.org/web/packages/actuar/vignettes/risk.pdf
  6. https://www.vaia.com/en-us/explanations/business-studies/actuarial-science-in-business/ruin-theory/
  7. https://www.scirp.org/reference/referencespapers?referenceid=3531394
  8. https://www.tandfonline.com/doi/pdf/10.1080/10920277.2016.1245624
  9. https://www.researchgate.net/publication/312081532_On_Cramer%27s_First_Contributions_to_Ruin_Theory
  10. https://www.cambridge.org/core/journals/annals-of-actuarial-science/article/moments-of-the-time-of-ruin-in-sparre-andersen-risk-models/C104EC1B002AC4309DC7CF513FACEED7
  11. https://www.researchgate.net/publication/233630549_On_the_Time_Value_of_Ruin
  12. https://books.google.com/books/about/Collective_Risk_Theory.html?id=i8LuAAAAMAAJ
  13. https://doi.org/10.1007/978-3-319-72005-0_5
  14. https://www.cambridge.org/core/journals/journal-of-the-staple-inn-actuarial-society/article/on-the-mathematical-theory-of-risk/CC05E74A2C8A60C905EEDDD2BA148120
  15. https://fiveable.me/actuarial-mathematics/unit-8/lundbergs-inequality-adjustment-coefficients/study-guide/UawmNDXVAdWvseuE
  16. http://www.mysmu.edu/faculty/yktse/NAM/NAM_S5.pdf
  17. https://onlinelibrary.wiley.com/doi/abs/10.1002/9780470012505.tal026
  18. https://arxiv.org/pdf/1110.2642
  19. https://www.researchgate.net/profile/Jan-Dhaene/publication/319801833_Ruin_theory/links/59c784ba0f7e9bd2c0143f00/Ruin-theory.pdf
  20. https://doi.org/10.1142/7731
  21. https://doi.org/10.1007/978-1-4613-0111-0
  22. https://doi.org/10.1002/9780470317044
  23. https://www.casact.org/sites/default/files/database/astin_vol19no1_71.pdf
  24. https://www.sciencedirect.com/science/article/pii/0167668788900959
  25. https://www.soa.org/globalassets/assets/library/research/transactions-of-society-of-actuaries/1969/january/tsa69v21pt1n59ab5.pdf
  26. https://www.tandfonline.com/doi/abs/10.1080/10920277.1998.10595671
  27. https://arxiv.org/pdf/2203.10680
  28. https://www.mdpi.com/2227-7390/7/10/891
  29. https://www.sciencedirect.com/science/article/abs/pii/S0167668797000279
  30. https://ir.library.oregonstate.edu/downloads/cz30pw80c
  31. https://www.sciencedirect.com/science/article/abs/pii/S0167668712001473
  32. https://www.aimsciences.org/article/doi/10.3934/puqr.2023023
  33. https://www.mdpi.com/2227-9091/9/6/109
  34. https://arxiv.org/abs/2108.10075
  35. https://doi.org/10.3390/e25040698
  36. https://www.researchgate.net/publication/353022049_Asymptotic_sum-ruin_probability_for_a_bidimensional_renewal_risk_model_with_subexponential_claims
  37. https://www.mdpi.com/2227-9091/10/6/129
  38. https://arxiv.org/abs/2408.06168
Add your contribution
Related Hubs
Contribute something
User Avatar
No comments yet.