Next Article in Journal
Bounds for the Ruin Probability in the Sparre–Andersen Model
Next Article in Special Issue
A Quantitative Comparison of Mortality Models with Jumps: Pre- and Post-COVID Insights on Insurance Pricing
Previous Article in Journal
Model for Technology Risk Assessment in Commercial Banks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LSTM-Based Coherent Mortality Forecasting for Developing Countries

1
Department of Mathematics and Statistics, Concordia University, Montreal, QC H3G 1M8, Canada
2
Department of Financial and Actuarial Mathematics, Xi’an Jiaotong–Liverpool University, Suzhou 215123, China
*
Author to whom correspondence should be addressed.
Risks 2024, 12(2), 27; https://doi.org/10.3390/risks12020027
Submission received: 13 December 2023 / Revised: 19 January 2024 / Accepted: 26 January 2024 / Published: 1 February 2024
(This article belongs to the Special Issue Extreme Events: Mortality Modelling and Insurance)

Abstract

:
This paper studies a long short-term memory (LSTM)-based coherent mortality forecasting method for developing countries or regions. Many of such developing countries have experienced a rapid mortality decline over the past few decades. However, their recent mortality development trend is not necessarily driven by the same factors as their long-term behavior. Hence, we propose a time-varying mortality forecasting model based on the life expectancy and lifespan disparity gap between these developing countries and a selected benchmark group. Here, the mortality improvement trend for developing countries is expected to converge gradually to that of the benchmark group during the projection phase. More specifically, we use a unified deep neural network model with LSTM architecture to project the life expectancy and lifespan disparity difference, which further controls the rotation of the time-varying weight parameters in the model. This approach is applied to three developing countries and three developing regions. The empirical results show that this LSTM-based coherent forecasting method outperforms classical methods, especially for the long-term projections of mortality rates in developing countries.

1. Introduction

In the last few decades, human mortality has improved significantly in several countries, especially in developing regions. These mortality reductions can generate important longevity risks for life insurance companies and pension schemes. The study of these longevity improvements is fundamental in life insurance and annuity research in actuarial science literature.
There are statistical techniques for the forecasting of future mortality. A popular method, proposed by Lee and Carter (1992), is the so-called Lee–Carter (LC) model, where the log force of mortality ln ( m x , t ) is represented as the sum of an age component a x plus the product of an age-specific function b x and a time component k t . Obviously, such a model cannot be fitted as a regular regression model because of the product of parameter terms. In work by Lee and Carter (1992), a two-step method was applied, where singular value decomposition (SVD) is used to fit the model, followed by an autoregressive component, or random walk, to forecast the time component k t .
In the last two decades, the LC method has been used in the practice frequently, and several papers propose various extensions of the method; interested readers can refer to Lee (2000), Pitacco (2004), Wong-Fupuy and Haberman (2004), as well as the Cairns–Blake–Dowd (CBD) model (Cairns et al. 2006), and references therein. This literature has shown how the original LC method can lack flexibility with regard to the effect of age. Renshaw and Haberman (2003) extended the LC method to a multi-factor version adding an age-specific enhancement. In addition, a single-factor model with cohort effects is proposed by Renshaw and Haberman (2006). For a comprehensive review of the early literature on these various forecasting methods, refer to Cairns et al. (2008) and Cairns et al. (2011a).
Note that apart from the LC method and its extensions, another approach developed in the literature for mortality modeling is based on generalized linear models (GLMs); for example, see Brouhns et al. (2002), Renshaw and Haberman (2006), and O’Hare and Li (2012). For a comprehensive survey on fitting GLMs to mortality data, refer to Currie (2016). In addition, the Bayesian approach appears in the literature on mortality modeling. For example, Czado et al. (2005) and Pedroza (2006) extend the LC model to Bayesian analyses using Markov chain Monte Carlo (MCMC) methods. Cairns et al. (2011b) further extended Bayesian stochastic mortality modeling to two populations; see Antonio et al. (2015) for an application of a Bayesian method under multiple populations. For more recent studies on mortality modeling with a Bayesian approach, refer to Li and Lu (2018), Li et al. (2019), and Wong et al. (2023), and references therein.
The LC method, and many of its early extensions, focuses on a single population (for example, modeling only one gender at a time or combined genders, and only one country). In particular, let m x , t denote the mortality rate of age x at time t for t = 0 , 1 , 2 , T and x = 0 , 1 , 2 ω for a given population. Then, let a x denote the mean log mortality rate over time, that is, a x = 1 T + 1 t = 0 T ln m x , t . Using the LC method, one obtains
ln m x , t = a x + b x k t + ϵ x , t ,
where ϵ x , t is the mean zero random noise, and k t is modeled using random walk with drift:
k t = d + k t 1 + δ t , δ t N ( 0 , σ 2 ) , E ( δ t δ s ) = 0 , for t s .
The literature shows that it is difficult for the LC method to forecast mortality rates for two genders at the same time in one population or in multiple populations and regions, where a certain divergence could be reached due to the differences in b x and d (and in turn, different k t ) in the model. For example, Carter and Lee (1992) suggested using the same k t , but gender-specific b x values to forecast male and female mortality rates separately in the U.S. Lee and Nault (1993) used the same k t and b x for mortality forecasting in each province of Canada, which works only when the b x values of different provinces, as obtained from historical data, are similar.
However, some early discussions on the convergence of life expectancy around the world (see, e.g., White 2002, Wilmoth 1998, and Vaupel and Schnabel 2004) show that there is convergence in long-term life expectancy, diverging forecasts of mortality rates for different populations in a group of countries is unrealistic. Therefore, Li and Lee (2005) introduced the so-called coherent extension of the LC method for mortality forecasting of a group of populations (we call it the Li–Lee method), where the log mortality rates for each member in the group are decomposed into three parts, namely, member-specific a i , x , common age and period effects B x and K t , and member-specific age and period effects, b i , x and k i , t . More precisely, the log mortality rates ln ( m i , x , t ) , for member i in the group at age x and time t, can be expressed as
ln m i , x , t = a i , x + B x K t + b i , x k i , t + ε i , x , t ,
with
K t = d 0 + K t 1 + ν t , k i , t = α 0 , i + α 1 , i k i , t 1 + ϵ i , t ,
where a i , x measures the average mortality level at age x in country i. K t is the common period effect for all countries and is modeled by a random walk with drift d 0 . B x is the common age effect, i.e., the common mortality sensitivity at age x, with respect to K t . In addition, k i , t and b i , x are the country-specific period and age effects, respectively, which measure the fluctuations around the common mortality patterns in the group for country i. Finally, ε i , x , t , ν t , and ϵ i , t are normally distributed i.i.d. errors. It turns out that the additional information provided by similar members/countries in the group can improve the forecast accuracy for individual countries.
On the other hand, according to Hanewald (2011), there is a strong long-term connection between the mortality dynamics and the gross domestic product (GDP) per capita and unemployment rate in a country, which points to the essential difference of the mortality improvements between developed countries and developing countries during the same period of time. Based on such an observation, Niu and Melenberg (2014) improved the LC model with an extra factor (namely, GDP) describing the economic growth; see Boonen and Li (2017) for the study under multiple populations. More recently, Ma and Boonen (2023) further argued that the consumer price index (CPI) is a more suitable factor, added to the LC model, to explain the mortality trends in a country. It explains the affordability of healthcare, food, and housing in that country. The above-mentioned studies verify, from a different point of view, that life expectancy (or any other similar indexes) can play an important role in mortality forecasting.
Moreover, as mentioned by Li et al. (2013), mortality decline decelerates in younger ages and accelerates at older ages in many developed countries. Such a “rotation" can generate problems in the results of long-term mortality projections using the LC method for developing countries that do not exhibit such a subtle rotation in their historical data, e.g., the projected mortality rates are low for younger ages. Hence, Li et al. (2013) developed a rotation-based LC method, where the out-of-sample b x was assumed to be converging to an ultimate structure based on the development of life expectancy. But, as argued by Li and Lu (2017), mortality rates should change smoothly and continuously across ages; such a problem, known as age–coherent mortality forecasting, is present in the above-mentioned Li–Lee model. Thus, Gao and Shi (2021) proposed two alternative extensions to the ordinary LC method: the LC-Geometric and LC-Hyperbolic models. The goal was to achieve long-term age coherence in mortality forecasts while retaining the short-term rotation-type forecasting adopted by Li et al. (2013). Here, Geometric and Hyperbolic refer to the type of decay allowed in the autoregressive (AR) model.
Recently, borrowing the concept of “rotation" proposed by Li et al. (2013), Li et al. (2021) developed a so-called rotation algorithm for the coherent mortality projections of less developed countries, which included all regions of Africa, Asia (except Japan), Latin America, the Caribbean, Melanesia, Micronesia, and Polynesia. Using Li–Lee’s model, where “rotation" refers to the effects of age and time components in the projection phase for developing countries; rotation may occur in b i , x and d i , based on their own historical data and the corresponding data from a group of developed countries used as the benchmark group. In their model, the rotation algorithm is controlled by a life expectancy gap function, between the target developing country and the benchmark group, where the gap function can be fitted by (double) logistic functions with some selected threshold levels for convergence of the gap.
Therefore, reliable mortality projection methods, especially the long-term projections of age-specific mortality rates, are crucial for developing countries. However, the recent fast decline in aggregate mortality might not be a long-term behavior (according to Müller and Krawinkel (2005), Austin and McKinney (2012), and Jeuland et al. (2013), the main factors contributing to recent mortality improvements in developing countries (especially for infants, the young, and the working-age population) are modernization, improved healthcare coverage, better nutrition, and prevention of infectious diseases, which can obviously only last for a short period of time); in the long run, the mortality patterns of developing countries could more closely resemble those of more developed countries (see, e.g., Li and Lee 2005). Therefore, predicting long-term, age-specific mortality rates in a developing country by simply extrapolating its historical patterns may lead to implausible results.
Hence, a method with the aforementioned “rotation" helps find a balance between the historical mortality pattern of the developing country and the average mortality patterns of a group of developed countries (the benchmark group). Note that for different target countries, the method proposed by Li et al. (2021) needs to use expert judgments and fit different gap functions with possible different gap thresholds. This can bring restrictions to the unified application of the method.
On the other hand, as mentioned by Aburto et al. (2020), populations with the same life expectancy level may experience substantial differences in the time of death. This indicates that the mortality pattern of a developing country can still be different from the benchmark group even if there is a convergence in the life expectancy gap. Hence, the lifespan disparity, which describes life expectancy lost due to death by an individual at different ages and times (see, e.g., Vaupel and Romo (2003) and Zhang and Vaupel (2009)), may provide additional information when examining the convergence of mortality development between the developing country and the benchmark group. Therefore, to continue the study of coherent mortality forecasting (especially long-term projections based on lifespan disparity) for developing countries, here we propose a unified coherent mortality forecasting method with time-dependent rotation weights based on a benchmark group. A deep neural network, in particular a long short-term memory (LSTM), is used for the projection of the life expectancy and lifespan disparity gaps between the target developing country and the corresponding benchmark group. The projected gaps are used in the control of the rotated time-varying weight parameters in the model during the projection phase for mortality forecasting of the developing country.
In the last decade, neural networks, especially deep neural networks (DNNs), have gained attention in the human mortality modeling and forecasting literature. The purpose of building neural network models for human mortality is to extend the modeling and forecasting ability passed by classical parametric models, such as the LC method and its many extensions. For example, Hainaut (2018) proposed a type of neural network analyzer, which uses an encoding and decoding network structure in the approximation of the non-linearity among ages, for each year in the data, for a single country. The model is essentially a feedforward neural network extension of the LC method, where a simple, fully connected feedforward neural network was used to learn the common nonlinearities in the lower dimensional structure of the log-forces of mortality, for different ages crossing the years.
Nigri et al. (2019) extended the classical LC method by introducing an LSTM model for the time series prediction in the forecasting phase of the LC method, where the related time series (i.e., k t ) were extracted by following the same SVD method by Lee and Carter (1992). Their results show the prediction power of LSTM, compared to the classical time series prediction method (e.g., ARIMA), especially in capturing the nonlinearities.
More recently, Lindholm and Palmborg (2022) discussed the procedures to efficiently use training data in mortality forecasting when applying an LSTM-based Poisson LC method. Marino et al. (2022) further confirmed that an LSTM model can improve the predictive power of the classical LC method by providing a rigorous analysis of the prediction interval for their so-called LC–LSTM model. Note that Nigri et al. (2021) also used the LSTM model for life expectancy and lifespan disparity forecasting. According to the above-mentioned literature, the LSTM model has shown great prediction power for the forecasting of period effects in the classical LC model, as well as the forecasting of life expectancy and lifespan disparity. It is interesting to further consider the question of whether such a powerful tool (LSTM) can help improve long-term mortality forecasting for developing countries.
On the other hand, deep feedforward neural networks (FNNs), as a different way of extending the LC method, may also be applied to mortality modeling; e.g., see Richman and Wüthrich (2021), where they treat human mortality modeling as a classical supervised learning problem. For the application of convolution neural networks (CNNs) in mortality modeling, refer to Wang et al. (2021) and Schnürch and Korn (2022). However, their neural networks and methods are fundamentally different from the LSTM-based neural network structure here; therefore, their results are not directly comparable. Furthermore, some non-neural network machine learning models have appeared in the literature to forecast human mortality rates; for example, see Deprez et al. (2017) and Levantesi and Pizzorusso (2019).
As a result, in this paper, we develop an LSTM-based coherent mortality forecasting method with time-varying rotation structures based on a benchmark group. It provides a unified and model-free mortality projection method for developing countries. The paper is organized as follows: Section 2 introduces some preliminaries, including neural networks, the LSTM, the classical LC, and the Li–Lee model, as well as the definitions of life expectancy and lifespan disparity. The LSTM-based coherent mortality forecasting model is presented in Section 3. Finally, the mortality data and the empirical results are presented in Section 4, followed by the conclusion and some remarks in Section 5.

