Actuarial Principles and Insurance Mathematics

Chapter 8
Aggregate Loss, Risk and Reinsurance Models

Post Graduate In Insurance Management
Actuarial Principles and Insurance Mathematics   •   Post Graduate

8. Aggregate Loss, Risk and Reinsurance Models

Introduction and Context

Insurance pricing and solvency decisions are rarely based on a single claim. An insurer normally receives a large portfolio of risks, each producing zero, one, or several claims of uncertain size. The central actuarial problem is therefore to understand the distribution of the total loss arising from a portfolio during a specified period. This total, commonly called aggregate loss, determines premium adequacy, capital requirements, reinsurance needs, profitability, and the probability that an insurer will be unable to meet its obligations.

Two complementary perspectives are used. The collective risk model treats the portfolio as a random number of claims with independently and identically distributed severities. The individual risk model treats each policy separately and aggregates the losses generated by all policies. Both models provide a foundation for pricing, reserving, capital modelling, catastrophe analysis, and enterprise risk management.

Reinsurance modifies the insurer’s retained aggregate loss. Proportional arrangements share premiums and losses according to fixed proportions, whereas non-proportional arrangements respond to claims or aggregate losses exceeding specified thresholds. The choice among these structures involves a trade-off between risk reduction, cost, moral hazard, claims volatility, counterparty exposure, and operational complexity.

Learning Objectives

After studying this chapter, the student should be able to:

1. Define aggregate loss and distinguish frequency, severity, and dependence risk.

2. Construct collective and individual risk models.

3. Derive moments and probability distributions of aggregate losses.

4. Explain convolution, probability generating functions, recursive methods, approximation, and simulation.

5. Interpret ruin probability and the classical compound Poisson risk model.

6. Analyse proportional, excess-of-loss, stop-loss, and catastrophe reinsurance.

7. Evaluate reinsurance alternatives using quantitative and managerial criteria.

8. Recognise model risk, dependence, data limitations, governance requirements, and emerging analytical practices.

1. The Architecture of an Aggregate-Loss Model

Let the time horizon be one year. For a portfolio, define:

  • N = number of claims during the year;
  • X_i = amount of the ith claim;
  • S = total or aggregate loss.

In the collective model,

S = X_1 + X_2 + ... + X_N,

where the number of claims N is random and the severities X_i are random variables, usually assumed initially to be independent of N and mutually independent. If N = 0, then S = 0.

The distribution of S is influenced by two separate sources of uncertainty. Frequency risk arises from uncertainty in N. Severity risk arises from uncertainty in the amount of each claim. A portfolio may have a stable expected claim count but highly volatile losses because of a heavy-tailed severity distribution. Conversely, claim amounts may be predictable while frequency fluctuates materially.

The expected aggregate loss follows from conditional expectation:

E(S) = E[E(S|N)] = E(N)E(X),

provided the claims are identically distributed and independent of N. The variance is

Var(S) = E(N)Var(X) + Var(N)[E(X)]².

The first term represents variation in claim sizes conditional on the number of claims. The second represents variation caused by the number of claims itself. This decomposition is useful in portfolio diagnostics: a frequency-control programme and a claims-severity programme address different components of aggregate volatility.

If the claims are not identically distributed, then

E(S) = E\left(\sum_{i=1}^{N}X_i\right),

and the simple formulas may need to be replaced by exposure-based or segmented calculations. Dependence among claims also requires care. Claims arising from the same flood, earthquake, cyber event, or economic shock are unlikely to be independent.

2. Collective and Individual Risk Models

2.1 The collective risk model

The collective model begins with a frequency distribution for N and a common severity distribution F_X. Typical frequency models include the Poisson, binomial, and negative binomial distributions.

For a Poisson frequency with parameter lambda,

P(N=n) = e^(-lambda) lambda^n/n!, n = 0,1,2,...

The aggregate loss is then a compound Poisson random variable. Its probability generating function, when claim amounts are integer-valued, is

G_S(z) = G_N[G_X(z)] = exp{lambda[G_X(z)-1]}.

For a compound Poisson model,

E(S) = lambda E(X),

Var(S) = lambda E(X²).

The compound Poisson structure is attractive because it is mathematically tractable and often a reasonable first model for independent claim events. It is not automatically appropriate for catastrophe portfolios, highly exposed commercial accounts, or portfolios with strong policyholder heterogeneity.

