Next Article in Journal
Effect of Variance Swap in Hedging Volatility Risk
Next Article in Special Issue
The Impact of Model Uncertainty on Index-Based Longevity Hedging and Measurement of Longevity Basis Risk
Previous Article in Journal
Neural Networks and Betting Strategies for Tennis
Previous Article in Special Issue
A Two-Population Extension of the Exponential Smoothing State Space Model with a Smoothing Penalisation Scheme
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Retiree Mortality Forecasting: A Partial Age-Range or a Full Age-Range Model?

1
Department of Actuarial Studies and Business Analytics, Macquarie University, Balaclava Rd, Macquarie Park, NSW 2109, Australia
2
Cass Business School, City, University of London, London EC1Y 8TZ, UK
*
Author to whom correspondence should be addressed.
Submission received: 2 June 2020 / Revised: 24 June 2020 / Accepted: 24 June 2020 / Published: 1 July 2020
(This article belongs to the Special Issue Mortality Forecasting and Applications)

Abstract

:
An essential input of annuity pricing is the future retiree mortality. From observed age-specific mortality data, modeling and forecasting can take place in two routes. On the one hand, we can first truncate the available data to retiree ages and then produce mortality forecasts based on a partial age-range model. On the other hand, with all available data, we can first apply a full age-range model to produce forecasts and then truncate the mortality forecasts to retiree ages. We investigate the difference in modeling the logarithmic transformation of the central mortality rates between a partial age-range and a full age-range model, using data from mainly developed countries in the Human Mortality Database (2020). By evaluating and comparing the short-term point and interval forecast accuracies, we recommend the first strategy by truncating all available data to retiree ages and then produce mortality forecasts. However, when considering the long-term forecasts, it is unclear which strategy is better since it is more difficult to find a model and parameters that are optimal. This is a disadvantage of using methods based on time-series extrapolation for long-term forecasting. Instead, an expectation approach, in which experts set a future target, could be considered, noting that this method has also had limited success in the past.

1. Introduction

Improving human survival probability contributes greatly to an aging population. To guarantee one individual’s financial income in retirement, a policyholder may purchase a fixed-term or lifetime annuity. A fixed-term or lifetime annuity is a contract offered by insurers guaranteeing regular payments in exchange for an initial premium. Since an annuity depends on survival probabilities and interest rates, pension funds and insurance companies are more likely to face a risk of longevity. Longevity risk is a potential systematic risk attached to the increasing life expectancy of policyholders, which can eventually result in a higher payout ratio than expected (Crawford et al. 2008). The concerns about longevity risk have led to a surge of interest in modeling and forecasting age-specific mortality rates.
Many models for forecasting age-specific mortality indicators have been proposed in demographic literature (see Booth and Tickle 2008, for reviews). Of these, Lee and Carter (1992) implemented a principal component method to model the logarithm of age-specific mortality rates ( m x , t ) and extracted a single time-varying index representing the trend in the level of mortality, from which the forecasts are obtained by a random walk with drift. Since then, the Lee–Carter (LC) method has been extended and modified (see Booth and Tickle 2008; Pitacco et al. 2009; Shang and Haberman 2018; Shang et al. 2011, for reviews). The LC method has been applied to many countries, including Belgium (Brouhns et al. 2002), Austria (Carter and Prkawetz 2001), England and Wales (Renshaw and Haberman 2003), and Spain (Debón et al. 2008; Guillen and Vidiella-i-Anguera 2005).
Many statistical models focus on time-series extrapolation of past trends exhibited in age-specific mortality rates (see, e.g., Booth and Tickle 2008). We consider two modeling strategies: On the one hand, we can first truncate all available data to retiree ages and then produce mortality forecasts (see, e.g., Cairns et al. 2006; 2009). On the other hand, we can first use the available data to produce forecasts and then truncate the mortality forecasts to retiree ages (see, e.g., Shang and Haberman 2017). In this paper, our contribution is to investigate the difference in modeling the logarithmic transformation of the central mortality rates between a partial age-range and a full age-range model, using data from mainly developed countries in the Human Mortality Database (2020). By evaluating and comparing the short-term point and interval forecast accuracies, we recommend the first strategy by truncating all available data to retiree ages and then produce mortality forecasts. However, when we consider the long-term forecasts, it is unclear which strategy is better since it is more difficult to find a model and parameters that are optimal. This is a disadvantage of using methods based on time-series extrapolation for long-term forecasting. Instead, an expectation approach, in which experts set a future target, could be considered, noting that this method has also had limited success in the past (Booth and Tickle 2008). Our recommendations could be useful to actuaries for choosing a better modeling strategy and more accurately pricing a range of annuity products.
The article is organized as follows: In Section 2, we describe the mortality data sets of 19 mainly developed countries. We revisit five time-series extrapolation models for forecasting age-specific mortality rates, which have been shown in the literature to work well across the full age range for some data sets (for more details, consult Shang 2012; Shang and Haberman 2018). Using these models as a testbed, we compare point and interval forecast accuracies between the two modeling strategies and provide our recommendations in Section 3. Conclusions are presented in Section 4.

