Lee–Carter model

The Lee–Carter model is a numerical algorithm used in

ordered monotonically by time, usually with ages in columns and years in rows. The output is a forecasted matrix of mortality rates in the same format as the input.

The model uses singular value decomposition (SVD) to find:

• A univariate time series vector ${\displaystyle \mathbf {k} _{t}}$ that captures 80–90% of the mortality trend (here the subscript ${\displaystyle t}$ refers to time),
• A vector ${\displaystyle \mathbf {b} _{x}}$ that describes the relative mortality at each age (here the subscript ${\displaystyle x}$ refers to age), and
• A scaling constant (referred to here as ${\displaystyle s_{1}}$ but unnamed in the literature).

Surprisingly, ${\displaystyle \mathbf {k} _{t}}$ is usually linear, implying that gains to life expectancy are fairly constant year after year in most populations. Prior to computing SVD, age specific mortality rates are first transformed into ${\displaystyle \mathbf {A} _{x,t}}$, by taking their

them by subtracting their age-specific means over time. The age-specific mean over time is denoted by ${\displaystyle \mathbf {a} _{x}}$. The subscript ${\displaystyle x,t}$ refers to the fact that ${\displaystyle \mathbf {A} _{x,t}}$ spans both age and time.

Many researchers adjust the ${\displaystyle \mathbf {k} _{t}}$ vector by fitting it to empirical life expectancies for each year, using the ${\displaystyle \mathbf {a} _{x}}$ and ${\displaystyle \mathbf {b} _{x}}$ generated with SVD. When adjusted using this approach, changes to ${\displaystyle \mathbf {k} _{t}}$ are usually small.

To forecast mortality, ${\displaystyle \mathbf {k} _{t}}$ (either adjusted or not) is projected into ${\displaystyle n}$ future years using an

model. The corresponding forecasted ${\displaystyle \mathbf {A} _{x,t+n}}$ is recovered by multiplying ${\displaystyle \mathbf {k} _{t+n}}$ by ${\displaystyle \mathbf {b} _{x}}$ and the first diagonal element of S (when ${\displaystyle \mathbf {U} \mathbf {S} \mathbf {V^{*}} ={\text{svd}}(\mathbf {A} _{x,t})}$). The actual mortality rates are recovered by taking exponentials of this vector.

Because of the linearity of ${\displaystyle \mathbf {k} _{t}}$, it is generally modeled as a random walk with trend. Life expectancy and other life table measures can be calculated from this forecasted matrix after adding back the means and taking exponentials to yield regular mortality rates.

In most implementations,

. A band of mortality between 5% and 95% percentiles of the simulated results is considered to be a valid forecast. These simulations are done by extending ${\displaystyle \mathbf {k} _{t}}$ into the future using randomization based on the standard error of ${\displaystyle \mathbf {k} _{t}}$ derived from the input data.

Algorithm

The algorithm seeks to find the least squares solution to the equation:

${\displaystyle \ln {(\mathbf {m} _{x,t})}=\mathbf {a} _{x}+\mathbf {k} _{t}\mathbf {b} _{x}+\epsilon _{x,t}}$

where ${\displaystyle \mathbf {m} _{x,t}}$ is a matrix of mortality rate for each age ${\displaystyle x}$ in each year ${\displaystyle t}$.

1. Compute ${\displaystyle \mathbf {a} _{x}}$ which is the average over time of ${\displaystyle \ln {(\mathbf {m} _{x,t})}}$ for each age:

   ${\displaystyle \mathbf {a} _{x}={\frac {\sum _{t=1}^{T}{\ln {(\mathbf {m} _{x,t})}}}{T}}}$

2. Compute ${\displaystyle \mathbf {A} _{x,t}}$ which will be used in SVD:

   ${\displaystyle \mathbf {A} _{x,t}=\ln {(\mathbf {m} _{x,t})}-\mathbf {a} _{x}}$

3. Compute the singular value decomposition of ${\displaystyle \mathbf {A} _{x,t}}$:

   ${\displaystyle \mathbf {U} \mathbf {S} \mathbf {V^{*}} ={\text{svd}}(\mathbf {A} _{x,t})}$

4. Derive ${\displaystyle \mathbf {k} _{t}}$, ${\displaystyle s_{1}}$ (the scaling eigenvalue), and ${\displaystyle \mathbf {b} _{x}}$ from ${\displaystyle \mathbf {U} }$, ${\displaystyle \mathbf {S} }$, and ${\displaystyle \mathbf {V^{*}} }$:

   ${\displaystyle \mathbf {k} _{t}=(u_{1,1},u_{2,1},...,u_{t,1})}$

   ${\displaystyle \mathbf {b} _{x}=(v_{1,1},v_{1,2},...,v_{1,x})}$

5. Forecast ${\displaystyle \mathbf {k} _{t}}$ using a standard univariate
   ARIMA
   model to ${\displaystyle n}$ additional years:

   ${\displaystyle \mathbf {k} _{t+n}={\text{ARIMA}}(\mathbf {k} _{t},n)}$

6. Use the forecasted ${\displaystyle \mathbf {k} _{t+n}}$, with the original ${\displaystyle \mathbf {b} _{x}}$, and ${\displaystyle \mathbf {a} _{x}}$ to calculate the forecasted mortality rate for each age:

   ${\displaystyle \mathbf {m} _{x,t+n}=\exp(\mathbf {a} _{x}+s_{1}\mathbf {k} _{t+n}\mathbf {b} _{x})}$

Discussion

Without applying SVD or some other method of

algorithm.

The Lee–Carter model was introduced by

There have been extensions to the Lee–Carter model, most notably to account for missing years, correlated male and female populations, and large scale coherency in populations that share a mortality regime (western Europe, for example). Many related papers can be found on Professor Ronald Lee's website.

Implementations

There are surprisingly few software packages for forecasting with the Lee–Carter model.

• LCFIT is a web-based package with interactive forms.
• Professor Rob J. Hyndman provides an R package for demography that includes routines for creating and forecasting a Lee–Carter model.
• Alternatives in R include the StMoMo package of Villegas, Millossovich and Kaishev (2015).
• Professor German Rodriguez provides code for the Lee–Carter Model using Stata.
• Using
  Matlab, Professor Eric Jondeau and Professor Michael Rockinger have put together the Longevity Toolbox
  for parameter estimation.

References