2. Preliminaries

2.1. RNN with LSTM Architecture

Recurrent neural networks (RNNs) are a class of artificial neural networks (ANNs), that can store representations of recent input data through their feedback connections. RNNs have many significant applications in the areas of the speech process, non-Markov control, or time series prediction (see, e.g., Mozer 1991). For instance, let { x 1 , x 2 , , x n } denote the time sequence of input vectors, and { h 1 , h 2 , , h n } denote the time sequence of output vectors; for the simple RNN, the output vector at time-step t is defined as follows:
h t = ϕ ( W h h h t 1 + W h x x t + b h ) ,
where ϕ is an activation function, W h h and W h x are the kernel weights for previous time step outputs and current inputs, respectively, and b h is the corresponding bias.
However, with the conventional gradient-based back-propagation through time (BPTT) algorithm (see Williams and Zipser (1995)), simple RNNs suffer from the problem of vanishing or exploding gradients (Pascanu et al. 2013). Then, RNNs with an LSTM architecture, or simply LSTMs, were introduced by Hochreiter and Schmidhuber (1997) in order to overcome such vanishing gradient problems. Instead of using all the memory dynamically when processing the data, the LSTM architecture relies both on the memory block and a few gates for controlling data elaborations.
LSTMs have shown great power in natural language processing and time series predictions. According to Marino et al. (2022), the LSTM can be expressed in the following mathematical form. Let N 0 denote the number of neurons within the input layer, N p denote the number of neurons of the p-th hidden layer with n { 1 , , P } , and N P + 1 denote the number of neurons of the output layer, where P, N 0 , and N p for p { 1 , , P } and N P + 1 N . Then, the activation of the p-th hidden layer may expressed as an affine mapping, A ( p ) : R N p 1 R N p , where R N p 1 refers to the output produced by the ( p 1 ) -th hidden layer. The output of an LSTM neuron at any time t in the p-th hidden layer can be expressed as follows:
h t ( p ) = o t ( p ) tanh ( c t ( p ) ) ,
where ⊙ denotes the element-wise product. The key to the LSTM lies in the following equations, which describe the outputs of four different gates in the architecture:
Forget gate : f t ( p ) = σ f A ( p ) = σ W f ( p ) , h t ( p 1 ) + U f ( p ) , h t 1 ( p ) + b f ( p ) , Input gate : i t ( p ) = σ i A ( p ) = σ W i ( p ) , h t ( p 1 ) + U i ( p ) , h t 1 ( p ) + b i ( p ) , Output gate : o t ( p ) = σ o A ( p ) = σ W o ( p ) , h t ( p 1 ) + U o ( p ) , h t 1 ( p ) + b o ( p ) , Memory state : c t ( p ) = f t ( p ) c t 1 ( p ) + i t ( p ) tanh W c ( p ) , h t ( p 1 ) + U c ( p ) , h t 1 ( p ) + b c ( p ) ,
where σ ( x ) = ( 1 + e x ) 1 is the sigmoid activation function, tanh ( x ) = e x e x e x + e x is the hyperbolic tangent activation function, W k ( p ) for k = f , i , o , c are, respectively, the weight matrices for the four different gates of feedforward connections in the structure, U k ( p ) for k = f , i , o , c are the corresponding weight matrices for the gates of recurrent connections, and b k ( p ) for k = f , i , o , c are the bias terms in the model.
Then let D = ( x t , y t ) , x t R N 0 , y t R N P + 1 be a dataset where x t and y t are the input variables and associated responses at time t, respectively. Hence, the LSTM is essentially a function, say g L S T M : R N 0 R N P + 1 , with
y t = g L S T M ( x t ; W ) + γ t = ψ h t ( p ) h t ( p 1 ) h t ( 1 ) ( x t ; W ) + γ t ,
where ψ : R N p R N P + 1 is the activation function at the output layer, W is the set of all weight parameters in the network, and γ t is a noise term, with zero mean and variance σ t 2 , independent of g L S T M .