2. Data Sets

The data sets used in this study were taken from the Human Mortality Database (2020). For each sex in a given calendar year, the mortality rates obtained by the ratio between “number of deaths" and “exposure to risk" are arranged in a matrix for age and calendar year. Nineteen countries, mainly developed countries, were selected, and thus 38 sub-populations of age- and sex-specific mortality rates were obtained for all analyses. The 19 countries selected all have reliable data series commencing at/before 1950. Due to possible structural breaks (i.e., two world wars), we truncate all data series from 1950 onwards. The omission of Germany is because the Human Mortality Database (2020) for a reunited Germany only dates back to 1990. The selected countries and their abbreviations are shown in Table 1, along with their last year of available data (recorded in April 2019). To avoid fluctuations at older ages, we consider ages from 0 to 99 in a single year of age and the last age group is from 100 onwards. Should we consider all ages from 0 to 110+, and we may encounter the missing-value issue and observe mortality rates outside the range of [ 0 , 1 ] for some years.

3. Results

3.1. Forecast Evaluation

We present 19 countries that begin in 1950 and end in the last year listed in Table 1. We keep the last 30 observations for forecasting evaluation, while the remaining observations are treated as initial fitting observations, from which we produce the one-step-ahead to 30-step-ahead forecasts. Via an expanding window approach (see also Zivot and Wang 2006), we re-estimate the parameter in the time-series forecasting models by increasing the fitted observations by one year and produce the one-step-ahead to 29-step-ahead forecasts. We iterate this process by increasing the sample size by one year until the end of the data period. The process produces 30 one-step-ahead forecasts, 29 two-step-ahead forecasts, ⋯, one 30-step-ahead forecast. We compare these forecasts with the holdout samples to determine the out-of-sample forecast accuracy.

3.2. Forecast Error Criteria

To evaluate the point forecast accuracy, we consider the mean absolute percentage error (MAPE) and root mean squared percentage error (RMSPE). The MAPE and RMSPE criteria measure how close the forecasts compare with the actual values of the variable being forecast, regardless of the error sign. The MAPE and RMSPE criteria can be expressed as:
MAPE h = 1 p × ( 31 h ) j = 1 ( 31 h ) x = 1 p m x , j m ^ x , j m x , j × 100 , x = 1 , , p , RMSPE h = 1 p × ( 31 h ) j = 1 ( 31 h ) x = 1 p m x , j m ^ x , j m x , j 2 × 100
where m x , j represents the actual holdout sample for age x in the forecasting year j, p denotes the total number of ages, and m ^ x , j represents the forecasts for the holdout sample.
To evaluate the pointwise interval forecast accuracy, we consider the interval score criterion of Gneiting and Raftery (2007). We consider the common case of the symmetric 100 ( 1 α ) % prediction intervals, with lower and upper bounds that were predictive quantiles at α / 2 and 1 α / 2 , denoted by m ^ x , j lb and m ^ x , j ub . As defined by Gneiting and Raftery (2007), a scoring rule for evaluating the pointwise interval forecast accuracy at time point j is
S α m ^ x , j lb , m ^ x , j ub ; m x , j = m ^ x , j ub m ^ x , j lb + 2 α m ^ x , j lb m x , j 𝟙 m x , j < m ^ x , j lb 2 α m x , j m ^ x , j ub 𝟙 m x , j > m ^ x , j ub ,
where 𝟙 { · } denotes the binary indicator function. The optimal interval score is achieved when m x , j lies between m ^ x , j lb and m ^ x , j ub , with the distance between the upper and lower bounds being minimal. To obtain summary statistics of the interval score, we take the mean interval score across different ages and forecasting years. The mean interval score can be expressed as:
S ¯ α , h = 1 p × ( 31 h ) j = 1 ( 31 h ) x = 1 p S α m ^ x , j lb , m ^ x , j ub ; m x , j .

3.3. Comparison of Point and Interval Forecast Errors

