Back to Mathematics
ActuarialFall 2025 - Directed Reading Program

Insurance Surplus Modeling with Poisson Processes

Exploring ruin probability through stochastic processes and heavy-tailed distributions

Stochastic ProcessesActuarial ScienceMonte Carlo

Overview

This project investigates how an insurance company's financial reserve evolves over time when claims arrive randomly—a fundamental problem in actuarial science and applied probability. By modeling claim arrivals as a Poisson process and analyzing the resulting jump process, we explore the critical question: What is the probability that the company eventually runs out of money?

Key Focus: Understanding when classical risk theory works, when it fails, and what that failure costs in terms of required capital.

Motivation

Real-World Context

Hurricane Katrina (2005)

Losses: $41 billion

Impact: Multiple insurers went bankrupt

Heavy-tailed risks can wipe out years of profit in one event

AIG Crisis (2008)

Losses: $85 billion bailout

Impact: Near-total collapse

Underestimating tail risk leads to systemic failure

These events expose a critical gap in classical insurance theory: traditional models assume claim sizes follow light-tailed distributions, but reality often produces heavy-tailed catastrophic events.

The Central Questions

  1. When does classical risk theory provide accurate predictions?
  2. When does it fail catastrophically?
  3. How much more capital is needed to account for heavy-tailed risks?

Part I: The Mathematical Model

The Cramér-Lundberg Model

The surplus (reserve) of an insurance company at time tt is modeled as:

U(t)=u+cti=1N(t)XiU(t) = u + ct - \sum_{i=1}^{N(t)} X_i

where:

  • uu = initial capital (starting reserve)
  • ctct = cumulative premium income at rate cc
  • N(t)N(t) = number of claims by time tt
  • XiX_i = size of the ii-th claim
  • i=1N(t)Xi=S(t)\sum_{i=1}^{N(t)} X_i = S(t) = total claims paid (compound Poisson process)

Model Assumptions

1. Claim arrivals: N(t)Poisson(λ)N(t) \sim \text{Poisson}(\lambda)

  • Claims arrive according to a Poisson process with rate λ\lambda
  • Why Poisson? It's the unique continuous-time process with independent, stationary increments and discrete jumps, plus the memoryless property makes calculations tractable

2. Claim sizes: XiX_i are i.i.d. with distribution FF and mean μ\mu

  • Common choices for FF: Exponential (light-tailed), Pareto (heavy-tailed), Gamma, Log-Normal

3. Independence: XiX_i and N(t)N(t) are independent

4. Net profit condition: c>λμc > \lambda\mu

  • Premium income exceeds expected claim outflow
  • Derived from E[U(t)]=u+(cλμ)tE[U(t)] = u + (c - \lambda\mu)t

Part II: Ruin Theory

Defining Ruin

Time of Ruin:

τ=inf{t>0:U(t)<0}\tau = \inf\{t > 0 : U(t) < 0\}

The first time the surplus becomes negative. If U(t)U(t) never goes negative, then τ=\tau = \infty.

Ultimate Ruin Probability:

ψ(u)=P(τ<U(0)=u)\psi(u) = P(\tau < \infty \mid U(0) = u)

The probability of eventual ruin starting with initial capital uu.

Key Properties

  • limuψ(u)=0\lim_{u \to \infty} \psi(u) = 0 (more capital = safer)
  • limu0ψ(u)=1\lim_{u \to 0} \psi(u) = 1 if cλμc \leq \lambda\mu (certain ruin without net profit condition)
  • Question: How fast does ψ(u)0\psi(u) \to 0 as uu \to \infty?

Part III: Lundberg's Inequality

The Classical Result

Theorem (Lundberg): Assuming the net profit condition c>λμc > \lambda\mu holds, if there exists R>0R > 0 satisfying:

λE[eRX]=λ+cR\lambda \cdot E[e^{RX}] = \lambda + cR

then:

ψ(u)CeRu\psi(u) \leq C \cdot e^{-Ru}

where CC is a constant depending on RR.

Interpretation

  • RR is called the adjustment coefficient
  • Quantifies company safety: higher premium rate → larger RR → safer
  • Exponential decay: Ruin probability decreases exponentially with initial capital
  • This assumes light-tailed distributions (exponential safety)

Proof Sketch: Martingale Method

The proof uses the Optional Stopping Theorem:

  1. Given the adjustment coefficient equation: λE[eRX]=λ+cR\lambda E[e^{RX}] = \lambda + cR
  2. Define the exponential martingale: M(t)=eRU(t)eλt(E[eRX]1)M(t) = e^{RU(t)} \cdot e^{-\lambda t(E[e^{RX}] - 1)}
  3. Apply Optional Stopping Theorem at time τ\tau: E[M(τ)]=M(0)=eRuE[M(\tau)] = M(0) = e^{Ru}
  4. At ruin: U(τ)0    eRU(τ)1U(\tau) \leq 0 \implies e^{RU(\tau)} \leq 1
  5. Therefore: eRuP(τ<)e^{Ru} \geq P(\tau < \infty), which gives ψ(u)eRu\psi(u) \leq e^{-Ru}

Part IV: When Theory Fails - Heavy Tails

The Critical Assumption

Lundberg's inequality requires:

E[eRX]<for some R>0E[e^{RX}] < \infty \quad \text{for some } R > 0