2.2. The Mortality Models

Here, we use two classical mortality projection methods, namely the LC models and Li–Lee model. Specifically, we use the LC method for a rough/first-step mortality forecast of the target developing countries. That is, for such single populations, the Lee–Carter method (Lee and Carter 1992) assumes that the logarithm of the crude death rates ( m x , t ) for each age x and year t satisfies (1) and (2), where a x summarizes the average level of mortality throughout the time at age x, k t provides the overall level of mortality at year t, and b x measures the age effect of mortality on different periods. Note that the following is also assumed in the LC method (for identification purposes):
a x = 1 T + 1 t = 0 T log m x , t , x = 0 ω b x = 1 .
On the other hand, the Li–Lee method introduced by Li and Lee (2005) is an extension of the classical LC method, which can generate coherent mortality projections for multiple countries. In this study, the Li–Lee method (given by (3) and (4)) is used as the first step in mortality forecasting of the benchmark countries (i.e., a selected group of developed countries). Note that K t , in general, can be fitted by a random walk with drift (i.e., non-stationary), whereas k i , t for all i is assumed to be stationary; that is, k i , t shall be fitted by a random walk without drift or first-order autoregressive model A R ( 1 ) with a coefficient that yields a bounded short-term trend; for more details, refer to Li and Lee (2005).
It is obvious that under the Li–Lee method, the long-term mortality trend is uniquely determined by the common period effect K t , which makes the mortality forecasts coherent for all member countries in the group. In the empirical application, if the k i , t of a country is non-stationary, then this country is considered non-coherent with other countries in the group, i.e., there is significant divergence between the historical mortality experience and the common mortality patterns B x and K t . Hence, we may need to exclude them from the selected benchmark group.
Note that in order to ensure comparability between the parameters in the LC and Li–Lee methods in our subsequent analysis, one should impose the same normalization constraints on the key parameters of the two methods (see, e.g., Li et al. 2018).

2.3. Life Expectancy and Lifespan Disparity

Most mortality forecasting methods aim to predict how many additional years of life people will gain in the future. Life expectancy at birth, which measures the central tendency, is frequently applied to evaluate the precision of mortality forecasting methods. However, as mentioned by Aburto et al. (2020), populations with the same life expectancy level might experience substantial differences in the time of death, i.e., life expectancy cannot detect distributional variations in lifespan. Therefore, lifespan disparity can serve as an additional indicator to evaluate mortality forecasting methods (see, e.g., Bohk-Ewald et al. 2017).
Let us introduce the notation and definitions of life expectancy and lifespan disparity. Let S ( x , t ) and μ ( x , t ) denote the survival function and the force of mortality for an individual age x at time t, respectively, for a given population. These are assumed to be two continuous functions with respect to x and t. Also, denote by e x , t the life expectancy for age x at time t as
e x , t = x S ( y , t ) d y S ( x , t ) ,
where
S ( x , t ) = exp 0 x μ ( a , t ) d a ,
and μ ( a , t ) is the corresponding force of mortality at age a and time t.
To measure lifespan disparity, we take the average number of life years lost at birth (see Vaupel and Romo 2003; Zhang and Vaupel 2009):
e 0 , t = 0 e y , t d ( y , t ) d y ,
where d ( y , t ) are the deaths at age y and time t, and e y , t is the remaining life expectancy at age y and time t. Obviously, (7) shows that lifespan disparity is an indicator representing the life expectancy lost due to death by an individual at age x at time t. Note that lifespan disparity can be described by other measures, such as the standard deviation, interquartile range, Gini coefficient, or prolate index. However, we shall use e 0 , t for lifespan disparity in our analysis. Demographically, apart from its interpretation as the average life years lost or lost living potential, it also provides information about the capacity for further increases in life expectancy (see Bohk-Ewald et al. 2017).

3. LSTM-Based Coherent Method