Modeling and forecasting mortality can be taken place in two routes. On the one hand, we can first truncate the available data to certain ages, such as from 60 to 99 in a single year of age and 100+ as the last age, and then produce mortality forecasts for these retiree ages. On the other hand, we can first use the available data, i.e., age-specific mortality from 0 to 99 in a single year of age and 100+ as the last age, to produce forecasts for these 101 ages and then truncate the mortality forecasts to certain ages, such as 60 to 99 in a single year of age and 100+ as the last age.
We study five time-series extrapolation models for forecasting age-specific mortality, which have been shown in the literature to work well across the full age range for some data sets. Please note that the Cairns-Blake-Dowd suite of models are not included in this paper, because they are designed just for ages 55 and over. The models that we have considered are subjective and far from extensive, but they suffice to serve as a testbed for comparing the forecast accuracy. These five models are: the Lee–Carter model with Poisson errors (see, e.g., Brouhns et al. 2002; Renshaw and Haberman 2003; 2006), the Lee–Carter model with Gaussian errors (see, e.g., Booth et al. 2002; Koissi et al. 2006; Renshaw and Haberman 2003), age-period-cohort model (see, e.g., Renshaw and Haberman 2006) and the Plat model (Plat 2009).
For the short-term forecast horizon (i.e., the one-step-ahead forecast horizon), we compute the mean of the MAPEs and RMSPEs to evaluate the point forecast accuracy. From Table 2, there is an advantage of directly modeling and forecasting the truncated series for the female mortality. For modeling the male mortality, the advantage of directly modeling and forecasting the truncated series intensifies. By comparing the mean errors of the 19 countries, Table 2 shows that the most accurate forecasting method is the Plat model for providing best estimates of the female and male mortality forecasts. The Lee–Carter model with Poisson errors produces smaller MAPEs and RMSPEs than the Lee–Carter model with Gaussian errors. With the Lee–Carter model with Gaussian errors, it is advantageous to consider two components rather than only the first component.
For the short-term forecast horizon (i.e., the one-step-ahead forecast horizon), we compute the S ¯ α = 0.2 to evaluate the interval forecast accuracy. From Table 3, there is an advantage of directly modeling and forecasting the truncated series for both female and male mortality. By comparing the mean errors of the 19 countries, Table 3 shows that the most accurate forecasting method is the Plat model for providing the interval forecasts of both female and male mortality rates.
For the long-term forecast horizon (i.e., the 30-step-ahead forecast horizon), we also compute the mean of the MAPEs and RMSPEs to evaluate the point forecast accuracy. From Table 4, it is unclear from the comparison of the point forecast errors if there is an advantage of modeling and forecasting the truncated series for the female mortality. In contrast, there is an advantage of modeling and forecasting the male mortality and forecasting the whole series and then truncating the mortality forecasts. By comparing the mean errors of the 19 countries, Table 4 shows that the most accurate forecasting method is the Lee–Carter model with Poisson errors for providing best estimates of the female mortality forecasts. The most accurate forecasting method is the APC model for providing the best estimates of the male mortality forecasts. The Lee–Carter model with Poisson errors produces smaller MAPEs and RMSPEs than the Lee–Carter model with Gaussian errors.
For the long-term forecast horizon (i.e., the 30-step-ahead forecast horizon), we also compute the S ¯ α = 0.2 to evaluate the interval forecast accuracy. From Table 5, there is a slight advantage of directly modeling and forecasting the truncated series for the female mortality. For modeling the male mortality, there is an advantage of modeling and forecasting the whole series and then truncating the mortality forecasts. By comparing the mean errors of the 19 countries, Table 5 shows that the most accurate forecasting method is the Lee–Carter model with Poisson errors for providing the interval forecasts of the female mortality rates. The most accurate forecasting method is the APC model for providing the interval forecasts of the male mortality rates. The Lee–Carter model with Poisson errors produces smaller MAPEs than the Lee–Carter model with Gaussian errors.

4. Conclusions

We consider forecasting retiree mortality using two modeling strategies. On the one hand, we can first truncate the available data to retiree ages and then produce mortality forecasts. On the other hand, we can first use the available data to produce forecasts and then truncate the mortality forecasts to retiree ages. Using the empirical data from Human Mortality Database (2020), we investigate the short-term and long-term point and interval forecast accuracies. Between the two model strategies, we recommend the first strategy by truncating all available data to retiree ages and then produce short-term mortality forecasts. Our recommendations could be useful to actuaries for choosing a better modeling strategy and more accurately pricing a range of annuity products.
For the long-term mortality forecasts, we recommend the first strategy for modeling female mortality but the second strategy for modeling male mortality. It is difficult to recommend a strategy when the model and its parameters may not be optimal for the long-term forecasts. This is a disadvantage of using methods based on time-series extrapolation for long-term forecasting. Instead, an expectation approach, in which experts set a future target, could be considered, noting that this method has also had limited success in the past (Booth and Tickle 2008).
There are several ways in which the present study can be further extended, and we briefly mention two: (1) We could consider some machine learning methods to model mortality forecasts (see, e.g., Perla et al. 2020; Richman and Wüthrich 2020). (2) The results depend on the age range considered as well as the selected countries. The R code for reproducing the results can be provided upon request from the corresponding author.

Author Contributions

