Ruin theory

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 c > 0 from customers and claims arrive according to a Poisson process ${\displaystyle N_{t}}$ with intensity λ and are independent and identically distributed non-negative random variables ${\displaystyle \xi _{i}}$ with distribution F and mean μ (they form a compound Poisson process). So for an insurer who starts with initial surplus x, the aggregate assets ${\displaystyle X_{t}}$ are given by:[5]

${\displaystyle X_{t}=x+ct-\sum _{i=1}^{N_{t}}\xi _{i}\quad {\text{ for t}}\geq 0.}$

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

${\displaystyle \psi (x)=\mathbb {P} ^{x}\{\tau <\infty \}}$

where the time of ruin is ${\displaystyle \tau =\inf\{t>0\,:\,X(t)<0\}}$ with the convention that ${\displaystyle \inf \varnothing =\infty }$. 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])

${\displaystyle \psi (x)=\left(1-{\frac {\lambda \mu }{c}}\right)\sum _{n=0}^{\infty }\left({\frac {\lambda \mu }{c}}\right)^{n}(1-F_{l}^{\ast n}(x))}$

where ${\displaystyle F_{l}}$ is the transform of the tail distribution of ${\displaystyle F}$,

${\displaystyle F_{l}(x)={\frac {1}{\mu }}\int _{0}^{x}\left(1-F(u)\right){\text{d}}u}$

and ${\displaystyle \cdot ^{\ast n}}$ denotes the ${\displaystyle n}$-fold convolution. In the case where the claim sizes are exponentially distributed, this simplifies to[7]

${\displaystyle \psi (x)={\frac {\lambda \mu }{c}}e^{-\left({\frac {1}{\mu }}-{\frac {\lambda }{c}}\right)x}.}$

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

${\displaystyle X_{t}=x+ct-\sum _{i=1}^{N_{t}}\xi _{i}\quad {\text{ for }}t\geq 0,}$

where the claim number process ${\displaystyle (N_{t})_{t\geq 0}}$ is a renewal process and ${\displaystyle (\xi _{i})_{i\in \mathbb {N} }}$ are independent and identically distributed random variables. The model furthermore assumes that ${\displaystyle \xi _{i}>0}$ almost surely and that ${\displaystyle (N_{t})_{t\geq 0}}$ and ${\displaystyle (\xi _{i})_{i\in \mathbb {N} }}$ are independent. The model is also known as the renewal risk model.

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

${\displaystyle m(x)=\mathbb {E} ^{x}[e^{-\delta \tau }K_{\tau }]}$,

where ${\displaystyle \delta }$ is the discounting force of interest, ${\displaystyle K_{\tau }}$ is a general penalty function reflecting the economic costs to the insurer at the time of ruin, and the expectation ${\displaystyle \mathbb {E} ^{x}}$ corresponds to the probability measure ${\displaystyle \mathbb {P} ^{x}}$. The function is called expected discounted cost of insolvency by Powers.[10]

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

${\displaystyle m(x)=\mathbb {E} ^{x}[e^{-\delta \tau }w(X_{\tau -},X_{\tau })\mathbb {I} (\tau <\infty )]}$,

where ${\displaystyle \delta }$ is the discounting force of interest and ${\displaystyle w(X_{\tau -},X_{\tau })}$ 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 ${\displaystyle X_{\tau -}}$ and the deficit at ruin ${\displaystyle X_{\tau }}$), and the expectation ${\displaystyle \mathbb {E} ^{x}}$ corresponds to the probability measure ${\displaystyle \mathbb {P} ^{x}}$. Here the indicator function ${\displaystyle \mathbb {I} (\tau <\infty )}$ 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 ${\displaystyle \tau }$, the penalty function is multiplied by the discounting factor ${\displaystyle e^{-\delta \tau }}$, and then averaged over the probability distribution of the waiting time to ${\displaystyle \tau }$. 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	${\displaystyle \mathbb {P} ^{x}\{\tau <\infty \}}$	${\displaystyle \delta =0,w(x_{1},x_{2})=1}$
Joint (defective) distribution of surplus and deficit	${\displaystyle \mathbb {P} ^{x}\{X_{\tau -}<x,X_{\tau }<y\}}$	${\displaystyle \delta =0,w(x_{1},x_{2})=\mathbb {I} (x_{1}<x,x_{2}<y)}$
Defective distribution of claim causing ruin	${\displaystyle \mathbb {P} ^{x}\{X_{\tau -}-X_{\tau }<z\}}$	${\displaystyle \delta =0,w(x_{1},x_{2})=\mathbb {I} (x_{1}+x_{2}<z)}$
Trivariate Laplace transform of time, surplus and deficit	${\displaystyle \mathbb {E} ^{x}[e^{-\delta \tau -sX_{\tau -}-zX_{\tau }}]}$	${\displaystyle w(x_{1},x_{2})=e^{-sx_{1}-zx_{2}}}$
Joint moments of surplus and deficit	${\displaystyle \mathbb {E} ^{x}[X_{\tau -}^{j}X_{\tau }^{k}]}$	${\displaystyle \delta =0,w(x_{1},x_{2})=x_{1}^{j}x_{2}^{k}}$

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.

• 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)

• Financial risk
• Volterra integral equation#Ruin theory

• Gerber, H.U. (1979). An Introduction to Mathematical Risk Theory. Philadelphia: S.S. Heubner Foundation Monograph Series 8.
• Asmussen S., Albrecher H. (2010). Ruin Probabilities, 2nd Edition. Singapore: World Scientific Publishing Co.