In this section, we introduce our LSTM-based coherent mortality forecasting method with an embedded rotation in the time-varying model weight parameters during the projection phase for developing countries. The key to our method is a time-varying LC model, see (8) below, with an LSTM-based component that controls the rotation of the weight parameters in the model. Here, the rotation is referred to as the gradual change in the mortality development pattern during the projection phase (see Li et al. 2021). The time-varying parameters are defined through time-dependent weighted averages of the corresponding projected parameters based on the historical data from the target developing country and the selected benchmark group, respectively, where the weights are rotated based on how close the target developing country is to the benchmark group, under given criteria.
More specifically, we consider the following projection method for the logarithm of central death rate m x , t j for a particular developing country, say j at age x and year t, such that
ln m x , t j = a x j + b x , t j k t j + ε x , t , k t j = d t j + k t 1 j + ϵ t ,
where a x j can be estimated by the average mortality level at age x for developing country j. Here, k t j is the period effect and ε x , t and ϵ t are two zero mean random noises.
The main difference between our method, in (8), and the classical LC method is the time-varying b x , t j that measures a time-dependent age effect on mortality at different periods, and d t j describes a time-dependent drift in the random walk model used to project k t j .
Then, we select a group of developed countries as the benchmark group, where the classical Li–Lee method is applied in the projection of the logarithm of the central death rates m i , x , t , for member i in the group, that is,
ln m i , x , t = a i , x + B x K t + b i , x k i , t + ε i , x , t ,
with
K t = d 0 + K t 1 + ν t , k i , t = α 0 , i + α 1 , i k i , t 1 + ϵ i , t .
B x measures the common age effect in the benchmark group, and d 0 gives the drift in the random walk model of the common period effect K t . B x and d 0 are the two key components to be extracted from the benchmark group, using the Li–Lee model, to then be used in the rotation during the projection phase.
Here is how the rotation works in our method. As explained in Section 2, first select the life expectancy and lifespan disparity, respectively, as the criteria to describe the gap in terms of mortality levels, between the target developing country and the benchmark group. The varying gap will control the rotation of age and period effects in the projection phase for developing countries.
However, instead of selecting and fitting various (double) logistic functions as well as tailored threshold levels for the gap (see Li et al. 2021), we propose using a unified LSTM model for the gap forecasting. In addition, we use life expectancy, given by (6), and lifespan disparity, given by (7), respectively, in the construction of the gap function that describes the mortality distance between the target developing country and the benchmark group.
In particular, for notation simplicity, we let e u i and e u i denote, respectively, the (projected) life expectancy and lifespan disparity at birth for the i-th member in the benchmark group, in year u, and define the corresponding average life expectancy and lifespan disparity at birth, in year u, for the whole benchmark group as follows:
e u b = 1 N i = 1 N e u i , for u = , T , , e u b = 1 N i = 1 N e u i , for u = , T , ,
where T is the number of years in the training/in-sample data and N is the total number of members in the benchmark group.
Let e u j and e u j denote the corresponding life expectancy and lifespan disparity at birth for the target developing country/region j, in year u, for u = , T , . A unified LSTM model is introduced to forecast the life expectancy and lifespan disparity for both the target developing country and the benchmark group, such that the projected gaps in mortality levels between the developing country and the benchmark group can be expressed as the forecast for life expectancy or lifespan disparity difference.
More specifically, we construct LSTM models for the projection of e t · and e t · for both target developing countries/regions (i.e., e t j and e x , t j ) and the benchmark group (i.e., e t b and e t b ) as follows:
e t · = g L S T M e ( e t 1 · ; W e ) + ϵ t e , e t · = g L S T M e ( e t 1 · ; W e ) + ϵ t e ,
where ϵ t e and ϵ t e are zero mean errors. g L S T M e and g L S T M e are given by (5), respectively, for life expectancy and lifespan disparity. And W e and W e are the weight parameters in the corresponding LSTM models (see for example Nigri et al. 2021). The parameters in the LSTM model (i.e., the functional form of g L S T M e and g L S T M e ) are optimized using an L 2 loss function, namely
min W e 1 2 t e t · g L S T M e ( e t 1 · ; W e ) 2 , min W e 1 2 t e t · g L S T M e ( e t 1 · ; W e ) 2 .
Note that, to show the long-term projection power of our method, we need a sufficient number of years of mortality rates in the out-of-sample data, which results in a limited size of the in-sample mortality data for the training of the neural network model. Therefore, in this paper, we only consider a first-order autoregressive approach in the LSTM model, as illustrated in (9), where the neural network learns at each time step the relationship between two consecutive values during the training period (i.e., one-to-one structure). In addition, the method can be extended to more complex LSTM models with the structure of many-to-one or many-to-many, given that the available data are sufficiently large.
Next, we illustrate in detail our LSTM-based coherent mortality forecasting method. Consider the method based on lifespan disparity; for the case with life expectancy, one simply replaces all e t · by e t · in the corresponding equations. As noted, at the core of the method is the time-varying LC model given in (8), where the term for the age effect, b x , t j , and the drift term, d t j , of the period effect, depend on time t. The time dependence is described through a set of time-varying weights, in terms of lifespan disparity gaps, linking the mortality improvements between the target developing country and the benchmark group. Let b ^ x j and d ^ j denote the estimated age effect term and the drift parameter of the period effect term for the target developing country j, based on the classical LC method.
Now, let B ^ x denote the estimated common age effect term and d ^ 0 denote the drift parameter of the common period effect term, obtained by using Li–Lee’s method for the benchmark group. These two parameters provide information on the common mortality improvements of the benchmark group. Then, the next step is to specify how the time-varying b x , t j and d t j in (8) are defined in the mortality projection phase. To be specific, at the beginning of the projection phase, one can simply rely on the historical mortality data of the developing country when forecasting the short-term mortality rates.
In addition, denote the lifespan disparity at birth at time t projected through the LSTM model given in (9) as e t j and e t b for the developing country/region j and the benchmark group, respectively. Then, define the lifespan disparity gap at time t between the target country/region and the benchmark group as follows:
g t : = e t b e t j .
Then, for intermediate or long-term projections, include data from the benchmark group so that the long-term mortality development of the developing country converges gradually to the common trend in this benchmark group. Hence, redefine the age effect and drift terms of the period effect in the LC model as
b x , t + 1 j : = ( 1 ω t ) b ^ x j + ω t B ^ x , d t + 1 j : = ( 1 ω t ) d ^ j + ω t d ^ 0 ,
where
ω t = 1 2 1 + sin π 2 2 max g T g t g T , 0 1 p ,
for each age x and t = T , T + 1 , , and g · is given by (10). ω t denotes the time-varying weights that link the projected time-dependent age and period effect parameters, in year t + 1 , to the weighted average of the estimated b ^ x j and d ^ j , respectively, with B ^ x and d ^ 0 in the first step (see, e.g., Li et al. 2013). To simplify the analysis, we apply here the same weight parameter for both b j and d j in (11). In addition, p [ 0 , 1 ] in (12) is a tuning parameter that controls the functional form of ω t . We set p = 1 in our analysis such that ω t has a considerably low rate of change when its value is close to zero or one. Note that when t = T , we have ω t = 0 , which means that at the beginning of the projection phase, the method relies only on the historical data from the developing country. For t > T , the lifespan disparity gap decreases, such that ω t increases smoothly to one if the lifespan disparity gap diminishes in the future projection phase. Note that if the projected lifespan disparity gap for a particular developing country diverges (e.g., g t > g T for t > T ), we simply forecast the mortality based on its own historical data (that is ω t = 1 for all t > T ).
Finally, we summarize the model structure schematically in Figure 1 below. To be specific, the LSTM-based coherent mortality forecasting model contains two parts: (1) on the left of Figure 1 is a neural network component that contains an input layer, two layers of the LSTM structure (both LSTM layers contain 128 neurons with accompanied dropout layers), and two (fully connected) dense layers that contain 64 and 32 neurons, respectively, with accompanied dropout layers for the output. Note that adding dense layers in the model can provide flexibility in the control of the non-linearity of the model. (2) On the right of Figure 1, the projected life expectancy or lifespan disparity is transferred into a component of the rotation algorithm for the calculating of time-varying weights, and then the projected weights are applied to the time-varying LC model for the forecasting of mortality rates.

4. Empirical Analysis