This assumes the moment generating function exists.

Heavy-Tailed Distributions

When E[eRX]=E[e^{RX}] = \infty for all R>0R > 0:

  • No adjustment coefficient exists
  • MGF does not exist
  • Tails decay polynomially (much slower than exponential)
  • Examples: Pareto, Log-Normal distributions

Exponential vs Polynomial Decay

DistributionTail BehaviorRuin Decay
ExponentialP(X>x)eλxP(X > x) \sim e^{-\lambda x}ψ(u)eRu\psi(u) \sim e^{-Ru}
Pareto(α\alpha)P(X>x)xαP(X > x) \sim x^{-\alpha}ψ(u)u(α1)\psi(u) \sim u^{-(\alpha-1)}

Key Insight: Exponential decay is dramatically faster than polynomial decay.

Part V: Simulation Results

Interactive Visualization

Ruin Probability vs Initial Capital

Note: Y-axis uses logarithmic scale to better show exponential decay

Detailed Results Table

Initial CapitalExponentialParetoRatio (P/E)
$081.0%79.8%1.0×
$50,00022.8%15.2%0.7×
$100,0006.3%4.8%0.8×
$150,0001.9%1.8%0.9×
$200,0000.530%0.930%1.8×
$250,0000.130%0.600%4.6×
$300,0000.066%0.530%8.8×

All values represent probability of eventual ruin over a 50-year horizon with 10,000 Monte Carlo simulations

Part VI: Analysis & Interpretation

Key Findings

Exponential Distribution

  • Ruin probability decays exponentially: ψ(u) ≤ Ce^(-Ru)
  • Simulation closely matches Lundberg theoretical bound
  • At $200K capital: 0.53% ruin probability
  • Safe regime: theory accurately predicts behavior

Implication: Classical risk theory works well for light-tailed claims

Pareto Distribution

  • Ruin probability decays polynomially (much slower)
  • No exponential bound exists
  • At $200K capital: 0.93% ruin probability
  • Ratio (Pareto/Exponential) increases with capital: 0.7× → 8.8×

Implication: Heavy tails require substantially more capital at high safety levels

Capital Gap

  • For 1% ruin target: Exponential needs $183K, Pareto needs $196K
  • 7% more capital required for Pareto
  • Gap widens dramatically at lower ruin probabilities
  • At 0.1% target: difference would be much larger

Implication: Tail risk premium increases with safety requirements

The "Tail Risk Premium"

The extra capital needed for heavy-tailed risks:

  • At 1% ruin: Only 7% more capital
  • At 0.1% ruin: Difference would be much larger
  • The gap widens at higher safety levels

Light-Tailed vs Heavy-Tailed: Complete Comparison

CharacteristicLight-Tailed (Exponential)Heavy-Tailed (Pareto)
Tail DecayExponential: P(X > x) ~ e^(-λx)Polynomial: P(X > x) ~ x^(-α)
MGFExists for all R > 0Does not exist (infinite for R > 0)
Ruin Probabilityψ(u) ~ e^(-Ru)ψ(u) ~ u^(-(α-1))
Lundberg BoundValid and tightDoes not exist
Typical RuinMany small claims accumulateOne giant catastrophic claim
Capital EfficiencyLogarithmic growth in safetyPolynomial growth needed
Real ExamplesAuto insurance, standard propertyHurricanes, earthquakes, terrorism

Limitations & Future Work

Current Limitations

  1. Simple assumptions: Pure compound process (no investment returns)
  2. No risk management: No reinsurance, capital injections, or dividends
  3. Single tail type: Most insurers face both light-tailed and sudden catastrophic claims
  4. Parameter certainty: In reality, limited data makes α\alpha, λ\lambda, cc hard to estimate

Proposed Extensions

1. Regime-Switching Models

  • Light-tailed in normal times
  • Transition to heavy-tailed during crises
  • Captures realistic dynamics

2. Optimal Reinsurance

  • Transfer heavy-tail risks to reinsurers
  • Capital injection strategies
  • Minimize ruin probability vs. cost

3. Multi-Line Insurance

  • Multiple correlated portfolios
  • Diversification effects
  • Contagion risks

4. Data-Driven Parameter Estimation

  • Bayesian estimation with limited data
  • Uncertainty quantification
  • Robust optimization under parameter uncertainty

Acknowledgments

This project was completed as part of the Directed Reading Program (DRP) in Fall 2025. Special thanks to Daniel Viray for his invaluable guidance and mentorship throughout this research.

Resources

References

  1. Asmussen, S., & Albrecher, H. (2010). Ruin Probabilities (2nd ed.). World Scientific.
  2. Cramér, H. (1930). On Collective Risk Theory. Skandia Jubilee Volume.
  3. Lundberg, F. (1903). Approximations of the Probability Function. Almqvist & Wiksell.
  4. Ross, S. M. (2014). Introduction to Probability Models (11th ed.). Academic Press.
  5. Barrera, M., Rojas, F., & Villaseñor, J. A. (2020). On the ruin probability of a generalized Cramér–Lundberg model driven by mixed Poisson processes.
  6. Aurzada, F. & Buck, M. (2020). Ruin probabilities with insurance and financial risks having temporarily negative capital.

This project was completed as part of the Directed Reading Program (DRP) in Fall 2025, exploring fundamental questions in actuarial mathematics and stochastic processes.