Insurance Surplus Modeling with Poisson Processes
Exploring ruin probability through stochastic processes and heavy-tailed distributions
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
- When does classical risk theory provide accurate predictions?
- When does it fail catastrophically?
- 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 is modeled as:
where:
- = initial capital (starting reserve)
- = cumulative premium income at rate
- = number of claims by time
- = size of the -th claim
- = total claims paid (compound Poisson process)
Model Assumptions
1. Claim arrivals:
- Claims arrive according to a Poisson process with rate
- 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: are i.i.d. with distribution and mean
- Common choices for : Exponential (light-tailed), Pareto (heavy-tailed), Gamma, Log-Normal
3. Independence: and are independent
4. Net profit condition:
- Premium income exceeds expected claim outflow
- Derived from
Part II: Ruin Theory
Defining Ruin
Time of Ruin:
The first time the surplus becomes negative. If never goes negative, then .
Ultimate Ruin Probability:
The probability of eventual ruin starting with initial capital .
Key Properties
- (more capital = safer)
- if (certain ruin without net profit condition)
- Question: How fast does as ?
Part III: Lundberg's Inequality
The Classical Result
Theorem (Lundberg): Assuming the net profit condition holds, if there exists satisfying:
then:
where is a constant depending on .
Interpretation
- is called the adjustment coefficient
- Quantifies company safety: higher premium rate → larger → 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:
- Given the adjustment coefficient equation:
- Define the exponential martingale:
- Apply Optional Stopping Theorem at time :
- At ruin:
- Therefore: , which gives
Part IV: When Theory Fails - Heavy Tails
The Critical Assumption
Lundberg's inequality requires:
This assumes the moment generating function exists.
Heavy-Tailed Distributions
When for all :
- No adjustment coefficient exists
- MGF does not exist
- Tails decay polynomially (much slower than exponential)
- Examples: Pareto, Log-Normal distributions
Exponential vs Polynomial Decay
| Distribution | Tail Behavior | Ruin Decay |
|---|---|---|
| Exponential | ||
| Pareto() |
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 Capital | Exponential | Pareto | Ratio (P/E) |
|---|---|---|---|
| $0 | 81.0% | 79.8% | 1.0× |
| $50,000 | 22.8% | 15.2% | 0.7× |
| $100,000 | 6.3% | 4.8% | 0.8× |
| $150,000 | 1.9% | 1.8% | 0.9× |
| $200,000 | 0.530% | 0.930% | 1.8× |
| $250,000 | 0.130% | 0.600% | 4.6× |
| $300,000 | 0.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
| Characteristic | Light-Tailed (Exponential) | Heavy-Tailed (Pareto) |
|---|---|---|
| Tail Decay | Exponential: P(X > x) ~ e^(-λx) | Polynomial: P(X > x) ~ x^(-α) |
| MGF | Exists for all R > 0 | Does not exist (infinite for R > 0) |
| Ruin Probability | ψ(u) ~ e^(-Ru) | ψ(u) ~ u^(-(α-1)) |
| Lundberg Bound | Valid and tight | Does not exist |
| Typical Ruin | Many small claims accumulate | One giant catastrophic claim |
| Capital Efficiency | Logarithmic growth in safety | Polynomial growth needed |
| Real Examples | Auto insurance, standard property | Hurricanes, earthquakes, terrorism |
Limitations & Future Work
Current Limitations
- Simple assumptions: Pure compound process (no investment returns)
- No risk management: No reinsurance, capital injections, or dividends
- Single tail type: Most insurers face both light-tailed and sudden catastrophic claims
- Parameter certainty: In reality, limited data makes , , 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
- Asmussen, S., & Albrecher, H. (2010). Ruin Probabilities (2nd ed.). World Scientific.
- Cramér, H. (1930). On Collective Risk Theory. Skandia Jubilee Volume.
- Lundberg, F. (1903). Approximations of the Probability Function. Almqvist & Wiksell.
- Ross, S. M. (2014). Introduction to Probability Models (11th ed.). Academic Press.
- 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.
- 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.