A comparison of the Lee-Carter model and AR-ARCH model for forecasting mortality rates

With the decline in the mortality level of populations, national social security systems and insurance companies of most developed countries are reconsidering their mortality tables taking into account the longevity risk. The Lee and Carter model is the first discrete-time stochastic model to consider the increased life expectancy trends in mortality rates and is still broadly used today. In this paper, we propose an alternative to the Lee-Carter model: an AR(1)-ARCH(1) model. More specifically, we compare the performance of these two models with respect to forecasting age-specific mortality in Italy. We fit the two models, with Gaussian and t-student innovations, for the matrix of Italian death rates from 1960 to 2003. We compare the forecast ability of the two approaches in out-of-sample analysis for the period 2004-2006 and find that the AR(1)-ARCH(1) model with t-student innovations provides the best fit among the models studied in this paper.     

Is the Home Equity Conversion Mortgage in the United States sustainable? Evidence from pricing mortgage insurance premiums and non-recourse provisions using the conditional Esscher transform

The purpose of this paper is to build a modeling and pricing framework to investigate the sustainability of the Home Equity Conversion Mortgage (HECM) program in the United States under realistic economic scenarios, i.e., whether the premium payments cover the fair premiums for the inherent risks in the HECM program. We note that earlier HECM models use static mortality tables, neglecting the dynamics of mortality rates and extreme mortality jumps. The earlier models also assume housing prices follow a geometric Brownian motion, which contradicts the fact that housing prices exhibit strong autocorrelation and varying volatility over time. To solve these problems, we propose a generalized Lee-Carter model with asymmetric jump effects to fit the mortality data, and model the house price index via an ARIMA-GARCH process. We then employ the conditional Esscher transform to price the non-recourse provision of reverse mortgages and compare it with the calculated mortgage insurance premiums. The HECM program turns out to be sustainable based on our model setup and parameter settings.     

Modeling mortality with a Bayesian vector autoregression

Parametric mortality models capture the cross section of mortality rates. These models fit the older ages better, because of the more complex cross section of mortality at younger and middle ages. Dynamic parametric mortality models fit a time series to the parameters, such as a Vector-auto-regression (VAR), in order to capture trends and uncertainty in mortality improvements. We consider the full age range using the Heligman and Pollard (1980) model, a cross-sectional mortality model with parameters that capture specific features of different age ranges. We make the Heligman–Pollard model dynamic using a Bayesian Vector Autoregressive (BVAR) model for the parameters and compare with more commonly used VAR models. We fit the models using Australian data, a country with similar mortality experience to many developed countries. We show how the Bayesian Vector Autoregressive (BVAR) models improve forecast accuracy compared to VAR models and quantify parameter risk which is shown to be significant.     

Mortality density forecasts: An analysis of six stochastic mortality models

This paper develops a framework for developing forecasts of future mortality rates. We discuss the suitability of six stochastic mortality models for forecasting future mortality and estimating the density of mortality rates at different ages. In particular, the models are assessed individually with reference to the following qualitative criteria that focus on the plausibility of their forecasts: biological reasonableness; the plausibility of predicted levels of uncertainty in forecasts at different ages; and the robustness of the forecasts relative to the sample period used to fit the model. An important, though unsurprising, conclusion is that a good fit to historical data does not guarantee sensible forecasts. We also discuss the issue of model risk, common to many modelling situations in demography and elsewhere. We find that even for those models satisfying our qualitative criteria, there are significant differences among central forecasts of mortality rates at different ages and among the distributions surrounding those central forecasts.     

Forecasting mortality rate improvements with a high-dimensional VAR

Forecasting mortality rates is a problem which involves the analysis of high-dimensional time series. Most of usual mortality models propose to decompose the mortality rates into several latent factors to reduce this complexity. These approaches, in particular those using cohort factors, have a good fit, but they are less reliable for forecasting purposes. One of the major challenges is to determine the spatial–temporal dependence structure between mortality rates given a relatively moderate sample size. This paper proposes a large vector autoregressive (VAR) model fitted on the differences in the log-mortality rates, ensuring the existence of long-run relationships between mortality rate improvements. Our contribution is threefold. First, sparsity, when fitting the model, is ensured by using high-dimensional variable selection techniques without imposing arbitrary constraints on the dependence structure. The main interest is that the structure of the model is directly driven by the data, in contrast to the main factor-based mortality forecasting models. Hence, this approach is more versatile and would provide good forecasting performance for any considered population. Additionally, our estimation allows a one-step procedure, as we do not need to estimate hyper-parameters. The variance–covariance matrix of residuals is then estimated through a parametric form. Secondly, our approach can be used to detect nonintuitive age dependence in the data, beyond the cohort and the period effects which are implicitly captured by our model. Third, our approach can be extended to model the several populations in long run perspectives, without raising issue in the estimation process. Finally, in an out-of-sample forecasting study for mortality rates, we obtain rather good performances and more relevant forecasts compared to classical mortality models using the French, US and UK data. We also show that our results enlighten the so-called cohort and period effects for these populations.     