Investigation, H.L.S.; Methodology, H.L.S.; Software, H.L.S.; Supervision, S.H. These authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are grateful to the comments and suggestions received from three reviewers, and the scientific committee of the Living to 100 symposium in 2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Booth, Heather, and Leonie Tickle. 2008. Mortality modelling and forecasting: A review of methods. Annals of Actuarial Science 3: 3–43. [Google Scholar] [CrossRef]
  2. Booth, Heather, John Maindonald, and Len Smith. 2002. Applying Lee-Carter under conditions of variable mortality decline. Population Studies 56: 325–36. [Google Scholar] [CrossRef] [PubMed]
  3. Brouhns, Natacha, Michel Denuit, and Jeroen K. Vermunt. 2002. A Poisson log-bilinear regression approach to the construction of projected life-tables. Insurance: Mathematics and Economics 31: 373–93. [Google Scholar]
  4. 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]
  5. Cairns, Andrew J. G., David Blake, Kevin Dowd, Guy D. Coughlan, David Epstein, Alen Ong, and Igor Balevich. 2009. A quantitative comparison of stochastic mortality models using data from England and Wales and the United States. North American Actuarial Journal 13: 1–35. [Google Scholar] [CrossRef]
  6. Carter, Lawrence R., and Alexia Prskawetz. 2001. Examining Structural Shifts in Mortality Using the Lee-Carter Model. Technical Report. Rostock: Center for Demography and Ecology Information, University of Wisconsin-Madison. [Google Scholar]
  7. Crawford, Thomas, Richard de Haan, and Chad Runchey. 2008. Longevity Risk Quantification and Management: A Review of rElevant Literature. Technical Report. Schaumburg: Society of Actuaries. [Google Scholar]
  8. Debón, Ana, F. Montes, and Francisco Puig. 2008. Modelling and forecasting mortality in Spain. European Journal of Operational Research 198: 624–37. [Google Scholar] [CrossRef] [Green Version]
  9. Gneiting, Tilmann, and Adrian E. Raftery. 2007. Strictly proper scoring rules, prediction and estimation. Journal of the American statistical Association 102: 359–78. [Google Scholar] [CrossRef]
  10. Guillen, Montserrat, and Antoni Vidiella-i-Anguera. 2005. Forecasting Spanish natural life expectancy. Risk Analysis: An International Journal 25: 1161–70. [Google Scholar] [CrossRef] [PubMed]
  11. Human Mortality Database. 2020. University of California, Berkeley (USA), and Max Planck Institute for Demographic Research (Germany). Available online: www.mortality.org or www.humanmortality.de (accessed on 19 December 2019).
  12. Koissi, Marie-Claire, Arnold F. Shapiro, and Göran Högnäs. 2006. Evaluating and extending the Lee-Carter model for mortality forecasting: Bootstrap confidence interval. Insurance: Mathematics and Economics 38: 1–20. [Google Scholar] [CrossRef]
  13. Lee, Ronald D., and Lawrence R. Carter. 1992. Modeling and forecasting U.S. mortality. Journal of the American Statistical Association 87: 659–71. [Google Scholar] [CrossRef]
  14. Perla, Francesca, Ronald Richman, Salvatore Scognamiglio, and Mario V. Wuthrich. 2020. Time-Series Forecasting of Mortality Rates Using Deep Learning. Technical Report. Amsterdam: SSRN. [Google Scholar]
  15. Pitacco, Ermanno, Michel Denuit, Steven Haberman, and Annamaria Olivieri. 2009. Modelling Longevity Dynamics for Pensions and Annuity Business. Oxford: Oxford University Press. [Google Scholar]
  16. Plat, Richard. 2009. On stochastic mortality modeling. Insurance: Mathematics and Economics 45: 393–404. [Google Scholar]
  17. 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]
  18. 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]
  19. Richman, Ronald, and Mario V. Wüthrich. 2020. A neural network extension of the Lee-Carter model to multiple populations. Annals of Actuarial Science. in press. [Google Scholar]
  20. Shang, Han Lin. 2012. Point and interval forecasts of age-specific life expectancies: A model averaging approach. Demographic Research 27: 593–644. [Google Scholar] [CrossRef] [Green Version]
  21. Shang, Han Lin, and Steven Haberman. 2017. Grouped multivariate and functional time series forecasting: An application to annuity pricing. Insurance: Mathematics and Economics 75: 166–79. [Google Scholar]
  22. Shang, Han Lin, and Steven Haberman. 2018. Model confidence sets and forecast combination: An application to age-specific mortality. Genus 74: 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Shang, Han Lin, Heather Booth, and Rob J. Hyndman. 2011. Point and interval forecasts of mortality rates and life expectancy: A comparison of ten principal component methods. Demographic Research 25: 173–214. [Google Scholar] [CrossRef] [Green Version]
  24. Zivot, Eric, and Jiahui Wang. 2006. Modeling Financial Time Series with S-PLUS. New York: Springer. [Google Scholar]
