---
title: Lee–Carter Mortality Forecasting
url: https://www.emergentmind.com/topics/lee-carter-framework
type: topic
---

# Lee–Carter Mortality Forecasting

The Lee–Carter framework is a foundational and widely adopted approach for modeling and forecasting mortality rates in demographic and actuarial contexts. It is based on the decomposition of the logarithm of age-specific mortality rates into fixed age effects and a time-varying period effect, with stochastic dynamics applied to capture systematic mortality improvement over time. The framework offers extensions for joint multi-population modeling, incorporation of cohort and volatility effects, and various statistical estimation strategies.

## 1. Mathematical Formulation of the Lee–Carter Model

The classical Lee–Carter (LC) model expresses the log central death rate at age $x$ and year $t$ as

$$
\log m_{x,t} = \alpha_x + \beta_x\,\kappa_t,
$$

where:
- $\alpha_x$ is the average log mortality rate at age $x$ across time,
- $\beta_x$ quantifies the sensitivity of age $x$ to the period factor (i.e., how improvements at age $x$ respond to changes in the period index),
- $\kappa_t$ is a time-varying mortality index reflecting period-specific mortality improvement.

For stochastic projection, $\kappa_t$ is assumed to follow a time series process, typically a random walk with drift:

$$
\kappa_t = \kappa_{t-1} + \theta + \omega_t,\quad \omega_t \sim N(0,\sigma^2)
$$

with identifiability constraints $\sum_x \beta_x = 1$ and $\sum_t \kappa_t = 0$. Extensions can incorporate cohort terms (e.g., Renshaw–Haberman: $\gamma_{t-x}$) and generalizations to Age–Period–Cohort (APC) models [1605.09484].

## 2. Estimation and Inference

Parameter estimation in the LC class is conducted via maximum likelihood, typically under a Poisson or negative binomial framework for observed death counts, using exposures $E_{x,t}$ as an offset:

$$
D_{x,t} \sim \text{Pois}(E_{x,t} m_{x,t}),
$$

or, for overdispersed settings,

$$
D_{x,t} \sim \text{NegBin}(E_{x,t} m_{x,t},\phi).
$$

For fully Bayesian variants and for the incorporation of smooth age, period, and cohort effects, this framework is extended using generalized additive models (GAM) with P-splines and Hamiltonian Monte Carlo for joint posterior sampling [1802.03242]. Model selection among competing cut-offs (e.g., age thresholds for model transition) is achieved using Pareto-smoothed importance sampling leave-one-out cross-validation (PSIS-LOOIC) and Bayesian stacking to form predictive mixtures of models.

In joint multi-population extensions, as in the Li–Lee model, the period index is decomposed into a "common" (global) trend and population-specific deviations, both modeled as stochastic processes subject to identifiability constraints. Cross-population period innovations may be modeled jointly as correlated random walks [2111.10164].

## 3. Extensions: Old-Age, Volatility, and State-Space Generalizations

The Lee–Carter framework struggles at extreme old ages, where data are sparser and the log-linear structure becomes inadequate. To address this, parametric “Kannisto-like” formulae for $x \geq x_\text{old}$ can be used, imposing a logistic asymptotic behavior:

$$
m_{x t} =
\frac{\exp(\beta_0^{\text{old}} + \beta_1^{\text{old}} x + \beta_2^{\text{old}} t + \beta_3^{\text{old}} x\,t)}{1 + \exp(\beta_0^{\text{old}} - \log \psi + \beta_1^{\text{old}} x + \beta_2^{\text{old}} t + \beta_3^{\text{old}} x\,t)}
\times \exp[s_\gamma(t-x) + \kappa_t],
$$

with $\psi$ governing the old-age mortality plateau. Multi-population models may share $\psi$ across populations to borrow statistical strength in old-age extrapolation [1802.03242].

For full dynamic modeling, the framework is often cast in a state-space or hidden Markov model setting, which enables:
- Unified likelihood or Bayesian inference,
- Flexible identification constraints suited to filtering,
- Inclusion of observation error heteroscedasticity,
- Stochastic volatility in period effects (via time-varying process innovations),
- Cohort effects with ARIMA dynamics [1605.09484].

Bayesian estimation is facilitated by Kalman filtering and smoothing, particle MCMC, or Hamiltonian Monte Carlo.

## 4. Forecasting and Uncertainty Quantification

Forecasts are generated by simulating the time series of the period index (and cohort effects, if included), then combining with age-$x$ coefficients to yield projected log-mortality, and converting to death probabilities:

$$
q_{x,t+h} \approx 1 - \exp(-m_{x,t+h}).
$$

Credible bands (e.g., 90%) are obtained by summarizing the empirical quantiles of projected mortality rates or derived quantities (e.g., period/cohort life expectancy) across posterior draws [1802.03242]. This inherently probabilistic approach enables comprehensive quantification of both parameter and process uncertainty, which is increasingly required in pension, insurance, and regulatory forecasting contexts.

Model averaging across competing transition-ages or structural choices is handled via Bayesian stacking, forming a weighted ensemble whose predictive distribution is:

$$
p(\hat y|y) = \sum_k w_k p(\hat y|y, M_k),
$$

with $w$ selected to maximize expected log predictive density.

## 5. Performance, Empirical Results, and Limitations

Empirical validation on UK data (1961–2013) demonstrates that stacking across models with different old-age cut-offs provides superior robust forecasts, with credible bands that largely conform to observed log-mortality in out-of-sample tests (88–92% of observations covered at 90% level) [1802.03242]. Some undercoverage can occur in rapidly changing subpopulations or ages where recent trend acceleration exceeds stochastic volatility captured in the model.

Joint-sex fitting with shared old-age asymptote and highly correlated period innovations substantially increases sampling efficiency and produces biologically plausible co-movement in male and female projections.

Life expectancy intervals under the Bayesian framework are wider than scenario-based official bands on short horizons and become narrower beyond 30 years, reflecting the proper propagation of parameter uncertainty.

Key limitations of the basic Lee–Carter model include:
- Inadequacy at very young and very old ages without explicit model extensions,
- Inability to accommodate one-off shocks or regime changes (e.g., pandemic years) except by scenario analysis,
- Implicit assumption of log-linearity, which may not always hold after structural breaks,
- Sensitivity to the handling of cohort effects or heteroscedasticity, addressed in state-space and Bayesian extended frameworks.

## 6. Practical and Regulatory Impact

The Lee–Carter class and its Bayesian GAM or state-space generalizations are core tools for life-insurance, social-security, and longevity-risk management. Their flexibility in handling joint-population inference, old-age smoothing, uncertainty quantification, and direct model averaging make them especially suited for modern stochastic projections required by Solvency II and similar regulatory regimes.

Key features enabling regulatory adoption include:
- Full probabilistic forecasting with credible bands,
- Direct scenario analysis for stress testing,
- Modular implementation in computational frameworks (e.g., Stan, R, Python) exploiting fast HMC or Kalman filtering,
- Empirical validation against holdout data and coverage properties.

Models in this class offer transparent, technically rigorous, and extensible methods for long-range mortality risk assessment, supporting actuarial pricing, reserving, and policy setting at scale [1802.03242, 1605.09484].

Source: https://www.emergentmind.com/topics/lee-carter-framework