The age pattern of transitory mortality jumps and its impact on the pricing of catastrophic mortality bonds

To value catastrophic mortality bonds, a number of stochastic mortality models with transitory jump effects have been proposed. Rather than modeling the age pattern of jump effects explicitly, most of the existing models assume that the distributions of jump effects and general mortality improvements across ages are identical. Nevertheless, this assumption does not seem to be in line with what we observe from historical data. In this paper, we address this problem by introducing a Lee–Carter variant that captures the age pattern of mortality jumps by a distinct collection of parameters. The model variant is then further generalized to permit the age pattern of jump effects to vary randomly. We illustrate the two proposed models with mortality data from the United States and English and Welsh populations, and use them to value hypothetical mortality bonds with similar specifications to the Atlas IX Capital Class B note that was launched in 2013. It is found that the features we consider have a significant impact on the estimated prices.     

Modelling life tables with advanced ages: An extreme value theory approach

We propose a new model – we call it a smoothed threshold life table (STLT) model – to generate life tables incorporating information on advanced ages. Our method allows a smooth mortality transition from non-extreme to extreme ages, and provides objectively determined highest attained ages with which to close the life table.We proceed by modifying the threshold life table (TLT) model developed by Li et al. (2008). In the TLT model, extreme value theory (EVT) is used to make optimal use of the relatively small number of observations at high ages, while the traditional Gompertz distribution is assumed for earlier ages. Our novel contribution is to constrain the hazard function of the two-part lifetime distribution to be continuous at the changeover point between the Gompertz and EVT models. This simple but far-reaching modification not only guarantees a smooth transition from non-extreme to extreme ages, but also provides a better and more robust fit than the TLT model when applied to a high quality Netherlands dataset. We show that the STLT model also compares favourably with other existing methods, including the Gompertz–Makeham model, logistic models, Heligman–Pollard model and Coale–Kisker method, and that a further generalisation, a time-dependent dynamic smooth threshold life table (DSTLT) model, generally has superior in-sample fitting as well as better out-of-sample forecasting performance, compared, for example, with the Cairns et al. (2006) model.     

Evaluating the goodness of fit of stochastic mortality models

This study sets out a framework to evaluate the goodness of fit of stochastic mortality models and applies it to six different models estimated using English & Welsh male mortality data over ages 64-89 and years 1961-2007. The methodology exploits the structure of each model to obtain various residual series that are predicted to be iid standard normal under the null hypothesis of model adequacy. Goodness of fit can then be assessed using conventional tests of the predictions of iid standard normality. The models considered are: Lee and Carter’s (1992) one-factor model, a version of Renshaw and Haberman’s (2006) extension of the Lee-Carter model to allow for a cohort-effect, the age-period-cohort model, which is a simplified version of the Renshaw-Haberman model, the 2006 Cairns-Blake-Dowd two-factor model and two generalized versions of the latter that allow for a cohort-effect. For the data set considered, there are some notable differences amongst the different models, but none of the models performs well in all tests and no model clearly dominates the others.     

Modelling stochastic mortality for dependent lives

Stochastic mortality, i.e. modelling death arrival via a jump process with stochastic intensity, is gaining an increasing reputation as a way to represent mortality risk. This paper is a first attempt to model the mortality risk of couples of individuals, according to the stochastic intensity approach. Dependence between the survival times of the members of a couple is captured by an Archimedean copula.We also provide a methodology for fitting the joint survival function by working separately on the (analytical) marginals and on the (analytical) copula. First, we provide a sample-based calibration for the intensity, using a time-homogeneous, non mean-reverting, affine process: this gives the marginal survival functions. Then we calibrate and select the best fit copula according to the Wang and Wells [Wang, W., Wells, M.T., 2000b. Model selection and semiparametric inference for bivariate failure-time data. J. Amer. Statis. Assoc. 95, 62-72] methodology for censored data. By coupling the calibrated marginals with the best fit copula, we obtain a joint survival function, which incorporates the stochastic nature of mortality improvements.We apply the methodology to a well known insurance data set, using a sample generation. The best fit copula turns out to be one listed in [Nelsen, R.B., 2006. An Introduction to Copulas, Second ed. In: Springer Series], which implies not only positive dependence, but dependence increasing with age.     

