Skip to main content

Disparities in length of life across developed countries: measuring and decomposing changes over time within and between country groups

Abstract

Background

Over the past half century the global tendency for improvements in longevity has been uneven across countries. This has resulted in widening of inter-country disparities in life expectancy. Moreover, the pattern of divergence appears to be driven in part by processes at the level of country groupings defined in geopolitical terms. A systematic quantitative analysis of this phenomenon has not been possible using demographic decomposition approaches as these have not been suitably adapted for this purpose. In this paper we present an elaboration of conventional decomposition techniques to provide a toolkit for analysis of the inter-country variance, and illustrate its use by analyzing trends in life expectancy in developed countries over a 40-year period.

Methods

We analyze trends in the population-weighted variance of life expectancy at birth across 36 developed countries and three country groups over the period 1970–2010. We have modified existing decomposition approaches using the stepwise replacement algorithm to compute age components of changes in the total variance as well as variance between and within groups of Established Market Economies (EME), Central and Eastern Europe (CEE), and the Former Soviet Union (FSU). The method is generally applicable to the decomposition of temporal changes in any aggregate index based on a set of populations.

Results

The divergence in life expectancy between developed countries has generally increased over the study period. This tendency dominated from the beginning of 1970s to the early 2000s, and reversed only after 2005. From 1970 to 2010, the total standard deviation of life expectancy increased from 2.0 to 5.6 years among men and from 1.0 to 3.6 years among women. This was determined by the between-group effects due to polarization between the EME and the FSU. The latter contrast was largely fueled by the long-term health crisis in Russia. With respect to age, the increase in the overall divergence was attributable to between-country differences in mortality changes at ages 15–64 years compared to those aged 65 and older. The within-group variance increased, especially among women. This change was mostly produced by growing mortality differences at ages 65 and older.

Conclusions

From the early 1970s to the mid-2000s, the strong divergence in life expectancy across developed countries was largely determined by the between-group variance and mortality polarization linked to the East–West geopolitical division.

Peer Review reports

Introduction

The lengthening of the human lifespan is one of the most remarkable achievements of human civilization [13]. Over the past 60 years, global life expectancy has greatly increased. However, life expectancy gains have been uneven across countries, and during the last decades of the 20th century, inter-country disparities expanded [4, 5]. While life expectancy increased steadily in most countries, some countries and regions had larger gains than others, and some even experienced mortality reversals [610]. Among developed countries, the most significant mortality reversals occurred in Eastern Europe. Starting in the late 1960s, this region experienced several decades of rising mortality in particular from cardiovascular diseases, alcohol-related and external causes of death (i.e., injuries, poisonings, and violence). However, at the very end of the 1980s, around the time of the fall of the Berlin Wall, a more favorable dynamic of decreasing mortality emerged in Central and Eastern Europe. Some years later, in the mid-1990s, a similar decline in mortality took place in the former Soviet Baltic states. Finally, in the mid-2000s, favorable changes started in the Slavic countries of the former Soviet Union [1117].

Although in the established market economies (EME) life expectancy has been rising continuously, there has also been variation in the size of these increases, as certain countries (Japan, Spain, Italy, France) have improved their relative positions, while others (Denmark, the United States, the Netherlands) have lost their former advantages.

A number of studies have carried out classic decomposition analyses of life expectancy changes in single countries of Europe and North America with the goal of identifying the age- and cause-specific changes in mortality that are most responsible for the changes in longevity. These have found that in Russia, Ukraine, Belarus, the Baltic countries, Hungary, Poland, and other countries of Eastern Europe, the declining or stagnating life expectancy levels over the 1970s–1980s were largely attributable to rising mortality among young and middle-aged adults, combined with a lack of mortality improvements among the elderly [6, 13, 1821]. It has also been shown also that the speed at which longevity was increasing was higher in some of the advanced countries and lower in others due to variation in mortality improvements among the elderly [2224].

Over the past 15 years studies looking at mortality and life expectancy changes in single countries or comparing trends among a few countries were complemented by analyses of mortality trends across all or most of the countries of the world [4, 8, 25, 26]. The latter studies by Moser et al. [4] and by Smits and Monden [8] used single summary measures of inequality to quantify changes in the amount of worldwide length of life disparity. These studies provided an instrument for an objective detection of transitions from mortality convergence to mortality divergence. This sort of transition has been conceptualized by Meslé and Vallin’s convergence-divergence theory and is attributed to the emergence of new epidemiological challenges which are addressed at first by a few vanguard countries and only later (after some years or even decades) by other countries [27, 28].

The conventional decomposition analyses discussed above have been used for investigating components and driving forces of changes in life expectancy in single countries. In the present study, we show how these methods can be adapted and elaborated in order to analyze changes in summary measures of inter-country disparity based on data from many countries, and then illustrate their use by looking at trends in mortality in 36 developed countries over the period 1970–2010.

The technical details of this method are provided in the next section, but for the sake of clarity we provide in the remainder of the Introduction a high-level summary of our approach. We use the population-weighted cross-country variance as a measure of the amount of inter-country disparity. Because this measure takes into account the population size of each country, it has the advantage of reflecting the overall public health burden produced by the disparity. The weighted disparity measure expresses not only the lifetime differences among countries and country groups, but also how many people in different countries experience longer or shorter lifetimes. However, to see how the weighting influences the results, we provide results of alternative calculations without population weighting.

Our empirical analysis includes decomposition of changes in the life expectancy dispersion across a set of 36 developed countries using the general stepwise replacement algorithm [29]. We further substantiate the converge-divergence theory by new information about components of the life expectancy divergence of the 1970s–1990s and convergence of the 2000s. In particular, for the first time we evaluate the importance of changes in the amount of length of life disparity: 1) of the former East–West geopolitical divide, and specifically of the country groups of established market economies (EME), Central and Eastern Europe (CEE), and the former Soviet Union (FSU); 2) of the three principal age groups (childhood, midlife, and old ages); 3) of mortality and population composition.

Data and methods

Data

We use the Human Mortality Database (HMD) data [30] on deaths and population exposure by sex and age (0, 1–4, 5–9, 10–14,…, 95+) for 36 developed countries and regions (see Appendix 1 for the list of countries) over the period 1970 to 2010.Footnote 1 The starting point approximates the beginning of a new phase of longevity divergence [28], while the end point is the most recent year for which data for the majority of the countries are available (see Appendix 2 for the country-specific life expectancy values).

Between 1970 and 2010, the total population of the countries under study increased from 1.0 to 1.2 billion. Over this period, the share of the total population in the EME countries rose from 71 to 75 %, the share of the total population in the CEE countries decreased from 9 to 8 %, and the share of the total population in the FSU countries declined from 20 to 17 % (see Appendix 3 for the country-specific population sizes).

The use of the HMD ensures that the mortality data are of high quality up to advanced ages, which is important due to the substantial contributions of old-age mortality to the changes in longevity during the period under study.

Methods

Measuring inter-country disparity

Consider a set of populations i (i = 1, 2,…, n) split by age x. In year t, this set can be described by the matrix of death rates M(t) = [m x,i (t)] and by the matrix of population exposures P(t) = [p x,i (t)]. In matrix M(t), each column M i (t) is a vector of age-specific death rates in country i at time t. The corresponding life expectancy at birth e 0, as a function of vector M i (t), is denoted as e 0(M i (t)). This function produces e 0 values through computing life tables from death rates m x,i (t).

The aggregate scalar index for expressing central tendency in a group of n populations at time t is the population-weighted average length of life:

$$ \overline{e_0(t)}={\displaystyle {\sum}_i{\pi}_i(t){e}_0\left({M}_i(t)\right),} $$
(1)