Table 1. The 19 countries with the initial year of 1950 and their ending year listed below.
Table 1. The 19 countries with the initial year of 1950 and their ending year listed below.
CountryAbbreviationLast YearCountryAbbreviationLast Year
AustraliaAUS2014NorwayNOR2014
BelgiumBEL2015PortugalPRT2015
CanadaCAN2011SpainSPA2016
DenmarkDEN2016SwedenSWE2016
FinlandFIN2015SwitzerlandSWI2016
FranceFRA2016ScotlandSCO2016
ItalyITA2014England & WalesEW2016
JapanJPN2016IrelandIRE2014
NetherlandNET2016United States of AmericaUSA2016
New ZealandNZ2013
Table 2. For the one-step-ahead forecast horizon h = 1 , we compute the MAPE and RMSPE for each country and each model. The most accurate model for each country is highlighted in bold. For each model, we consider modeling the data with either a partial age range (termed as Partial) or a full age range (termed as Full).
Table 2. For the one-step-ahead forecast horizon h = 1 , we compute the MAPE and RMSPE for each country and each model. The most accurate model for each country is highlighted in bold. For each model, we consider modeling the data with either a partial age range (termed as Partial) or a full age range (termed as Full).
ErrorCountryLC (Poisson)LC (Gaussian)LC 2 (Gaussian)APCPlat
FullPartialFullPartialFullPartialFullPartialFullPartial
MAPEFemale
AUS4.454.386.014.974.694.427.326.384.834.66
BEL4.924.906.085.015.845.047.856.244.444.74
CAN3.293.223.943.573.743.034.984.293.202.73
DEN7.057.158.256.907.746.367.587.155.314.79
FIN6.436.1311.566.808.266.809.317.945.635.73
FRA4.284.474.974.473.963.807.676.164.223.65
ITA3.253.158.363.713.123.176.495.563.683.43
JPN5.996.0223.087.033.843.236.005.024.043.31
NET3.623.614.583.674.393.657.536.294.193.78
NZ7.707.758.657.908.337.8510.029.387.757.74
NOR5.255.286.925.546.415.418.187.585.665.56
PRT6.065.4211.685.847.825.788.757.306.136.58
SPA5.155.0412.995.736.034.197.926.244.805.21
SWE4.034.034.314.094.554.118.197.004.914.31
SWI4.694.725.994.655.894.848.947.395.285.03
SCO6.796.838.786.557.626.497.216.975.485.15
EW4.794.726.274.495.203.516.325.803.923.15
IRE8.457.3315.018.279.878.019.529.529.029.19
USA3.333.223.843.203.812.374.504.283.432.13
Mean5.245.128.495.395.854.857.596.665.054.78
Male
AUS5.394.847.855.395.745.587.715.985.155.27
BEL7.866.1110.539.748.437.436.295.445.314.86
CAN5.544.287.585.055.854.765.194.254.123.31
DEN8.196.6710.2910.517.778.337.776.535.785.52
FIN8.127.4911.848.5710.028.438.117.106.616.59
FRA3.713.174.423.964.113.667.546.024.163.32
ITA5.854.2111.295.964.633.735.164.963.742.87
JPN3.953.769.584.313.784.134.924.613.593.27
NET8.025.5812.657.757.235.375.834.795.244.03
NZ9.218.1712.6310.869.8710.909.788.798.128.38
NOR8.666.3713.616.847.216.857.426.596.076.02
PRT6.125.578.707.268.466.307.446.905.855.88
SPA4.393.888.485.447.104.386.525.804.114.29
SWE5.054.519.825.386.015.288.076.024.534.08
SWI5.625.376.515.766.865.997.716.555.275.30
SCO7.066.2010.638.568.427.517.516.716.295.80
EW4.784.048.525.155.994.286.734.933.493.43
IRE13.249.3417.2812.4312.5910.5510.079.849.869.71
USA4.653.755.904.335.623.315.794.263.712.42
Mean6.605.449.907.017.146.147.146.115.264.97
RMSPEFemale
AUS5.845.777.626.496.115.828.867.876.206.19
BEL6.406.437.556.397.376.479.467.645.746.30
CAN4.174.114.964.594.743.986.065.313.963.50
DEN9.369.7010.059.239.558.419.328.836.676.10
FIN8.518.3214.548.8510.889.1411.6910.127.468.05
FRA5.575.876.225.515.034.859.157.535.334.94
ITA4.144.019.804.624.054.197.906.924.754.59
JPN7.137.5425.038.264.784.077.326.305.004.51
NET4.624.605.744.735.504.688.867.555.234.88
NZ10.4110.4411.4810.5511.1610.6712.7512.1310.2910.49
NOR6.856.928.637.208.137.1310.159.437.237.35
PRT7.857.0313.557.479.857.5110.939.168.288.71
SPA6.536.5114.707.187.855.269.597.606.226.86
SWE5.275.265.545.315.955.389.758.486.135.56
SWI5.996.037.555.917.486.1210.638.986.636.54
SCO8.968.9911.578.589.868.569.078.897.116.77
EW6.336.258.295.946.744.427.697.154.863.96
IRE11.109.6318.7510.8212.8810.5812.2212.4311.7612.26
USA4.194.134.794.134.733.275.915.654.292.77
Mean6.806.7110.346.947.516.349.338.316.486.33
Male
AUS7.216.789.727.297.427.249.677.837.127.48
BEL9.988.4713.1413.5010.4810.028.867.837.397.19
CAN6.805.489.356.227.376.026.775.675.194.42
DEN11.059.2313.5014.3310.4911.5910.208.757.977.86
FIN13.1613.0715.2712.4915.4513.3911.8410.8810.5611.04
FRA4.724.115.605.015.124.759.217.685.414.45
ITA7.175.3513.727.286.314.846.836.504.884.04
JPN4.914.7810.985.454.775.176.205.934.454.38
NET9.687.0515.479.828.966.807.616.365.605.54
NZ13.3612.0816.8615.0413.5914.7212.9411.9011.7212.27
NOR11.328.9817.029.459.829.739.898.988.428.62
PRT8.117.5910.559.5910.568.609.729.078.068.46
SPA5.705.319.977.019.175.738.237.715.475.64
SWE6.566.0411.716.997.757.049.727.665.895.60
SWI7.567.308.507.779.248.1910.028.607.167.34
SCO9.218.4312.9810.7311.1010.2810.049.058.507.97
EW5.935.1610.386.347.425.418.396.374.504.59
IRE19.0514.7823.1717.4518.2415.4415.4914.9015.1915.23
USA5.594.826.955.306.664.087.876.144.783.26
Mean8.797.6212.369.329.478.379.458.317.287.13
Table 3. For the one-step-ahead forecast horizon h = 1 , we compute the mean interval score ( × 100 ) for each country and each model.
Table 3. For the one-step-ahead forecast horizon h = 1 , we compute the mean interval score ( × 100 ) for each country and each model.
CountryLC (Poisson)LC (Gaussian)LC 2 (Gaussian)APCPlat
FullPartialFullPartialFullPartialFullPartialFullPartial
Female
AUS3.513.442.702.682.592.786.605.262.702.53
BEL3.583.283.343.113.283.196.094.682.612.64
CAN2.912.612.071.891.851.855.174.121.581.28
DEN3.813.694.624.274.684.548.467.473.412.60
FIN5.024.936.175.375.706.198.737.843.783.42
FRA2.071.851.991.931.861.874.903.842.021.75
ITA2.331.893.212.041.992.244.383.952.102.14
JPN3.222.137.242.741.821.702.963.021.561.40
NET3.203.082.762.432.712.527.585.992.492.22
NZ5.875.754.864.684.905.639.188.195.014.63
NOR4.364.293.383.203.483.236.926.193.232.94
PRT5.344.335.544.214.584.205.574.694.243.92
SPA4.003.323.852.592.912.114.243.392.622.69
SWE3.002.932.832.762.943.188.196.442.462.14
SWI3.093.024.304.144.344.577.125.112.622.64
SCO4.794.624.013.804.124.457.677.514.123.40
EW2.522.351.961.841.861.796.495.852.141.80
IRE6.575.746.544.935.325.286.546.425.725.55
USA3.182.922.031.781.911.275.895.512.671.19
Mean3.813.483.863.183.313.296.465.553.002.68
Male
AUS5.905.534.754.444.354.5711.277.634.604.57
BEL8.797.717.738.246.877.819.878.855.995.49
CAN5.314.783.863.393.623.318.686.423.552.65
DEN8.237.047.488.337.087.9711.809.966.375.49
FIN10.369.959.138.639.659.1513.6712.305.207.39
FRA3.042.662.912.852.993.0510.628.203.322.73
ITA4.753.454.163.443.053.336.686.353.153.04
JPN3.322.524.973.623.313.854.724.762.382.22
NET7.466.265.855.125.254.358.306.695.203.79
NZ10.519.459.539.118.479.3914.3712.248.627.90
NOR8.367.576.955.825.926.1610.219.116.395.88
PRT7.456.407.326.977.567.549.249.165.875.76
SPA5.354.324.543.804.403.457.627.213.483.37
SWE5.365.055.214.444.834.6212.019.024.113.65
SWI6.236.006.766.637.017.0011.778.965.154.68
SCO8.067.357.967.378.058.1411.0210.036.735.33
EW4.293.714.033.263.633.199.376.862.812.79
IRE15.3712.2711.7810.1010.4310.1313.8413.3211.7811.12
USA5.264.753.252.983.131.9111.198.334.121.85
Mean7.026.156.225.715.775.7310.338.715.204.72
Table 4. For the 30-step-ahead forecast horizon h = 30 , we compute the MAPE and RMSPE for each country and each model. The most accurate model for each country is highlighted in bold. For each model, we consider modeling the data with either a partial age range (termed as Partial) or a full age range (termed as Full).
Table 4. For the 30-step-ahead forecast horizon h = 30 , we compute the MAPE and RMSPE for each country and each model. The most accurate model for each country is highlighted in bold. For each model, we consider modeling the data with either a partial age range (termed as Partial) or a full age range (termed as Full).
ErrorCountryLC (Poisson)LC (Gaussian)LC 2 (Gaussian)APCPlat
FullPartialFullPartialFullPartialFullPartialFullPartial
MAPEFemale
AUS18.7820.0519.0920.4618.6221.5815.2026.1120.5326.99
BEL15.1915.5017.0215.9920.2216.6115.7721.6526.2729.78
CAN13.2112.7212.569.939.0710.6817.9616.2922.7610.32
DEN11.9112.3611.4410.7114.6411.4123.5922.0334.1519.92
FIN16.3518.2816.7918.3721.5920.9223.3827.2726.1029.23
FRA16.8117.3516.7417.3020.0018.2313.8018.5119.0230.02
ITA18.5521.1726.0024.2224.1622.9015.0632.5423.4148.43
JPN22.8323.8940.7127.4024.6523.0352.4819.8252.9630.09
NET15.0414.4220.2514.4111.7012.8214.9116.9831.1915.17
NZ34.6135.9138.2340.1138.5636.2128.8441.5935.9444.18
NOR10.179.5212.3211.3910.4410.4414.2219.6835.7312.53
PRT41.8934.6153.0140.1552.4145.1617.4046.5447.4654.60
SPA25.0924.3352.6627.0327.0028.8535.9627.8235.9240.34
SWE14.5413.7315.8512.5013.5012.5516.0216.2624.5512.05
SWI11.2710.9310.9910.8312.2512.3017.5619.7637.2518.94
SCO16.9118.2917.5519.7219.7520.5610.7725.8631.5320.71
EW17.2419.4155.2519.2517.4620.119.5734.1225.0334.45
IRE38.9338.0968.2043.5945.5841.7719.4744.0938.1941.49
USA8.907.5111.135.987.755.5516.6513.2412.567.83
Mean19.3819.3727.1520.4921.5420.6219.9325.8030.5627.74
Male
AUS50.1352.4756.6855.4149.3355.9538.2152.5453.8172.66
BEL58.3454.8865.6992.4758.8385.9731.6346.4287.2446.11
CAN45.4642.2647.8544.4647.8044.7833.7043.2679.9449.49
DEN63.4563.0163.3667.6766.8767.4466.0157.7173.5756.30
FIN43.3948.1144.3496.5847.7851.7720.8531.2975.9652.26
FRA32.8133.8335.8134.7737.0035.3639.0935.8660.4939.45
ITA64.4162.3570.3867.2253.6362.0332.1161.0581.6460.60
JPN12.3012.6820.9812.7912.9112.8239.5320.9432.1215.61
NET79.6876.1478.6779.2965.1971.0845.6251.8979.8847.61
NZ99.3282.86112.42117.61104.32114.2863.1076.3992.0594.38
NOR87.5073.9788.6271.4365.6467.3545.0966.8570.0060.19
PRT42.7938.1149.0183.7958.2449.5918.0135.2932.2647.62
SPA25.4323.1358.4529.8527.1530.8415.2125.0226.2928.47
SWE51.8947.5755.0157.8953.0551.9145.4643.2876.6542.27
SWI36.0037.4641.7342.2739.0145.4830.6446.4250.8255.93
SCO52.6655.6355.8692.8661.3686.4138.5157.0265.6354.61
EW52.4455.9652.8160.0350.0655.5235.4248.1647.2654.41
IRE110.2091.02112.71111.84111.82104.5243.1097.2663.0793.34
USA24.0724.9926.0726.8126.1228.1321.5623.2070.2131.07
Mean54.3351.3959.8165.5354.5359.0136.9948.4164.1552.76
RMSPEFemale
AUS21.7423.1122.1723.5121.8224.6019.0928.3424.4130.76
BEL18.2518.9220.8219.4424.7020.3818.6724.8931.5435.64
CAN15.3714.9514.5312.1811.3113.0122.7621.3828.3714.13
DEN14.0814.4413.3313.4616.6213.5025.5825.1140.8423.08
FIN19.8222.3720.1422.6726.1825.8027.3730.7833.0235.66
FRA19.0119.6019.2720.1723.4721.1817.9721.6024.5534.67
ITA21.1823.4528.2526.3626.4425.1718.2836.1528.1957.57
JPN25.7326.9147.7230.0627.2725.7153.7527.3755.5033.61
NET16.3315.7721.1315.6713.9314.4320.1521.3237.4317.97
NZ44.3745.8446.5549.7146.6845.7033.1048.7942.2950.98
NOR12.6512.4514.5313.9012.8013.1317.9422.6141.0814.97
PRT48.0940.4061.1746.2860.6652.0620.9654.0952.8062.86
SPA30.7329.5159.1733.9633.8936.1638.2430.7442.1045.40
SWE16.2315.5017.5414.6117.1814.7322.0522.2131.0413.67
SWI13.9713.5213.6013.4215.1717.2723.5124.3842.6221.26
SCO22.0024.5120.5625.9522.1926.9814.0529.2935.8224.81
EW21.8424.0861.1923.8921.6224.4411.7838.7929.5939.16
IRE44.0643.5973.2748.9250.7547.7224.3357.1442.5547.34
USA10.108.7012.137.029.116.7718.4516.3215.298.77
Mean22.9223.0330.9024.2725.3624.6723.5830.6035.7432.23
Male
AUS58.8261.4866.4364.8057.7366.1245.9264.3564.1982.67
BEL72.3563.8578.04103.4071.5295.8637.4554.3690.5653.55
CAN57.5752.4261.4055.5460.6055.8140.5152.4084.4957.41
DEN81.1180.8978.7175.3882.5382.0377.0066.7591.5268.67
FIN51.7157.5252.22111.4955.4061.7823.3936.3078.5160.90
FRA38.1538.8441.5340.3840.5840.9143.6840.6468.4844.68
ITA80.7077.3287.7984.1567.3678.1536.0171.6285.6575.05
JPN15.4616.4825.7715.2516.6615.9140.6226.8237.0618.33
NET102.5895.23100.24102.7681.6989.7453.6460.8187.9764.18
NZ112.8398.95128.69136.49117.87132.7574.8291.95102.96107.61
NOR108.6285.81110.9682.2580.6279.4350.6077.0479.8576.25
PRT49.9344.0356.6395.1764.8955.5720.4041.5238.6755.07
SPA32.8230.0463.1333.1432.1534.6019.1727.5130.2333.06
SWE62.1956.6566.3768.9862.4162.3952.8651.0184.6350.20
SWI43.7145.2748.8450.1745.9353.2236.9155.6661.5260.67
SCO64.2566.3768.29105.0272.1897.9044.1267.6171.8166.18
EW62.4565.9863.7271.7458.8065.5642.0758.1156.8963.73
IRE122.60104.66125.27128.85124.03118.3449.25116.0280.35105.59
USA28.4429.3330.2631.4130.4132.0525.8526.0972.7034.40
Mean65.5961.6471.2876.6564.3969.3742.8657.1972.0062.01
Table 5. For the 30-step-ahead forecast horizon h = 30 , we compute the mean interval score ( × 100 ) for each country and each model.
Table 5. For the 30-step-ahead forecast horizon h = 30 , we compute the mean interval score ( × 100 ) for each country and each model.
CountryLC (Poisson)LC (Gaussian)LC 2 (Gaussian)APCPlat
FullPartialFullPartialFullPartialFullPartialFullPartial
Female
AUS8.257.209.308.947.758.8013.4210.4612.319.29
BEL5.994.996.856.187.8310.108.2411.8019.5312.75
CAN12.6410.767.548.095.117.0923.1521.9323.714.47
DEN4.815.235.375.076.3518.669.728.8528.918.85
FIN12.809.579.239.2912.38162.2325.8520.5727.4916.37
FRA5.394.515.497.335.455.6911.3711.8214.189.99
ITA8.814.5412.6210.216.168.466.167.097.2513.39
JPN6.734.3012.4710.435.406.0627.2420.4239.637.77
NET9.949.0910.619.666.968.1122.0623.5942.298.94
NZ8.758.819.6110.3916.1582.217.788.999.1313.52
NOR6.817.825.396.415.217.6814.4713.5544.897.80
PRT20.6310.5926.5813.2123.6418.538.518.4953.8522.47
SPA16.7013.3518.0816.1212.879.2915.9512.6832.888.74
SWE9.668.889.477.856.1613.7125.9726.3235.746.55
SWI3.673.395.305.086.5842.7422.6118.8431.9414.23
SCO8.366.089.657.6510.3730.247.136.6426.0211.29
EW3.844.0411.424.454.816.205.097.4413.978.55
IRE17.2512.3920.5716.6418.8921.438.069.8115.2620.33
USA4.112.503.111.802.522.7015.6613.934.474.05
Mean9.227.2710.468.678.9824.7314.6513.8525.4511.02
Male
AUS15.9712.7818.9115.6913.6223.8510.1510.1816.5121.01
BEL22.2818.1829.4347.28197.68274.409.2512.3039.3517.42
CAN12.9911.4312.5311.8012.9413.6313.5013.6335.6414.97
DEN19.7419.2821.7230.66441.92463.8122.0916.4023.1715.61
FIN16.4516.9815.5433.3626.1143.948.548.8034.8223.53
FRA10.588.2012.8512.7414.7114.9310.699.3816.7014.70
ITA17.3416.2419.6417.4712.5617.477.8915.9227.7617.58
JPN3.583.698.704.737.2723.3730.4125.3238.0210.07
NET23.7421.6324.0420.7818.9123.4314.9610.4228.7815.71
NZ45.0026.1354.3941.5052.5995.2918.2316.6128.1042.92
NOR26.4525.7925.8627.4016.8923.4914.7713.6322.2016.68
PRT19.8713.6824.0137.2630.92118.328.5012.3212.2622.79
SPA27.4414.8035.0519.8815.8114.268.9210.8011.0910.42
SWE17.7113.4618.5317.6019.8720.0114.9611.3425.1910.65
SWI9.379.0514.4813.3228.5728.639.0411.1914.0022.49
SCO20.8519.0121.6143.0384.5196.6810.2212.7725.0718.27
EW15.2414.2016.6821.0617.3517.397.4811.1714.0817.50
IRE70.2637.2673.4048.7969.77109.0112.7125.2128.2236.51
USA10.959.6510.609.2610.4710.3412.5911.0942.6411.69
Mean21.3616.3924.1024.9357.5075.3812.8913.6025.4518.97

Share and Cite

MDPI and ACS Style

Shang, H.L.; Haberman, S. Retiree Mortality Forecasting: A Partial Age-Range or a Full Age-Range Model? Risks 2020, 8, 69. https://0-doi-org.brum.beds.ac.uk/10.3390/risks8030069

AMA Style

Shang HL, Haberman S. Retiree Mortality Forecasting: A Partial Age-Range or a Full Age-Range Model? Risks. 2020; 8(3):69. https://0-doi-org.brum.beds.ac.uk/10.3390/risks8030069

Chicago/Turabian Style

Shang, Han Lin, and Steven Haberman. 2020. "Retiree Mortality Forecasting: A Partial Age-Range or a Full Age-Range Model?" Risks 8, no. 3: 69. https://0-doi-org.brum.beds.ac.uk/10.3390/risks8030069

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