A semiparametric panel approach to mortality modeling

During the past twenty years, there has been a rapid growth in life expectancy and an increased attention on funding for old age. Attempts to forecast improving life expectancy have been boosted by the development of stochastic mortality modeling, for example the Cairns–Blake–Dowd (CBD) 2006 model. The most common optimization method for these models is maximum likelihood estimation (MLE) which relies on the assumption that the number of deaths follows a Poisson distribution. However, several recent studies have found that the true underlying distribution of death data is overdispersed in nature (see Cairns et al. 2009 and Dowd et al. 2010). Semiparametric models have been applied to many areas in economics but there are very few applications of such models in mortality modeling. In this paper we propose a local linear panel fitting methodology to the CBD model which would free the Poisson assumption on number of deaths. The parameters in the CBD model will be considered as smooth functions of time instead of being treated as a bivariate random walk with drift process in the current literature. Using the mortality data of several developed countries, we find that the proposed estimation methods provide comparable fitting results with the MLE method but without the need of additional assumptions on number of deaths. Further, the 5-year-ahead forecasting results show that our method significantly improves the accuracy of the forecast.     

A continuous-time stochastic model for the mortality surface of multiple populations

We formulate, study and calibrate a continuous-time model for the joint evolution of the mortality surface of multiple populations. We model the mortality intensity by age and population as a mixture of stochastic latent factors, that can be either population-specific or common to all populations. These factors are described by affine time-(in)homogeneous stochastic processes. Traditional, deterministic mortality laws can be extended to multi-population stochastic counterparts within our framework. We detail the calibration procedure when factors are Gaussian, using centralized data-fusion Kalman filter. We provide an application based on the joint mortality of UK and Dutch males and females. Although parsimonious, the specification we calibrate provides a good fit of the observed mortality surface (ages 0–89) of both sexes and populations between 1960 and 2013.     

Modeling multi-country mortality dependence and its application in pricing survivor index swaps—A dynamic copula approach

This paper introduces mortality dependence in multi-country mortality modeling using a dynamic copula approach. Specifically, we use time-varying copula models to capture the mortality dependence structure across countries, examining both symmetric and asymmetric dependence structures. In addition, to capture the phenomenon of a heavy tail for the multi-country mortality index, we consider not only the setting of Gaussian innovations but also non-Gaussian innovations under the Lee–Carter framework model. As tests of the goodness of fit of different dynamic copula models, the pattern of mortality dependence, and the distribution of the innovations, we used empirical mortality data from Finland, France, the Netherlands, and Sweden. To understand the effect of mortality dependence on longevity derivatives, we also built a valuation framework for pricing a survivor index swap, then investigated the fair swap rates of a survivor swap numerically. We demonstrate that failing to consider the dynamic copula mortality model and non-Gaussian innovations would lead to serious underestimations of the swap rates and loss reserves.     

A common age effect model for the mortality of multiple populations

We introduce a model for the mortality rates of multiple populations. To build the proposed model we investigate to what extent a common age effect can be found among the mortality experiences of several countries and use a common principal component analysis to estimate a common age effect in an age–period model for multiple populations. The fit of the proposed model is then compared to age–period models fitted to each country individually, and to the fit of the model proposed by Li and Lee (2005).Although we do not consider stochastic mortality projections in this paper, we argue that the proposed common age effect model can be extended to a stochastic mortality model for multiple populations, which allows to generate mortality scenarios simultaneously for all considered populations. This is particularly relevant when mortality derivatives are used to hedge the longevity risk in an annuity portfolio as this often means that the underlying population for the derivatives is not the same as the population in the annuity portfolio.     

On stochastic mortality modeling

In the last decennium a vast literature on stochastic mortality models has been developed. All well-known models have nice features but also disadvantages. In this paper a stochastic mortality model is proposed that aims at combining the nice features from the existing models, while eliminating the disadvantages. More specifically, the model fits historical data very well, is applicable to a full age range, captures the cohort effect, has a non-trivial (but not too complex) correlation structure and has no robustness problems, while the structure of the model remains relatively simple. Also, the paper describes how to incorporate parameter uncertainty in the model. Furthermore, a risk neutral version of the model is given, that can be used for pricing.     

Indifference pricing of pure endowments and life annuities under stochastic hazard and interest rates

We study indifference pricing of mortality contingent claims in a fully stochastic model. We assume both stochastic interest rates and stochastic hazard rates governing the population mortality. In this setting we compute the indifference price charged by an insurer that uses exponential utility and sells k contingent claims to k independent but homogeneous individuals. Throughout we focus on the examples of pure endowments and temporary life annuities. We begin with a continuous-time model where we derive the linear pdes satisfied by the indifference prices and carry out extensive comparative statics. In particular, we show that the price-per-risk grows as more contracts are sold. We then also provide a more flexible discrete-time analog that permits general hazard rate dynamics. In the latter case we construct a simulation-based algorithm for pricing general mortality-contingent claims and illustrate with a numerical example.     

Longevity risk management for life and variable annuities: The effectiveness of static hedging using longevity bonds and derivatives

For many years, the longevity risk of individuals has been underestimated, as survival probabilities have improved across the developed world. The uncertainty and volatility of future longevity has posed significant risk issues for both individuals and product providers of annuities and pensions. This paper investigates the effectiveness of static hedging strategies for longevity risk management using longevity bonds and derivatives (q-forwards) for the retail products: life annuity, deferred life annuity, indexed life annuity, and variable annuity with guaranteed lifetime benefits. Improved market and mortality models are developed for the underlying risks in annuities. The market model is a regime-switching vector error correction model for GDP, inflation, interest rates, and share prices. The mortality model is a discrete-time logit model for mortality rates with age dependence. Models were estimated using Australian data. The basis risk between annuitant portfolios and population mortality was based on UK experience. Results show that static hedging using q-forwards or longevity bonds reduces the longevity risk substantially for life annuities, but significantly less for deferred annuities. For inflation-indexed annuities, static hedging of longevity is less effective because of the inflation risk. Variable annuities provide limited longevity protection compared to life annuities and indexed annuities, and as a result longevity risk hedging adds little value for these products.     

Classification Trees With Bivariate Linear Discriminant Node Models

This article introduces a classification tree algorithm that can simultaneously reduce tree size, improve class prediction, and enhance data visualization. We accomplish this by fitting a bivariate linear discriminant model to the data in each node. Standard algorithms can produce fairly large tree structures because they employ a very simple node model, wherein the entire partition associated with a node is assigned to one class. We reduce the size of our trees by letting the discriminant models share part of the data complexity. Being themselves classifiers, the discriminant models can also help to improve prediction accuracy. Finally, because the discriminant models use only two predictor variables at a time, their effects are easily visualized by means of two-dimensional plots. Our algorithm does not simply fit discriminant models to the terminal nodes of a pruned tree, as this does not reduce the size of the tree. Instead, discriminant modeling is carried out in all phases of tree growth and the misclassification costs of the node models are explicitly used to prune the tree. Our algorithm is also distinct from the “linear combination split” algorithms that partition the data space with arbitrarily oriented hyperplanes. We use axis-orthogonal splits to preserve the interpretability of the tree structures. An extensive empirical study with real datasets shows that, in general, our algorithm has better prediction power than many other tree or nontree algorithms.     

One-year Value-at-Risk for longevity and mortality

Upcoming new regulation on regulatory required solvency capital for insurers will be predominantly based on a one-year Value-at-Risk measure. This measure aims at covering the risk of the variation in the projection year as well as the risk of changes in the best estimate projection for future years. This paper addresses the issue how to determine this Value-at-Risk for longevity and mortality risk. Naturally, this requires stochastic mortality rates. In the past decennium, a vast literature on stochastic mortality models has been developed. However, very few of them are suitable for determining the one-year Value-at-Risk. This requires a model for mortality trends instead of mortality rates. Therefore, we will introduce a stochastic mortality trend model that fits this purpose. The model is transparent, easy to interpret and based on well known concepts in stochastic mortality modeling. Additionally, we introduce an approximation method based on duration and convexity concepts to apply the stochastic mortality rates to specific insurance portfolios.     

A comparative study of parametric mortality projection models

The relative merits of different parametric models for making life expectancy and annuity value predictions at both pensioner and adult ages are investigated. This study builds on current published research and considers recent model enhancements and the extent to which these enhancements address the deficiencies that have been identified of some of the models. The England & Wales male mortality experience is used to conduct detailed comparisons at pensioner ages, having first established a common basis for comparison across all models. The model comparison is then extended to include the England & Wales female experience and both the male and female USA mortality experiences over a wider age range, encompassing also the working ages.