2.2 The individual risk model

Suppose a portfolio contains m policies. Let Y_j denote the loss from policy j during the period. The aggregate loss is

S = Y_1 + Y_2 + ... + Y_m.

A simple binary individual risk model assumes that policy j produces a claim amount c_j with probability q_j and zero otherwise. Thus,

Y_j = c_j I_j,

where I_j is a Bernoulli indicator. If the risks are independent,

E(S) = sum from j=1 to m of q_j c_j,

Var(S) = sum from j=1 to m of q_j(1-q_j)c_j².

More generally, Y_j may have a full distribution, including several possible claim amounts. The individual model is particularly useful when policy characteristics differ materially, when policy limits are important, or when the number of policies is manageable. It permits explicit recognition of heterogeneity, although computation becomes more difficult as the portfolio grows.

The collective model aggregates similar risks before modelling total loss. The individual model preserves policy-level information. Neither is universally superior. The collective model is often efficient for large homogeneous portfolios, while the individual model can be more transparent for underwriting, capital allocation, and portfolio concentration analysis.

3. Aggregate-Loss Distributions and Computation

3.1 Convolution

If N is fixed at n and the claims are independent, the distribution of S conditional on N=n is the n-fold convolution of the severity distribution. For continuous severities, the density of the sum of two claims is

f_{X_1+X_2}(s) = integral from 0 to s of f_X(x)f_X(s-x) dx.

The unconditional aggregate distribution is a mixture:

P(S in A) = sum over n of P(N=n)P(X_1+...+X_n in A).

Convolution is conceptually direct but may be computationally expensive, especially for continuous distributions or large claim counts.

3.2 Recursive calculation and generating functions

When losses are discretised into monetary units, recursive algorithms can calculate aggregate probabilities efficiently. If the severity probability generating function is G_X(z), the compound Poisson generating function given above provides the aggregate probabilities through its coefficients.

A practical actuarial implementation normally rounds or groups claim amounts into bands. Discretisation introduces approximation error, but a sufficiently fine grid can provide accurate estimates of deductibles, quantiles, and tail probabilities. Recursive methods are especially valuable when a full distribution is required rather than only its mean and variance.

3.3 Approximation methods

For a large portfolio of relatively small, independent risks, the Central Limit Theorem suggests

S approximately follows a Normal distribution with mean E(S) and variance Var(S).

The normal approximation is useful for central probabilities but can be unreliable for skewed insurance losses. Aggregate losses are frequently right-skewed, and a normal distribution assigns positive probability to negative losses and may underestimate extreme quantiles.

A translated normal or normal-power approximation may improve the fit by incorporating skewness. A gamma approximation is also often suitable for positive, right-skewed aggregates. For extreme tail analysis, analysts may use heavy-tailed severity models, saddlepoint approximations, extreme-value methods, or a fitted parametric distribution. The choice must be validated against observed data and the intended use of the model.

3.4 Simulation

Monte Carlo simulation is flexible when analytical calculation is difficult. A basic aggregate-loss simulation consists of the following steps:

1. Generate a claim count N from the selected frequency distribution.

2. Generate N claim severities from the selected severity distribution.

3. Sum the simulated severities to obtain one aggregate loss.

4. Repeat the process many times.

5. Estimate the mean, variance, quantiles, tail probabilities, and reinsurance recoveries from the simulated sample.

Simulation can incorporate policy limits, deductibles, inflation, seasonal effects, dependence, multiple lines of business, and reinsurance terms. Its principal disadvantages are sampling error, computational cost, and the possibility of producing precise-looking results from an inappropriate model. Variance-reduction methods, such as importance sampling, are useful when estimating very small ruin probabilities.

4. Risk Measures for Aggregate Loss

The mean aggregate loss describes expected cost but not uncertainty. Important risk measures include:

  • Variance and standard deviation, which measure dispersion around the mean.
  • Value at Risk (VaR) at level alpha, defined as the alpha-quantile of S.
  • Tail Value at Risk (TVaR), also called Expected Shortfall, defined as the expected loss conditional on exceeding the selected VaR threshold, subject to the convention used at atoms.
  • Probability of exceeding a retention or capital level.

VaR is easy to communicate but does not describe the average severity of losses beyond the threshold and may not be subadditive in general. TVaR provides more information about tail severity and is often better suited to capital and reinsurance analysis.

A risk measure should be interpreted in relation to the time horizon, confidence level, valuation basis, and treatment of dependence. A 99.5 per cent one-year capital measure is not directly comparable with a 99 per cent ten-year measure.

5. Ruin and the Classical Risk Process

In the classical insurance risk model, the insurer’s surplus at time t is

U(t) = u + ct - S(t),

where u is initial surplus, c is the premium income rate, and S(t) is cumulative aggregate claims by time t. If claims arrive according to a Poisson process and have independent severities, this is the Cramér–Lundberg model.

Ruin occurs when U(t) becomes negative. The ultimate ruin probability is

psi(u) = P[U(t)<0 for some t greater than or equal to 0 | U(0)=u].

The safety loading condition is essential. If c is no greater than the expected claim rate, long-run ruin is unavoidable under standard assumptions. If c exceeds expected claims, the insurer has a positive expected surplus drift, although ruin remains possible because of random large claims.

The adjustment coefficient R is a positive solution, when it exists, of

E(e^{R X}) = 1 + (c/lambda)R.

Under suitable assumptions, the ultimate ruin probability has an exponential upper bound of the form

psi(u) less than or equal to e^(-Ru).

This bound illustrates the effect of capital: greater initial surplus reduces ruin probability, but the rate of reduction depends on the frequency, severity, and premium structure. In practice, finite-horizon ruin probabilities, investment returns, inflation, claim settlement delays, dependence, and operational risks must also be considered. The classical model is a benchmark rather than a complete solvency model.

6. Reinsurance as Transformation of Aggregate Risk

Reinsurance transfers part of an insurer’s risk to another insurer or reinsurer in return for premium. It can reduce earnings volatility, protect capital, support underwriting capacity, and provide catastrophe protection. It also introduces ceded premium costs, counterparty credit risk, collateral issues, basis risk, and contract interpretation risk.

6.1 Proportional reinsurance

Under quota-share reinsurance, the reinsurer accepts a fixed proportion q of every covered premium and loss. If the original claim is X, the ceded claim is qX and the insurer retains (1-q)X. Ignoring expenses,

E(S_retained) = (1-q)E(S),

and

Var(S_retained) = (1-q)² Var(S).

Quota share is simple and provides automatic protection for frequency and severity. It may be useful for a new line of business, a capital-constrained insurer, or a portfolio for which the cedant wants underwriting support. However, it cedes profitable as well as unprofitable business and may provide less efficient protection against a small number of extreme claims.

Under surplus-share reinsurance, the cedant retains a fixed monetary line and cedes amounts above that line, usually subject to a maximum number of lines. It is more responsive to policy-size variation than quota share.

6.2 Per-risk excess of loss

A per-risk excess-of-loss treaty responds when an individual claim exceeds a retention d. The reinsurer’s recovery for claim X is

R(X) = (X-d)+ = max(X-d,0),

subject to any treaty limit. The insurer retains min(X,d), before considering aggregate limits or reinstatements.

The expected recovery is

E[(X-d)+] = integral from d to infinity of (1-F_X(x)) dx.

This identity, often called the stop-loss premium formula, shows that the cost of protection depends on the severity tail beyond the retention. Per-risk excess of loss is appropriate when the main concern is unusually large individual claims, such as major property, liability, or engineering losses.

6.3 Catastrophe excess of loss

Catastrophe excess of loss responds to the aggregate losses arising from one event, such as a windstorm or earthquake. If event loss is Z, a layer with attachment a and limit l pays

R(Z) = min[(Z-a)+, l].

The attachment point and limit should reflect the insurer’s event-loss profile and solvency tolerance. Catastrophe covers often include event definitions, hours clauses, geographical restrictions, occurrence aggregation rules, exclusions, reinstatement provisions, and reporting requirements. Modelling must account for event frequency, event severity, geographical accumulation, and multiple policies affected by the same event.

6.4 Aggregate excess of loss and stop-loss reinsurance

Aggregate excess of loss protects the insurer when total losses over a period exceed an attachment point a. The recovery is commonly

R(S) = min[(S-a)+, l].

Stop-loss reinsurance is conceptually similar but is often expressed as protection for the loss ratio or aggregate loss above a specified threshold. It is designed to stabilise overall underwriting results rather than protect only individual claims.

For an unlimited stop-loss layer,

E[(S-a)+] = integral from a to infinity of P(S>x) dx.

This formula links reinsurance pricing directly to the aggregate-loss survival function. Stop-loss cover may be valuable for volatile portfolios, but it can be expensive and may create disputes concerning earned premium, loss development, accounting periods, and catastrophe aggregation.

7. Pricing and Evaluation of Reinsurance

A theoretical pure premium for a reinsurance layer is the expected recovery. A commercial premium normally adds a loading for expenses, capital, profit, uncertainty, brokerage, collateral, and the reinsurer’s cost of capital. Common approaches include expected-value pricing, variance or standard-deviation loading, risk-measure pricing, burning-cost pricing, exposure rating, and catastrophe-model output.

Burning-cost pricing uses historical ceded losses, adjusted for exposure changes, inflation, claims development, and large-loss distortions. Exposure rating estimates losses from an external severity curve and the cedant’s exposure profile. Catastrophe models combine hazard, vulnerability, financial terms, and event simulation.

A sound evaluation compares the retained and ceded distributions rather than considering price alone. Relevant measures include expected retained loss, standard deviation, VaR, TVaR, probability of breaching capital, return on capital, and the effect on earnings volatility. The analysis should also test alternative assumptions about inflation, claims development, parameter uncertainty, reinstatement premiums, and reinsurer default.

8. Applied Case: Reinsurance Design for a Regional Property Portfolio

Background

Consider an illustrative insurer writing commercial property risks in a region exposed to ordinary fire losses and occasional severe storms. The insurer has adequate capital for normal claims but is concerned that one storm may generate losses across many policies. Management is considering quota share, per-risk excess of loss, and catastrophe excess of loss.

Situation and problem

Historical data suggest that most claims are small or moderate, while storm losses are highly clustered. A quota-share treaty would reduce all losses but would also transfer a substantial portion of profitable ordinary business. A per-risk treaty would protect against large individual fires but would not fully address a storm affecting hundreds of properties. The central problem is therefore one of matching the reinsurance structure to the dominant risk mechanism.

Analysis

The actuary first separates ordinary claims from catastrophe events and models them differently. Ordinary frequency is modelled using segmented policy data, while storm losses are generated from event-level simulations incorporating geographical accumulation. The retained-loss distribution is calculated under each alternative.

The quota-share option produces a substantial reduction in both expected loss and variance, but the insurer gives up premium and underwriting margin on every policy. The per-risk excess-of-loss option sharply reduces the effect of large individual claims but leaves considerable event aggregation. The catastrophe layer has little effect on ordinary annual volatility but materially reduces the upper tail caused by a severe storm.

Strategy and solution

Management selects a layered programme: a modest quota share for capital support during portfolio growth, per-risk excess of loss for large individual claims, and catastrophe excess of loss for event accumulation. The programme includes clear event definitions, an appropriate attachment point, a limit aligned with modelled probable maximum loss, reinstatement terms, and counterparty security standards.

Outcome

The combined structure retains ordinary underwriting margin while reducing the probability that a single event will threaten solvency. The insurer also establishes monitoring triggers for geographical accumulation, inflation, changes in construction quality, and reinsurer credit quality. The result is not simply a lower expected loss; it is a more controlled distribution of retained risk.

Learning points

The case demonstrates that reinsurance should be designed around the source and shape of risk. A treaty that is effective for independent large claims may be ineffective for correlated event losses. Quantitative optimisation must therefore be integrated with underwriting, capital, legal, and counterparty analysis.

9. Limitations, Model Risk, and Governance

Aggregate-loss models can fail because of inadequate data, incorrect distributional assumptions, parameter uncertainty, structural change, and dependence that is not represented. A fitted Poisson distribution may understate frequency variation if risks are heterogeneous. A light-tailed severity distribution may seriously underestimate extreme losses. Historical data may contain no observation comparable with a future catastrophe.

Common mistakes include confusing expected loss with required capital, using a normal approximation for a strongly skewed portfolio, ignoring policy limits and deductibles, treating catastrophe claims as independent, and evaluating reinsurance solely by its quoted premium. Another error is to model gross losses carefully but apply reinsurance terms incorrectly, especially aggregate definitions, hours clauses, reinstatements, exclusions, and claims-development provisions.

Model governance should include documented assumptions, independent validation, sensitivity testing, back-testing where possible, scenario analysis, version control, and clear ownership. Senior management should understand the uncertainty around outputs rather than treating a single estimated percentile as an objective fact. Regulatory solvency frameworks generally require disciplined assessment of underwriting, catastrophe, market, credit, and operational risks, although the precise requirements vary by jurisdiction.

Ethically, actuaries and managers must avoid deliberately optimistic assumptions designed to reduce apparent capital needs or justify inadequate reinsurance. The interests of policyholders, shareholders, cedants, reinsurers, and broader financial stability should be considered. Contract clarity is also an ethical and governance issue: ambiguous terms can create disputes precisely when protection is most needed.

10. Technology, Analytics, and Emerging Practice

Modern aggregate-loss modelling relies increasingly on simulation platforms, distributed computing, cloud infrastructure, geospatial data, and automated data pipelines. Machine-learning methods can improve frequency segmentation, anomaly detection, claims triage, and exposure classification. However, predictive accuracy does not remove the need for actuarial judgement. Black-box models may be difficult to explain, unstable under changing conditions, or unfair to particular groups.

Digital risk creates new challenges. Cyber losses may be correlated across firms through common software, cloud providers, or supply chains, making historical independence assumptions unreliable. Climate change may alter hazard frequency, severity, and geographical patterns, weakening the assumption that historical catastrophe data represent the future. Scenario analysis, stress testing, expert judgement, and model ensembles are consequently becoming more important.

Best practice combines several elements: sound exposure data; transparent frequency and severity assumptions; explicit treatment of dependence; validation against internal and external evidence; comparison of analytical and simulation results; reinsurance contract testing; and regular review after material portfolio or environmental changes. The most useful model is not necessarily the most complex one, but the one whose assumptions are understood, whose limitations are disclosed, and whose outputs support better decisions.

Chapter Summary

Aggregate-loss modelling combines claim frequency and claim severity to estimate the distribution of total insurance losses. The collective model treats claims as a random number of random severities, while the individual model aggregates policy-level losses. Expected aggregate loss and variance can often be derived using conditional expectation and variance, but full distributions may require convolution, recursive methods, approximation, or simulation.

Ruin models connect aggregate losses with premium income and capital, showing how adverse claims can exhaust surplus even when expected premiums exceed expected claims. Reinsurance transforms the retained-loss distribution. Quota share and surplus share are proportional structures; per-risk excess of loss protects against large individual claims; catastrophe excess of loss addresses event accumulation; and aggregate excess of loss or stop-loss protects against poor overall results. Reinsurance evaluation should consider tail risk, capital, price, counterparty quality, contract terms, and model uncertainty. Effective actuarial practice requires quantitative analysis combined with governance, professional judgement, and continual validation.

Key Terms

Aggregate loss: The total amount of claims incurred during a specified period.

Frequency: The number of claims occurring in a period.

Severity: The monetary amount of an individual claim.

Collective risk model: A model representing aggregate loss through a random claim count and claim severities.

Individual risk model: A model that aggregates losses separately for each policy or exposure unit.

Compound Poisson model: An aggregate model with Poisson claim frequency and independent claim severities.

Convolution: The operation used to derive the distribution of a sum of independent random variables.

Probability generating function: A function that represents the probabilities of an integer-valued random variable.

Monte Carlo simulation: Repeated random generation used to approximate a loss distribution or risk measure.

Value at Risk: A selected quantile of a loss distribution.

Tail Value at Risk: The expected loss beyond a specified tail quantile.

Ruin probability: The probability that an insurer’s surplus becomes negative.

Adjustment coefficient: A parameter used to obtain exponential bounds for ruin probability.

Quota-share reinsurance: Proportional reinsurance under which a fixed percentage of premiums and losses is ceded.

Surplus-share reinsurance: Proportional reinsurance based on retaining a fixed monetary line and ceding excess policy amounts.

Excess-of-loss reinsurance: Non-proportional cover responding to losses above a specified retention.

Catastrophe excess of loss: Cover responding to aggregate losses caused by a defined catastrophic event.

Stop-loss reinsurance: Protection against aggregate losses exceeding a specified attachment point.

Burning-cost pricing: Reinsurance pricing based primarily on adjusted historical ceded losses.

Model risk: The possibility that model assumptions, structure, data, or implementation produce materially misleading results.

Review Questions

Short Answer Questions

1. Distinguish frequency risk from severity risk.

2. State the expected value and variance formulas for a compound Poisson aggregate loss.

3. What is the principal difference between collective and individual risk models?

4. Define the adjustment coefficient in the classical risk model.

5. Explain why catastrophe losses cannot normally be treated as independent ordinary claims.

Long Answer Questions

1. Explain the construction of a collective risk model and discuss the assumptions underlying the compound Poisson model.

2. Compare convolution, recursive methods, approximation, and simulation for aggregate-loss calculation.

3. Derive and interpret the mean and variance of aggregate loss using conditional expectation and conditional variance.

4. Describe proportional, per-risk excess-of-loss, catastrophe excess-of-loss, and stop-loss reinsurance.

5. Discuss model risk and governance requirements in actuarial aggregate-loss and reinsurance modelling.

Analytical / Application Questions

1. An insurer has claim count N with mean 500 and variance 700. Individual claims have mean 8,000 and variance 144,000,000. Calculate the expected aggregate loss and its variance, assuming independence.

2. An excess-of-loss treaty has retention 100,000 and limit 400,000. Explain how the recovery is calculated for individual claims of 60,000, 250,000, and 700,000.

3. A portfolio is exposed to many moderate independent claims and a small probability of a severe regional storm affecting many policies. Recommend a reinsurance structure and justify it using frequency, severity, and dependence considerations.

22. Multiple Choice Questions with Answers

Question 1

Under the collective risk model, aggregate loss is represented as:

A. A fixed claim amount multiplied by a fixed claim count

B. A random sum of claim severities with a random claim count

C. The maximum loss among all policies only

D. Premium income less operating expenses

Correct Answer: B. A random sum of claim severities with a random claim count

Question 2

If N and X are independent, which expression gives the expected aggregate loss when S is the sum of N claims?

A. E(N)+E(X)

B. E(N)/E(X)

C. E(N)E(X)

D. Var(N)Var(X)

Correct Answer: C. E(N)E(X)

Question 3

For a compound Poisson model with frequency parameter lambda, the variance of aggregate loss is:

A. lambda Var(X)

B. lambda[E(X)]²

C. lambda E(X²)

D. E(X²)/lambda

Correct Answer: C. lambda E(X²)

Question 4

Which method is generally most flexible for incorporating complex dependence, policy limits, and reinsurance terms?

A. Normal approximation

B. Monte Carlo simulation

C. Simple variance calculation

D. A fixed-frequency model

Correct Answer: B. Monte Carlo simulation

Question 5

A normal approximation may be unsuitable for an insurance aggregate loss primarily because the distribution is often:

A. Symmetric and bounded below by zero

B. Right-skewed and potentially heavy-tailed

C. Always uniform

D. Independent of claim severity

Correct Answer: B. Right-skewed and potentially heavy-tailed

Question 6

In the classical risk model, ruin occurs when:

A. Premium income exceeds expected claims

B. The insurer’s surplus becomes negative

C. The number of claims is exactly zero

D. The reinsurer earns a profit

Correct Answer: B. The insurer’s surplus becomes negative

Question 7

Under quota-share reinsurance with ceded proportion q, the insurer generally retains:

A. q of every covered loss

B. 1/q of every covered loss

C. (1-q) of every covered loss

D. Only losses exceeding q

Correct Answer: C. (1-q) of every covered loss

Question 8

A per-risk excess-of-loss treaty with retention d pays, before any limit:

A. min(X,d)

B. max(X-d,0)

C. d-X

D. qX for a fixed q

Correct Answer: B. max(X-d,0)

Question 9

Catastrophe excess-of-loss reinsurance is particularly designed to address:

A. Only ordinary small claims

B. Losses from one event affecting multiple risks

C. Administrative expenses

D. Every claim in exact proportion to premium

Correct Answer: B. Losses from one event affecting multiple risks

Question 10

If an insurer’s portfolio has many correlated claims from one storm, which modelling assumption is most important to reconsider?

A. That claims are independent

B. That premiums are received annually

C. That currency has a numerical value

D. That a policy has an insured party

Correct Answer: A. That claims are independent

Actuarial Principles and Insurance Mathematics
Chapter 8 of 15