Future life expectancy in Europe taking into account the impact of smoking, obesity, and alcohol
Abstract
Introduction: In Europe, women can expect to live on average 82 years and men 75 years. Forecasting how life expectancy will develop in the future is essential for society. Most forecasts rely on a mechanical extrapolation of past mortality trends, which leads to unreliable outcomes because of temporal fluctuations in the past trends due to lifestyle ‘epidemics’. Methods: We project life expectancy for 18 European countries by taking into account the impact of smoking, obesity, and alcohol on mortality, and the mortality experiences of forerunner populations. Results: We project that life expectancy in these 18 countries will increase from, on average, 83.4 years for women and 78.3 years for men in 2014 to 92.8 years for women and 90.5 years for men in 2065. Compared to others (Lee–Carter, Eurostat, United Nations), we project higher future life expectancy values and more realistic differences between countries and sexes. Conclusions: Our results imply longer individual lifespans, and more elderly in society. Funding: Netherlands Organisation for Scientific Research (NWO) (grant no. 452-13-001).
eLife digest
On average, in Europe, men can currently expect to live till the age of 75 and women until they are 82. But what will their lifespans be in the next decades? Reliable answers to this question are essential to help governments plan for future health care and social security costs.
While medical improvements are likely to further extend lifespans, lifestyle factors can result in temporal distortions of this trend. Yet, most estimates of future life expectancy fail to consider changing lifestyles, as they only use past mortality trends in their calculations. This can make these projections unreliable: for example, increases in smoking rates among Northern and Western European men led to stagnating male life expectancies in the 1950s and 1960s, but these picked up again after smoking declined. The same pattern is showing for women, except it is lagging as they took up smoking later than men. Based simply on the extrapolation of past mortality trends, current projection models fail to consider the past and predicted modifications of life expectancy trends prompted by changing rates of health behaviours – such as increases followed by (anticipated) declines in alcohol consumption and obesity rates, similar to what was observed with smoking.
To produce a more reliable forecast, Janssen et al. incorporated trends in smoking, obesity, and alcohol use into life expectancy projections for 18 European countries. The predictions suggest that life expectancy for women in these countries will increase from 83.4 years in 2014 to 92.8 years in 2065. For men, it will also go up, from 78.3 to 90.5 years.
In the future, this integrative approach may help to track the effects of health-behaviour related prevention policies on life expectancy, and allow scientists to account for changes caused by the COVID-19 pandemic. In the meantime, these estimates are higher than those obtained using more traditional methods; they suggest that communities should start to adjust to the possibility of longer individual lifespans, and of larger numbers of elderly people in society.
Introduction
Mortality projections are essential for forecasting how long people will live on average, for predicting the future extent of population ageing, for determining the sustainability of pension schemes and social security systems, for setting life insurance premiums, and for helping governments plan for the demand for services such as healthcare (European Commission, 2009). Recently, the societal and academic relevance of obtaining reliable and robust mortality forecasts has increased due to the revived debate on the existence and proximity of a limit to human life expectancy (Fries, 1980; Manton et al., 1991; Couzin, 2005; de Beer et al., 2017; Olshansky and Carnes, 2019; Vaupel et al., 2021); the linkage between life expectancy developments and retirement age in many European countries (Carone et al., 2016; OECD, 2016); the stagnation of increases in life expectancy in the United States, the United Kingdom, and other European countries since 2011 (Murphy et al., 2019; Raleigh, 2019); and the recent COVID-19 pandemic.
The growing relevance of obtaining reliable mortality forecasts has led to important advances in mortality forecasting by professionals and researchers from – among others – demography, actuarial sciences, and public health (Janssen, 2018). In the past, mortality forecasts were often based on expectations fuelled by the debate on the limit to life expectancy (Oeppen and Vaupel, 2002; Booth and Tickle, 2008). The majority of the currently existing mortality forecasting methods, including those used by statistical offices in Europe, are mainly based on extrapolations of past trends in age-specific mortality (Booth and Tickle, 2008; Stoeldraijer et al., 2013a). These extrapolative methods rely on the regularities observed in mortality trends over time and in the age patterns (Booth and Tickle, 2008). An extrapolative approach to mortality forecasting is considered more objective, easier to apply, and more likely to result in reliable forecasts than the previously employed expectation approaches and explanation approaches (mortality forecasting by cause of death or with an explanatory model) (Booth and Tickle, 2008). The stochastic Lee–Carter mortality projection methodology (Lee and Carter, 1992) has long been the benchmark extrapolative mortality forecasting method (see Janssen, 2018).
It is, however, increasingly being acknowledged that mechanically extrapolating past national mortality trends can result in unreliable outcomes (e.g., Janssen and Kunst, 2007; Janssen et al., 2013; Bongaarts, 2014; Stoeldraijer, 2019; Li and Raftery, 2021; Blakely, 2018). First, when past trends in mortality are non-linear, the projection outcome will depend to a large extent on the period on which the extrapolation is based (Janssen and Kunst, 2007; Janssen et al., 2013; Stoeldraijer, 2019). Second, extrapolative mortality forecasts made for individual populations often project unrealistic crossovers or large future differences in life expectancy between populations (Li and Lee, 2005; Janssen et al., 2013; Hyndman et al., 2013). For example, extrapolating the large mortality declines observed in Eastern European countries – which have historically had relatively low life expectancy, but have experienced improvements since 2005 (Trias-Llimós et al., 2018) – would result in long-term future life expectancy values being higher in Eastern than in Western European countries with historically high life expectancy values and more moderate mortality declines.
We argue that to obtain a reliable and robust mortality forecast, it is essential to distinguish between (1) the general and gradual long-term mortality decline due to socio-economic and medical progress that can be extrapolated into the future, thereby taking into account the mortality experiences of other countries, and (2) remaining factors that cause deviations from this general mortality decline as well as country and sex differences in this trend.
While life expectancy has increased overall as a result of socio-economic and medical progress (Omran, 1998; Oeppen and Vaupel, 2002; Mackenbach, 2013), past trends in mortality and life expectancy were characterised by periods of stagnation or even decline and exhibited large differences between countries and sexes (e.g., Vallin and Meslé, 2004; Leon, 2011). Lifestyle factors, particularly smoking, alcohol abuse, and obesity, contributed greatly to these deviations and differences (e.g., Lhachimi et al., 2016; WHO, 2018a; Janssen et al., 2021). For example, the decelerating life expectancy improvements among men in many north-western European countries in the 1950s and 1960s, and, in later decades, in other European countries and among women, can be largely explained by the smoking epidemic (Rostron and Wilmoth, 2011; Janssen et al., 2015; Lindahl-Jacobsen et al., 2016). The pattern of a rapid increase and a subsequent decline in smoking prevalence occurred first among men in Anglo-Saxon countries, and was followed later in other countries and by women (Lopez et al., 1994). These developments resulted in similar patterns in smoking-attributable mortality about 30–40 years later (Lopez et al., 1994; Janssen, 2020). The large increases in alcohol consumption and alcohol-attributable mortality between 1990 and 2005 in Eastern Europe contributed to the stalling of life expectancy improvements in Eastern Europe in this period (Trias-Llimós et al., 2018). The more favourable alcohol-related trends in Eastern Europe since 2005 have contributed to the convergence in life expectancy levels between Eastern and Western Europe. The stagnation in life expectancy improvements since approximately 2011 in the United States and the United Kingdom, but also in other selected European countries (Murphy et al., 2019), can be partly attributed to the high obesity prevalence and obesity-attributable mortality in these countries (Raleigh, 2019; Janssen et al., 2020a; Janssen et al., 2020b ) and to the increasing impact on life expectancy levels of past large increases in obesity prevalence (Vidra et al., 2019).
Previous mortality forecasts that incorporated lifestyle information (Bongaarts, 2006; Pampel, 2005; Wang and Preston, 2009; King and Soneji, 2011; Preston et al., 2014; French and O'Hare, 2013; Foreman et al., 2018; Janssen and Kunst, 2007; Janssen et al., 2013; Li and Raftery, 2021) generally included information on only one lifestyle factor, usually smoking. Moreover, these forecasts employed very different techniques and ignored the mortality experiences of other countries. The previous so-called coherent or multi-population forecasting methods that take into account the mortality experiences of other countries (Li and Lee, 2005; Hyndman et al., 2013; Enchev et al., 2017; Eurostat, 2020b; Bergeron-Boucher et al., 2018) have, despite their growing popularity, so far been applied almost exclusively to all-cause mortality. The two previous projections that incorporated both elements simultaneously only included the effect of smoking (Janssen et al., 2013; Li and Raftery, 2021). The approach developed for this purpose by Janssen et al., 2013 consists of the combination of a coherent projection of non-smoking-related mortality with an age-period-cohort projection of smoking-attributable mortality. This approach was applied to the Netherlands (Janssen et al., 2013) and was adopted by Statistics Netherlands as part of their official population forecast (Stoeldraijer et al., 2013b). Li and Raftery, 2021 applied an amended version of this approach to 69 countries.
Our objective is to project future life expectancy in 18 European countries, while simultaneously taking into account the time-varying impact of smoking, obesity, and alcohol on mortality, and the mortality experiences of forerunner populations.
Materials and methods
Our approach to mortality forecasting relied on an analysis of past mortality trends and their determinants. In line with the evidence gathered, we distinguished between (1) the general and gradual long-term mortality decline not affected by the three lifestyle factors that could be extrapolated into the future, while taking into account the mortality experiences of other countries, and (2) deviations from and differences in this general mortality decline caused predominantly by the time-varying impact of smoking, obesity, and alcohol on mortality, which required the use of more advanced projection techniques.
The projection involved four steps. First, we determined the long-term decline in mortality and life expectancy without the combined effect of smoking, obesity, and alcohol. To this end, we used existing recent age- and sex-specific estimates of mortality after excluding the combined impact smoking, obesity, and alcohol (Janssen et al., 2021), which we refer to as non-lifestyle-attributable mortality. Second, we projected this long-term underlying mortality decline into the future, while taking into account the mortality experiences of other populations. Specifically, we performed a Li–Lee coherent projection (Li and Lee, 2005) of non-lifestyle-attributable mortality in which we regarded the high life expectancy values and more favourable long-term trends among women in France, Spain, and Italy as the values and trends towards which the other populations will converge, since we regard these three populations as forerunners. Third, we obtained future estimates of mortality attributable to smoking, obesity, and alcohol, which we refer to as lifestyle-attributable mortality. We did so by utilizing recently published projections of smoking-, obesity-, and alcohol-attributable mortality that were both data and theory driven (Janssen et al., 2020c; Janssen et al., 2020b; Janssen et al., 2020d). Fourth, we combined the projections of non-lifestyle-attributable mortality and lifestyle-attributable mortality by extending the approach that Janssen et al. developed for combining the separate projections of non-smoking- and smoking-attributable mortality (Janssen et al., 2013).
We applied our projection approach to the national populations of 18 European countries (see Table 1 for the included countries) using country, sex, and age-specific all-cause mortality and exposure data for the 1990–2014 period from the Human Mortality Database (HMD and University of California, Berkeley (USA), and Max Planck Institute for Demographic Research (Germany), 2019), and similar data for lifestyle-attributable mortality and non-lifestyle-attributable mortality from a recent study (Janssen et al., 2021). We were restricted in both the starting year and the end year because lifestyle-attributable mortality data was available only for this period.
Comparison of observed gains in life expectancy at birth (e0) between 1990 and 2014 with the gains when the effects of smoking, obesity, and alcohol are removed ( = for non-lifestyle-attributable mortality), 18 European countries, by country and sex.
Country | Gain e0 1990–2014 | |||
---|---|---|---|---|
Observed | Non-lifestyle-attributable mortality | |||
Men | Women | Men | Women | |
Austria | 6.69 | 4.87 | 5.05 | 5.29 |
Belgium | 5.88 | 4.18 | 4.22 | 4.90 |
Czech Republic | 8.19 | 6.33 | 5.50 | 6.62 |
Denmark | 6.54 | 4.94 | 4.98 | 5.01 |
Finland | 7.19 | 4.96 | 5.19 | 5.56 |
France | 6.56 | 4.47 | 4.95 | 4.84 |
Germany | 6.52 | 4.92 | 4.99 | 5.45 |
Greece | 3.80 | 4.18 | 3.44 | 4.23 |
Hungary | 7.12 | 5.47 | 5.61 | 6.16 |
Ireland | 7.03 | 5.52 | 5.97 | 5.90 |
Italy | 6.92 | 4.91 | 4.67 | 4.97 |
Netherlands | 6.05 | 3.20 | 4.20 | 4.50 |
Norway | 6.58 | 4.29 | 6.10 | 5.21 |
Poland | 7.41 | 6.13 | 6.56 | 7.04 |
Slovenia | 8.20 | 6.00 | 5.62 | 6.39 |
Sweden | 5.54 | 3.66 | 5.09 | 4.35 |
Switzerland | 6.98 | 4.40 | 5.25 | 4.76 |
United Kingdom | 6.41 | 4.49 | 4.89 | 4.67 |
Average | 6.64 | 4.83 | 5.13 | 5.33 |
Min | 3.80 | 3.20 | 3.44 | 4.23 |
Max | 8.20 | 6.33 | 6.56 | 7.04 |
Variance | 0.92 | 0.66 | 0.51 | 0.62 |
Annual change | 0.28 | 0.20 | 0.21 | 0.22 |
Forerunners* | 4.81 | 4.81 | 5.00 | 5.00 |
-
*French, Italian, and Spanish women (unweighted average). For Spanish women, the gain in e0 (1990–2014) was 5.05 years for all-cause mortality and 5.20 years for non-lifestyle-attributable mortality.
-
Table 1—source data 1
Data behind Table 1.
- https://cdn.elifesciences.org/articles/66590/elife-66590-table1-data1-v2.xlsx
We obtained estimates of life expectancy at birth (e0), including 95% projection intervals up to 2065, by applying standard life table techniques to the projected mortality rates for ages 0–130 (Preston et al., 2001). We focussed on life expectancy at birth because it is a very common summary measure of health, and the most common output measure of mortality forecasts. We have chosen a relatively large projection horizon (2015–2065) given the comparatively short historical time series available (1990–2014) to illustrate that our approach is able to generate reliable outcomes for the long-term future.
We compared our projection outcomes with the outcomes of the benchmark Lee–Carter extrapolation (Lee and Carter, 1992) applied to all-cause mortality and examined the separate effects of incorporating lifestyle factors, and of including the mortality experiences of other countries. We report the outcomes of additional projections and comparisons in the Results section ‘The differences explained’. We also compared our outcomes to the official forecasts by Eurostat and United Nations (see the Discussion section ‘Comparison with other projections’).
More detailed information on the data and methods can be found in Appendix 1.
Results
Underlying long-term increase in life expectancy
Over the 1990–2014 period, life expectancy at birth (e0) increased, on average, across the 18 European countries studied by 6.6 years for men (from 71.8 to 78.3 years) and by 4.8 years for women (from 78.6 to 83.4 years) (see Table 1). These values translate into a yearly increase in e0 of 0.28 years for men and 0.20 years for women. This difference between men and women can largely be attributed to lifestyle factors. Without the impact of smoking, obesity, and alcohol, the increase in e0 over the 1990–2014 period was more similar between men (5.1 years; 0.21 annually) and women (5.3 years; 0.22 annually). In addition, the increase in e0 was more similar between countries for non-lifestyle-attributable mortality (variance: 0.5 for men, 0.6 for women) than for all-cause mortality (variance: 0.9 for men, 0.7 for women).
The greater increase in e0 for all-cause mortality than for non-lifestyle-attributable mortality among men can be explained by the significant declines in lifestyle-attributable mortality men experienced (see Figure 1). These declines stemmed predominantly from large declines in smoking-attributable mortality (Figure 1—figure supplement 1), after a period of sharp increases. In contrast, the past increases in lifestyle-attributable mortality among women were driven by increases in all three factors (except alcohol in Eastern Europe) (Figure 1—figure supplement 2) and resulted in smaller increases in e0 for all-cause mortality than for non-lifestyle-attributable mortality.

Past and projected age-standardised lifestyle-attributable mortality fractions (LAMF) (20–84), 1990–2065, by sex, for 18 European countries.
-
Figure 1—source data 1
Data behind Figure 1.
- https://cdn.elifesciences.org/articles/66590/elife-66590-fig1-data1-v2.xlsx
Figure 2 and Figure 2—figure supplement 1 show that whereas trends in e0 since 1950 have been rather unstable, the trends in e0 without lifestyle factors can be regarded as more stable. This is indicated by (1) the more stable trend in e0 for non-smoking-related mortality for men in the 1950s and 1960s, which seems to be in line with the more recent trends in e0 for non-lifestyle-attributable mortality and (2) the close correspondence among women between the recent trends in e0 for non-lifestyle-attributable mortality and the trends in e0 for all-cause mortality in the 1950s–1960s that were probably affected only minimally by lifestyle factors.

Comparison of trends in observed life expectancy at birth (e0) ( = all-cause) with trends in e0 when the effects of smoking, obesity and/or alcohol are removed ( = non-smoking, non-obesity, non-alcohol, non-lifestyle combined), 1950–2014*, men.
*Based on the available information. Reproduced from various panels in Supplementary Figure 3a from Janssen et al., 2021, under the terms of a Creative Commons Attribution License (CC-BY 4.0; https://creativecommons.org/licenses/by/4.0/).
-
Figure 2—source data 1
Data behind Figure 2.
- https://cdn.elifesciences.org/articles/66590/elife-66590-fig2-data1-v2.xlsx
Projected future level of life expectancy
We extrapolated the more universal and more stable mortality trends that we observed for non-lifestyle-attributable mortality, while assuming that the trends for men and women in the individual countries will eventually move towards the more favourable long-term trends for women in France, Spain, and Italy. We selected these populations as the forerunner populations in terms of life expectancy in Europe because they exhibit both very high recent e0 values and very favourable long-lasting past trends in e0 (Stoeldraijer, 2019).
We added to these projections the projected levels of lifestyle-attributable mortality (see Figure 1). For men, the past declines in age-standardised lifestyle-attributable mortality fractions (LAMF) over the 1990–2014 period are projected to further decline until 2065, albeit at a different pace than in the past. Among Eastern European men in particular, the decline is projected to decelerate. For women, we project that the past increases in LAMF will (eventually) turn into declines. Figure 1—figure supplements 1 and 2 show the projections of smoking-, obesity-, and alcohol-attributable mortality (Janssen et al., 2020c; Janssen et al., 2020b; Janssen et al., 2020d). The projections follow the observed wave-shaped dynamic of the smoking epidemic (Lopez et al., 1994; Janssen, 2020) and the hypothesised wave shape of the obesity epidemic (Xu and Lam, 2018; Jaacks et al., 2019). For alcohol-attributable mortality, unrealistic future differences between countries were avoided by assuming that the current increases observed for selected countries will eventually turn into declines (Janssen et al., 2020d).
We project that in the 18 European countries studied, life expectancy at birth will increase from, on average, 78.3 years for men and 83.4 years for women in 2014 to 90.5 years for men and 92.8 years for women in 2065 (Table 2; Figure 3). This represents an average increase per year of 0.24 years for men and 0.18 years for women. The projected increase is greater for men than for women because all-cause mortality levels are currently further away from non-lifestyle-attributable mortality levels for men than for women and because lifestyle-attributable mortality is higher among men than among women. There is, therefore, more room for improvement for men than for women.

Observed and projected life expectancy at birth (including 95% projection intervals) using our projection approach, which takes into account the impact of smoking, obesity, and alcohol on past and future mortality trends, and the mortality experiences of forerunner countries, 18 European countries, by country and sex, 1990–2065.
-
Figure 3—source data 1
Data behind Figure 3.
- https://cdn.elifesciences.org/articles/66590/elife-66590-fig3-data1-v2.xlsx
Observed (2014) and projected (2065) life expectancy at birth (e0), by country and sex, according to our projection methodology, which takes into account the impact of smoking, obesity, and alcohol and the mortality experiences of forerunner countries (‘lifestyle and coherent’); the benchmark Lee–Carter extrapolative mortality projection applied to all-cause mortality (‘Lee–Carter’); and when purely accounting for smoking, obesity, and alcohol (‘adding lifestyle’).
Country | Men | Women | ||||||
---|---|---|---|---|---|---|---|---|
e0 2014 | e0 2065 | e0 2014 | e0 2065 | |||||
Lee–Carter | Adding lifestyle | Lifestyle and coherent | Lee–Carter | Adding lifestyle | Lifestyle and coherent | |||
Austria | 78.91 | 88.91 | 89.02 | 90.38 | 83.73 | 91.42 | 92.57 | 92.54 |
Belgium | 78.57 | 87.76 | 88.20 | 90.93 | 83.52 | 90.09 | 92.43 | 92.98 |
Czechia | 75.72 | 87.77 | 87.23 | 88.86 | 81.73 | 90.87 | 92.28 | 91.37 |
Denmark | 78.56 | 87.72 | 88.48 | 90.50 | 82.67 | 89.66 | 93.33 | 93.08 |
Finland | 78.13 | 88.25 | 88.30 | 90.22 | 83.85 | 90.90 | 92.65 | 92.69 |
France | 79.28 | 89.11 | 90.03 | 91.40 | 85.44 | 92.77 | 94.24 | 94.01 |
Germany | 78.43 | 88.36 | 88.67 | 90.18 | 83.35 | 90.78 | 92.60 | 92.62 |
Greece | 78.46 | 85.07 | 86.32 | 90.97 | 83.83 | 89.75 | 90.51 | 92.44 |
Hungary | 72.26 | 81.74 | 85.89 | 88.64 | 79.24 | 87.09 | 90.21 | 90.96 |
Ireland | 79.15 | 89.59 | 90.57 | 90.95 | 83.23 | 91.67 | 94.38 | 93.40 |
Italy | 80.55 | 89.25 | 89.14 | 91.43 | 85.16 | 92.30 | 93.34 | 93.53 |
Netherlands | 79.87 | 88.67 | 88.73 | 91.11 | 83.29 | 88.80 | 92.59 | 93.18 |
Norway | 80.03 | 88.86 | 90.23 | 90.87 | 84.09 | 90.72 | 93.48 | 93.35 |
Poland | 73.66 | 85.21 | 88.28 | 89.37 | 81.41 | 90.44 | 92.74 | 91.88 |
Slovenia | 77.96 | 88.68 | 88.95 | 89.91 | 83.69 | 92.20 | 93.79 | 92.60 |
Sweden | 80.35 | 88.37 | 89.18 | 90.64 | 84.05 | 89.92 | 92.27 | 93.01 |
Switzerland | 80.93 | 89.77 | 89.31 | 90.88 | 85.11 | 91.41 | 92.70 | 93.16 |
United Kingdom | 79.25 | 88.99 | 89.35 | 91.04 | 82.99 | 89.93 | 92.70 | 93.34 |
Average | 78.34 | 87.89 | 88.66 | 90.46 | 83.35 | 90.59 | 92.71 | 92.78 |
Min | 72.26 | 81.74 | 85.89 | 88.64 | 79.24 | 87.09 | 90.21 | 90.96 |
Max | 80.93 | 89.77 | 90.57 | 91.43 | 85.44 | 92.77 | 94.38 | 94.01 |
Variance | 4.98 | 3.77 | 1.41 | 0.62 | 2.02 | 1.73 | 1.07 | 0.55 |
-
Table 2—source data 1
Data behind Table 2.
- https://cdn.elifesciences.org/articles/66590/elife-66590-table2-data1-v2.xlsx
Compared to the average past increases (1990–2014) in e0 (0.28 years for men; 0.20 years for women), these projected increases are lower, particularly for men. This is because – as mentioned – the increase in e0 for men over the 1990–2014 period is actually an acceleration compared to the underlying trend in non-lifestyle-attributable mortality because of the decline in lifestyle-attributable mortality, predominantly in smoking-attributable mortality.
The highest e0 in 2065 is projected for France among women (94.0) and for Italy among men (91.4). Hungary is the country with the lowest projected e0 in 2065, at 88.6 years among men and 91.0 years among women. These countries also had in 2014, respectively, the (second) highest and the lowest values. We project the largest gains in life expectancy for Hungary and Poland, particularly among men, in line with past large increases in non-lifestyle-attributable mortality due to a process of catch-up (Table 1), and large projected declines in lifestyle-attributable mortality in these countries (Figure 1). Among women, we project large gains as well for Denmark and the United Kingdom, which can be attributed to their high current levels of lifestyle-attributable mortality, predominantly due to their high levels of smoking-attributable mortality, and the large projected (continued) declines therein (Figure 1—figure supplement 2). The smallest gains are projected for Switzerland, which could be related to the relatively small projected declines in lifestyle-attributable mortality resulting from a long projected continuation of increases in obesity prevalence.
The sex difference in e0 is projected to decline from, on average, 5.0 years in 2014 to 2.3 years in 2065 (Supplementary file 1A; Figure 3). The gap is projected to be smallest in Greece, at 1.5 years; and largest in Slovenia, at 2.7 years. Among the explanations for the projected convergence of the sex gap in e0 is that the sex differences are smaller for non-lifestyle-attributable mortality than for all-cause mortality; the sex differences in lifestyle-attributable mortality are smaller in the long run; and our assumption that the trends in non-lifestyle-attributable mortality will converge.
Comparison with the benchmark Lee–Carter mortality projection
Compared to the benchmark Lee–Carter (LC) extrapolative mortality projection, our projection results in higher future e0 in 2065 for all populations under study (Table 2). On average, the projected e0 values in 2065 are 2.6 years higher for men and 2.1 years higher for women in our projection than in the individual Lee–Carter mortality projection, which projected e0 values of 87.9 years for men and 90.6 years for women in 2065 (Table 2). For men, particularly in Hungary, the differences between our projected values and those of the LC projection are considerable, mainly because of large projected declines in lifestyle-attributable mortality. Also for Dutch women a considerable difference between our projection and the LC projection can be observed. For Dutch women not only large declines in lifestyle-attributable mortality are projected, but they also experienced a more favourable trend in non-lifestyle-attributable mortality than in all-cause mortality.
In addition, our projection results in smaller differences in future e0 between countries and sexes. That is, the LC projection projects that the variance between countries in e0 for 2065 is 3.8 years for men and 1.7 years for women, whereas in our projection, the expected variance is 0.6 years for both men and women (Table 2). In addition, Figure 3—figure supplement 1 illustrates that the LC projection, unlike our projection, can result in unlikely crossovers and divergences between countries (e.g., Greece and Hungary for men; the Netherlands and Slovenia for women).
In addition, the sex difference in e0 in 2065 is projected to be 2.7 years in the LC projection, but 2.3 years in our projection (Table 2). Furthermore, a closer look at the sex differences projected by the LC projection (Supplementary file 1A) reveals large differences in the projected sex differences between countries ranging from 0.1 to 5.4 years. In our projection, by contrast, the range is between 1.5 and 2.7 years. The very small difference in e0 between women and men projected for 2065 in the Netherlands (0.1 years) and in the United Kingdom (0.9 years) indicates that the individual LC projection of all-cause mortality might result in an unlikely crossover between men and women in e0 in the long-term future due to the more favourable past trends in all-cause mortality among men than among women in the 1990–2014 period. Thus, our projection avoids an unlikely crossover in future life expectancy values between men and women.
Moreover, unlike the LC mortality projection, our projection results in non-linear future increases in e0 (Figure 3; Figure 3—figure supplement 2). Particularly for men in Belgium, Greece, the Czech Republic, Hungary, and Poland, but also for women in Hungary, Denmark, and the United Kingdom, we see faster projected increases in life expectancy in the early decades, but slower increases in life expectancy further in the future, in line with the projections of lifestyle-attributable mortality. The LC projection, by contrast, anticipates more constant future increases in life expectancy at birth.
The differences explained
The observed differences between our projection and the Lee–Carter projection can stem from either the lifestyle dimension, the coherent dimension, or both. Within the lifestyle dimension, two mechanisms apply: first, the projection of non-lifestyle mortality instead of all-cause mortality, and, second, the effect of projected lifestyle-attributable mortality itself.
The higher values estimated by our projection than by the LC projection stem from both dimensions (Table 2; Supplementary file 1B). The coherent dimension led to higher values because the French, Italian, and Spanish women combined exhibited more favourable trends in non-lifestyle-attributable mortality than the individual populations combined, after controlling for the underlying age pattern of mortality. The lifestyle dimension resulted in an additional effect, which was particularly large for women. For women, unlike for men, the past trends in non-lifestyle-attributable mortality were more favourable than the past trends in all-cause mortality. For both men and women, an additional positive effect emerges from the projected (eventual) decline in lifestyle-attributable mortality.
The smaller projected differences between countries in our projection than in the LC projection also stem from both elements, although for men, the impact of integrating lifestyle is more pronounced (Table 2; Supplementary file 1B). Particularly for men, the past trends in non-lifestyle-attributable mortality were much more similar between countries than the past trends in all-cause mortality. In addition, for both men and women, country differences in lifestyle-attributable mortality in 2065 are projected to be smaller than they are currently.
The smaller sex differences in projected e0 values in our projection stem purely from the convergence dimension, whereas the nonlinearity of future trends in our projection stems completely from the lifestyle dimension. Figure 3—figure supplement 3 illustrates that the inclusion of the lifestyle dimension in the Lee–Carter projection (thus ignoring the coherent dimension) leads to a future e0 that is moving back upwards towards the future trend in e0 for non-lifestyle-attributable mortality.
Discussion
Appraisal of our projection methodology
An important element of our projection approach is the identification of the underlying long-term mortality trend that can be extrapolated into the future. We illustrated that past trends in non-lifestyle-attributable mortality serve this purpose better than past trends in all-cause mortality. For non-lifestyle-attributable mortality, the differences in past trends between countries and between men and women are smaller, and the trends are more robust. In contrast, the trends in all-cause mortality exhibit a wave-shaped fluctuation due to initial increases in smoking-, obesity-, and alcohol-attributable mortality, followed (eventually) by declines. This pattern leads, first, to a deceleration of increases in e0, and a move away from the underlying increase in e0 in non-lifestyle attributable mortality, followed by an acceleration of increases in e0, and a move back towards the underlying increase in e0 in non-lifestyle-attributable mortality. Thus, extrapolating past declines in all-cause mortality will lead to non-reliable outcomes: i.e., the outcomes will be either too low when decelerating increases are extrapolated or too high when accelerating increases are extrapolated. Similarly, the extrapolation of past trends in all-cause mortality will result in non-robust outcomes, when updating the projection after a few years. Because the past trends in non-lifestyle-attributable mortality are more linear than the past trends in all-cause mortality, especially from a historical point of view, our projection is less dependent on the data period used for the extrapolation.
Our finding of an average annual past increase in e0 for non-lifestyle-attributable mortality of 0.21 years for men and 0.22 years for women across the 18 European countries over the 1990–2014 period is lower than the continuous upward trend in record life expectancy of 0.243 years that Oeppen and Vaupel, 2002 observed over the 1840–2000 period. There are two main reasons. First, whereas record values always increase, the levels in individual countries may not, and since the 1990s, the record life expectancy has been observed for a non-European country (Japan). Second, more recent increases in life expectancy are driven more by declines in mortality from chronic diseases at older ages with a less strong effect on life expectancy than from declines in mortality from infectious diseases at younger ages, which drove increases in life expectancy up to 1950 (Omran, 1998; Vallin and Meslé, 2009). As a result, a deceleration of increases in life expectancy over time is inevitable. Nonetheless, these large past increases in e0 for non-lifestyle-attributable mortality indicate that considerable increases in e0 are likely in the future.
We applied a coherent mortality projection approach to past trends in non-lifestyle-attributable mortality. This resulted in more realistic estimates of the differences in future e0 values between countries and sexes than those of the individual Lee–Carter projection of past trends. Another added value of our inclusion of the mortality experiences of other countries in our projection is that it provided a broader empirical basis for the determination of the most likely future developments (Janssen and Kunst, 2007), which, in turn, improved the robustness of our projection outcomes (Stoeldraijer, 2019). However, next to a potential effect on the outcomes of the exact coherent methodology that is applied (Stoeldraijer, 2019), it is important to take into account the choice of the reference population (Stoeldraijer, 2019; Booth and Tickle, 2008). Using forerunner populations as the reference population, as advocated by among others Oeppen and Vaupel, 2002 and Bergeron-Boucher et al., 2018, instead of an average across many populations (Li and Lee, 2005; Eurostat, 2020a) will generally result in higher future life expectancy values (see the next section as well). However, given that the past trend in e0 for non-lifestyle-attributable mortality is more equal between countries and sexes than the past trend in e0 for all-cause mortality (see Table 1), which reference populations are chosen will likely have a larger impact when coherently forecasting all-cause mortality than when coherently forecasting non-lifestyle-attributable mortality (see also Janssen and Kunst, 2007 and Janssen et al., 2013). In addition, our convergence assumption relies on – at least – a continuation of improvements in socio-economic developments and medical care. Such continued improvements depend on continued investments in these areas, including diminishing socio-economic inequalities, better preparations for additional possible outbreaks of infectious diseases, and for the potential harmful effects of climate change. Although this might be considered optimistic, our assumption is supported by evidence that life expectancy has long followed a continuously increasing trend (Oeppen and Vaupel, 2002), and that this trend is even stronger from a historical perspective if the effects of lifestyle factors are excluded (see the description of Figure 2 and its supplement in the Results section).
To the projected long-term trend in non-lifestyle-attributable mortality, we added estimates of how smoking-, obesity-, and alcohol-attributable mortality are likely to develop in the near future – and, consequently, how long it will take before all-cause mortality trends again follow more closely the underlying trends in non-lifestyle-attributable mortality. We generated these estimates based on previously projected age-specific smoking-, obesity-, and alcohol-attributable mortality (Janssen et al., 2020c; Janssen et al., 2020b; Janssen et al., 2020d). To avoid generating unrealistic outcomes for the (long-term) future (e.g., continuing increases in smoking-attributable mortality for women, in obesity-attributable mortality, and in alcohol-attributable mortality in the few Eastern European countries where it is still increasing), the projections were based on advanced approaches that did not simply extrapolate past trends, but also took into account knowledge about the future progression of smoking-, alcohol-, and obesity-attributable mortality. Nonetheless, the projection of lifestyle-attributable mortality involves more subjectivity than the extrapolation of non-lifestyle-attributable mortality (see Janssen et al., 2020c; Janssen et al., 2020b; Janssen et al., 2020d for a critical evaluation). For example, in the separate projections of smoking-, obesity-, and alcohol-attributable mortality, population-specific lower bounds were implemented to avoid unlikely future estimates of zero smoking, obesity, and alcohol prevalence, in line with the theory (Lopez et al., 1994; Xu and Lam, 2018), and to prevent unlikely crossovers between sexes and between countries in alcohol-attributable mortality. This approach could be considered conservative. On the other hand, however, the use of a wave-shaped pattern for obesity, and the large declines in alcohol-attributable mortality we project for Eastern European countries, could be considered too optimistic.
However, in our view, the added value of integrating lifestyle into mortality projections outweighs the uncertainties that come with it. First, trends in non-lifestyle-attributable mortality provide a better approximation of the underlying mortality trend. Second, our projection is driven not only by data, but also by theory, and is preceded by a careful study of past trends. Third, by distinguishing between the underlying mortality trend and the additional effect of the three lifestyle factors combined, these projections provide much greater transparency and insight than a mere mechanical extrapolation of past trends could. Fourth, our approach resulted in more realistic differences between countries and sexes than those projected by other approaches (see the next section as well). The abovementioned four added values of our mortality projection approach can be clearly linked to important criteria for evaluating the performance of mortality forecasts in addition to accuracy: namely, robustness, plausibility, reasonableness, and transparency (Cairns et al., 2009). Moreover, the outcomes of our projection approach – which rely above all on the coherent projection of non-lifestyle-attributable mortality – are less dependent on the explicit choices that underlie general mortality forecasts (e.g., the choice of the forerunner populations and the historical time period).
Comparison with other projections
Compared to the benchmark Lee–Carter mortality extrapolation, our projection approach results in higher e0 in the long run; smaller differences in future e0 between countries and sexes; and non-linear future increases in e0, in line with the shift towards the higher e0 values for non-lifestyle-attributable mortality.
Moreover, compared to the most recent official mortality forecast by Eurostat, 2020a and the United Nations, 2020, our mortality projection resulted in substantially higher projected e0 values in 2065 (Supplementary file 1C). Our values were, on average, 5.0 years higher for men and 3.2 years (Eurostat) and 4.3 years (UN) higher for women. The UN projects very large and unrealistic differences in future life expectancy between countries and sexes, including higher projected life expectancy for men than for women in a number of countries. These outcomes clearly show the added value of adopting a coherent mortality projection approach. While Eurostat has performed a coherent mortality projection, it used a different methodology. Eurostat made a partial convergence assumption, applied to men and women separately, and considered the aggregate all-cause mortality of 12 countries (Belgium, Denmark, Germany, Spain, France, Italy, the Netherlands, Austria, Portugal, Finland, Sweden, and the United Kingdom) as the target (Lanzieri, 2009; Eurostat, 2020b). In addition, Eurostat did not integrate the role of lifestyle factors.
The added value of our projection approach over pure coherent mortality forecasting is that we projected the more robust past trends in non-lifestyle-attributable mortality rather than in all-cause mortality, and we projected non-linear future trends in line with the projection of lifestyle-attributable mortality. The reason why coherent mortality projections have to be applied to non-lifestyle-attributable mortality rather than to all-cause mortality is that there are important differences in smoking-, obesity-, and alcohol-attributable mortality trends between countries and sexes (Janssen, 2020; Vidra et al., 2019; Janssen et al., 2020d).
Among the advantages of taking the effects of obesity and alcohol in addition to those of smoking into account is that e0 values for non-lifestyle-attributable mortality are higher than e0 values for non-smoking-attributable mortality (Figure 2, Figure 2—figure supplement 1), and thus have the potential for larger future increases. In addition, as was illustrated, the effects of both obesity-attributable mortality and alcohol-attributable mortality (particularly among Eastern European men) on past and future life expectancy trends are large.
Overall conclusion and implications
All in all, the outcomes of our projection approach, which simultaneously takes into account the lifestyle factors smoking, obesity and alcohol, and the mortality experiences of forerunner populations, can be considered more realistic, more robust, and more insightful than previous extrapolative mortality projections, including those by the United Nations and Eurostat.
Our approach, which distinguishes between the underlying long-term mortality decline and the remaining factors that cause deviations from this mortality decline, could also be adopted to account for the impact of the COVID-19 pandemic on mortality and life expectancy. First, it would be illogical to extrapolate the mortality trends up to 2020 to obtain life expectancy estimates for 2021 and beyond, because life expectancy in 2020 has been importantly negatively affected by COVID-19 (Marois et al., 2020), and will likely to continue to do so for at least 2021. This effect cannot be predicted using a simple extrapolation approach. Second, the underlying long-term mortality decline is argued not to be affected by COVID-19 (e.g., Stoeldraijer, 2020), in line with the predominant temporal dip in life expectancy as a result of the Spanish flu in 1918 (see Roser, 2020). Although, as mentioned, the continuation of the underlying long-term mortality decline does depend on important investments to achieve continued improvements in socio-economic developments and medical care.
Our findings of higher projected life expectancy values than those estimated by the benchmark Lee–Carter mortality projection and the official mortality forecasts by Eurostat and the United Nations have important implications for society. First, individuals could consider planning their life course (e.g., education, work) differently given their longer expected lifespan. Second, a higher projected life expectancy suggests that there will be more older people in the population, which will likely bring both opportunities (e.g., an experienced workforce) and challenges in terms of healthcare planning and social security. Our results imply that either the pension age will have to be raised or, if it is not, that the capital requirements of pension providers will need to be increased. For determining the pension age based on projected life expectancy, a projection that is more robust will have significant added value. More generally, having a more stable projection of mortality will allow for more robust population forecasts.
In addition, our findings of (1) higher projected life expectancy values than those currently used in most mortality projection methodologies and (2) of a rather robust long-term increase in e0 after the effects of smoking, obesity, and alcohol are controlled for, indicate that a limit to human life expectancy (average lifespan) is not within reach – an issue that continues to spark scientific debate (see the first paragraph of the introduction).
Data sharing
Some of our original data regarding lifestyle-attributable mortality were based on previous publications, which, in turn, used data that are openly available. The all-cause mortality data and the exposure data can be obtained through the Human Mortality Database. We have provided source data files for all our tables and figures. These comprise the numerical data that are represented in the different figures, and the output on which the different tables are based. In addition, one excel file with all the final numerical/output data that were used as input for the tables and figures will be made available at the Open Science Framework: https://osf.io/ghu45/. In addition, we will upload there the underlying observed age-specific mortality rates (all-cause mortality, non-lifestyle-attributable mortality, lifestyle-attributable mortality) as well as the adjusted and projected age-specific mortality rates (medians and 90% and 95% projection intervals). The different R codes used for the different steps of the data analyses will be shared – as well – through the Open Science Framework link above.
Appendix 1
Supplementary data and methods
Setting
We included data for the national populations, by sex and age, for 18 European countries, over the period 1990 up to 2014: Denmark, Finland, Norway, Sweden, Austria, Belgium, France, Germany, Ireland, the Netherlands, Switzerland, the United Kingdom, Greece, Italy, the Czech Republic, Hungary, Poland, and Slovenia.
We were restricted to the period 1990–2014 because data on alcohol-attributable mortality were only available from 1990 onwards, which meant, in turn, that we only had data on lifestyle-attributable mortality from 1990 onwards (see Janssen et al., 2021). We included the abovementioned European countries because realistic projections for smoking-, obesity-, and alcohol-attributable mortality were generated for these countries (see Janssen et al., 2020c; Janssen et al., 2020b; Janssen et al., 2020d).
Data
For the abovementioned 18 European countries, for the years 1990 up to 2014, we obtained all-cause mortality (both death numbers and rates) and exposure data by country, sex and single year of age (0–100) from the HMD and University of California, Berkeley (USA), and Max Planck Institute for Demographic Research (Germany), 2019 (update May 1 2019). In addition, we used estimated lifestyle-attributable mortality fractions by country, sex, and single year of age (0–100) from a recent study (Janssen et al., 2021). These lifestyle-attributable mortality fractions represent the share of mortality that can be attributed to smoking, obesity, and alcohol. Multiplying these fractions to the age- and sex-specific all-cause mortality rates, we obtained the age- and sex-specific lifestyle-attributable mortality rates. Subtracting these lifestyle-attributable mortality rates from all-cause mortality rates resulted in non-lifestyle-attributable mortality rates, which – consequently – refers to mortality that is not attributable to smoking, obesity, and alcohol.
In addition, we used for our analysis already published projections of age-specific smoking-, obesity-, and alcohol-attributable mortality fractions for ages 0–84 (Janssen et al., 2020c; Janssen et al., 2020b; Janssen et al., 2020d).
Overall projection approach
Our approach to mortality forecasting relied on an analysis of past mortality trends and their determinants. In line with the evidence gathered, we distinguished between (1) the general and gradual long-term mortality decline not affected by the three lifestyle factors that could be extrapolated into the future, while taking into account the mortality experiences of other countries, and (2) deviations from and differences in this general mortality decline caused predominantly by the time-varying impact of smoking, obesity, and alcohol on mortality, which required the use of more advanced projection techniques.
More specifically, our mortality approach consists of the combination of two projections:
A projection of non-lifestyle-attributable mortality using the coherent Li–Lee (LL) projection methodology, which takes into account the mortality experiences of other populations; and.
A projection of lifestyle-attributable mortality fractions based on already published projections of smoking-, obesity-, and alcohol-attributable mortality fractions.
Main outcome measure
Our main outcome measure is life expectancy at birth (e0) by country and sex, which we estimated for the period 1990 up until 2014, and which we projected up to 2065. Life expectancy at birth is a very important summary measure of health. In addition, it is the most common summary output of mortality forecasts, although the reporting of other summary output measures, such as lifespan disparity, has recently attracted attention (Bohk-Ewald et al., 2017). We have chosen a relatively long projection horizon (2015–2065) given the comparatively short historical time series available (1990–2014) to illustrate that our approach is able to generate reliable outcomes for the long-term future.
To obtain life expectancy at birth, we applied standard life table techniques (Preston et al., 2001) to either the observed or the projected age-specific mortality rates (0–130) for the different national populations under study. In doing so, we used the exponential method to transfer the age-specific mortality rates (mx) into age-specific probabilities (qx), by using the formula qx=1-exp(-mx).
To obtain age- and sex-specific all-cause and non-lifestyle-attributable mortality rates for ages 0–130 to be used as input in our life table calculations, we obtained the respective mortality rates for ages 101–130 by applying the Kannisto model of old-age mortality (Thatcher et al., 1998) to mortality for ages 80–100 for all-cause mortality and non-lifestyle-attributable mortality, respectively.
We compared our projection outcomes foremost with the outcomes of a Lee–Carter extrapolation (Lee and Carter, 1992) applied to all-cause mortality (see next section). In addition, however, we also performed additional projections, to illustrate the effects of the different steps in our projection approach (see the prefinal section of this Appendix). All our projections are based on 50,000 simulations, generated from the forecasts of the time series model. For future e0 we used the median values, and the 95% projection intervals.
Lee–Carter mortality projection
The Lee–Carter (LC) model decomposes logged age-specific mortality rates over time [ln(mx,t)] into an average age effect across time a(x), an overall time trend across all ages k(t), and an age-specific adjustment b(x) of this overall time trend (Lee and Carter, 1992). Projection into the future involves extrapolating the overall time trend k(t). The formula for the Lee–Carter model reads as:
We generated the parameter estimates based on Poisson-likelihood (Brouhns et al., 2002; Renshaw and Haberman, 2006). For the projection of the period parameter in the Lee–Carter model, we applied the common assumption of a random walk with drift ( = ARIMA(0,1,0) with drift) (Lee and Carter, 1992; Renshaw and Haberman, 2006).
Steps involved in our projection
Our projection involved four steps.
Determining the long-term decline in mortality and life expectancy without the combined effect of smoking, obesity, and alcohol.
Projecting non-lifestyle-attributable mortality into the future, thereby taking into account the mortality experience of other populations.
Obtaining estimates of future mortality that is attributable to smoking, obesity, and alcohol.
Combining the projected non-lifestyle attributable mortality rates with the projected lifestyle-attributable mortality fractions.
Step (1) Determining the long-term trend in mortality and life expectancy without the effect of smoking, obesity, and alcohol
As a first step in our projection, we assessed the trends over time in mortality and life expectancy without the effect of smoking, obesity, and alcohol. More specifically we assessed the gain in life expectancy at birth (e0) over the period 1990–2014 for non-lifestyle-attributable mortality and compared this with the observed gain for all-cause mortality. To this end, we used existing recent age- and sex-specific estimates of the share of mortality due to smoking, obesity and alcohol (Janssen et al., 2021) for ages 20–100, which we refer to as lifestyle-attributable mortality fractions (LAMF).
These age- and sex-specific lifestyle-attributable mortality fractions were estimated based on age- and sex-specific smoking-, obesity-, and alcohol-attributable mortality fractions from earlier studies (Janssen, 2020; Janssen et al., 2020b; Janssen et al., 2020d). Smoking-attributable mortality fractions (1950–2014; 35–100 M; 40–100 F) were indirectly estimated by applying a simplified version of the commonly applied Peto–Lopez methodology (Peto et al., 1992; Janssen et al., 2013; Janssen, 2020) to lung cancer mortality data from the WHO Mortality Database (update April 11 2018) (WHO, 2018b), and to additional epidemiological information regarding smoking from the American Cancer Society’s Cancer Prevention Study II (ACS-CPSII) (Peto et al., 1992; Thun et al., 1997). Obesity-attributable mortality fractions (1975–2016; 20–100) were estimated by applying the population-attributable fraction (PAF) formula to age- and sex-specific obesity prevalence data from the NCD Risk Factor Collaboration study 2017 (Abarca-Gómez et al., 2017) and age group- and sex-specific relative risks of dying from obesity from a meta-review (DYNAMO-HIA Consortium, 2010; Janssen et al., 2020a). Alcohol-attributable mortality fractions (1990–2016; 20–100) were estimated based on alcohol-attributable mortality rates obtained from the Global Burden of Disease Study 2017 (Stanaway et al., 2018; GBD, 2018). To ages 65+, we applied the age pattern observed for the main group of causes of death wholly attributable to alcohol (Semyonova et al., 2014), using data from the WHO Mortality Database (WHO, 2018b) and the Human Cause-of-Death Database, 2017. In the estimation of smoking-, obesity-, and alcohol-attributable mortality fractions, Loess smoothing was used to convert estimates by five-year age groups into estimates by single years of age. By setting the smoking-attributable mortality fractions for ages 20–34 (men) and ages 20–39 (women) to zero, and by restricting the data to the same time period, we ended up with estimated age- and sex-specific smoking-, obesity-, and alcohol-attributable mortality fractions for the 1990–2014 period for ages 20–100.
We, consequently, estimated the past age-and sex-specific lifestyle-attributable mortality fractions for the three lifestyles combined (LAMF) (1990-2014; 20-100), by multiplicatively aggregating the population-attributable mortality fractions for the individual risk factors () using the formula:
(Ezzati et al., 2003; Janssen et al., 2021). For ages 0–20 we assumed these fractions to be zero.
By applying these age- and sex-specific lifestyle-attributable mortality fractions to age-and sex-specific all-cause mortality rates, we obtained age- and sex-specific non-lifestyle-attributable mortality rates (1990–2014; 0–100).
To obtain the age- and sex-specific all-cause and non-lifestyle-attributable mortality rates for ages 0–130 to be used as input in our life table calculations, we obtained the respective mortality rates for ages 101–130 by applying the Kannisto model of old-age mortality (Thatcher et al., 1998) to mortality for ages 80–100 for all-cause mortality and non-lifestyle-attributable mortality, respectively.
Step (2) Projecting non-lifestyle-attributable mortality into the future, thereby taking into account the mortality experience of other populations
To project age- and sex-specific non-lifestyle-attributable mortality rates (1990–2014; 0–100) into the future, thereby taking into account the mortality experience of other populations, we applied the often used Li and Lee coherent projection methodology (Li and Lee, 2005; Stoeldraijer, 2019). In doing so, we selected women in France, Spain and Italy as the forerunner populations in terms of life expectancy in Europe because they exhibit both very high recent life expectancy values and very favourable long-lasting past trends in life expectancy (measured by strong increases in e0 over the 1970–2011 period among non-Eastern European countries) (Stoeldraijer, 2019, p. 106). Our decision to use forerunner populations as the reference populations is in line with the seminal work by James Vaupel (e.g., Oeppen and Vaupel, 2002), and our decision to use female populations is in line with the recent work by Bergeron-Boucher et al., 2018. In addition, using three populations instead of one is regarded more robust.
The Li–Lee (LL) coherent projection methodology, in essence, applies the Lee–Carter (LC) methodology twice. First, the LC model is applied to the forerunner populations, resulting in a common time trend (=a common improvement over time across all ages) K(t), a common age effect that adjusts the rate of mortality decline for each age B(x) and a common average age-specific death rate over time A(x). Second, the LC model is applied to the residuals, resulting in population-specific time trends k(t,i), population-specific age effects b(x,i), and population-specific average age-specific death rates a(x,i). These parameters are subsequently combined into separate models for each population, by the formula below:
Similar as to the Lee–Carter model, we generated the parameter estimates based on Poisson-likelihood (Brouhns et al., 2002; Renshaw and Haberman, 2006).
The past trend in the common time trend K(t) is extrapolated using the ARIMA(1,1,0) model with drift, which provided the best fit (=minimum corrected Akaike Information Criterion (AICc)) for both all-cause mortality and non-lifestyle-attributable mortality. For the remaining population-specific time trends k(t,i), we used a random walk with no drift, and thus assumed non-stationarity, which means that we used the last observed difference between the population and the common (instead of mean reverting). This assumption is in line with the results of the statistical test of it in most cases (see Appendix 1—table 1). For comparative purposes, we applied the assumption to all populations. Our assumption implies that the population-specific time trends k(t,i) will rely on the past trend in K(t) among the forerunner populations. However this will not always result in parallel age-specific mortality trends between the respective population and the forerunner populations, because of differences between populations in the average age-specific death rate a(x,i) and the age-specific adjustments of the overall time trend b(x,i).
Test for stationarity of k(t,i) within the coherent Li-Lee mortality projections, for all-cause mortality and non-lifestyle-attributable mortality in which French, Italian, and Spanish women are selected as the forerunner populations, for the 18 European countries, by sex.
All-cause mortality | Non-lifestyle-attributable mortality | |||
---|---|---|---|---|
Men | Women | Men | Women | |
Austria | 1 | 1 | 1 | 1 |
Belgium | 1 | 1 | 1 | 1 |
Czech Republic | 1 | 1 | 1 | 1 |
Denmark | 1 | 1 | 1 | 1 |
Finland | 1 | 1 | 1 | 1 |
France | 0 | 1 | 1 | 1 |
Germany | 1 | 1 | 1 | 0 |
Greece | 1 | 1 | 1 | 1 |
Hungary | 2 | 1 | 1 | 1 |
Ireland | 1 | 1 | 1 | 1 |
Italy | 1 | 0 | 1 | 0 |
Netherlands | 1 | 1 | 1 | 0 |
Norway | 1 | 1 | 1 | 1 |
Poland | 1 | 1 | 1 | 1 |
Slovenia | 2 | 1 | 1 | 1 |
Sweden | 1 | 1 | 1 | 1 |
Switzerland | 1 | 1 | 1 | 1 |
United Kingdom | 1 | 1 | 1 | 1 |
-
0 = stationary time series; i.e., the time series does not exhibit a trend.
1 or 2 = non-stationary time series.
We obtained, for each sex-specific population, 50,000 simulated matrices of age (x) times period (t) of non-lifestyle-attributable mortality rates.
To obtain projected mortality rates for ages 101–130, for each sex-specific population, we extrapolated the age pattern for ages 80–100 up until 130 for each forecast year, by applying the Kannisto model of old-age mortality (Thatcher et al., 1998).
Step (3) Obtaining estimates of future mortality that is attributable to smoking, obesity and alcohol
We obtained projected lifestyle-attributable mortality fractions (0–130) based on previously published data- and theory-driven projections of smoking-, obesity-, and alcohol-attributable mortality (20–84) (Janssen et al., 2020d; Janssen et al., 2020b; Janssen et al., 2020a). See Appendix 1—box 1 for some essential information on these projections.
We obtained the projected fractions by lifestyle factor and sex-specific population, by means of 50,000 simulation matrices of age (x) times period (t), for ages 20 up until 84, and for 2015 (or the first available year in the projection) up to 2065. See Appendix 1—table 2 for the available projection years per country and lifestyle factor. The simulation matrices were generated from the forecasts of the different underlying time series models.
Overview of the data used for the projections per lifestyle factor and combined.
Age | Smoking-attributable mortality fractions (SAMF) | Alcohol-attributable mortality fractions (AAMF) | Obesity-attributable mortality fractions (OAMF)* | Lifestyle-attributable mortality fractions (LAMF)† | ||||
---|---|---|---|---|---|---|---|---|
35–100 M; 40–100 F | 20–100 | 20–100 | 20–100‡ | |||||
Years by country | ||||||||
Austria | 1955 | 2014 | 1990 | 2014 | 1975 | 2016 | 1990 | 2014 |
Belgium | 1954 | 2015 | 1990 | 2015 | 1975 | 2016 | 1990 | 2015 |
Czech Republic | 1953 | 2016 | 1990 | 2016 | 1975 | 2016 | 1990 | 2016 |
Denmark | 1951 | 2015 | 1990 | 2016 | 1975 | 2016 | 1990 | 2015 |
Finland | 1952 | 2015 | 1990 | 2015 | 1975 | 2016 | 1990 | 2015 |
France | 1950 | 2015 | 1990 | 2015 | 1975 | 2016 | 1990 | 2015 |
Germany | 1970 | 2015 | 1990 | 2015 | 1975 | 2016 | 1990 | 2015 |
Greece | 1961 | 2013 | 1990 | 2013 | 1981 | 2016 | 1990 | 2014§ |
Hungary | 1955 | 2014 | 1990 | 2014 | 1975 | 2016 | 1990 | 2014 |
Ireland | 1950 | 2014 | 1990 | 2014 | 1975 | 2016 | 1990 | 2014 |
Italy | 1951 | 2014 | 1990 | 2014 | 1975 | 2016 | 1990 | 2014 |
Netherlands | 1950 | 2016 | 1990 | 2016 | 1975 | 2016 | 1990 | 2016 |
Norway | 1951 | 2014 | 1990 | 2014 | 1975 | 2016 | 1990 | 2014 |
Poland | 1959 | 2014 | 1990 | 2014 | 1975 | 2016 | 1990 | 2014 |
Slovenia | 1985 | 2014 | 1990 | 2014 | 1983 | 2016 | 1990 | 2014 |
Sweden | 1951 | 2016 | 1990 | 2016 | 1975 | 2016 | 1990 | 2016 |
Switzerland | 1951 | 2015 | 1990 | 2016 | 1975 | 2016 | 1990 | 2015 |
United Kingdom | 1950 | 2015 | 1990 | 2016 | 1975 | 2016 | 1990 | 2015 |
-
*We projected obesity prevalence data (1975–2016) into the future. To obtain the OAMF, we only needed the RRs of dying from obesity. However, to obtain age-standardised fractions, we also needed all-cause mortality data.
†See projection step three for details on how we dealt with projections not being available for all the three lifestyle factors for the years 2015–2016.
-
‡By setting the SAMF for ages 20–34 (M) and ages 20–39 (F) to zero.
§We implemented the median projected AAMF and SAMF for Greece as the observed AAMF and SAMF in 2014 so that we could calculate the lifestyle-attributable mortality up until 2014, and so that we could apply it to the all-cause mortality rates to obtain non-lifestyle-attributable mortality rates. For this purpose, we used all-cause mortality from HMD from 26 February 2020.
Subsequently, we combined the projections (50,000 simulations) of the smoking-, obesity-, and alcohol-attributable mortality fractions to obtain 50,000 simulations of lifestyle-attributable mortality fractions (LAMF) (20–84) by again aggregating them multiplicatively (see Equation 2). That is, for b = 1, …, 50,000 simulations we applied Equation 2 to the matrices SAMF_x,t[b], OAMF_x,t[b], and AAMF_x,t[b].
Because the projections for smoking-, obesity-, and alcohol-attributable mortality are not available for the same years (for obesity from 2017 onwards and for smoking and alcohol in some cases only from 2015 onwards) (see Appendix 1—table 1), we adopted a slightly different procedure for the years 2015 and 2016. That is, for these years, we combined the available observed values for the specific lifestyle factor(s) with the projected values (50,000 simulations) for the specific lifestyle factor(s) for which we already had projections. In practice, we did so by turning the observed values into matrices so that we obtained equally large matrices everywhere.
As a next step, we extrapolated each of the 50,000 simulated age patterns of LAMF (20–84) to age 100 by applying Loess smoothing (for each year-simulation combination separately). More specifically, we applied Loess smoothing with degree = 2 and span = 0.75 on log(LAMF) for ages 20–84 and 130, where the LAMF at age 130 is set to 0.000001 so that a decline is enforced. To obtain LAMF for ages 85–100, we extrapolated using the fitted Loess curve, and subsequently took the exponent. This generated 50,000 simulated matrices of LAMF_x,t for x = 20–100 and t. For ages 0–19 and 101–130, we set the LAMF_x,t to zero.
Appendix 1—box 1.
Projections of smoking-, obesity-, and alcohol-attributable mortality (20–84).
For our estimate of future mortality that is attributable to smoking, obesity and alcohol, we used previously published data- and theory-driven projections of smoking-, obesity-, and alcohol-attributable mortality fractions (20–84) (Janssen et al., 2013; Janssen et al., 2020b; Janssen et al., 2020d). These projections follow the observed wave-shaped dynamic of the smoking epidemic (Lopez et al., 1994; Janssen, 2020) and the hypothesised wave shape of the obesity epidemic (Xu and Lam, 2018; Jaacks et al., 2019). For alcohol-attributable mortality, unrealistic future differences between countries were avoided by assuming that the current increases observed for selected countries will eventually turn into declines (Janssen et al., 2020d). All projections were based on the integration of these wave patterns into either a conventional age-period projection (for obesity prevalence, and consequently obesity-attributable mortality) or age-period-cohort projections (for smoking- and alcohol-attributable mortality). This was done by utilising the fact that a wave pattern occurs when the logistically transformed outcome measure has a quadratic shape.
The projections of obesity-attributable mortality were based on projections of obesity prevalence (Janssen et al., 2020b; Janssen et al., 2020a). For obesity prevalence, a wave shape was obtained by linearly extrapolating the speed of change ( = the first-order difference of the time trend) in the logit of prevalence, from a positive speed (increase) to a zero speed (maximum level of the quadratic shape) to a negative speed (decline). This is implemented by applying the age-period Lee–Carter mortality model to logistically transformed age-specific obesity prevalence data (1975–2016). Subsequently, the first-order difference of the general time trend k(t) is linearly extrapolated by applying a second-order random walk with negative drift. This linear extrapolation is applied from 2000 onwards in the majority of countries, but from 1980 onwards for Central and Eastern European women.
Smoking-attributable mortality fractions (1950–2014/LAY) and alcohol-attributable mortality fractions (1990–2016/LAY) are projected by integrating the wave pattern of epidemics into age-period-cohort mortality modelling and projection (Janssen et al., 2020c; Janssen et al., 2020d). That is, the age-period-cohort model by Cairns et al., 2009 is applied to the respective fractions using a generalised logit function as a link function. Subsequently, the (decelerating) increases in the period parameter (for smoking-attributable mortality among women, and for alcohol-attributable mortality in most Central and Eastern European countries and selected North-Western European countries) are projected by a quadratic curve with correlated errors to obtain future declines. Declines in the period parameter as well as the recent trend in the cohort parameter (after removing unreliable outer cohorts) are extrapolated using time-series forecasting (best ARIMA model under some restrictions).
In performing these projections, upper and lower bounds are implemented to ensure that the projections will not lead to unlikely crossovers between countries or between men and women, and will not lead to zero, in line with the theory (Lopez et al., 1994; Xu and Lam, 2018; Jaacks et al., 2019). More specifically, for obesity prevalence, the population- and age-specific prevalence in 1975 is assumed as the lower bound, except for Eastern European women (other than women in Estonia and Hungary), for whom the age-specific prevalence for men in 1975 is assumed. For alcohol-attributable mortality, different lower bounds by country group and sex were implemented based on their past trends and their past (peak) levels of age-standardised alcohol-attributable mortality fractions (AAMF). For the smoking-attributable mortality fractions (SAMF), a 5% smoking prevalence was assumed as the lower bound for both men and women. In addition, for women, the age-specific SAMF observed for Danish women in the year in which the all-age SAMF was highest is set as the upper bound.
To illustrate our projection of lifestyle-attributable mortality fractions, we calculated past and future all-age estimates of lifestyle-attributable mortality fractions. To avoid dependence on future death or population numbers, we calculated age-standardised mortality fractions by applying the sex- and country-specific age composition of all-cause mortality in 2010. The age-standardised lifestyle-attributable mortality fractions were calculated by directly applying the multiplicative aggregation to the past and future (median) age-standardised smoking-, obesity-, and alcohol-attributable mortality (ages 20–84).
Step (4) Combining the projected non-lifestyle attributable mortality rates with the projected lifestyle-attributable mortality fractions
In the final step of our projection approach, we combined the projected non-lifestyle mortality rates with the projected lifestyle-attributable mortality fraction to obtain our projected all-cause mortality rates. For this purpose, we extended the approach that Janssen et al., 2013 developed for combining the projection of non-smoking-attributable mortality with smoking-attributable mortality fractions.
We ended up with the following formula:
Where stands for age- and year-specific all-cause mortality rates, stands for age- and year-specific non-lifestyle-related mortality rates, and stands for the age- and year-specific lifestyle-attributable mortality fraction.
We applied this formula to the 50,000 simulation matrices of coherently projected non-lifestyle-attributable mortality rates (0–130) and the 50,000 simulation matrices of projected LAMF (0–130) to obtain the 50,000 simulation matrices for all-cause-mortality based on our projection.
Additional projections
As mentioned, we compared our projection outcomes foremost with the outcomes of a Lee–Carter extrapolation (Lee and Carter, 1992) applied to all-cause mortality. In addition, however, we also performed additional projections, to illustrate the effects of the different steps in our projection approach. That is, we also performed a Lee–Carter projection of non-lifestyle-attributable mortality, a Li–Lee coherent projection of all-cause mortality, and a Li–Lee coherent projection of non-lifestyle-attributable mortality. In addition, we assessed the effect of purely incorporating lifestyle, without employing a coherent mortality projection, which constitutes of applying the Lee–Carter projection to non-lifestyle-attributable mortality and adding the estimated future lifestyle-attributable mortality. For these projections, we adopted similar specifications as for our projection approach, and as for the Lee–Carter extrapolation of all-cause mortality. Supplementary file 1C compares e0 in 2065 for all the different projections we employed.
Software
We used the R software for all elements of our analysis and projection. We implemented the Lee–Carter and Li–Lee models using the StMoMo package in R (Millossovich, 2018).
Data availability
Some of our original data regarding lifestyle-attributable mortality were based on previous publications, which, in turn, used data that are openly available. The all-cause mortality data and the exposure data can be obtained through the Human Mortality Database. We have provided source data files for all our tables and figures. These comprise the numerical data that are represented in the different figures, and the output on which the different tables are based. In addition, at the Open Science Framework (https://osf.io/ghu45/) we uploaded (i) one excel file with all the final numerical / output data that were used as input for the tables and figures, (ii) the underlying observed age-specific mortality rates (all-cause mortality, non-lifestyle-attributable mortality, lifestyle-attributable mortality) as well as the adjusted and projected age-specific mortality rates (medians and 90% and 95% projection intervals), and (iii) the different R codes used for the different steps of the projection.
-
Open Science FrameworkID ghu45. Future life expectancy in Europe taking into account the impact of smoking, obesity and alcohol.
References
-
Major strides in forecasting future healthThe Lancet 392:e14–e15.https://doi.org/10.1016/S0140-6736(18)31861-0
-
How long will we live?Population and Development Review 32:605–628.https://doi.org/10.1111/j.1728-4457.2006.00144.x
-
Trends in causes of death in low-mortality countries: implications for mortality projectionsPopulation and Development Review 40:189–212.https://doi.org/10.1111/j.1728-4457.2014.00670.x
-
Mortality modelling and forecasting: a review of methodsAnnals of Actuarial Science 3:3–43.https://doi.org/10.1017/S1748499500000440
-
A poisson log-bilinear regression approach to the construction of projected lifetablesInsurance: Mathematics and Economics 31:373–393.https://doi.org/10.1016/S0167-6687(02)00185-3
-
A quantitative comparison of stochastic mortality models using data from England and Wales and the United StatesNorth American Actuarial Journal 13:1–35.https://doi.org/10.1080/10920277.2009.10597538
-
Pension reforms in the EU since the early 2000's: Achievements and Challenges AheadEuropean Economy Discussion Papers 42:.https://doi.org/10.2765/620267
-
ReportWorkpackage 7: Overweight and Obesity. Report on Data Collection for Overweight and Obesity Prevalence and Related Relative RisksInternational Association for the Study of Obesity.
-
Multi-population mortality models: fitting, forecasting and comparisonsScandinavian Actuarial Journal 2017:319–342.https://doi.org/10.1080/03461238.2015.1133450
-
ReportDealing with the Impact of an Ageing Population in the EU (2009 Ageing Report)Brussels: Commission of the European Communities.
-
ReportMethodology of the Eurostat population projections 2019-based (EUROPOP2019)European Commission/Eurostat.
-
A dynamic factor approach to mortality modelingJournal of Forecasting 32:587–599.https://doi.org/10.1002/for.2254
-
Aging, natural death, and the compression of morbidityNew England Journal of Medicine 303:130–135.https://doi.org/10.1056/NEJM198007173030304
-
ReportGlobal Burden of Disease Study 2017 - ResultsSeattle, United States: Institute for health metrics and evaluation.
-
WebsiteHuman Cause-of death database (HCD)French Institute for Demographic Studies (France) and Max Planck Institute for Demographic Research (Germany) . Accessed June 30, 2017.
-
The obesity transition: stages of the global epidemicThe Lancet Diabetes & Endocrinology 7:231–240.https://doi.org/10.1016/S2213-8587(19)30026-9
-
ReportObesity-Attributable Mortality in the Long-Term Future in EuropeNetherlands Interdisciplinary Demographic Institute NIDI Working Paper 2020-06.
-
Smoking epidemic in Europe in the 21st centuryTobacco Control 22:tobaccocontrol-2020-055658.https://doi.org/10.1136/tobaccocontrol-2020-055658
-
Past and future alcohol-attributable mortality in EuropeInternational Journal of Environmental Research and Public Health 17:9024.https://doi.org/10.3390/ijerph17239024
-
The combined impact of smoking, obesity and alcohol on life-expectancy trends in EuropeInternational Journal of Epidemiology 11:dyaa273.https://doi.org/10.1093/ije/dyaa273
-
The future of death in AmericaDemographic Research 25:1–38.https://doi.org/10.4054/DemRes.2011.25.1
-
ConferenceEUROPOP2008: a set of population projections for the European UnionConference: XXVI IUSSP International Population Conference.
-
Modeling and Forecasting U.S. MortalityJournal of the American Statistical Association 87:659–671.https://doi.org/10.1080/01621459.1992.10475265
-
Trends in European life expectancy: a salutary viewInternational Journal of Epidemiology 40:271–277.https://doi.org/10.1093/ije/dyr061
-
Accounting for smoking in forecasting mortality and life expectancyThe Annals of Applied Statistics 15:437–459.https://doi.org/10.1214/20-AOAS1381
-
Why did Danish women's life expectancy stagnate? The influence of interwar generations' smoking behaviourEuropean Journal of Epidemiology 31:1207–1211.https://doi.org/10.1007/s10654-016-0198-7
-
A descriptive model of the cigarette epidemic in developed countriesTobacco Control 3:242–247.https://doi.org/10.1136/tc.3.3.242
-
Convergence and divergence of life expectancy in Europe: a centennial viewEuropean Journal of Epidemiology 28:229–240.https://doi.org/10.1007/s10654-012-9747-x
-
Limits to human life expectancy: evidence, prospects, and implicationsPopulation and Development Review 17:603–637.https://doi.org/10.2307/1973599
-
StMoMo: an R package for stochastic mortality modellingJournal of Statistical Software 84:i03.https://doi.org/10.18637/jss.v084.i03
-
ConferenceStalling of mortality in the United Kingdom and Europe: an analytical review of the evidenceLSE Social Policy Working Paper. pp. 11–19.
-
BookFragmentation of retirement markets due to differences in life expectancyIn: Ayuso M, editors. OECD Business and Finance Outlook. Paris: OECD Publishing. pp. 177–205.https://doi.org/10.1787/9789264257573-en
-
Inconvenient truths about human longevityThe Journals of Gerontology: Series A 74:S7–S12.https://doi.org/10.1093/gerona/glz098
-
The Epidemiologic transition theory revisited thirty years laterWorld Health Statistics Quarterly 52:99–119.
-
BookDemography - Measuring and Modeling Population ProcessesOxford: Blackwell Publishers.
-
BookTrends in Life Expectancy in EU and Other OECD Countries: Why Are Improvements Slowing?Paris: OECD.
-
A cohort-based extension to the Lee–Carter model for mortality reduction factorsInsurance: Mathematics and Economics 38:556–570.https://doi.org/10.1016/j.insmatheco.2005.12.001
-
BookApproaches to the Assessment of Alcohol-Related Losses in the Russian PopulationIn: Anson J, Luy M, editors. Mortality in an International Perspective. London: Springer. pp. 137–168.https://doi.org/10.1007/978-3-319-03029-6
-
Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017The Lancet 392:1923–1994.https://doi.org/10.1016/S0140-6736(18)32225-6
-
Bevolkingsprognose 2012-2060: model en veronderstellingen betreffende de sterfteBevolkingstrends 61:1–27.
-
ThesisMortality Forecasting in the Context of Non-Linear Past Mortality Trends: An EvaluationPhD Thesis, University of Groningen.
-
BookTrends in tobacco smoking and mortality from cigarette use in Cancer Prevention Studies I (1959 through 1965) and II (1982 through 1988)In: Burns D. M, Garfinkel L, Samet J. M, editors. Changes in Cigarette-Related Disease Risks and Their Implications for Prevention and Control. Bethesda, Md: National Cancer Institute. Smoking and Tobacco Control Monograph. pp. 305–382.
-
The contribution of alcohol to the East-West life expectancy gap in Europe from 1990 onwardInternational Journal of Epidemiology 47:731–739.https://doi.org/10.1093/ije/dyx244
-
WebsiteWorld population prospects 2019Department of Economic and Social Affairs - Population Dynamics. Accessed December 22, 2020.
-
The segmented trend line of highest life expectanciesPopulation and Development Review 35:159–187.https://doi.org/10.1111/j.1728-4457.2009.00264.x
-
WebsiteWHO mortality database. health statistics and health information systemsAccessed April 11, 2018.
Decision letter
-
Belinda NicolauReviewing Editor; McGill University, Canada
-
Eduardo FrancoSenior Editor; McGill University, Canada
-
Ugofilippo BaselliniReviewer
-
Collin PayneReviewer
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
Combining a few analytical strategies, this manuscript proposes a new approach to forecast mortality that distinguishes future changes in lifestyle- and non-lifestyle attributable mortality on life expectancy projection. The authors use data from 18 European countries to test their model and the results yield more optimistic forecasts than other well-established forecasting methods; future generations will have an increased life expectancy than previously expected. The proposed methodology could bring significant benefits when applied to other contexts and represents an important contribution to the field of life expectancy forecasting.
Decision letter after peer review:
Thank you for submitting your article "Future life expectancy in Europe taking into account the impact of smoking, obesity and alcohol" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and a Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Ugofilippo Basellini (Reviewer #1); Collin Payne (Reviewer #2).
As is customary in eLife, the reviewers have discussed their critiques with one another. What follows below is the Reviewing Editor's edited compilation of the essential and ancillary points provided by reviewers in their critiques and in their interaction post-review. Please submit a revised version that addresses these concerns directly. Although we expect that you will address these comments in your response letter, we also need to see the corresponding revision clearly marked in the text of the manuscript. Some of the reviewers' comments may seem to be simple queries or challenges that do not prompt revisions to the text. Please keep in mind, however, that readers may have the same perspective as the reviewers. Therefore, it is essential that you attempt to amend or expand the text to clarify the narrative accordingly.
Essential Revisions (for the authors):
1. Line 206: "France, Spain, and Italy". I am unsure why you included Spain in your forerunner population, since you do not have Spain in your results. This seems not very consistent, and it should be better justified and clarified. I suppose that you do not have lifestyle-attributable mortality for Spain? How do your results change if you exclude Spain from the forerunners (I imagine very little)?
2. Line 392-393: "However, in our view, the added value of integrating lifestyle into mortality projections outweighs the uncertainties that come with it." This should, and can, be empirically tested with out-of-sample exercises.
3. SM Methodology: "All our projections are based on 50,000 simulations" and "We obtained, for each sex-specific population, 50,000 simulated matrices". I was not able to understand how you performed these simulations. Are these generated from the forecasts of the time series model? Or do you employ a bootstrap on the deviance residuals of the fitted models? Or both? More details on how you compute variability of your forecasts are needed.
4. The following text is an expanded explanation of the Public Evaluation of Reviewer 2: "My main issue with the paper as written concerns the assumptions built into the Li-Lee projection of non-lifestyle attributable mortality. I wonder how sensitive the results are to the choice of reference populations here--that is, I'm not sure that female populations of France, Spain, and Italy are the best populations to use as reference here.
These populations currently have high life expectancies, but also have histories of substantial mortality shocks and had much higher mortality rates in the near past. My issue here comes down to the essential difference between current mortality rates and current mortality conditions. The low age-specific mortality rates of these populations are a function of both current mortality conditions (e.g. period effects/current conditions in terms of development, medical care, etc) and of the history of cohort mortality selection pressures among those currently in the population. The chosen reference populations all have older cohorts who have experienced fairly severe mortality selection in early life, which may lead to smaller groups of more selected, "hardier" individuals surviving to later life. I worry that choosing these populations may provide an over-optimistic reference point for convergence (indeed, these three countries have among the largest gaps between period e0 and currently attained cohort e0--see https://doi.org/10.1080/00324728.2019.1618480). So I'd be interested in seeing whether results differed by using a set of low-mortality countries that have faced less severe life-course mortality pressures--e.g. Sweden, Norway, and Switzerland. My hunch is that your life expectancy projections will decline a little, though not massively. But I think it's worth pushing on this assumption a bit to see how sensitive it is."
Reviewer #1 (Recommendations for the authors):
The manuscript introduces a new approach to forecast mortality that explicitly considers the role of lifestyle epidemics (smoking, obesity and alcohol) on the dynamics of mortality. Specifically, a four-step projection methodology is proposed, which distinguishes between non-lifestyle-attributable and lifestyle-attributable mortality. The former is forecast using an extrapolative model, the latter is forecast using a data- and theory-driven approach, and the two components are combined to obtain all-cause mortality forecasts.
The model is applied to 18 European countries and compared to three other well-established forecasting methods (Lee-Carter, United Nations and Eurostat). The proposed approach results in more optimistic forecasts than the other methods, suggesting that future individuals will live longer lives than previously expected. Moreover, forecasts of the proposed methodology are more realistic as they do not display implausible crossovers between sexes and countries which characterize the other methods.
For these reasons, the proposed methodology appears to be an important contribution to the mortality forecasting literature. However, the manuscript has some weaknesses that the authors should address to further improve their work.
Weaknesses
While the forecasts of the proposed approach appear to be more realistic than those of the benchmark model of Lee and Carter (1992), the forecast accuracy of the proposed model is not empirically evaluated. In the recent mortality forecasting literature, great attention has been devoted to out-of-sample validation exercises, which are employed to measure the accuracy of the forecasts (see, e.g.,Shang et al. 2011,Bergeron-Boucher et al. 2017). A point and interval forecast accuracy evaluation of the proposed model as compared to the Lee-Carter one would provide concrete evidence of the new model's strengths. Clearly, given the limited time series available to the authors, only a short-term forecast evaluation can be performed (e.g. 5 and 10 years forecasts).
Moreover, the results and discussions of the authors are exclusively based on life expectancy, which is certainly a very important summary measure of mortality, but it does not capture everything. No results are shown or discussed in terms of forecast mortality rates and lifespan inequality. The former are very relevant to understand the plausibility of the forecast mortality age-pattern. The latter is another important summary measure of mortality that complements life expectancy (van Raalte et al., 2018), and it should be considered in mortality forecasting (Bohk-Ewald et al., 2017).
The data employed to fit the model is limited to 25 years (1990--2014) because lifestyle-attributable mortality data is available only for this period of time. Based on this (rather short) fitting period, the authors forecast mortality for a period of time twice as long (up to 2065, that is a 51 year forecast). This seems rather unbalanced. Moreover, forecasts of lifestyle-attributable mortality seem to be driven more by theory/expert opinion rather that by observed time series, especially for females. I think these potential limitations should be more clearly specified in the manuscript.
Finally, the principal added value of the manuscript – the new approach to forecast mortality – is currently not available to the general public: "Because of the large size of the many different underlying simulation matrices, we are restricted in sharing these data. … The R codes can be requested from the author." One simple way to address this problem is to share a much smaller set of simulations (say 100 instead of 50,000) and warn the users to increase the number of simulations when running the codes on their machines. This would directly provide the research community with the routines to employ the proposed methodology.
References
Bergeron-Boucher, M.-P., V. Canudas-Romo, J. Oeppen, and J. W. Vaupel (2017). Coherent forecasts of mortality with compositional data analysis. Demographic Research 37, 527-566.
Bohk-Ewald, C., M. Ebeling, and R. Rau (2017). Lifespan Disparity as an Additional Indicator for Evaluating Mortality Forecasts. Demography 54 (4), 1559-1577.
Lee, R. D. and L. R. Carter (1992). Modeling and forecasting US mortality. Journal of the American Statistical Association 87 (419), 659-671.
Shang, H. L., H. Booth, and R. Hyndman (2011). Point and interval forecasts of mortality rates and life expectancy: A comparison of ten principal component methods. Demographic Research 25, 173-214.
van Raalte, A. A., I. Sasson, and P. Martikainen (2018). The case for monitoring life-span inequality. Science 362 (6418), 1002-1004.
Reviewer #2 (Recommendations for the authors):
The authors have conducted an ambitious projection exercise, drawing together a number of methods and approaches to seek to understand the impacts of future changes in lifestyle- and non-lifestyle attributable mortality on prospects for life expectancy increase in a set of European countries. This focus on the underlying dynamics of mortality change--that is, disaggregating the long-term changes in non-lifestyle attributable mortality patterns from the much more variable lifestyle-attributable mortality patterns--represents a real step forward in life expectancy projection modeling. These methods could have considerable benefits for projections in other contexts.
My main issue with the paper as written concerns the assumptions built into the Li-Lee projection of non-lifestyle attributable mortality, where I would encourage the authors to test the sensitivity of their projections to alternative choices of reference countries (looking beyond Italy, Spain, and France).
https://doi.org/10.7554/eLife.66590.sa1Author response
Essential Revisions (for the authors):
1. Line 206: "France, Spain, and Italy". I am unsure why you included Spain in your forerunner population, since you do not have Spain in your results. This seems not very consistent, and it should be better justified and clarified. I suppose that you do not have lifestyle-attributable mortality for Spain? How do your results change if you exclude Spain from the forerunners (I imagine very little)?
We selected women in France, Spain, and Italy as the forerunner populations in terms of life expectancy in Europe because they exhibit both very high recent life expectancy values and very favourable long-lasting past trends in e0 (measured by strong increases in e0 over the 1970-2011 period among non-Eastern European countries) (Stoeldraijer 2019) (see page 4 of the Supplementary Data and Methods (now Appendix 1)). Our decision to use forerunner populations as the reference populations is in line with the seminal work by James Vaupel (e.g., Oeppen and Vaupel 2002), and our decision to use female populations is in line with the recent work by Bergeron-Boucher et al. 2018. In addition, using three populations instead of one is regarded more robust.
While we had lifestyle-attributable mortality data for Spain, we ended up excluding the projection outcomes for Spain because we were not able to obtain realistic long-term projections of smoking-attributable mortality for Spanish women. Specifically, we found that the projections of smoking-attributable mortality for Spanish women were not robust to the selection of different historical periods (see Janssen et al. 2020 Tobacco Control). Nonetheless, our main argument for why we selected women in France, Spain, and Italy as the forerunner populations in Europe in terms of life expectancy holds.
Given that the gain in e0 over the 1990-2014 period was slightly higher for Spanish women than it was for women in Italy and France (see Table 1), but that the population size was smaller in Spain than in Italy and France, the exclusion of Spanish women from the forerunner populations will most likely result in only slightly lower future life expectancy values. This expectation proves indeed in line with the results of our sensitivity analysis regarding the selection of the forerunner populations that we present in Author response table 1. That is, the coherent projection of non-lifestyle-attributable mortality in which we did not include Spanish women as forerunner population resulted in life expectancy values in 2065 that were – on average – 0.21 years lower for men, and 0.14 years lower for women. For all-cause mortality we see a similar effect (slightly larger difference).
Effect of the inclusion of alternative reference populations within coherent mortality projections* on the coherently projected e0 values in 2065 for all-cause mortality and for non-lifestyle-attributable mortality.
All-cause mortality | Non-lifestyle attributable mortality | |||||
---|---|---|---|---|---|---|
Country | Coh. Li-Lee: sensivity I | Coh. Li-Lee: sensivity II | Coh. Li-Lee: sensivity III | Coh. Li-Lee: sensivity I | Coh. Li-Lee: sensivity II | Coh. Li-Lee: sensivity III |
Men | ||||||
Austria | -0.20 | -2.19 | -0.76 | -0.13 | -0.77 | -0.94 |
Belgium | -0.17 | -2.08 | -0.61 | -0.16 | -0.81 | -0.90 |
Czech Republic | -0.24 | -2.33 | -0.71 | -0.23 | -0.76 | -1.09 |
Denmark | -0.17 | -2.11 | -0.60 | -0.14 | -0.71 | -0.91 |
Finland | -0.26 | -2.13 | -0.83 | -0.22 | -0.73 | -1.00 |
France | -0.21 | -2.12 | -0.75 | -0.09 | -0.78 | -0.78 |
Germany | -0.23 | -2.09 | -0.85 | -0.22 | -0.70 | -0.95 |
Greece | -0.27 | -2.05 | -0.80 | -0.31 | -0.82 | -1.01 |
Hungary | -0.19 | -2.11 | -0.65 | -0.24 | -0.83 | -1.14 |
Ireland | -0.17 | -2.05 | -0.80 | -0.13 | -0.82 | -1.06 |
Italy | -0.28 | -1.86 | -0.74 | -0.18 | -0.73 | -0.92 |
Netherlands | -0.16 | -1.94 | -0.67 | -0.17 | -0.76 | -0.84 |
Norway | -0.35 | -1.83 | -0.77 | -0.31 | -0.82 | -1.03 |
Poland | -0.22 | -2.14 | -0.67 | -0.28 | -0.79 | -1.01 |
Slovenia | -0.22 | -2.03 | -0.68 | -0.15 | -0.67 | -1.00 |
Sweden | -0.27 | -1.79 | -0.71 | -0.29 | -0.64 | -0.89 |
Switzerland | -0.23 | -1.82 | -0.67 | -0.19 | -0.69 | -0.89 |
United Kingdom | -0.26 | -2.02 | -0.70 | -0.23 | -0.80 | -0.91 |
Average | -0.23 | -2.04 | -0.72 | -0.21 | -0.76 | -0.96 |
Women | ||||||
Austria | -0.19 | -1.61 | -0.56 | -0.16 | -0.73 | -0.81 |
Belgium | -0.21 | -1.67 | -0.72 | -0.09 | -0.76 | -0.83 |
Czech Republic | -0.24 | -1.70 | -0.64 | -0.10 | -0.68 | -0.87 |
Denmark | -0.19 | -1.77 | -0.66 | -0.08 | -0.53 | -0.58 |
Finland | -0.22 | -1.69 | -0.69 | -0.14 | -0.73 | -0.86 |
France | -0.23 | -1.37 | -0.54 | -0.12 | -0.71 | -0.77 |
Germany | -0.21 | -1.55 | -0.66 | -0.16 | -0.66 | -0.82 |
Greece | -0.27 | -1.52 | -0.72 | -0.13 | -0.69 | -0.85 |
Hungary | -0.18 | -1.79 | -0.65 | -0.16 | -0.79 | -1.03 |
Ireland | -0.28 | -1.83 | -0.77 | -0.18 | -0.81 | -0.80 |
Italy | -0.25 | -1.59 | -0.65 | -0.16 | -0.70 | -0.85 |
Netherlands | -0.25 | -1.57 | -0.68 | -0.15 | -0.66 | -0.75 |
Norway | -0.24 | -1.67 | -0.66 | -0.09 | -0.66 | -0.85 |
Poland | -0.18 | -1.85 | -0.69 | -0.19 | -0.71 | -0.90 |
Slovenia | -0.20 | -1.62 | -0.68 | -0.14 | -0.60 | -0.85 |
Sweden | -0.23 | -1.56 | -0.65 | -0.14 | -0.59 | -0.74 |
Switzerland | -0.19 | -1.45 | -0.59 | -0.16 | -0.66 | -0.73 |
United Kingdom | -0.23 | -1.65 | -0.68 | -0.18 | -0.56 | -0.70 |
Average | -0.22 | -1.64 | -0.66 | -0.14 | -0.68 | -0.81 |
-
*original = using women in France, Spain and Italy as the reference population.
sensitivity 1 = using only women in France and Italy as the reference population.
-
sensitivity 2 = using women in Sweden, Norway and Switzerland as the reference population.
sensitivity 3 = using the average across the 18 countries included in the analysis as the reference population.
-
Note: the sensitivity analyses were all based on 1000 simulations, and we assumed non-stationarity for the remaining population-specific time trends k(t,i).
In our revised manuscript, we now provide some additional details regarding why we selected women in France, Spain, and Italy as the forerunner populations in Europe (Results section; lines 219-222; Appendix 1 – Supplementary Data and Methods – page 4). We have also included a discussion regarding the selection of forerunner populations in our discussion (see as well our reply to essential revision nr. 4).
2. Line 392-393: "However, in our view, the added value of integrating lifestyle into mortality projections outweighs the uncertainties that come with it." This should, and can, be empirically tested with out-of-sample exercises.
Although out-of-sample exercises can indeed be used in assessing the accuracy of forecasts, we do not think that such an exercise would be able to test whether the added value of integrating lifestyle into mortality projections outweighs the uncertainties that come with doing so. That is, as well as the accuracy of forecasts, which out-of-sample exercises can indeed help to determine, there are many additional evaluation criteria for assessing the performance of forecasts (see Cairns et al. 2009). These include, for example, biological reasonableness, the plausibility of outcomes, and the robustness of the forecast outcomes (Cairns et al. 2009). Because these latter issues are more difficult to quantify, they are often ignored. Instead, accuracy alone is tested by means of out-of-sample (=backtesting) exercises. We believe that our mortality projection approach is particularly promising because its outcomes can be considered more realistic, more robust, and more insightful than previous mortality projections that were purely based on extrapolation.
These points are highlighted in the remainder of the respective paragraph:
“However, in our view, the added value of integrating lifestyle into mortality projections outweighs the uncertainties that come with it. First, it provides a better approximation of the underlying mortality trend. Second, our projection is driven not only by data, but by theory, and is preceded by a careful study of past trends. Third, by distinguishing between the underlying mortality trend and the additional effect of the three lifestyle factors combined, these projections provide much greater transparency and insight than a mere mechanical extrapolation of past trends could.”
Moreover, in response to the comment by reviewer 1 on this issue, we have now added at the end of this paragraph an explicit reference to the different evaluation criteria used for mortality projections, and have also added a reference to our more realistic outcomes (lines 438-446):
“Fourth, our approach resulted in more realistic differences between countries and sexes than those projected by other approaches (see the next section as well). The above mentioned four added values of our mortality projection approach can be clearly linked to important criteria for evaluating the performance of mortality forecasts in addition to accuracy: namely, robustness, plausibility, reasonableness, and transparency (Cairns et al. 2009). Moreover, the outcomes of our projection approach – which rely above all on the coherent projection of non-lifestyle-attributable mortality – are less dependent on the explicit choices that underlie general mortality forecasts (e.g., the choice of the forerunner populations and the historical time period)”.
In addition, in our text, we avoided using the words “accurate” and “inaccurate” as much as possible and replaced these terms with “reliable” and “unreliable”.
3. SM Methodology: "All our projections are based on 50,000 simulations" and "We obtained, for each sex-specific population, 50,000 simulated matrices". I was not able to understand how you performed these simulations. Are these generated from the forecasts of the time series model? Or do you employ a bootstrap on the deviance residuals of the fitted models? Or both? More details on how you compute variability of your forecasts are needed.
Yes, the simulations are generated from the forecasts of the time series model. We did not employ a bootstrap on the deviance residuals of the fitted models. Thus, we did not take into account parameter uncertainty. We have added these details to the Supplementary Data and Methods (page 2 and page 5)(=now Appendix 1).
4. The following text is an expanded explanation of the Public Evaluation of Reviewer 2: "My main issue with the paper as written concerns the assumptions built into the Li-Lee projection of non-lifestyle attributable mortality. I wonder how sensitive the results are to the choice of reference populations here--that is, I'm not sure that female populations of France, Spain, and Italy are the best populations to use as reference here.These populations currently have high life expectancies, but also have histories of substantial mortality shocks and had much higher mortality rates in the near past. My issue here comes down to the essential difference between current mortality rates and current mortality conditions. The low age-specific mortality rates of these populations are a function of both current mortality conditions (e.g. period effects/current conditions in terms of development, medical care, etc) and of the history of cohort mortality selection pressures among those currently in the population. The chosen reference populations all have older cohorts who have experienced fairly severe mortality selection in early life, which may lead to smaller groups of more selected, "hardier" individuals surviving to later life. I worry that choosing these populations may provide an over-optimistic reference point for convergence (indeed, these three countries have among the largest gaps between period e0 and currently attained cohort e0--see https://doi.org/10.1080/00324728.2019.1618480). So I'd be interested in seeing whether results differed by using a set of low-mortality countries that have faced less severe life-course mortality pressures--e.g. Sweden, Norway, and Switzerland. My hunch is that your life expectancy projections will decline a little, though not massively. But I think it's worth pushing on this assumption a bit to see how sensitive it is."
Yes, the choice of the reference or forerunner populations is indeed an important aspect of multi-population or coherent mortality forecasting (Stoeldraijer 2019; Kjærgaard et al. 2016). However, given that the trend in non-lifestyle-attributable mortality is more equal between countries and sexes than the trend in all-cause mortality (see Table 1), this choice will likely have a larger impact when coherently forecasting all-cause mortality than when coherently forecasting non-lifestyle-attributable mortality (see also Janssen and Kunst 2007 and Janssen et al. 2013). A comparison of the average (unweighted) gains in e0 over the 1990-2014 period for women in France, Italy, and Spain (4.8 for all-cause mortality; 5.0 for non-lifestyle-attributable mortality) with the average (unweighted) gains in e0 for women in Sweden, Norway, and Switzerland (4.1 for all-cause mortality; 4.8 for non-lifestyle-attributable mortality) shows that the selection of these latter populations as the forerunner populations would likely result in lower projection outcomes, particularly for all-cause mortality, and to a much lesser extent for non-lifestyle-attributable mortality. This expectation proves indeed in line with the results of our sensitivity analysis regarding the selection of the forerunner populations that we present in Author response table 1.
In response to the remark by the reviewer, we have now added a discussion on the choice of the reference populations in the “Appraisal of our projection methodology” section (lines 384 – 396): “However, as well as the potential effects on the outcomes of the exact coherent methodology that is applied (Stoeldraijer, 2019), it is important to take into account the choice of the reference population (Stoeldraijer 2019; Booth 2020). Using forerunner populations as the reference population, as advocated by among others Oeppen and Vaupel 2002 and Bergeron-Boucher et al. 2018, instead of an average across many populations (Li and Lee 2005; Eurostat 2020) will generally result in higher future life expectancy values (see the next section as well). However, given that the past trend in e0 for non-lifestyle-attributable mortality is more equal between countries and sexes than the past trend in e0 for all-cause mortality (see Table 1), which reference populations are chosen will likely have a larger impact when coherently forecasting all-cause mortality than when coherently forecasting non-lifestyle-attributable mortality (see also Janssen and Kunst 2007 and Janssen et al. 2013)”.
Reviewer #1 (Recommendations for the authors):
The manuscript introduces a new approach to forecast mortality that explicitly considers the role of lifestyle epidemics (smoking, obesity and alcohol) on the dynamics of mortality. Specifically, a four-step projection methodology is proposed, which distinguishes between non-lifestyle-attributable and lifestyle-attributable mortality. The former is forecast using an extrapolative model, the latter is forecast using a data- and theory-driven approach, and the two components are combined to obtain all-cause mortality forecasts.The model is applied to 18 European countries and compared to three other well-established forecasting methods (Lee-Carter, United Nations and Eurostat). The proposed approach results in more optimistic forecasts than the other methods, suggesting that future individuals will live longer lives than previously expected. Moreover, forecasts of the proposed methodology are more realistic as they do not display implausible crossovers between sexes and countries which characterize the other methods. For these reasons, the proposed methodology appears to be an important contribution to the mortality forecasting literature. However, the manuscript has some weaknesses that the authors should address to further improve their work.
Weaknesses
While the forecasts of the proposed approach appear to be more realistic than those of the benchmark model of Lee and Carter (1992), the forecast accuracy of the proposed model is not empirically evaluated. In the recent mortality forecasting literature, great attention has been devoted to out-of-sample validation exercises, which are employed to measure the accuracy of the forecasts (see, e.g.,Shang et al., 2011,Bergeron-Boucher et al. 2017). A point and interval forecast accuracy evaluation of the proposed model as compared to the Lee-Carter one would provide concrete evidence of the new model's strengths. Clearly, given the limited time series available to the authors, only a short-term forecast evaluation can be performed (e.g. 5 and 10 years forecasts).
It is indeed the case that we did not empirically evaluate the forecast accuracy of the proposed model. There are three reasons why we made this choice. First, accuracy is only one of the many criteria for assessing the performance of forecasts (see Cairns et al. 2009). Among the other evaluation criteria are biological reasonableness, the plausibility of outcomes, the robustness of the forecast outcomes, and transparency (Cairns et al. 2009). Because these latter issues are more difficult to quantify, they are often ignored. Instead, only accuracy is tested by means of out-of-sample exercises (=backtesting). However, as we indicated in our discussion, we consider our mortality projection approach to be particularly promising because its outcomes can be considered more realistic, more robust, more theoretically informed, and more transparent than the outcomes of other approaches. Second, given the limited time series for lifestyle-attributable morality, we could only perform a short-term forecast evaluation, whereas our main aim is to produce reliable long-term forecasts. Third, given the profile of the readers of eLife, we believe that applying additional technical exercises would not necessarily be an improvement.
Related to our choice to avoid empirically evaluating the forecast accuracy of our model, we have avoided making any statements in our discussion that would suggest that our model is more accurate than the Lee-Carter model. Moreover, we now avoid the use of “accurate” and “inaccurate” as much as possible in the remainder of our text. In addition, we now explicitly refer to the different criteria for evaluating mortality projections in the last paragraph of the “Appraisal of our projection methodology” section.
Moreover, the results and discussions of the authors are exclusively based on life expectancy, which is certainly a very important summary measure of mortality, but it does not capture everything. No results are shown or discussed in terms of forecast mortality rates and lifespan inequality. The former are very relevant to understand the plausibility of the forecast mortality age-pattern. The latter is another important summary measure of mortality that complements life expectancy (van Raalte et al. 2018), and it should be considered in mortality forecasting (Bohk-Ewald et al., 2017).
We completely agree that in addition to life expectancy, other outcome measures of mortality forecasts are relevant and interesting. However, we decided to focus on life expectancy at birth in our manuscript because it is a very common summary measure of health; it is the most common output measure of mortality forecasts; it is easier to grasp than lifespan variation; and it summarizes the age pattern of mortality using a single measure. Given the more general audience of eLife, we do not consider it particularly helpful to include additional outcome measures as part of our submission. However, in response to the comment of the reviewer, we have clarified why we have chosen to focus on life expectancy, thereby referring to other outcome measures that are important as well (see lines 168 – 170; see Supplementary Data and Methods (now Appendix 1) page 2 – first paragraph).
In addition, we have added figures with the projected age-specific mortality rates (Author response image 1 and Author response image 2) and a table (Author response table 2) with the observed (2014) and projected (2065) standard deviation (SD) of the age-at-death distribution from age zero onwards ( = the dx column of the lifetable). This measure is often used to denote lifespan variation (see Edwards and Tuljapurkar, 2005; Garcia and Aburto, 2019).
Furthermore, we will make available the projected age-specific mortality rates (median and 95% and 90% projection intervals) on Open Science Framework (OSF).

Age-specific mortality pattern in 1990*, 2014*, and for different projections in 2065, by country, ages 0-130, men. * based on the adjusted rates up to 130 (because of zeros in the observed rates).

Age-specific mortality pattern in 1990*, 2014*, and for different projections in 2065, by country, ages 0-130, women. * based on the adjusted rates up to 130 (because of zeros in the observed rates).
Observed (2014) and projected (2065) lifespan variation from age zero onwards, by country and sex, according to our projection methodology, which takes into account the impact of smoking, obesity, and alcohol and the mortality experiences of forerunner countries (“lifestyle and coherent”); the benchmark Lee-Carter extrapolative mortality projection applied to all-cause mortality (“Lee-Carter”); and when purely accounting for smoking, obesity, and alcohol (“adding lifestyle”).
MEN | WOMEN | |||||||
---|---|---|---|---|---|---|---|---|
Projected SD* 2065 | Projected SD* 2065 | |||||||
Country | SD* 2014 | Lee-Carter | Adding lifestyle | Lifestyle and coherent | SD* 2014 | Lee-Carter | Adding lifestyle | Lifestyle and coherent |
Austria | 14.33 | 10.88 | 11.11 | 10.31 | 13.05 | 9.69 | 8.73 | 8.91 |
Belgium | 14.70 | 11.57 | 11.27 | 9.78 | 13.36 | 11.09 | 9.53 | 9.18 |
Czechia | 14.36 | 10.96 | 12.02 | 10.69 | 12.61 | 9.23 | 8.53 | 9.22 |
Denmark | 14.25 | 10.51 | 10.37 | 10.09 | 13.26 | 10.15 | 10.10 | 9.27 |
Finland | 14.69 | 11.61 | 11.18 | 10.35 | 12.66 | 10.90 | 9.65 | 9.03 |
France | 15.32 | 11.75 | 11.35 | 10.04 | 13.59 | 11.11 | 9.92 | 9.27 |
Germany | 14.25 | 10.66 | 11.31 | 10.39 | 12.93 | 9.61 | 8.85 | 9.23 |
Greece | 15.10 | 14.31 | 14.26 | 10.56 | 12.67 | 10.20 | 9.98 | 9.25 |
Hungary | 15.11 | 12.77 | 11.23 | 11.21 | 13.90 | 11.38 | 9.33 | 10.01 |
Ireland | 14.52 | 11.43 | 10.84 | 9.78 | 13.18 | 10.59 | 9.98 | 9.25 |
Italy | 13.53 | 9.11 | 10.95 | 9.76 | 12.41 | 9.06 | 9.01 | 8.94 |
Netherlands | 13.62 | 9.82 | 10.16 | 9.60 | 13.19 | 11.12 | 9.68 | 8.93 |
Norway | 13.74 | 10.17 | 9.93 | 9.68 | 12.50 | 10.22 | 9.39 | 8.79 |
Poland | 15.93 | 13.49 | 12.52 | 11.50 | 13.92 | 11.08 | 9.65 | 9.86 |
Slovenia | 13.89 | 10.05 | 10.93 | 10.33 | 12.28 | 9.00 | 8.39 | 8.96 |
Sweden | 13.60 | 10.31 | 10.13 | 9.47 | 12.49 | 9.88 | 9.15 | 8.69 |
Switzerland | 13.84 | 9.45 | 10.29 | 9.81 | 12.59 | 9.37 | 8.79 | 8.66 |
United Kingdom | 14.64 | 12.24 | 11.87 | 10.11 | 13.62 | 11.06 | 10.58 | 9.37 |
Average | 14.41 | 11.17 | 11.21 | 10.19 | 13.01 | 10.26 | 9.40 | 9.16 |
Min | 13.53 | 9.11 | 9.93 | 9.47 | 12.28 | 9.00 | 8.39 | 8.66 |
Max | 15.93 | 14.31 | 14.26 | 11.50 | 13.92 | 11.38 | 10.58 | 10.01 |
Variance | 0.42 | 1.81 | 0.99 | 0.28 | 0.25 | 0.61 | 0.34 | 0.12 |
-
* SD = Standard deviation of the age at death distribution (age 0-100).
The data employed to fit the model is limited to 25 years (1990-2014) because lifestyle-attributable mortality data is available only for this period of time. Based on this (rather short) fitting period, the authors forecast mortality for a period of time twice as long (up to 2065, that is a 51 year forecast). This seems rather unbalanced. Moreover, forecasts of lifestyle-attributable mortality seem to be driven more by theory/expert opinion rather that by observed time series, especially for females. I think these potential limitations should be more clearly specified in the manuscript.
In general, it is indeed considered optimal to base projections – and, in particular, extrapolations – on a historical period of about 40-50 years when projecting 40-50 years into the future (Technical Panel on Assumptions and Methods, 2003). Thus, it is a pity that because data on alcohol-attributable mortality were only available from 1990 onwards, we could base our forecast only on data from 1990 to 2014/LAY. However, we decided against restricting our projection horizon to 25 years in the future because we considered it essential to show that our approach is able to generate reliable predictions for the long-term future. It should also be noted that the use of a relatively long historical period stems from the desire to select the most reliable long-term past trend on which the extrapolations could be based. However, because of the substantial non-linearity in the past trends in life expectancy in recent decades, many of the approaches currently used rely on a shorter past period (Stoeldraijer et al. 2013a), or rely on the trends for women that were more linear (e.g. Bergeron-Boucher et al. 2018). Moreover, it is in cases in which the past trends were non-linear that the projection outcomes rely most heavily on the historical time period chosen (Stoeldraijer et al. 2013a). Because our approach relies on the extrapolation of non-lifestyle-attributable mortality, which exhibits more linear past trends than all-cause mortality (especially from a historical point of view)(see Figure 2 and Figure 2 —figure supplement 1)(=previous Figure S2), the effect of the choice of the historical time period is smaller than it is in extrapolations of all-cause mortality. However, including only the trends from 1990 onwards has likely resulted in an underestimation of the effects of the separate extrapolation of non-lifestyle-attributable mortality and of lifestyle-attributable mortality, as for men in particular, the differences between the trends in e0 for all-cause and for non-lifestyle-related mortality would have been much larger if a longer past trend had been included. This is primarily because of the large increases in smoking-attributable mortality among men in many Western European countries, which resulted in stagnating increases in life expectancy in the 1960s and 1970s (see Figure 2); and because of the large increases in alcohol-attributable mortality among Eastern European men since 1975 (e.g. Trias-Llimós et al. 2017), which contributed as well to the past stagnating trends in life expectancy in these regions (Meslé et al. 2002).
We have now added a justification for our decision to use a relatively long projection horizon given the comparatively short historical period (see lines 170-173; appendix 1 – supplementary data and methods – page 2 paragraph 1). In addition, in the discussion, we explained that our approach – because it is focused on non-lifestyle-attributable mortality – is less dependent on the choice of the historical period (lines 360 – 363 and lines 443- 446).
The forecasting of lifestyle-attributable mortality indeed involves more subjectivity than the extrapolation of non-lifestyle-attributable mortality, as we acknowledge in our paper (see lines 419-421). To avoid generating unrealistic outcomes for the (long-term) future (e.g., continuing increases in smoking-attributable mortality for women, in obesity-attributable mortality, and in alcohol-attributable mortality in some Eastern European countries), we had to adopt approaches that did not simply extrapolate past trends, but that also took into account knowledge about the likely future progression of smoking-, alcohol-, and obesity-attributable mortality. In doing so, however, we still made as much use as possible of the past trends (for example, by making use of past trends for men in the case of smoking and of past trends in other European countries in the case of alcohol; and by adopting an extrapolation approach, but then applying it to the speed of change in the logit of prevalence for obesity to ensure a wave-shaped pattern, and by applying a quadratic curve to the past trends for populations that were still experiencing increases in smoking and obesity), and implemented our approach in frequently used models (e.g., the Lee-Carter model and age-period-cohort models). See Appendix 1 – Box 1 for a more extended summary of the different approaches used to project obesity-, smoking-, and alcohol-attributable mortality; and see Janssen, van der Broek et al., 2020; Janssen, El Gewily, Bardoutsos, 2020; and Janssen, El Gewily et al., 2020 for the actual projections, including a critical appraisal of each of them.
We have now extended the discussion of the projection of lifestyle-attributable mortality (see lines 413 – 422):
“To avoid generating unrealistic outcomes for the (long-term) future (e.g. continuing increases in smoking-attributable mortality for women, in obesity-attributable mortality, and in alcohol-attributable mortality in the few Eastern European countries where it is still increasing), the projections were based on advanced approaches that did not simply extrapolate past trends, but also took into account knowledge about the future progression of smoking-, alcohol-, and obesity-attributable mortality. Nonetheless, the projection of lifestyle-attributable mortality involves more subjectivity than the extrapolation of non-lifestyle-attributable mortality (see Janssen, van der Broek et al., 2020; Janssen et al., 2020; Janssen, El Gewily and Bardoutsos, 2020 for a critical evaluation)”.
Finally, the principal added value of the manuscript – the new approach to forecast mortality – is currently not available to the general public: "Because of the large size of the many different underlying simulation matrices, we are restricted in sharing these data. … The R codes can be requested from the author." One simple way to address this problem is to share a much smaller set of simulations (say 100 instead of 50,000) and warn the users to increase the number of simulations when running the codes on their machines. This would directly provide the research community with the routines to employ the proposed methodology.
We carefully considered what data we could share in order to make the new approach to forecasting mortality available to the general public. Although much of the information is already provided on Open Science Framework (https://osf.io/ghu45/) in the Excel file with the final data used as input for the tables and figures, we will now also provide information on the underlying observed age-specific mortality rates (all-cause mortality, non-lifestyle-attributable mortality, lifestyle-attributable mortality), as well as the fitted and projected age-specific mortality rates (by means of medians and 90% and 95% projection intervals). In addition, we will make the different R codes used for the different steps of the data analyses public.
In line with this change, we rewrote our data sharing statement: “Some of our original data regarding lifestyle-attributable mortality were based on previous publications, which, in turn, used data that are openly available. The all-cause mortality data and the exposure data can be obtained through the Human Mortality Database. We have provided source data files for all our tables and figures. These comprise the numerical data that are represented in the different figures, and the output on which the different tables are based. In addition, one excel file with all the final numerical / output data that were used as input for the tables and figures will be made available at the Open Science Framework: https://osf.io/ghu45/. In addition, we will upload there the underlying observed age-specific mortality rates (all-cause mortality, non-lifestyle-attributable mortality, lifestyle-attributable mortality) as well as the adjusted and projected age-specific mortality rates (medians and 90% and 95% projection intervals). The different R codes used for the different steps of the data analyses will be shared – as well – through the Open Science Framework link.
References
Bergeron-Boucher, M.-P., V. Canudas-Romo, J. Oeppen, and J. W. Vaupel (2017). Coherent forecasts of mortality with compositional data analysis. Demographic Research 37, 527-566.
Bohk-Ewald, C., M. Ebeling, and R. Rau (2017). Lifespan Disparity as an Additional Indicator for Evaluating Mortality Forecasts. Demography 54 (4), 1559-1577.
Lee, R. D. and L. R. Carter (1992). Modeling and forecasting US mortality. Journal of the American Statistical Association 87 (419), 659-671.
Shang, H. L., H. Booth, and R. Hyndman (2011). Point and interval forecasts of mortality rates and life expectancy: A comparison of ten principal component methods. Demographic Research 25, 173-214.
van Raalte, A. A., I. Sasson, and P. Martikainen (2018). The case for monitoring life-span inequality. Science 362 (6418), 1002-1004.
Reviewer #2 (Recommendations for the authors):
The authors have conducted an ambitious projection exercise, drawing together a number of methods and approaches to seek to understand the impacts of future changes in lifestyle- and non-lifestyle attributable mortality on prospects for life expectancy increase in a set of European countries. This focus on the underlying dynamics of mortality change--that is, disaggregating the long-term changes in non-lifestyle attributable mortality patterns from the much more variable lifestyle-attributable mortality patterns--represents a real step forward in life expectancy projection modeling. These methods could have considerable benefits for projections in other contexts.
My main issue with the paper as written concerns the assumptions built into the Li-Lee projection of non-lifestyle attributable mortality, where I would encourage the authors to test the sensitivity of their projections to alternative choices of reference countries (looking beyond Italy, Spain, and France).
Yes, the choice of the reference or forerunner populations is indeed an important aspect of multi-population or coherent mortality forecasting (Stoeldraijer 2019; Booth 2020). However, given that the trend in non-lifestyle-attributable mortality is more equal between countries and sexes than the trend in all-cause mortality (see Table 1), this choice will likely have a larger impact when coherently forecasting all-cause mortality than when coherently forecasting non-lifestyle-attributable mortality (see also Janssen and Kunst 2007 and Janssen et al. 2013). We selected women in France, Spain, and Italy as the forerunner populations in terms of life expectancy in Europe because they exhibit both very high recent life expectancy values and very favourable long-lasting past trends in e0 (measured by large increases in e0 over the 1970-2011 period among non-Eastern European countries) (Stoeldraijer 2019). Our decision to use forerunner populations as the reference populations is in line with the seminal work by James Vaupel (e.g., Oeppen and Vaupel 2002), and our decision to use female populations is in line with the recent work by Bergeron-Boucher et al. 2018. In addition, using three populations instead of one is regarded more robust.
An alternative approach would be to use the average mortality experience across many countries (Eurostat 2020; Stoeldraijer 2019), or – as indicated by reviewer 2 in the extended (non-public) explanation of this remark – to use the mortality experience of a set of low-mortality countries that have faced less severe life-course mortality pressures--e.g. Sweden, Norway, and Switzerland.
Author response table 1 reports on the outcomes of a sensitivity analysis in which we employed these two alternative selections of the reference population in coherent mortality projection. Our sensitivity analysis revealed that the use of the average mortality experience across all the included countries resulted – in line with our expectations – in lower projected e0 values in 2065, for both non-lifestyle attributable mortality and all-cause mortality. Furthermore our sensitivity analysis revealed that the use of women in Sweden, Norway and Switzerland as the reference population led to lower projected e0 values in 2065, but particularly so for all-cause mortality and much less so for non-lifestyle-attributable mortality. This latter outcome can be largely explained from the difference in the average (unweighted) gains in e0 over the 1990-2014 period for women in France, Italy, and Spain (4.8 for all-cause mortality; 5.0 for non-lifestyle-attributable mortality) compared to women in Sweden, Norway, and Switzerland (4.1 for all-cause mortality; 4.8 for non-lifestyle-attributable mortality).
The results of our sensitivity analysis, confirm that the use of forerunner populations instead of an average across many populations results in lower projected e0 values in 2065, and that the choice of the reference population will likely have a larger impact when coherently forecasting all-cause mortality than when coherently forecasting non-lifestyle-attributable mortality.
We have added a short discussion on the choice of the reference populations in the “Appraisal of our projection methodology” section of our revised manuscript (lines 384 – 396):
“However, as well as the potential effects on the outcomes of the exact coherent methodology that is applied (Stoeldraijer, 2019), it is important to take into account the choice of the reference population (Stoeldraijer 2019; Booth 2020). Using forerunner populations as the reference population, as advocated by among others Oeppen and Vaupel 2002 and Bergeron-Boucher et al., 2018, instead of an average across many populations (Li and Lee 2005; Eurostat 2020) will generally result in higher future life expectancy values (see the next section as well). However, given that the past trend in e0 for non-lifestyle-attributable mortality is more equal between countries and sexes than the past trend in e0 for all-cause mortality (see Table 1), which reference populations are chosen will likely have a larger impact when coherently forecasting all-cause mortality than when coherently forecasting non-lifestyle-attributable mortality (see also Janssen and Kunst 2007 and Janssen et al., 2013)”.
References
Edwards, R. D. and Tuljapurkar, S. (2005). Inequality in life spans and a new perspective on mortality convergence across industrialized countries. Population and Development Review, 31(4):645–674
Garcia, J., and Aburto, J.M. (2019). The impact of violence on Venezuelan life expectancy and lifespan inequality. International Journal of Epidemiology, 48(5), 1593-1601.
Meslé, F. and J. Vallin (2002). Mortality in Europe: The divergence between East and West. Population, 57: 157-197.
Technical Panel on Assumptions and Methods (2003). Report to the Social Security Advisory Board. Washington, D.C.
Trias-Llimós, S., Bijlsma, M. And Janssen, F. (2017). The role of birth cohorts in long-term trends in liver cirrhosis mortality across eight European countries. Addiction ,112(2), 250-258.
https://doi.org/10.7554/eLife.66590.sa2Article and author information
Author details
Funding
Netherlands Organisation for Scientific Research (Innovational Research Incentives Scheme Vidi/452-13-001)
- Fanny Janssen
- Anastasios Bardoutsos
- Shady El Gewily
Netherlands Organisation for Scientific Research (Innovational Research Incentives Scheme Vici/VI.C.191.019)
- Fanny Janssen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Mark van der Broek (Research Master Economics, University of Groningen) for creating the final figures in R.
Senior Editor
- Eduardo Franco, McGill University, Canada
Reviewing Editor
- Belinda Nicolau, McGill University, Canada
Reviewers
- Ugofilippo Basellini
- Collin Payne
Version history
- Received: January 15, 2021
- Accepted: June 28, 2021
- Accepted Manuscript published: July 6, 2021 (version 1)
- Version of Record published: August 4, 2021 (version 2)
Copyright
© 2021, Janssen et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,460
- Page views
-
- 209
- Downloads
-
- 13
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
Background:
In most of the world, the mammography screening programmes were paused at the start of the pandemic, whilst mammography screening continued in Denmark. We examined the mammography screening participation during the COVID-19 pandemic in Denmark.
Methods:
The study population comprised all women aged 50–69 years old invited to participate in mammography screening from 2016 to 2021 in Denmark based on data from the Danish Quality Database for Mammography Screening in combination with population-based registries. Using a generalised linear model, we estimated prevalence ratios (PRs) and 95% confidence intervals (CIs) of mammography screening participation within 90, 180, and 365 d since invitation during the pandemic in comparison with the previous years adjusting for age, year and month of invitation.
Results:
The study comprised 1,828,791 invitations among 847,766 women. Before the pandemic, 80.2% of invitations resulted in participation in mammography screening within 90 d, 82.7% within 180 d, and 83.1% within 365 d. At the start of the pandemic, the participation in screening within 90 d was reduced to 69.9% for those invited in pre-lockdown and to 76.5% for those invited in first lockdown. Extending the length of follow-up time to 365 d only a minor overall reduction was observed (PR = 0.94; 95% CI: 0.93–0.95 in pre-lockdown and PR = 0.97; 95% CI: 0.96–0.97 in first lockdown). A lower participation was, however, seen among immigrants and among women with a low income.
Conclusions:
The short-term participation in mammography screening was reduced at the start of the pandemic, whilst only a minor reduction in the overall participation was observed with longer follow-up time, indicating that women postponed screening. Some groups of women, nonetheless, had a lower participation, indicating that the social inequity in screening participation was exacerbated during the pandemic.
Funding:
The study was funded by the Danish Cancer Society Scientific Committee (grant number R321-A17417) and the Danish regions.
-
- Epidemiology and Global Health
- Genetics and Genomics
Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is partially maintained among repeated serial samples from the same host, it can transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control.