with the population weights being \( {\pi}_i(t)=\frac{{\displaystyle {\sum}_x}{p}_{x,i}(t)}{{\displaystyle {\sum}_i}{\displaystyle {\sum}_x}{p}_{x,i}(t)}. \)

Following Edwards [7], the population-weighted cross-country variance is used as an aggregate index of the lifetime disparity across countries:

$$ Var(t)={\displaystyle {\sum}_i{\pi}_i(t){\left[{e}_0\left({M}_i(t)\right)-\overline{e_0(t)}\right]}^2}. $$
(2)

This dispersion measure closely correlates with other measures of inequality, and has the advantage of being analytically decomposable into within- and between-group partitions. The split of the total variance into variance within and between country groups will be used in our study. While the total variance measures the amount of dispersion due to all potential factors, the between-group variance measures the contribution of the factors used to define the groups. The within-group variance measures dispersion caused by factors acting within the groups.

The units of variance are years squared, unlike life expectancy which is in units of years. For this reason, we also calculate the standard deviation as a measure of disparity measured in units of years:

$$ StD\ (t)={\left(Var(t)\right)}^{1/2}. $$
(3)

Decomposition problem and its solution

Consider an aggregated measure F (equal to \( \overline{e_0} \) or Var or StD or another index computed from the length of life distribution across countries) defined (according to equations (1) to (3)) as a function of matrices M and P, with its values changing between times t 0 and T. The decomposition task is to compute additive components of the total change F(T)-F(t 0) produced by age-specific changes in countries’ mortality rates (M-effects) and in countries’ population weights (P-effects).

The conventional decomposition equations which were (independently from each other) deduced in the 1980s by Andreev [31], Arriaga [32], and Pressat [33] cannot be used for completing this decomposition task, since these equations decompose a change only in the life expectancy and only in a single population between two different time points or between a pair of populations at one time point. However, the general stepwise replacement algorithm can be employed for completing the task. This method can be used for decompositions involving various output indexes calculated from data on more than one population (see Appendix 4 for a summary of the algorithm).

In our earlier work, we used the stepwise replacement algorithm for decomposition of changes in life expectancy of the total population between two time points into age-specific contributions of: 1) mortality within educational groups and 2) educational structure of the population [29]. The same method was subsequently applied to decompose changes in the total population’s length of life by occupational and marital status groups [3436]. In all these earlier decompositions, the replacement was running in the ascending order across ages, but within each age all possible replacement sequences (related to educational, occupational, or marital status groups) were realized followed by averaging components corresponding to these sequences.

In the present study, the stepwise replacement algorithm is used for decomposition of temporal change in measure F by: 1) age, mortality, and population composition; 2) by age and country group. Applying the method to the set of populations, the age components and their mortality and population-composition parts (M- and P-effects) can be obtained by running a sequence of replacements of the elements m x,i (t 0) and p x,i (t 0) by the elements m x,i (T) and p x,i (T). In matrices M and P, the replacement progresses from the youngest age zero to the oldest age 95+. Let M [x](t 0, T) and P [x](t 0, T) be matrices containing elements m y,i (T) and p y,i (T) at ages 0 ≤ y < x and elements m y,i (t 0) and p y,i (t 0) at ages ≥ x:

$$ {\mathbf{M}}^{\left[\boldsymbol{x}\right]}\left({t}_0,T\right)=\left[\begin{array}{cccccc}\hfill {}_1{m}_{0,1}(T)\hfill & \hfill {}_1{m}_{0,2}(T)\hfill & \hfill \cdots \hfill & \hfill {}_1{m}_{0,i}(T)\hfill & \hfill \cdots \hfill & \hfill {}_1{m}_{0,n}(T)\hfill \\ {}\hfill \underset{\cdots }{{}_4{m}_{1,1}(T)}\hfill & \hfill \underset{\cdots }{{}_4{m}_{1,2}(T)}\hfill & \hfill \cdots \hfill & \hfill \underset{\cdots }{{}_4{m}_{1,i}(T)}\hfill & \hfill \cdots \hfill & \hfill \underset{\cdots }{{}_4{m}_{1,n}(T)}\hfill \\ {}\hfill {}_5{m}_{x-5,1}(T)\hfill & \hfill {}_5{m}_{x-5,2}(T)\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill {}_5{m}_{x-5,i}(T)\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill {}_5{m}_{x-5,n}(T)\hfill \\ {}\hfill \underset{\cdots }{{}_5{m}_{x,1}\left({t}_0\right)}\hfill & \hfill \underset{\cdots }{{}_5{m}_{x,2}\left({t}_0\right)}\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill \underset{\cdots }{{}_5{m}_{x,i}\left({t}_0\right)}\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill \underset{\cdots }{{}_5{m}_{x,n}\left({t}_0\right)}\hfill \\ {}\hfill {}_{\infty }{m}_{95,1}\left({t}_0\right)\hfill & \hfill {}_{\infty }{m}_{95,2}\left({t}_0\right)\hfill & \hfill \cdots \hfill & \hfill {}_{\infty }{m}_{95,i}\left({t}_0\right)\hfill & \hfill \cdots \hfill & \hfill {}_{\infty }{m}_{95,n}\left({t}_0\right)\hfill \end{array}\right], $$
(4)
$$ {\mathbf{M}}^{\left[0\right]}\left({t}_0,T\right)=\mathbf{M}\left({t}_0\right), $$
$$ {\mathbf{P}}^{\left[\boldsymbol{x}\right]}\left({t}_0,T\right)=\left[\begin{array}{cccccc}\hfill {}_1{p}_{0,1}(T)\hfill & \hfill {}_1{p}_{0,2}(T)\hfill & \hfill \cdots \hfill & \hfill {}_1{p}_{0,i}(T)\hfill & \hfill \cdots \hfill & \hfill {}_1{p}_{0,n}(T)\hfill \\ {}\hfill \underset{\cdots }{{}_4{p}_{1,1}(T)}\hfill & \hfill \underset{\cdots }{{}_4{p}_{1,2}(T)}\hfill & \hfill \cdots \hfill & \hfill \underset{\cdots }{{}_4{p}_{1,i}(T)}\hfill & \hfill \cdots \hfill & \hfill \underset{\cdots }{{}_4{p}_{1,n}(T)}\hfill \\ {}\hfill {}_5{p}_{x-5,1}(T)\hfill & \hfill {}_5{p}_{x-5,2}(T)\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill {}_5{p}_{x-5,i}(T)\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill {}_5{p}_{x-5,n}(T)\hfill \\ {}\hfill \underset{\cdots }{{}_5{p}_{x,1}\left({t}_0\right)}\hfill & \hfill \underset{\cdots }{{}_5{p}_{x,2}\left({t}_0\right)}\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill \underset{\cdots }{{}_5{p}_{x,i}\left({t}_0\right)}\hfill & \hfill \begin{array}{l}\cdots \\ {}\cdots \end{array}\hfill & \hfill \underset{\cdots }{{}_5{p}_{x,n}\left({t}_0\right)}\hfill \\ {}\hfill {}_{\infty }{p}_{95,1}\left({t}_0\right)\hfill & \hfill {}_{\infty }{p}_{95,2}\left({t}_0\right)\hfill & \hfill \cdots \hfill & \hfill {}_{\infty }{p}_{95,i}\left({t}_0\right)\hfill & \hfill \cdots \hfill & \hfill {}_{\infty }{p}_{95,n}\left({t}_0\right)\hfill \end{array}\right]. $$
(5)
$$ {\mathbf{P}}^{\left[0\right]}\left({t}_0,T\right)=\mathbf{P}\left({t}_0\right)\ . $$

One step in the replacement sequence includes replacement pertaining to one elementary age group [x, x + a). The rows corresponding to this age group with elements m x,i (t 0) and p x,i (t 0) in the two matrices M and P should be replaced by respective elements m x,i (T) and p x,i (T) and the effect of these replacement on the value of F should be computed. Using the notation given in equations (4) and (5), the respective M- and P-effects are:

$$ {\varDelta}_M^x\left({t}_0,T\right)=\frac{1}{2}\left\{\left[F\left({\mathbf{M}}^{\left[x+a\right]},{\mathbf{P}}^{\left[x\right]}\right)-F\left({\mathbf{M}}^{\left[x\right]},{\mathbf{P}}^{\left[x\right]}\right)\right]+\left[F\left({\mathbf{M}}^{\left[x+a\right]},{\mathbf{P}}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}^{\left[x\right]},{\mathbf{P}}^{\left[x+a\right]}\right)\right]\right\}, $$
(6)
$$ {\varDelta}_P^x\left({t}_0,T\right)=\frac{1}{2}\left\{\left[F\left({\mathbf{M}}^{\left[x\right]},{\mathbf{P}}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}^{\left[x\right]},{\mathbf{P}}^{\left[x\right]}\right)\right]+\left[F\left({\mathbf{M}}^{\left[x+a\right]},{\mathbf{P}}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}^{\left[x+a\right]},{\mathbf{P}}^{\left[x\right]}\right)\right]\right\}. $$
(7)

The component corresponding to change in the elementary age group [x, x + a] is a sum of the M- and P-effects:

$$ {\varDelta}^x\left({t}_0,T\right)={\varDelta}_M^x\left({t}_0,T\right)+{\varDelta}_P^x\left({t}_0,T\right). $$
(8)

Finally, the total change in the function F between times t 0 and T is:

$$ F(T)-F\left({t}_0\right)={\displaystyle {\sum}_x{\varDelta}^x}\left({t}_0,T\right). $$
(9)

The approach to the decomposition by age and country (or country group) is similar to the one used for decomposition by age, mortality, and population composition (equations (4)-(9)). Again, one has to carry out a sequence of replacements of rows (ages) in matrices M and P. However, instead of replacing entire rows in the two matrices, it would be necessary to replace parts of these rows corresponding to certain countries or country groups.

Let us consider a super-simple case with two populations and one age group only. Accordingly, the function F(t) depends on four elements: m 1(t), m 2(t) and p 1(t), p 2(t). Each elementary age component ∆ has to be presented as a sum of two country components:

$$ F(T)-F\left({t}_0\right)={\varDelta}_1+{\varDelta}_2. $$
(10)

The component produced by population 1 is to be computed as an effect of replacements m 1(T) → m 1(t) , p 1(T) → p 1(t). Here, we should take into account two possible replacement sequences:

$$ {\varDelta}_{1(1)}=F\left({m}_1(T),\ {m}_2(t),\ {p}_1(T),\ {p}_2(t)\right)-F\left({m}_1(t),\ {m}_2(t),\ {p}_1(t),\ {p}_2(t)\right); $$
(11)
$$ {\varDelta}_{1(2)}=F\left({m}_1(T),\ {m}_2(T),\ {p}_1(T),\ {p}_2(T)\right)-F\left({m}_1(t),\ {m}_2(T),\ {p}_1(t),\ {p}_2(T)\right). $$
(12)

The final component produced by country 1 is the average of components produced by the two sequences:

$$ {\varDelta}_1=\frac{1}{2}\left({\varDelta}_{1(1)}+{\varDelta}_{1(2)}\right) $$
(13)

Accordingly, the component produced by population 2 is calculated as an effect of replacements m 2(T) →m 2(t) , p 2(T) → p 2(t) on function F. These effects are calculated similarly to equations (10)-(13).

For the three country-groups within each age, the group-specific components have to be computed for all replacement sequences and the results have to be averaged over these sequences. A more systematic description of the method for the decomposition by age and country or country group is given in Appendix 5.

Results

Life expectancy trends and cross-country disparity

Figure 1 shows the time trends in life expectancy at birth for the 36 developed countries included, the three country groups, and the aggregate of all 36 countries combined. Table 1 summarizes these trends at a level of country groups. During the period 1970–2010, the overall population-weighted average life expectancy increased at an annual average rate of 0.21 and 0.17 years for males and females, respectively. However, there was substantial variation across the three groups. Robust and sustained improvements were seen in the EMEs (0.24 and 0.20 years annually for males and females, respectively), while less favorable changes were observed in the CEE (0.15 and 0.18 years) and the FSU (−0.02 and 0.03 years) groups. In fact, the group-level averages for the CEE and (some of) the FSU countries obscures the fact that for much of the period, life expectancy was deteriorating in the FSU group, and was stagnating in the CEE group (see Appendix 2 and Appendix 3 for country-specific life expectancies and population weights in 1970, 1984, 1994, 2004, and 2010, respectively). Between 1970 and 2010, the range (maximum minus minimum) across the entire set of countries increased by 8.0 years for males and 4.4 years for females, and the StD nearly tripled for both sexes.

Fig. 1
figure 1

Life expectancy at birth for countries and country groups by sex in 1970–2010. Note: max-min range in life expectancy is shown above the trends; population-weighted standard deviation in life expectancy is shown below the trends

Table 1 Life expectancy at birth and measures of variance for the entire set of countries and for country groups in selected years (in years)

As the direction and the magnitude of the variance changes varied with time, the observation period was divided into sub-periods according to the character of the changes. While the StD increased continuously and gradually between 1970 and 1984, it fluctuated between 1985 and 2004, decreasing briefly around 1986 and then rising to the highest levels observed in 1994. Between 1995 and 2004, the StD continued to fluctuate, albeit to a lesser degree, and ended the period slightly higher. Finally, during the last period of 2005–2010, the StD decreased somewhat.

These changes in the StD were largely attributable to the high degree of variability in life expectancy in the FSU, and especially in Russia, the largest of the FSU countries. In Russia, the first period of 1970–1984 was characterized by a gradual deterioration in health [37]. The second period of 1985–1994 began with a sudden decrease in mortality associated with Gorbachev’s anti-alcohol campaign followed by a sharp rise in the early 1990s, which coincided with the resumption of the widespread availability of alcohol and the implementation of painful political and economic reforms [38]. The third period of 1995–2004 included another episode of recovery followed by further deterioration in life expectancy after the economic crisis of 1998. In 2005–2010, life expectancy in Russia and the other FSU countries increased substantially, a development which has been attributed to economic growth, the reduction of harm from alcohol, and improvements in the health care system [39].

In the CEE countries, life expectancy stagnation during the communist era in the 1970s and 1980s was followed by significant improvement between 1994 and 2004. In 2005–2010, the progress in life expectancy in the CEE continued but at a somewhat slower pace.

Life expectancy improved far more in the EME group than in the FSU and the CEE groups. Between 1970 and 2010, the gap between the EME and FSU countries increased from 3.9 to 14.4 years for males and from 0.9 years to 8.2 years for females. However, a more detailed look at the life expectancy dynamics across time and individual countries within the EME group reveals some variation in the magnitude of lifetime gains, with slower progress in the United States, the largest country in the group.

Table 1 suggests that the country groups experienced quite different patterns of dispersion change. Within the EME group, there was an important difference between males and females. While for males StD did not change much, for females it doubled and in the 2000s substantially exceeded the corresponding male values. Within the CEE group, StD increased steeply and continuously, with values for males always exceeding those for females. From 1970 to the 2000s it doubled for males and nearly quadrupled for females. Looking at country data in Appendix 2 and Appendix 3, it is apparent that there was a growing contrast between the group leader (East Germany), the group mainstream (Czech Republic, Poland, and Slovakia), and the countries who lagged behind (Bulgaria and Hungary). Finally, males of the FSU group experienced the greatest StD variation across time with a steep rise from 1970 to the mid-1990s and a sharp drop in the 2000s. Among the FSU females, StD followed a broadly similar trajectory. Although FSU showed the largest gap between male and female StDs, females were broadly similar to males with respect to temporal changes. It is surprising that in spite of the StD elevation in the 1970s–1990s, in 2010 the male StD was lower and female StD did not differ much from the starting levels of 1970. This is attributable to the fact that Ukraine and Belarus have been losing their starting life expectancy advantage relative to Russia.

Importantly, StD values within the three country groups were very much lower than the overall StD across all countries. This suggests a great role of the between-group lifetime variance.

Between- and within-group variance

Figure 2 further highlights the sharp increase in inter-country disparities from 1970 to 2004, which was interrupted twice by short episodes of increasing length of life in the FSU countries in the mid-1980s and the mid-1990s. It is also clear that this increase was almost entirely determined by the between-group variance. Whereas in the early 1970s there were moderate differences in the mortality levels of the three groups, in later decades these differences became much more pronounced, especially among men. Over the same period, however, the within-group variance was relatively stable. This suggests that mortality polarization linked to the former East–West geopolitical divide played an important role. The degree of polarization weakened slightly during the last five years of observation, but in 2010 it was still very much higher than in 1970.

Fig. 2
figure 2

Between- and within-group components of life expectancy variance by sex in 1970–2010

Table 2 provides further insights into the patterns of the between- and within-group variance. Among males the FSU group accounted for most (70–80 %) of the between-group variance over the whole observation period. Among females the FSU contribution to the between-group variance rose from around 20 % up to 70 % or more by the 1990s. Due to female life expectancy in the CEE countries being particularly low at the start of the period, the CEE contribution to the between-variance exceeded 50 % in the early 1970s but declined to about 2 % in 2010.

Table 2 Between- and within-group variance and its distribution by country groups in selected years

When we look at the total within-group disparity in life expectancy (second half of Table 2), a few interesting features stand out. First, the values of the within-group variance did not differ much between males and females. In fact, the female within-group variance exceeded the corresponding male values from the mid-1990s onward. Second, the EME group accounted for a large part of the total within-group variance, especially among females.

Components of the variance change

This section presents the results of decompositions of the change in the population-weighted StD. As explained in the Methods, a change in an aggregate measure in a set of countries depends on population (P-effects) as well as mortality (M-effects). While changes in countries’ population weights were usually minor, mortality changes were more substantial and also varied considerably across countries. Thus one can expect that the M-effects would be much greater than P-effects, as confirmed in Table 3. Only after the mid-1990s is it possible to see some P-effects of the overall and of the between-group StD changes among males. These effects were predominantly negative as a result of the redistribution of the population in favor of countries with lower mortality (the EME countries). In the within-group StD, P-effects of changes were negligible.

Table 3 Contributions of population composition (P-effects) and mortality (M-effects) to changes in standard deviation by time periods (in years)

We now focus on the decompositions of the StD changes by age group and by country group (Tables 4 and 5). The age components are aggregated into three broad intervals: childhood ages 0–14; working ages 15–64; and older ages 65+. Table 4 presents the age components of the changes in the StD across the whole set of countries, as well as in the between-group StD and the within-group StD.

Table 4 Age components of the total, between-group, and within-group standard deviation change by time periods (in years)
Table 5 Country-group components of changes in standard deviation by time period (in years)

Decompositions of the overall and the between-group disparity measures in Table 4 suggest that among males the sharp rise in the degree of divergence between 1970 and 1994 was largely determined by increases in the between-group mortality differences at working ages (15–64) followed by older ages (65+). Among females, the older age group contributed more than the working age group to the between-group mortality differences. Among both males and females, contributions of the mortality changes were much smaller for the childhood age group, but they were also positive (pro-divergence). In 1995–2004, the between-group mortality differences and the overall increases in the StD were produced by the older age group and (to a lesser extent) by the working age group. The childhood age group made small and negative (pro-convergence) contributions to the between-group mortality differences and to the overall StD.

After 2004, the overall and the between-group StD values declined. The decompositions show that in this period the StD decrease was produced by the working age group and (to a lesser extent) the childhood age group. Among males, the mortality changes at older ages worked against convergence; while among females, the respective components were low.

Compared to the changes in the between-group disparity, the changes in the within-group disparity were smaller, especially in 1970–1994, when they were 10 times smaller. The increases in the between-group StD coincided with decreases in the within-group StD in 1970–1984 among males and females, and in 1984–1994 among males. In 1970–1984, the within-group convergence was produced by the childhood age group (females and males) and the working age group (males only). In 1985–1994, the within-group convergence among males was determined by the working age group. Between 1984 and 2004 the within-group StD was nearly unchanged among males, but the within-group disparities increased considerably among females. After 2004, the within-group StD declined somewhat due to mortality changes at working and older ages.

Table 5 shows an unusual decomposition of changes in the overall, between-group, and within group disparities by country group. It suggests that the overall and the between-group dispersion changes (upper part of the table) were largely driven by EME and FSU groups. The part of the CEE group was much smaller due to its lower population weight and smaller temporal changes in life expectancy. Steep life expectancy increase in the EME group leading to widening of the gap between this group and FSU was contributing to the life expectancy divergence throughout the entire period. In the 1970s and especially in the early 1990s, deterioration in FSU largely contributed to the life expectancy divergence. In 2005–2010, life expectancy gains in the FSU produced negative (pro-convergence) contributions to the overall and the between-group StD change.

Although changes in the within-group dispersion were moderate, there was still a substantial rise in the female within-group disparity from 1985 to 2004, which was largely caused by the growing differences within the EME group.

Discussion

Although worldwide the average life expectancy has been increasing steadily for many decades, this positive trend has shown heterogeneity across countries. This variation can be investigated using the convergence-divergence framework, which helps to explain interchangeable stages of mortality divergence and convergence [28]. The present study looked solely at developed countries to avoid mixing populations with dramatically different mortality patterns [40].

Our aim was to examine changes in the amount of inter-country lifetime disparity over the last four decades and to assess relative importance of the former East–West geopolitical divide and of temporal changes in the mortality distribution by age and country-group for changes in the amount of disparity. To do this we adapted the general stepwise replacement algorithm for linking changes in countries’ age-specific mortality rates with changes in cross-country variance. This method is appropriate in an analysis of temporal changes in an aggregate demographic or public health measure based on a set of populations, rather than on a single population. While the conventional decomposition method permits one to assess the effects of age-specific mortality changes on the total change in life expectancy at birth in a single population, the method used here allowed us to assess the effects of simultaneous changes in numerous country-age-specific mortality rates on the total change in an aggregate scalar measure (the standard deviation of the length of life). This elaborated method is generally applicable to the decomposition of temporal changes in any aggregate index based on a set of populations.

Limitations

The set of countries under study did not include all developed countries as defined by the UN classification. Although 12 developed countries (Andorra, Bosnia and Herzegovina, Faeroe Islands, Greece, Israel, Macedonia, Lichtenstein, Moldova, Romania, Slovenia, San Marino, Serbia) were missing, there are two reasons why this omission is unlikely to have biased our findings. First, the life expectancy estimates for every country of the world by the UN Population Division signal that the life expectancy values in all of the missing countries are close to the group-specific average life expectancy levels either in the EME or the CEE groups [41]. Second, in 2010 the missing countries had a total population of 62 million, which constitutes only 4.8 % of the total population of all developed countries.

Although we recognize that the usefulness of decomposition analysis is higher if it includes decomposition by causes of death, inclusion of this additional dimension would lead to necessity to carry out first another study to address the cause-of-death inconsistency across time (accounting for changes in ICD revisions) and other changes in diagnostics and coding practices [42, 43]. Therefore, as a first step in using these elaborated methods it was decided to focus on the decomposition for all causes combined.

Principal findings

Our substantive results show the striking growth in life expectancy variance across the developed countries from the early 1970s to the early 2000s with the variance being much greater for males than for females. Despite the convergence in the period 2005–2010, life expectancy disparities across the developed countries are still enormous with the cross-country StD values being three times higher in 2010 than in 1970.

The rising between-group variance constitutes a core part of this phenomenon. Even 25 years after the fall of the Berlin Wall, a large fraction of the variance of the length of life is attributable to the East–West divide. This historical geopolitical gap between the two parts of the developed world grossly determines the life expectancy of the people living in respective countries. The influence of the politics of the 20th century on survival and longevity [12] has thus persisted far longer than was anticipated 25 years ago.

The overall and the between-group lifetime variance are much higher among males than among females. Rise in between-group disparity was mainly fueled by widening of the East–West mortality gap at ages 15 to 64 years among males and at ages 15+ years among females. Both the continuous progress in EME and the inconsistent and largely unfavorable changes in FSU contributed to the disparity. The remarkable mortality excess in adult age mortality (especially among males) in FSU and CEE compared to EME is in line with findings of earlier research on Russia and other ex-communist countries [10, 44].

While a large part of the between-group variance was determined by the FSU countries, the within-group variance was dominated by the EME. Unlike the between-group variance, the within-group variance and its changes were greater among females than they were among males starting in the late 1990s. The onset of this pattern coincided with the time when nearly all female mortality was already concentrated at old ages and reduction of old age mortality became an ultimate condition of further longevity progress. It appears that even across the EME countries there are substantial differences with respect to success in reduction of female mortality at old ages. This was examined by Meslé and Vallin, [24] who provided detailed comparison of components and possible drivers of old-age mortality trends in several advanced countries. According to their convergence-divergence theory, the growing disparity in longevity among women should be related to emergence of a new health challenge. This time, the challenge is connected to a hard transition from reduction of cardiovascular death at younger old ages to reduction of death from multiple pathologies at advanced ages.

Among males, the moderate changes in the within-group variance were mostly related to the FSU group and to mortality at working ages.

In the second half of the 2000s, both the between-group and the within-group parts of the inter-country variance declined due to mortality convergence at ages 15 to 64 in the FSU group and (to a much smaller extent) to mortality convergence at ages 65+ in the EME group. It is worth noting that the life expectancy gap between the CEE and the EME countries, which was shrinking over the 1990s as mortality in the CEE countries steeply declined, stabilized in the 2000s due to a slowdown in the positive trends in the CEE group.

Sensitivity analyses

The population-weighted variance metric used in this study is focused on individuals since it accounts for numbers of humans who are exposed to higher or lower death hazard in different places. There is, however, a disadvantage due to a low sensitivity of the metrics to mortality and mortality changes in countries with small populations. In addition, one may be particularly interested in to what extent results of our analysis depend on the USA, a country with the largest population size, about 300 million, which is somewhat lagging behind within the group of EME countries over the last two decades.

To evaluate how influential the very large contribution of the USA was to our findings we repeated all computations with this country excluded (outcomes not shown here but can be provided upon request). All in all, this change produces a minor impact on the results pertaining to the overall lifetime disparity and its changes. Average life expectancy in the EME group becomes slightly higher and its population weight becomes slightly lower, which results in slightly lower between-group and overall standard deviations. More important impacts are seen in the within-group standard deviation which becomes lower (especially for females) by about one-third in the 1990s and the 2000s. This happens mostly due to a diminished weight of the EME group in the within-group disparity. The temporal increase in the within-group disparity among females becomes somewhat smaller, contribution of ages 65+ diminishes by about 40 %, and contribution of EME to this increase diminishes by about one-fourth.

Use of the population-unweighted data produces a more visible change in the calculation outcomes (see Appendix 6 for detailed tables). This metric expresses the amount of inter-country difference in length of life among countries irrespective of their population sizes, counting each country as one unit.

With unweighted data, life expectancy of EME slightly increases and life expectancy of FSU increases substantially – by 2.5–3.5 years in the 1990s–2000s. The overall StD becomes substantially lower and its increase from 1970 to 2004 and decrease from 2004 to 2010 are becoming much less pronounced. All age components of the disparity changes are becoming substantially smaller. Contributions of ages 15 to 64 are still decisive for the overall and the between-group variance changes, especially among males, as well as contributions of ages 65+ to the within-group variance changes among females. Although contributions of the CEE group to the disparity changes increase, they are still substantially lower compared to FSU and EME.

Conclusions

Over most of the period 1970–2010 there was an overall tendency for life expectancy to diverge across the 36 developed countries examined. This was driven primarily by the growing gap between the EME and FSU groups due to diverse changes in mortality at working ages (especially among males) and at older ages. Within the EME group, divergence occurred in the 1990s and the early 2000s due to uneven progress with respect to female mortality at ages 65+. The within-group lifetime disparity rise among women was substantially dependent on the USA. Old-age and working-age mortality rates are still substantially higher in some places than in others, generating important variations in length of life across developed countries. This signals that many countries have significant scope for further health improvement, even at the existing levels of medical technology and economic wealth.

The inter-country mortality disparities point to the major health challenges, which even many developed countries continue to face. These disparities still persist in large part because four decades ago countries had moderately differing levels of length of life but have subsequently shown very variable capacity to address major health challenges [37]. Further monitoring and analysis of cross-country variation in longevity and survival and identification of the factors associated with success or failure in the fight to extend longevity should be a priority.

Notes

  1. We selected all HMD countries classified as developed by the UN for which data were available at least since 1970. Data series for Slovenia and Israel were excluded since they begin after 1970. For thirteen populations which had shorter data series, 2009 (twelve populations) or 2008 (one population) was used as the last data point. For more details, see Appendices 2 and 3.

References

  1. Oeppen J, Vaupel J. Broken limits to life expectancy. Science. 2002;296(5570):1029–31.

    Article  CAS  PubMed  Google Scholar 

  2. Wilmoth JR. Increase of human longevity: past, present, and future. Jpn J Popul. 2004;9:1.

    Google Scholar 

  3. Christensen K, Doblhammer G, Rau R, Vaupel JW. Ageing populations: the challenges ahead. Lancet. 2009;374:1196–208.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Moser K, Shkolnikov V, Leon DA. World mortality 1950–2000: divergence replaces convergence from the late 1980s. Bull World Health Organ. 2005;83(3):202–9.

    PubMed  PubMed Central  Google Scholar 

  5. Vallin J, Meslé F, Valkonen T. Trends in mortality and differential mortality. Popul Stud. 2001;36:332.

    Google Scholar 

  6. Meslé F, Vallin J. Historical trends in mortality. In: Rogers RG, Crimmins EM, eds. International Handbook of Adult Mortality. London New York: Springer Dordrecht Heidelberg; 2011:9-48. doi:10.1007/978-90-481-9996-9_2.

  7. Edwards R. Changes in world inequality in length of life: 1970–2000. Popul Dev Rev. 2011;37(3):499–528.

    Article  PubMed  Google Scholar 

  8. Smits J, Monden C. Length of life inequality around the globe. Soc Sci Med. 2009;68:1114–23.

    Article  PubMed  Google Scholar 

  9. McMichael AJ, McKee M, Shkolnikov V, Valkonen T. Mortality trends and setbacks: global convergence or divergence? Lancet. 2004;363(9415):1155–9.

    Article  PubMed  Google Scholar 

  10. Meslé F. Mortality in Eastern and Western Europe: a widening gap. In: Coleman D, editor. Europe’s Population in the 1990s. Oxford: Oxford University Press; 1996. p. 127–43.

    Google Scholar 

  11. Shkolnikov VM, Andreev EM, McKee M, Leon DA. Components and possible determinants of the decrease in Russian mortality in 2004–2010. Demogr Res. 2013;28(32):917–50.

    Article  Google Scholar 

  12. Mackenbach JP. Convergence and divergence of life expectancy in Europe: a centennial view. Eur J Epidemiol. 2013;28:229–40.

    Article  PubMed  Google Scholar 

  13. Meslé F. Mortality in central and eastern Europe: long-term trends and recent upturns. Demogr Res. 2004;3:45–70. Special Collection 2.

    Article  Google Scholar 

  14. Luy M, Wegner C, Lutz W. Adult mortality in Europe. In: Rogers RG, Crimmins EM, eds. International Handbook of Adult Mortality. London New York: Springer Dordrecht Heidelberg; 2011:49–81. doi:10.1007/978-90-481-9996-9_2.

  15. Leon DA. Trends in European life expectancy: a salutary view. Int J Epidemiol. 2011;40(2):271–7.

    Article  PubMed  Google Scholar 

  16. Tuljapurkar S, Li N, Boe C. A universal pattern of mortality decline in the G7 countries. Nature. 2000;405:789–92.

    Article  CAS  PubMed  Google Scholar 

  17. Horiuchi S. Epidemiological transition in developed countries: past, present, and future. The mortality, epidemiological and health transitions: their relevance for the study of health and mortality. In: Health and mortality issues of global concern. Proceedings of the Symposium on Health and Mortality. Brussels;1997: 54–71.

  18. Murphy M. Adult mortality in the former Soviet Union. In: Rogers RG, Crimmins EM, eds. International Handbook of Adult Mortality. London New York: Springer Dordrecht Heidelberg; 2011:83–100. doi:10.1007/978-90-481-9996-9_2.

  19. Shkolnikov VM, Andreev EM, Zhang Z, Oeppen J, Vaupel JW. Losses of expected lifetime in the United States and other developed countries: methods and empirical analyses. Demography. 2011;48:211–39.

    Article  PubMed  Google Scholar 

  20. Santow G. The mortality, epidemiological and health transitions: their relevance for the study of health and mortality. In: Health and Mortality Issues of Global Concern. Proceedings of the Symposium on Health and Mortality. Brussels; 1997: 39–53

  21. Lopez AD, Gaselli G, Valkonen T. Adult Mortality in Developed Countries: From Description to Explanation. Oxford: Clarendon; 1995. p. 361.

    Google Scholar 

  22. Canudas-Romo V. Three measures of longevity: time trends and record values. Demography. 2010;47(2):299–312.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Staetsky L. Diverging trends in female old-age mortality: a reappraisal. Demogr Res. 2009;21:885–914.

    Article  Google Scholar 

  24. Meslé F, Vallin J. Diverging trends in female old-age mortality: the United States and the Netherlands versus France and Japan. Popul Dev Rev. 2006;32(1):123–45.

    Article  Google Scholar 

  25. Wilson C. On the scale of global demographic convergence 1950–2000. Popul Demographic Rev. 2001;27(1):155–71.

    Article  CAS  Google Scholar 

  26. Caselli G, Meslé F, Vallin J. Epidemiologic transition theory exceptions. Genus. 2002;58(1):9–52.

    Google Scholar 

  27. Meslé F, Vallin J. The health transition: trends and prospects. In: Caselli G, Vallin J, Wunsch G, editors. Demography: Analysis and Synthesis. A Treatise in Population Studies, vol. II. New York: Academic; 2006. p. 29–44.

    Google Scholar 

  28. Vallin J, Meslé F. Convergences and divergences in mortality. A new approach to health transition. Demogr Res. 2004;2:1–44. Special Collection.

    Google Scholar 

  29. Shkolnikov VM, Begun AZ. Algorithm for decomposition of differences between aggregate demographic measures and its application to life expectancies, healthy life expectancies, parity-progression ratios and total fertility rates. Demogr Res. 2002;7(14):499–522.

    Google Scholar 

  30. Human Mortality Database. 2015. http://www.mortality.org. Accessed 21 Feb 2015.

  31. Andreev EM. Metod komponent v analize prodoljitelnosty zhizni [The method of components in the analysis of length of life]. Vestnik Statistiki. 1982;9:42–7.

    Google Scholar 

  32. Arriaga E. Measuring and explaining the change in life expectancies. Demography. 1984;21(1):83–96.

    Article  CAS  PubMed  Google Scholar 

  33. Pressat R. Contribution des écarts de mortalité par âge à la différance des vies moyennes. Population. 1985;4–5:766–70.

    Article  Google Scholar 

  34. Shkolnikov VM, Andreev EM, Jasilionis D, Leinsalu M, Antonova OI, McKee M. The changing relation between education and life expectancy in central and eastern Europe in the 1990s. J Epidemiol Community Health. 2006;60:875–81.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Andreev E, Hoffmann R, Carlson E, Shkolnikov V, Kharkova TL. Concentration of working-age male mortality among manual workers in urban Latvia and Russia, 1970–1989. Eur Soc. 2009;11(1):161–85.

    Article  Google Scholar 

  36. Jasilionis D, Andreev EM, Kharkova TL, Kingkade WW. Change in marital status structure as an obstacle for health improvement: evidence from six developed countries. Eur J Pub Health. 2011;22(4):602–4.

    Article  Google Scholar 

  37. Shkolnikov VM, Leon DA. Commentary: N Eberstadt’s ‘The health crisis in the USSR’ and sustainable mortality reversal in the post-Soviet space during communism and after. Int J Epidemiol. 2006;35(6):1406–9.

    Article  PubMed  Google Scholar 

  38. Shkolnikov VM, Cornia G, Leon D, Meslé F. Causes of the Russian mortality crisis: evidence and interpretations. World Dev. 1998;26:1995–2011.

    Article  Google Scholar 

  39. Grigoriev P, Meslé F, Shkolnikov VM, Andreev E, Fihel A, Pechholdova M, Vallin J. The recent mortality decline in Russia: beginning of the cardiovascular revolution? Popul Dev Rev. 2014;40(1):107–29.

    Article  Google Scholar 

  40. Murray CJL, Lopez A. Evidence-based health policy. Lessons from the global burden of disease study. Science. 1996;274:740–3.

    Article  CAS  PubMed  Google Scholar 

  41. United Nations. World Population Prospects. 2015. URL: https://esa.un.org/unpd/wpp/. Accessed 21 Feb 2015.

  42. Meslé F, Vallin J. The effect of ICD-10 on continuity in cause-of-death statistics. The example of France. Population. 2008;63(2):383–96.

    Article  Google Scholar 

  43. Meslé F, Vallin J. Reconstructing long-term series of causes of death. Historical Methods. 1996;29(2):72–8.

    Article  Google Scholar 

  44. Shkolnikov VM, Andreev EM, Leon DA, McKee M, Mesle F, Vallin J. Mortality reversal in Russia: the story so far. Hygiea Internationalis. 2004;4(4):29–80.

    Article  Google Scholar 

  45. Das Gupta P. Standardisation and decomposing of rates from cross-classified data. Genus. 1994;50(3–4):171–96.

    CAS  PubMed  Google Scholar 

  46. Horiuchi S, Wilmoth JR, Pletcher SD. A decomposition method based on a model of continuous change. Demography. 2008;45(4):785–801.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

This study was carried out within the framework of the projects “From Disparities in Mortality Trends to Future Health Challenges” (DIMOCHA) supported by the Agence national de la recherche (ANR, France) and the Deutsche Forschungsgemeinschaft (DFG, Germany) and “Mortality Divergence and Causes of Death” (MODICOD) supported by the AXA Foundation. Participation of VMS, and DAJ was additionally supported by the Dynasty Foundation (Russia).

Authors’ contributions

ST carried out the analysis including data processing and calculations, and drafted the manuscript. VMS conceived the idea, proposed the analytical technique, took part in interpreting the results and editing the manuscript. DJ aided in conceiving the study and interpreting the results and editing. PG contributed to the idea of the study and aided in drafting and revising the manuscript. DAJ aided in developing the R software for the decomposition and in drafting the methods part. DAL provided important inputs to the idea of the study and critically reviewed draft versions of the text. All authors have read and approved the manuscript.

Competing interests

The authors declare they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Financial disclosure

No external support was received for this study.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sergey Timonin.

Appendices

Appendix 1

Table 6 Countries and country groups under study

Appendix 2

Table 7 Life expectancy at birth in selected years (in years)

Appendix 3

Table 8 Population size in selected years (in millions)

Appendix 4

Brief summary of the general stepwise replacement algorithm

The idea of this method is to present decomposition of the total change in an index f(.) as a sequence of replacement of its arguments. Let f(.) depend on (say) three covariates a, b, and c which are defined at time points 1 and 2. The difference f(a 2,b 2,c 2)-f(a 1,b 1,c 1) can be presented as a sum of effects of a sequence of replacements:

$$ f\left({a}_2,{b}_2,{c}_2\right)-f\left({a}_1,{b}_1,{c}_1\right)=\varDelta f\left({a}_2\to {a}_1,{b}_1,{c}_1\right)+\varDelta f\left({a}_2,{b}_2\to {b}_1,{c}_2\right)+\varDelta f\left({a}_2,{b}_2,{c}_2\to {c}_1\right) $$

with the delta components

$$ \varDelta f\left({a}_2\to {a}_1,{b}_1,{c}_1\right)=f\left({a}_2,{b}_1,{c}_1\right)-f\left({a}_1,{b}_1,{c}_1\right); $$
$$ \varDelta f\left({a}_2,{b}_2\to {b}_1,{c}_1\right)=f\left({a}_2,{b}_2,{c}_1\right)-f\left({a}_2,{b}_1,{c}_1\right); $$
$$ \varDelta f\left({a}_2,{b}_2,{c}_2\to {c}_1\right)=f\left({a}_2,{b}_2,{c}_2\right)-f\left({a}_2,{b}_2,{c}_1\right). $$

Each of these three components is equal to the contribution of the shift of the value of respective covariate from point 1 to point 2. Note that once the value of the independent variable (say a) is shifted from a 1 to a 2, it remains equal to a 2 when delta-components of covariates that come later in the replacement sequence are calculated.

The three components given above correspond to the replacement sequence a-b-c. It is possible, however, to move from a 1 b 1 c 1 to a 2 b 2 c 2 using other pathways. For example, it is possible to carry out the replacement sequence b-a-c (b 2 → b 1, a 2 → a 1, c 2 → c 1). The number of all possible replacement sequences equals 6 (all possible permutations among three elements).

For a non-linear function f(a, b, c), values of the component produced by a movement from point 1 to point 2 for the same covariate in different replacement sequences are not exactly equal to each other. For example, the a-component Δf(a 2 → a 1, b 2, c 1) = f(a 2, b 2, c 1) − f(a 1, b 2, c 1) in sequence b-a-c is not exactly the same compared to the a-component in sequence a-b-c Δf(a 2 → a 1, b 1, c 1) = f(a 2, b 1, c 1) − f(a 1, b 1, c 1).

Thus, the stepwise replacement has to be carried out for all possible replacement sequences (permutations) and the final components are to be computed as averages of the components’ values over all these sequences [29, 45, 46]. In our example it means that every component (a, b, or c) should be calculated as average of six respective sequence-specific component values.

The stepwise replacement algorithm with the full run across all replacement sequences can be applied for decompositions concerning countries or regions as we have done here, or many important socio-demographic variables with small numbers of categories such as sex, cause of death by broad diagnostic groups, education or socioeconomic status by aggregate categories, birth order from 1 to 5+ and other.

However, there is an obvious difficulty when it comes to the age variable, which has high number of categories: about 20 and 100 in abridged and complete life tables, respectively, and about 30 in a fertility tables. Therefore, it was suggested to compute the age components by using the replacement sequence that is running in ascending order of ages [29]. Such an approach guarantees that the stepwise replacement algorithm’s results are exactly equal to results of the most used analytical decomposition formulae by Andreev (1982), Arriaga (1984), and Pressat (1985) in the case of life expectancy.

Appendix 5

Stepwise replacement algorithm for decomposition by age and country group

Similar to decomposition by age, mortality and population composition (M- and P-effects), one can estimate the contribution of a country (country group) to the total change. Following the logic of M-P decomposition, we estimate country-specific contributions within age components. In other words, every age-specific contribution has to be split by country effects.

Below we provide a general formal description for the case of three populations which can be extended to any number of countries. It is important to note here that the calculation of country effects is a resource-consuming task which requires in case of K countries 2K! calculations of the index function F(.) for every age group. Thus, even for the case of K = 10 this task might be not affordable.

Let the index function F(.) depends on matrices M and P. Following the general rule of the age decomposition the replacement in matrices M and P progresses from the youngest age 0 to the oldest age 95+. Let M [x] k (t 0, T) and P [x] k (t 0, T) be matrices containing elements m y,i (T) and p y,i (T) at ages 0 ≤ y < x − a and at age x - a for countries k, elements m y,i (t 0) and p y,i (t 0) at ages y ≥ x and at age x - a for counties other than k. For example, for k = 1,2:

$$ {\mathbf{M}}_{1,2}^{\left[\boldsymbol{x}\right]}\left({t}_0,T\right)=\left[\begin{array}{ccc}\hfill {}_1{m}_{0,1}(T)\hfill & \hfill {}_1{m}_{0,2}(T)\hfill & \hfill {}_1{m}_{0,3}(T)\hfill \\ {}\hfill {}_4{m}_{1,1}(T)\hfill & \hfill {}_4{m}_{1,2}(T)\hfill & \hfill {}_4{m}_{1,3}(T)\hfill \\ {}\hfill \cdots \hfill & \hfill \cdots \hfill & \hfill \cdots \hfill \\ {}\hfill {}_5{m}_{x-5,1}(T)\hfill & \hfill {}_5{m}_{x-5,2}(T)\hfill & \hfill {}_5{m}_{x-5,3}\left({t}_0\right)\hfill \\ {}\hfill {}_5{m}_{x,1}\left({t}_0\right)\hfill & \hfill {}_5{m}_{x,2}\left({t}_0\right)\hfill & \hfill {}_5{m}_{x,3}\left({t}_0\right)\hfill \\ {}\hfill \cdots \hfill & \hfill \cdots \hfill & \hfill \cdots \hfill \\ {}\hfill {}_{\infty }{m}_{95,1}\left({t}_0\right)\hfill & \hfill {}_{\infty }{m}_{95,2}\left({t}_0\right)\hfill & \hfill {}_{\infty }{m}_{95,3}\left({t}_0\right)\hfill \end{array}\right], $$
(14)
$$ {\mathbf{M}}_0^{\left[x\right]}\left({t}_0,T\right)={\mathbf{M}}_{1,2,3}^{\left[x-a\right]}\left({t}_0,T\right)\;\mathrm{and}\ {\mathbf{M}}_0^{\left[0\right]}\left({t}_0,T\right)=\mathbf{M}\left({t}_0\right), $$
$$ {\mathbf{P}}_{1,2}^{\left[\boldsymbol{x}\right]}\left({t}_0,T\right)=\left[\begin{array}{ccc}\hfill \begin{array}{l}{}_1{p}_{0,1}(T)\\ {}{}_4{p}_{1,1}(T)\end{array}\hfill & \hfill \begin{array}{l}{}_1{p}_{0,2}(T)\\ {}{}_4{p}_{1,2}(T)\end{array}\hfill & \hfill \begin{array}{l}{}_1{p}_{0,3}(T)\\ {}{}_4{p}_{1,3}(T)\end{array}\hfill \\ {}\hfill \cdots \hfill & \hfill \cdots \hfill & \hfill \cdots \hfill \\ {}\hfill {}_5{p}_{x-5,1}(T)\hfill & \hfill {}_5{p}_{x-5,2}(T)\hfill & \hfill {}_5{p}_{x-5,3}\left({t}_0\right)\hfill \\ {}\hfill {}_5{p}_{x,1}\left({t}_0\right)\hfill & \hfill {}_5{p}_{x,2}\left({t}_0\right)\hfill & \hfill {}_5{p}_{x,3}\left({t}_0\right)\hfill \\ {}\hfill \cdots \hfill & \hfill \cdots \hfill & \hfill \cdots \hfill \\ {}\hfill {}_{\infty }{p}_{95,2}\left({t}_0\right)\hfill & \hfill {}_{\infty }{p}_{95,2}\left({t}_0\right)\hfill & \hfill {}_{\infty }{p}_{95,3}\left({t}_0\right)\hfill \end{array}\right], $$
(15)
$$ {\mathbf{P}}_0^{\left[x\right]}\left({t}_0,T\right)={\mathbf{P}}_{1,2,3}^{\left[x-a\right]}\left({t}_0,T\right)\;\mathrm{and}\kern0.37em {\mathbf{P}}_0^{\left[0\right]}\left({t}_0,T\right)=\mathbf{P}\left({t}_0\right), $$

One step in the replacement sequence includes replacement pertaining a single age group [x,x + a). The elements corresponding to this age group in the two matrices M and P should be replaced simultaneously for every country k. The effect of this replacement on the value of F should be computed as an average effect of the replacement a country-specific element by all possible permutations. Using the notation given in (14) and (15), the effect of country 1 within the age component x is:

$$ \begin{array}{l}{\varDelta}_1^x\left({t}_0,T\right)=\frac{1}{4}\left\{\left[F\left({\mathbf{M}}_1^{\left[x+a\right]},{\mathbf{P}}_1^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_0^{\left[x+a\right]},{\mathbf{P}}_0^{\left[x+a\right]}\right)\right]\right.\\ {}\kern3.48em +\left[F\left({\mathbf{M}}_{1,2}^{\left[x+a\right]},{\mathbf{P}}_{1,2}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_2^{\left[x+a\right]},{\mathbf{P}}_2^{\left[x+a\right]}\right)\right]\\ {}\kern3.48em +\left[F\left({\mathbf{M}}_{1,3}^{\left[x+a\right]},{\mathbf{P}}_{1,3}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_3^{\left[x+a\right]},{\mathbf{P}}_3^{\left[x+a\right]}\right)\right]\\ {}\kern3.36em \left.+\left[F\left({\mathbf{M}}_{1,2,3,}^{\left[x+a\right]},{\mathbf{P}}_{1,2,3,}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_{2,3,}^{\left[x+n\right]},{\mathbf{P}}_{2,3,}^{\left[x+a\right]}\right)\right]\right\}\end{array} $$
(16)

The effects of second and third countries are:

$$ \begin{array}{l}{\Delta}_2^x\left({t}_0,T\right)=\frac{1}{4\ }\left\{\left[F\left({\mathbf{M}}_2^{\left[x+a\right]},{\mathbf{P}}_2^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_0^{\left[x+a\right]},{\mathbf{P}}_0^{\left[x+a\right]}\right)\right]\right.\\ {}\kern3.72em +\left[F\left({\mathbf{M}}_{1,2}^{\left[x+a\right]},{\mathbf{P}}_{1,2}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_1^{\left[x+a\right]},{\mathbf{P}}_1^{\left[x+a\right]}\right)\right]\\ {}\kern3.72em +\left[F\left({\mathbf{M}}_{2,3}^{\left[x+a\right]},{\mathbf{P}}_{2,3}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_3^{\left[x+a\right]},{\mathbf{P}}_3^{\left[x+a\right]}\right)\right]\;\\ {}\kern3.72em \left.+\left[F\left({\mathbf{M}}_{1,2,3}^{\left[x+a\right]},{\mathbf{P}}_{1,2,3}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_{2,3}^{\left[x+a\right]},{\mathbf{P}}_{2,3}^{\left[x+a\right]}\right)\right]\right\}\end{array} $$
(17)
$$ \begin{array}{l}{\Delta}_3^x\left({t}_0,T\right)=\frac{1}{4\ }\left\{\left[F\left({\mathbf{M}}_3^{\left[x+a\right]},{\mathbf{P}}_3^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_0^{\left[x+a\right]},{\mathbf{P}}_0^{\left[x+a\right]}\right)\right]\right.\\ {}\kern3.6em +\left[F\left({\mathbf{M}}_{2,3}^{\left[x+a\right]},{\mathbf{P}}_{2,3}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_2^{\left[x+a\right]},{\mathbf{P}}_2^{\left[x+a\right]}\right)\right]\\ {}\kern3.6em +\left[F\left({\mathbf{M}}_{1,3}^{\left[x+a\right]},{\mathbf{P}}_{1,3}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_1^{\left[x+a\right]},{\mathbf{P}}_1^{\left[x+a\right]}\right)\right]\\ {}\kern3.6em \left.+\left[F\left({\mathbf{M}}_{1,2,3}^{\left[x+a\right]},{\mathbf{P}}_{1,2,3}^{\left[x+a\right]}\right)-F\left({\mathbf{M}}_{1,2}^{\left[x+a\right]},{\mathbf{P}}_{1,2}^{\left[x+a\right]}\right)\right]\right\}\end{array} $$
(18)

The component corresponding to change in the elementary age group [xx + a) is equal to the sum of the respective country contributions:

$$ {\Delta}^x={\displaystyle \sum_{k=1}^3}{\Delta}_k^x. $$

The total change in the function F between times t 0 and T is:

$$ F(T)-F\left({t}_0\right)={\displaystyle \sum_x}{\Delta}^x. $$

Appendix 6

Alternative tables of results with the population-unweighted measures

Tables 9–13 below correspond to Tables 15 in the main text.

Table 9 Life expectancy at birth and measures of variance for the entire set of countries and for country groups in selected years (in years)
Table 10 Between- and within-group variance and its distribution by country groups in selected years
Table 11 Contributions of population composition (P-effects) and mortality (M-effects) to changes in standard deviation by time periods (in years)
Table 12 Age components of the total, between-group, and within-group standard deviation change by time periods (in years)
Table 13 Country-group components of changes in standard deviation by time period (in years)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Timonin, S., Shkolnikov, V.M., Jasilionis, D. et al. Disparities in length of life across developed countries: measuring and decomposing changes over time within and between country groups. Popul Health Metrics 14, 29 (2016). https://doi.org/10.1186/s12963-016-0094-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12963-016-0094-0

Keywords