Ruin Theory and the Lundberg Bound

The classical risk model treats an insurer’s surplus as a stochastic process. Initial capital u, premium income at constant rate c per unit time, and aggregate claim outflow S(t) drive U(t) forward in time. Claim arrivals are Poisson and claim sizes are i.i.d., independent of the arrival process.

Hardy QERM, 2nd ed., Ch. 5 develops this model as the original quantitative risk model. ASTAM, GI 301, and CFE 101 each treat it as the foundation for solvency, capital adequacy, and reinsurance threshold thinking.

Ruin is the event that the surplus ever falls below zero. The probability of ultimate ruin starting from initial capital u is ψ(u). The complement φ(u) = 1 - ψ(u) is the survival probability.

A finite-horizon variant ψ(u, t) restricts ruin to occur on or before time t. Numerically ψ(u, t) is computed by inverting Laplace transforms or by Monte Carlo; analytically only ψ(u) = lim_{t → ∞} ψ(u, t) has a clean structural theory.

A necessary condition for finite ruin probability is the net profit condition c > λ E[X], equivalent to a positive premium loading θ > 0. If this fails, ruin is certain regardless of u.

The adjustment coefficient R is the unique positive root of the Lundberg equation, when one exists. R measures how fast the upper bound on ψ(u) decays as initial capital grows.

R exists only when the severity has a moment generating function in a neighborhood of zero. Heavy-tailed severities (Pareto, lognormal) do not have an MGF, so the classical Lundberg theory does not apply and asymptotic ruin probabilities decay subexponentially. For light-tailed severity (exponential, gamma, Weibull with shape > 1), R exists and is found by solving a single nonlinear equation.

The headline result is a clean exponential upper bound on the probability of ultimate ruin. The bound holds for all u ≥ 0 and depends on the severity distribution only through R.

Lundberg gives an upper bound, not the exact value. The Cramer-Lundberg asymptotic refines it: ψ(u) is asymptotically equivalent to a constant C times e^{-Ru}, where C ∈ (0, 1) depends on severity. Together, the inequality and the asymptotic mean that for light-tailed severities, doubling capital roughly squares the surviving margin against ruin.

When X is exponential with mean 1/μ, the ruin probability has a closed form. The adjustment coefficient is R = μ - λ/c, and the constant in the asymptotic equals λ/(cμ), which gives the exact answer for all u.

This is the only severity family with such a clean exact formula. It is the canonical sanity check for the Lundberg inequality (the bound holds with equality for u = 0 only when C = 1, which fails here for θ > 0) and the canonical opening exam question.

An insurer faces Poisson claims at rate λ = 2 per year with exponential severity having mean 1,000 (so μ = 0.001). Premium loading θ = 0.20, so c = 1.2 · λ · E[X] = 1.2 · 2 · 1,000 = 2,400 per year.

By the closed form, R = μ · θ/(1 + θ) = 0.001 · 0.20/1.20 = 1.667 × 10^{-4}. The Lundberg bound on ruin starting from u = 10,000 is ψ(u) ≤ exp(-1.667) ≈ 0.189, and the exact probability is ψ(u) = (1/1.20) · exp(-1.667) ≈ 0.157.

Interpretation: with $10,000 of capital and a 20% loading, the long-run probability of ever going bankrupt is about 16%. To get ψ(u) below 1%, solve (1/1.20) e^{-Ru} = 0.01 for u: e^{-Ru} = 0.012, so u ≥ ln(0.012^{-1})/R = ln(83.33)/R ≈ 26,540.

A reinsurer is willing to write a treaty if the probability of ruin over the long run is at most 0.5%. Severity is exponential with mean 50,000, claim rate λ = 10 per year, premium loading θ = 0.15.

Compute R: μ = 1/50,000 = 2 × 10^{-5}, R = μ · θ/(1 + θ) = 2 × 10^{-5} · 0.15/1.15 = 2.609 × 10^{-6}. Required u solves e^{-Ru} ≤ 0.005 (using the Lundberg bound, the conservative direction). u ≥ ln(200)/R = 5.298/(2.609 × 10^{-6}) ≈ 2.03 million.

If the reinsurer instead uses the Cramer-Lundberg exact, the required capital drops to ln(200/1.15)/R ≈ ln(173.9)/R ≈ 5.158/2.609 × 10^{-6} ≈ 1.98 million. The bound is conservative by about 2.5% here. For tail-heavy severities the gap widens; for exponential severity it stays small.

Severity is Gamma(α = 2, θ = 500), so E[X] = αθ = 1,000 and the MGF is M_X(r) = (1 - θ r)^{-α} = (1 - 500r)^{-2} for r < 1/500. Claim rate λ = 2 per year, loading θ_load = 0.20, so c = 2,400.

Solve 2 · ((1 - 500r)^{-2} - 1) = 2,400 r. Rewrite: (1 - 500r)^{-2} = 1 + 1,200r. At r = 0.0006: LHS = (1 - 0.3)^{-2} = (0.7)^{-2} = 2.041; RHS = 1 + 0.72 = 1.720. LHS > RHS, so try smaller r. At r = 0.0004: LHS = (1 - 0.2)^{-2} = 1.5625; RHS = 1.48. Still LHS > RHS. At r = 0.0003: LHS = (1 - 0.15)^{-2} = 1.384; RHS = 1.36. Closer. At r = 0.00028: LHS = 1.353; RHS = 1.336. At r = 0.00025: LHS = (0.875)^{-2} = 1.306; RHS = 1.30. The root is approximately R ≈ 2.45 × 10^{-4}.

Capital required for ψ(u) ≤ 1% via Lundberg bound: u ≥ ln(100)/R = 4.605/(2.45 × 10^{-4}) ≈ 18,800. Compare to exponential severity with the same mean, same loading: R = 0.001 · 0.20/1.20 = 1.667 × 10^{-4}, requiring u ≥ ln(100)/1.667 × 10^{-4} ≈ 27,600. Gamma’s lighter tail (CV = 1/√2 < 1) gives a larger R and a smaller capital requirement.

Trap 1: Applying Lundberg to a heavy-tailed severity. Pareto and lognormal severities have no MGF in any neighborhood of zero, so R does not exist and the exponential decay structure breaks down. Ruin probabilities for heavy-tailed severities decay polynomially, not exponentially, and require subexponential ruin theory (Embrechts et al., Modelling Extremal Events) instead.

Trap 2: Forgetting the net profit condition c > λ E[X]. If premium income does not exceed expected claim outflow, ruin is certain (ψ(u) = 1 for all u), and there is no positive R to solve for. The Lundberg equation has only the trivial root r = 0 in that regime.

Trap 3: Confusing the Lundberg upper bound with the exact ruin probability. The exact ψ(u) is strictly less than e^{-Ru} for u > 0 except in degenerate cases. Using the bound where the question asks for the exact (or vice versa) is a routine ASTAM grading deduction.

The Lundberg inequality and the closed-form exponential result are infinite-horizon statements. They tell you the probability of ever being ruined.

Many regulatory and reinsurance applications use a finite horizon (one year, five years, runoff to maturity). ψ(u, t) requires either a numerical Laplace inversion of the Pollaczek-Khinchine transform or a Monte Carlo simulation of the surplus path. For typical insurance horizons (one to five years) the finite-horizon ruin probability is often a small fraction of the ultimate ruin probability, since ruin is a tail-time event for well-loaded portfolios.

The Lundberg adjustment coefficient is the cleanest single number summarizing how well an insurer’s premium loading offsets its severity tail. It feeds capital adequacy work in CFE 101, ERM frameworks in CERA, and ruin and reserving questions on ASTAM and GI 301.

Hardy, QERM, 2nd ed., Ch. 5 derives the surplus process, the Lundberg inequality, and the Cramer-Lundberg asymptotic; Klugman, Loss Models, 5th ed., Ch. 11 covers the same material from a long-term-actuarial-mathematics angle. Both are core ASTAM references.