This section presents the application of our LSTM-based coherent mortality forecasting method to three developing countries, namely China, Brazil, and Nigeria, which are the most populous countries in their respective continents and also belong to the emerging/emerged markets in the world. According to BBVA (2014), China and Brazil are classified as EAGLEs, i.e., emerging and growth-leading economies that are expected to have GDP increments larger than the average of G7 economies, excluding the US, in the next ten years. Nigeria is classified as NEST, i.e., an emerging country that is expected to have GDP increments lower than the average of the G7—excluding the US but higher than Italy’s—in the next ten years. In addition, we apply our method to three developing regions, namely less developed region(s) (LDR), less developed regions excluding China (LDRexChina), and less developed regions excluding the least developed countries (LDRexLDC). The United Nations defines the less developed countries/regions as all regions of Africa, Asia (except Japan), Latin America, and the Caribbean, plus Melanesia, Micronesia, and Polynesia, and categorizes 45 countries as the least developed countries (UN Source: https://unctad.org/topic/least-developed-countries/list (accessed on 21 November 2023)), including 33 countries in Africa, 8 countries in Asia, 1 in the Caribbean, and 3 in the Pacific. To proceed with the empirical results, we first introduce the mortality data used in the analysis.

4.1. Mortality Data

In this study, the benchmark group is made up of nine selected developed countries, namely Denmark, Finland, France, the Netherlands, Switzerland, Sweden, the UK, the US, and Japan. The mortality rates of these countries are obtained from the Human Mortality Database. In particular, we use the central death rates in the one-age and one-year blocks, i.e., ages equal to 0 , 1 , 2 , 3 , , 97 , 98 , 99 , and years ranging from 1950 to 2019.
The mortality data for the six target developing countries/regions mentioned above are not included in the Human Mortality Database. Hence, the corresponding data are obtained from the population division of the United Nations (UN Source: https://population.un.org/wpp/Download/Standard/Mortality/ (accessed on 21 November 2023)). Note that, a necessary condition for the application of our method is that the life expectancy or lifespan disparities of the target countries/regions converge to the ones of the benchmark group. Hence, a preliminary study is needed to select developed countries that can form a proper benchmark group. Figure 2 illustrates the convergence of life expectancy and lifespan disparity at birth between China and the benchmark group. The life expectancy and lifespan disparity at birth in the Ukraine do not converge. According to Figure 2, one can recognize a spike around the year 1960 in both the life expectancy and lifespan disparity in China. Such mortality outliers are due to the so-called Great Chinese Famine of 1959 to 1961. Hence, in order to reduce the effects of such extreme outliers from China, in the following analysis, we use only the data from the year 1962 to 2019 whenever the data of China are involved (i.e., the cases with China, LDR, and LDRexLDC).

4.2. LSTM for Life Expectancy and Lifespan Disparity

As discussed above, in order to construct time-varying weights that depend on the convergence of the life expectancy and lifespan disparity of a developing country to those of the benchmark group, one needs to develop projection models for the corresponding life expectancy or lifespan disparity gap.
Note that in the literature, (see Li et al. 2021), the forecasting of the life expectancy gap based on statistical methods uses different functions (logistic or double logistic) for developing countries/regions. Also, exogenous thresholds need to be introduced to test convergence in the model. The situation is even more complex if the lifespan disparity gap is also introduced in the method.
Here, instead of fitting various functions with thresholds to the life expectancy and lifespan disparity gaps, we use a unified LSTM model (see, e.g., Nigri et al. 2021) for the forecasting and identification of the gaps, compared with the benchmark group, for both the life expectancy and lifespan disparity of all six countries/regions. The projected life expectancy and lifespan disparity using the LSTM model for the six target countries/regions are presented in Appendix A.1. For similar results regarding the projection of life expectancy and lifespan disparity of other countries selected from the Human Mortality Database, refer to Nigri et al. (2021).

4.3. Empirical Results

In this empirical study, we carry out an out-of-sample test when training the model with the above-mentioned mortality dataset. For the purpose of long-term predictions, the data are divided into two parts, where the first 35 years of data, from 1950 to 1984 (from 1962 to 1984 for China, LDR, and LDRexLDC), are used for training, and the rest of the data, from 1985 to 2019, are set aside as the test data. To avoid possible overfitting to the training dataset, 20% of the training data are selected randomly as a validation part at each epoch.
To assess our models’ projections accuracy, the criteria used are mean square error, root mean square error, and mean absolute error for the projected log-mortality rates in the test data. The forecasting results for the LSTM-based, time-varying LC method, which includes life expectancy and lifespan disparity, respectively, are compared with the traditional LC and Li–Lee methods. All the experiments were performed using Keras with TensorFlow in Python, for the LSTM model, the R package “StMoMo" for LC, and the Li–Lee method for the initial mortality data processing. Note that in the following tables, we use LSTM–ex to denote our forecasting model based on life expectancy, and use LSTM–disp to denote our model based on lifespan disparity.

4.4. Results for China, Brazil, and Nigeria

The first empirical results are for the application of our model to the mortality data of China, Brazil, and Nigeria. These selected target developing countries represent demographic trends in their continents. This strategy removes the effect of different ethnic groups on life expectancy or lifespan disparity and demonstrates the generality of the model. The six-year average projection errors are listed in Table 1, Table 2 and Table 3; more detailed results are presented in Appendix A.2, in terms of projection errors for each year, for males and females, respectively.
Table 1, Table 2 and Table 3 (see Figure A5, Figure A6 and Figure A7) show a clear cumulative error in long-term forecasts, which reveals the difficulty of long-term mortality forecasting, especially for the classical LC method. In our method, especially that based on lifespan disparity, this accumulation error is reduced to some extent, making long-term forecasting more reliable.
The results clearly show that the classical LC method underperforms, as it is based on only the historical mortality data of the target country or region. If the current mortality development trend in a developing country is not sustainable, mortality rates will gradually approach those of developed countries (like the benchmark group selected here). Hence, in such a setting, projections based solely on national mortality data are not reasonable. From the above results, the LSTM-based, time-varying LC method with lifespan disparity controlling the rotation in the time-dependent weights, is the most accurate one among the four methods examined here, especially for long-term projections. On the other hand, it is interesting to observe that Nigeria has the most significant projection error reduction when transferring from the classical LC method to our LSTM-based time-varying LC method.

4.5. Results for LDR, LDRexChina, LDRexLDC

Finally, in order to demonstrate the projection accuracy of our method, the following illustration is for the mortality data of three developing regions, denoted as less developed region(s) (LDR), less developed regions excluding China (LDRexChina), and less developed regions excluding least developed country (LDRexLDC); see Table 4, Table 5 and Table 6 and also Figure A8, Figure A9 and Figure A10 in Appendix A.3).
Note that the error fluctuations in the prediction results for most of these less developed regions are reduced significantly (see Figure A8, Figure A9 and Figure A10), which is reasonable since the less developed regions contain larger populations (i.e., more stable) compared to individual countries. Overall, the LSTM-based, time-varying LC method with lifespan disparity as the control of the rotation in time-varying weights provides the most accurate projections within the four methods examined here.
It is worth mentioning that both our LSTM-based–time-varying LC method and the Li–Lee method incorporate mortality trend corrections based on a benchmark group. However, the empirical results show that, for developing countries or regions, such corrections are better modeled through the projection of life expectancy or lifespan disparity difference with an LSTM model, especially for long-term forecasts.
To end this section, we draw heatmaps that show the relative prediction errors (i.e., (predicted value − actual value)/actual value) across all ages and years for the out-of-sample data. The results are presented in Figure 3 and Figure 4. For most ages and years in the six target developing countries/regions (except young males in China and females in Brazil), our model performs well. However, we also observe some cohort effects in the results, especially for the data from China. A possible improvement could be to extend (11) and (12) to include age dependency in our model. This is a non-trivial extension that will be left for future studies.

5. Conclusions

Mortality improvements are linked to social progress, for instance, in terms of health, nutrition, education, hygiene, and access to medical assistance. It is difficult to accurately predict mortality development trends, especially over a long-term period. For developing countries or regions, it is particularly important to provide accurate long-term predictions of mortality rates for each age in the population, given that the current mortality data might not reveal sustainable development trends in the long–run.
The proposal here is an LSTM-based coherent mortality forecasting method for developing countries, where the life expectancy and lifespan disparity gaps between the target developing country and the selected benchmark group are used for long-term projections. In particular, we allow the mortality development pattern of a developing country to be a weighted average of trends generated by its own historical data and the selected benchmark group. And the rotation in the time-varying weights is controlled by the projected life expectancy and lifespan disparity gaps between the developing country and the benchmark group. In addition, we introduce a unified deep neural network model with an LSTM architecture for the long-term forecasting of the gaps in life expectancy and lifespan disparity for all six developing countries and regions in our analysis.
We apply this LSTM-based coherent mortality forecasting method to three developing countries, China, Brazil, and Nigeria, and three developing regions defined by the United Nations, namely LDR, LDRexChina, and LDRexLDC. The empirical results show that the LSTM-based coherent forecasting method with lifespan disparity outperforms the classical LC and Li–Lee methods, as well as the one with life expectancy, especially for long-term projections.

Author Contributions

Conceptualization, J.G., Y.S. and R.X.; methodology, Y.S., R.X.; software, Y.S. and R.X.; validation, J.G., Y.S. and R.X.; formal analysis, Y.S. and R.X.; investigation, Y.S. and R.X.; resources, J.G., Y.S. and R.X.; data curation, Y.S. and R.X.; writing—original draft preparation, Y.S. and R.X.; writing—review and editing, J.G., Y.S. and R.X.; visualization, Y.S. and R.X.; supervision, J.G. and R.X.; project administration, J.G. and R.X.; funding acquisition, J.G. and R.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Sciences and Engineering Research Council (NSERC) of Canada grant numbers RGPIN–2017–06643 and DGDND–2017–00096 (Garrido and Xu, during his stay at Concordia University), by the National Natural Science Foundation of China grant numbers 12201506 and 72171055, and by the XJTLU Research Development Funding grant number RDF-20-01-02 (Shang and Xu at XJTLU).

Data Availability Statement

The mortality data for the benchmark group is publicly available on the Human Mortality Database (https://www.mortality.org/, accessed on 21 November 2023); the mortality data for the developing countries and regions are publicly available on United Nations data source (https://population.un.org/wpp/Download/Standard/Mortality/, accessed on 21 November 2023).

Acknowledgments

The authors thank the three anonymous referees for their constructive comments, which helped improve the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Appendix A.1. Life Expectancy and Lifespan Disparity Forecasts

Figure A1. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target countries.
Figure A1. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target countries.
Risks 12 00027 g0a1
Figure A2. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target regions.
Figure A2. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target regions.
Risks 12 00027 g0a2
Figure A3. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target countries.
Figure A3. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target countries.
Risks 12 00027 g0a3
Figure A4. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target regions.
Figure A4. Historical (dotted lines) and forecast (blue for ARIMA, red for LSTM) values of e 0 , t for three target regions.
Risks 12 00027 g0a4

Appendix A.2. Projection Errors for China, Brazil, and Nigeria

Figure A5. Forecasting errors by year for China (females at the top, males at the bottom).
Figure A5. Forecasting errors by year for China (females at the top, males at the bottom).
Risks 12 00027 g0a5
Figure A6. Forecasting errors by year for Brazil (females at the top, males at the bottom).
Figure A6. Forecasting errors by year for Brazil (females at the top, males at the bottom).
Risks 12 00027 g0a6
Figure A7. Forecasting errors by year for Nigeria (females at the top, males at the bottom).
Figure A7. Forecasting errors by year for Nigeria (females at the top, males at the bottom).
Risks 12 00027 g0a7

Appendix A.3. Projection Errors for LDR, LDRexChina, and LDRexLDC

Figure A8. Forecasting errors by year for less developed regions (LDR) (females at the top, males at the bottom).
Figure A8. Forecasting errors by year for less developed regions (LDR) (females at the top, males at the bottom).
Risks 12 00027 g0a8
Figure A9. Forecasting errors by year for less developed regions, excluding China (LDRexChina) (females at the top, males at the bottom).
Figure A9. Forecasting errors by year for less developed regions, excluding China (LDRexChina) (females at the top, males at the bottom).
Risks 12 00027 g0a9
Figure A10. Forecasting errors by year for less developed regions, excluding least developed country (LDRexLDC) (females at the top, males on the bottom).
Figure A10. Forecasting errors by year for less developed regions, excluding least developed country (LDRexLDC) (females at the top, males on the bottom).
Risks 12 00027 g0a10

References

  1. Aburto, José Manuel, Francisco Villavicencio, Ugofilippo Basellini, Søren Kjærgaard, and James W. Vaupel. 2020. Dynamics of life expectancy and life span equality. Proceedings of the National Academy of Sciences 117: 5250–59. [Google Scholar] [CrossRef]
  2. Antonio, Katrien, Anastasios Bardoutsos, and Wilbert Ouburg. 2015. Bayesian Poisson log-bilinear models for mortality projections with multiple populations. European Actuarial Journal 5: 245–81. [Google Scholar] [CrossRef]
  3. Austin, Kelly F., and Laura A. McKinney. 2012. Disease, war, hunger, and deprivation: A cross-national investigation of the determinants of life expectancy in less-developed and sub-Saharan African nations. Sociological Perspectives 55: 421–47. [Google Scholar] [CrossRef]
  4. BBVA. 2014. Eagles Economic Outlook Annual Report. Available online: https://www.bbvaresearch.com/wp-content/uploads/2014/05/2014_EAGLEs_Economic_Outllok-Annual.pdf (accessed on 13 January 2024).
  5. Bohk-Ewald, Christina, Marcus Ebeling, and Roland Rau. 2017. Lifespan disparity as an additional indicator for evaluating mortality forecasts. Demography 54: 1559–77. [Google Scholar] [CrossRef]
  6. Boonen, Tim J., and Hong Li. 2017. Modeling and forecasting mortality with economic growth: A multipopulation approach. Demography 54: 1921–46. [Google Scholar] [CrossRef]
  7. Brouhns, Natacha, Michel Denuit, and Jeroen K Vermunt. 2002. A Poisson log-bilinear regression approach to the construction of projected lifetables. Insurance: Mathematics and Economics 31: 373–93. [Google Scholar] [CrossRef]
  8. Cairns, Andrew J. G., David Blake, and Kevin Dowd. 2006. A two-factor model for stochastic mortality with parameter uncertainty: Theory and calibration. Journal of Risk and Insurance 73: 687–718. [Google Scholar] [CrossRef]
  9. Cairns, Andrew J. G., David Blake, and Kevin Dowd. 2008. Modelling and management of mortality risk: A review. Scandinavian Actuarial Journal 2008: 79–113. [Google Scholar] [CrossRef]
  10. Cairns, Andrew J. G., David Blake, Kevin Dowd, Guy D. Coughlan, David Epstein, and Marwa Khalaf-Allah. 2011a. Mortality density forecasts: An analysis of six stochastic mortality models. Insurance: Mathematics and Economics 48: 355–67. [Google Scholar] [CrossRef]
  11. Cairns, Andrew J. G., David Blake, Kevin Dowd, Guy D. Coughlan, and Marwa Khalaf-Allah. 2011b. Bayesian stochastic mortality modelling for two populations. ASTIN Bulletin: The Journal of the IAA 41: 29–59. [Google Scholar]
  12. Carter, Lawrence R., and Ronald D. Lee. 1992. Modeling and forecasting US sex differentials in mortality. International Journal of Forecasting 8: 393–411. [Google Scholar] [CrossRef]
  13. Currie, Iain D. 2016. On fitting generalized linear and non-linear models of mortality. Scandinavian Actuarial Journal 2016: 356–83. [Google Scholar] [CrossRef]
  14. Czado, Claudia, Antoine Delwarde, and Michel Denuit. 2005. Bayesian Poisson log-bilinear mortality projections. Insurance: Mathematics and Economics 36: 260–84. [Google Scholar] [CrossRef]
  15. Deprez, Philippe, Pavel V. Shevchenko, and Mario V. Wüthrich. 2017. Machine learning techniques for mortality modeling. European Actuarial Journal 7: 337–52. [Google Scholar] [CrossRef]
  16. Gao, Guangyuan, and Yanlin Shi. 2021. Age-coherent extensions of the Lee–Carter model. Scandinavian Actuarial Journal 2021: 998–1016. [Google Scholar] [CrossRef]
  17. Hainaut, Donatien. 2018. A neural-network analyzer for mortality forecast. ASTIN Bulletin: The Journal of the IAA 48: 481–508. [Google Scholar] [CrossRef]
  18. Hanewald, Katja. 2011. Explaining mortality dynamics: The role of macroeconomic fluctuations and cause of death trends. North American Actuarial Journal 15: 290–314. [Google Scholar] [CrossRef]
  19. Hochreiter, Sepp, and Jürgen Schmidhuber. 1997. Long short-term memory. Neural Computation 9: 1735–80. [Google Scholar] [CrossRef]
  20. Jeuland, Marc A., David E. Fuente, Semra Ozdemir, Maura C. Allaire, and Dale Whittington. 2013. The long-term dynamics of mortality benefits from improved water and sanitation in less developed countries. PLoS ONE 8: e74804. [Google Scholar] [CrossRef] [PubMed]
  21. Lee, Ronald D. 2000. The Lee–Carter method for forecasting mortality, with various extensions and applications. North American Actuarial Journal 4: 80–91. [Google Scholar] [CrossRef]
  22. Lee, Ronald D., and Lawrence R. Carter. 1992. Modeling and forecasting us mortality. Journal of the American Statistical Association 87: 659–71. [Google Scholar]
  23. Lee, Ronald D., and Francois Nault. 1993. Modeling and forecasting provincial mortality in Canada. Paper presented at the World Congress of the IUSSP, Montreal, QC, Canada, August 24–September 1. [Google Scholar]
  24. Levantesi, Susanna, and Virginia Pizzorusso. 2019. Application of machine learning to mortality modeling and forecasting. Risks 7: 26. [Google Scholar] [CrossRef]
  25. Li, Hong, and Yang Lu. 2017. Coherent forecasting of mortality rates: A sparse vector-autoregression approach. ASTIN Bulletin: The Journal of the IAA 47: 563–600. [Google Scholar] [CrossRef]
  26. Li, Hong, and Yang Lu. 2018. A Bayesian non–parametric model for small population mortality. Scandinavian Actuarial Journal 2018: 605–28. [Google Scholar] [CrossRef]
  27. Li, Hong, Yang Lu, and Pintao Lyu. 2018. Modeling and Forecasting Chinese Population Dynamics in a Multi-Population Context. SOA Research Reports. Schaumburg: Society of Actuaries. [Google Scholar]
  28. Li, Hong, Yang Lu, and Pintao Lyu. 2021. Coherent mortality forecasting for less developed countries. Risks 9: 151. [Google Scholar] [CrossRef]
  29. Li, Johnny Siu-Hang, Kenneth Q. Zhou, Xiaobai Zhu, Wai-Sum Chan, and Felix Wai-Hon Chan. 2019. A Bayesian approach to developing a stochastic mortality model for China. Journal of the Royal Statistical Society Series A: Statistics in Society 182: 1523–60. [Google Scholar] [CrossRef]
  30. Li, Nan, and Ronald D. Lee. 2005. Coherent mortality forecasts for a group of populations: An extension of the Lee–Carter method. Demography 42: 575–94. [Google Scholar] [CrossRef] [PubMed]
  31. Li, Nan, Ronald D. Lee, and Patrick Gerland. 2013. Extending the lee-Carter method to model the rotation of age patterns of mortality decline for long-term projections. Demography 50: 2037–51. [Google Scholar] [CrossRef]
  32. Lindholm, Mathias, and Lina Palmborg. 2022. Efficient use of data for LSTM mortality forecasting. European Actuarial Journal 12: 749–778. [Google Scholar] [CrossRef]
  33. Ma, Qingxiao, and Tim J. Boonen. 2023. Longevity risk modeling with the consumer price index. North American Actuarial Journal, 1–18. [Google Scholar] [CrossRef]
  34. Marino, Mario, Susanna Levantesi, and Andrea Nigri. 2023. A neural approach to improve the Lee–Carter mortality density forecasts. North American Actuarial Journal 27: 148–165. [Google Scholar] [CrossRef]
  35. Mozer, Michael C. 1991. Induction of multiscale temporal structure. Paper presented at the Advances in Neural Information Processing Systems 4, NIPS Conference, Denver, CO, USA, December 2–5. [Google Scholar]
  36. Müller, Olaf, and Michael Krawinkel. 2005. Malnutrition and health in developing countries. Cmaj 173: 279–86. [Google Scholar] [CrossRef]
  37. Nigri, Andrea, Susanna Levantesi, and Mario Marino. 2021. Life expectancy and lifespan disparity forecasting: A long short-term memory approach. Scandinavian Actuarial Journal 2021: 110–33. [Google Scholar] [CrossRef]
  38. Nigri, Andrea, Susanna Levantesi, Mario Marino, Salvatore Scognamiglio, and Francesca Perla. 2019. A deep learning integrated Lee–Carter model. Risks 7: 33. [Google Scholar] [CrossRef]
  39. Niu, Geng, and Bertrand Melenberg. 2014. Trends in mortality decrease and economic growth. Demography 51: 1755–73. [Google Scholar] [CrossRef]
  40. O’Hare, Colin, and Youwei Li. 2012. Explaining young mortality. Insurance: Mathematics and Economics 50: 12–25. [Google Scholar] [CrossRef]
  41. Pascanu, Razvan, Tomas Mikolov, and Yoshua Bengio. 2013. On the difficulty of training recurrent neural networks. Paper presented at International Conference on Machine Learning, Atlanta, GA, USA, June 16–21; pp. 1310–18. [Google Scholar]
  42. Pedroza, Claudia. 2006. A Bayesian forecasting model: Predicting US male mortality. Biostatistics 7: 530–50. [Google Scholar] [CrossRef] [PubMed]
  43. Pitacco, Ermanno. 2004. Survival models in a dynamic context: A survey. Insurance: Mathematics and Economics 35: 279–98. [Google Scholar] [CrossRef]
  44. Renshaw, Arthur E., and Steven Haberman. 2003. Lee–Carter mortality forecasting with age-specific enhancement. Insurance: Mathematics and Economics 33: 255–72. [Google Scholar] [CrossRef]
  45. Renshaw, Arthur E., and Steven Haberman. 2006. A cohort-based extension to the Lee–Carter model for mortality reduction factors. Insurance: Mathematics and Economics 38: 556–70. [Google Scholar] [CrossRef]
  46. Richman, Ronald, and Mario V. Wüthrich. 2021. A neural network extension of the Lee–Carter model to multiple populations. Annals of Actuarial Science 15: 346–66. [Google Scholar] [CrossRef]
  47. Schnürch, Simo, and Ralf Korn. 2022. Point and interval forecasts of death rates using neural networks. ASTIN Bulletin: The Journal of the IAA 52: 333–60. [Google Scholar] [CrossRef]
  48. Vaupel, James W., and Vladimir Canudas Romo. 2003. Decomposing change in life expectancy: A bouquet of formulas in honor of Nathan Keyfitz’s 90th birthday. Demography 40: 201–16. [Google Scholar] [CrossRef]
  49. Vaupel, James W., and Sabine Schnabel. 2004. Forecasting best-practice life expectancy to forecast national life expectancy. Paper presented at the 2004 Annual Meeting of the Population Association of America, Boston, MA, USA, April 1–3. [Google Scholar]
  50. Wang, Chou-Wen, Jinggong Zhang, and Wenjun Zhu. 2021. Neighbouring prediction for mortality. ASTIN Bulletin: The Journal of the IAA 51: 689–718. [Google Scholar] [CrossRef]
  51. White, Kevin M. 2002. Longevity advances in high-income countries, 1955–96. Population and Development Review 28: 59–76. [Google Scholar] [CrossRef]
  52. Williams, Ronald J., and David Zipser. 1995. Gradient-based learning algorithms for recurrent networks and their computational complexity. In Back-Propagation: Theory, Architectures, and Applications, 1st ed. Edited by Yves Chauvin and David E Rumelhart. London: Psychology Press, Taylor & Francis Group, chp. 13. [Google Scholar] [CrossRef]
  53. Wilmoth, John R. 1998. Is the pace of Japanese mortality decline converging toward international trends? Population and Development Review 24: 593–600. [Google Scholar] [CrossRef]
  54. Wong, Jackie S. T., Jonathan J. Forster, and Peter W. F. Smith. 2023. Bayesian model comparison for mortality forecasting. Journal of the Royal Statistical Society Series C: Applied Statistics 72: 566–86. [Google Scholar] [CrossRef]
  55. Wong-Fupuy, Carlos, and Steven Haberman. 2004. Projecting mortality trends: Recent developments in the United Kingdom and the United States. North American Actuarial Journal 8: 56–83. [Google Scholar] [CrossRef]
  56. Zhang, Zhen, and James W. Vaupel. 2009. The age separating early deaths from late deaths. Demographic Research 20: 721–30. [Google Scholar] [CrossRef]
Figure 1. LSTM-based coherent mortality forecasting model.
Figure 1. LSTM-based coherent mortality forecasting model.
Risks 12 00027 g001
Figure 2. China and Ukraine vs. benchmark group.
Figure 2. China and Ukraine vs. benchmark group.
Risks 12 00027 g002
Figure 3. Relative prediction errors for three developing countries (males at the top, females at the bottom).
Figure 3. Relative prediction errors for three developing countries (males at the top, females at the bottom).
Risks 12 00027 g003
Figure 4. Relative prediction errors for three developing regions (males at the top, females at the bottom).
Figure 4. Relative prediction errors for three developing regions (males at the top, females at the bottom).
Risks 12 00027 g004
Table 1. Six-year (average) prediction errors for China.
Table 1. Six-year (average) prediction errors for China.
1985–19911992–19981999–20052006–20122013–2019Total
MSELC0.08810.09650.08560.07990.09120.0883
LL0.02040.02840.03170.03630.05130.0336
LSTM–ex0.00920.02340.02450.03510.05010.0284
LSTM–disp0.0073 *0.0166 *0.0119 *0.0083 *0.0104 *0.0109 *
MAELC0.21050.22940.21290.2150.23310.2202
LL0.11480.1320.12940.13080.15430.1323
LSTM–ex0.05940.09420.10940.13830.17380.115
LSTM–disp0.0586 *0.0849 *0.0677 *0.0569 *0.0681 *0.0672 *
RMSELC0.24830.29060.28940.28180.2890.2798
LL0.13880.16730.17330.17930.20140.172
LSTM–ex0.09590.15310.15650.18730.22380.1685
LSTM–disp0.0854 *0.1288 *0.1091 *0.0911 *0.1021 *0.1044 *
* indicates the smallest value of MSE/MAE/RMSE.
Table 2. Six-year (average) prediction errors for Brazil.
Table 2. Six-year (average) prediction errors for Brazil.
1985–19911992–19981999–20052006–20122013–2019Total
MSELC0.00940.02590.08480.16610.26910.1111
Li–Lee0.00610.01980.06470.12360.19990.0829
LSTM–ex0.00490.01950.06240.11250.17010.0739
LSTM–disp0.0038 *0.0024 *0.0157 *0.0205 *0.0566 *0.0198 *
MAELC0.07620.11910.20810.30840.41080.2245
Li–Lee0.06420.10330.17290.24420.32120.1811
LSTM–ex0.05060.10750.19960.25630.31610.1861
LSTM–disp0.0501 *0.0324 *0.0839 *0.0837 *0.1536 *0.0807 *
RMSELC0.09550.14610.26630.39240.51190.3239
Li–Lee0.07730.13350.24190.34570.44550.2845
LSTM–ex0.06970.13970.24960.33540.41210.2718
LSTM–disp0.0614 *0.0488 *0.1165 *0.1239 *0.2281 *0.1332 *
* indicates the smallest value of MSE/MAE/RMSE.
Table 3. Six-year (average) prediction errors for Nigeria.
Table 3. Six-year (average) prediction errors for Nigeria.
1985–19911992–19981999–20052006–20122013–2019Total
MSELC0.00350.00820.03620.09610.17260.0633
Li–Lee0.00070.00220.01020.02860.05420.0192
LSTM–ex0.00450.02910.03340.04220.05210.0323
LSTM–disp0.0011 *0.0008 *0.0002 *0.0006 *0.0013 *0.0008 *
MAELC0.04320.06890.16310.27710.37780.1861
Li–Lee0.02090.03690.08350.14690.20590.0988
LSTM–ex0.05830.15820.16860.19550.21640.1594
LSTM–disp0.0247 *0.0209 *0.0102 *0.0144 *0.0203 *0.0181 *
RMSELC0.05070.07590.18250.30590.41160.2481
Li–Lee0.02240.03910.09590.16630.23020.1361
LSTM–ex0.06450.16810.18220.20450.22620.1796
LSTM–disp0.0306 *0.0278 *0.0121 *0.0221 *0.0331 *0.0261 *
* indicates the smallest value of MSE/MAE/RMSE.
Table 4. Six-year (average) prediction errors for LDR.
Table 4. Six-year (average) prediction errors for LDR.
1985–19911992–19981999–20052006–20122013–2019Total
MSELC0.00940.03070.06750.10330.15780.0738
LL0.00320.01140.02930.05240.08290.0358
LSTM–ex0.0011 *0.009050.024250.04020.050950.0251
LSTM–disp0.00190.00135 *0.0008 *0.0031 *0.00335 *0.0021 *
MAELC0.08080.15370.22690.27720.34720.2172
LL0.04410.09060.14080.18470.24430.1409
LSTM–ex0.0253 *0.0710.11610.15760.18070.1101
LSTM–disp0.03440.0275 *0.0216 *0.0427 *0.0372 *0.033 *
RMSELC0.09690.17360.25850.32230.3970.2497
LL0.05290.10780.16880.22730.28760.1689
LSTM–ex0.0322 *0.08970.15110.19750.22290.1387
LSTM–disp0.04330.0349 *0.0285 *0.0554 *0.0549 *0.0434 *
* indicates the smallest value of MSE/MAE/RMSE.
Table 5. Six-year (average) prediction errors for LDRexChina.
Table 5. Six-year (average) prediction errors for LDRexChina.
1985–19911992–19981999–20052006–20122013–2019Total
MSELC0.00910.02710.07110.13580.22410.0934
Li–Lee0.00370.01260.03650.07380.12490.0504
LSTM–ex0.0008 *0.00560.01460.02020.02050.0123
LSTM–disp0.00230.0008 *0.0011 *0.0031 *0.0032 *0.0021 *
MAELC0.05560.11150.18610.25270.34040.1893
Li–Lee0.07320.13780.21910.30670.40220.2278
LSTM–ex0.0233 *0.05620.08980.11360.11150.0789
LSTM–disp0.03680.0221 *0.0222 *0.0411 *0.0369 *0.0317 *
RMSELC0.09380.16430.26640.36770.47270.3051
Li–Lee0.06040.11120.18890.26920.35110.2227
LSTM–ex0.0283 *0.07470.12030.14220.14280.1109
LSTM–disp0.04760.0283 *0.0332 *0.0541 *0.0516 *0.0447 *
* indicates the smallest value of MSE/MAE/RMSE.
Table 6. Six-year (average) prediction errors for LDRexLDC.
Table 6. Six-year (average) prediction errors for LDRexLDC.
1985–19911992–19981999–20052006–20122013–2019Total
MSELC0.01990.03280.050.07330.10250.0557
LL0.00740.01260.01820.03060.04540.0228
LSTM–ex0.0014 *0.00650.02540.04930.0640.0293
LSTM–disp0.00440.0032 *0.0009 *0.0061 *0.0022 *0.0034 *
MAELC0.11830.16750.2030.24060.29030.204
LL0.07220.09830.11650.14850.18970.1251
LSTM–ex0.0301 *0.0590.11760.17260.19560.115
LSTM–disp0.05510.045 *0.023 *0.0639 *0.037 *0.0448 *
RMSELC0.13990.1810.22240.26890.320.2264
LL0.08550.10730.13130.17140.21170.1414
LSTM–ex0.0368 *0.07720.15820.22170.25260.1493
LSTM–disp0.06610.0541 *0.0297 *0.0782 *0.04670.055 *
* indicates the smallest value of MSE/MAE/RMSE.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Garrido, J.; Shang, Y.; Xu, R. LSTM-Based Coherent Mortality Forecasting for Developing Countries. Risks 2024, 12, 27. https://doi.org/10.3390/risks12020027

AMA Style

Garrido J, Shang Y, Xu R. LSTM-Based Coherent Mortality Forecasting for Developing Countries. Risks. 2024; 12(2):27. https://doi.org/10.3390/risks12020027

Chicago/Turabian Style

Garrido, Jose, Yuxiang Shang, and Ran Xu. 2024. "LSTM-Based Coherent Mortality Forecasting for Developing Countries" Risks 12, no. 2: 27. https://doi.org/10.3390/risks12020027

APA Style

Garrido, J., Shang, Y., & Xu, R. (2024). LSTM-Based Coherent Mortality Forecasting for Developing Countries. Risks, 12(2), 27. https://doi.org/10.3390/risks12020027

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop