AgeSpecific Mortality and Fertility Rates for Probabilistic Population Projections
Abstract
The UN released official probabilistic population projections (PPP) for all
countries for the first time in July 2014. These were obtained by projecting
the period total fertility rate (TFR) and life expectancy at birth ()
using Bayesian hierarchical models, yielding a large set of future
trajectories of TFR and for all countries and future time periods
to 2100, sampled from their joint predictive distribution.
Each trajectory was then converted to agespecific mortality
and fertility rates, and population was projected using the cohortcomponent
method. This yielded a large set of trajectories of future age and
sexspecific population counts and vital rates for all countries.
In this paper we describe the methodology used for deriving the
agespecific mortality and fertility rates in the 2014 PPP,
we identify limitations of these methods, and we propose several
methodological improvements to overcome them.
The methods presented in this paper are implemented in the publicly
available bayesPop R package.
Keywords: Bayesian hierarchical model; Cohortcomponent method; Life expectancy at birth; Markov chain Monte Carlo; Total fertility rate; United Nations; World Population Prospects.
Contents
 1 Introduction
 2 AgeSpecific Mortality Rates for Probabilistic Population Projections
 3 AgeSpecific Fertility Rates for Probabilistic Population Projections
 4 Discussion
 5 Acknowledgements
List of Figures
 1 Original and Coherent Kannisto methods compared: Male and female mortality rates for Brazil and Lithuania
 2 Projected female mortality rates for Bangladesh and Pakistan in 20952100 projected using three different methods for computing
 3 Mortality rates for Bangladesh by time for three different age groups
 4 Rotating the parameter over time: Example from Japan
 5 Probabilistic projection of agespecific mortality rates for Japan and Kazakhstan in 20952100
 6 Joint predictive distribution of mortality rates for females and males for Kazakhstan in 20952100
 7 Joint distribution of mortality rates for females and males for Kazakhstan in 20952100 for individual age groups
 8 Projected agespecific mortality rates for Botswana, a country with a generalized HIV/AIDS epidemic.
 9 Example of projected Mean Age of Childbearing (MAC) for countries in Eastern Asia in WPP 2012.
 10 Trends in Mean Age at Childbearing in countries with the start of Phase III of fertility decline before 2000. Dots mark the time period when the country entered Phase III.
 11 Proportionate agespecific fertility rates (PASFR) by time for age groups 1519, 2529 and 3539 in Niger, Bangladesh and the Czech Republic
 12 Probabilistic projection of agespecific fertility rates for Ethiopia, Nepal and Japan in the time period 20952100
 13 PASFR by age over time for selected countries.
 14 Example of projected MAC for countries in Eastern Asia after applying the proposed methodology.
1 Introduction
The United Nations released official probabilistic population projections for all countries for the first time in July 2014 (Gerland et al., 2014). They were produced by probabilistically projecting the period total fertility rates (TFR) and life expectancies () for all countries using Bayesian hierarchical models (Alkema et al., 2011; Raftery et al., 2013). These probabilistic projections took the form of a large set of trajectories, each of which was sampled from the joint predictive distribution of TFR and female and male for all countries and all future time periods to 2100 using Markov chain Monte Carlo (MCMC) methods ^{1}^{1}1This general approach applies to countries experiencing normal mortality trends. For countries having ever experienced 2 per cent or more adult HIV prevalence during the period 1980 to 2010, all projected trajectories of life expectancy by sex for each of these countries were adjusted in such a way as to ensure that the median trajectory for each country was consistent with the 2012 Revision of the World Population Prospects deterministic projection that incorporates the impact of HIV/AIDS on mortality, as well as assumptions about future potential improvements both in the reduction of the epidemic and survival due to treatment.
For each trajectory, the life expectancies were converted to age and sexspecific mortality rates, and the total fertility rates were converted to agespecific fertility rates. The population was then projected forward using the cohortcomponent method. This yielded a large set of trajectories of population by age and sex, and agespecific fertility and mortality rates, for all countries and future time periods jointly. These were summarized by predictive medians and 80% and 95% prediction intervals for a wide range of population quantities of interest, for all countries and a wide range of regional and other aggregates. They were published as the UN’s Probabilistic Population Projections (PPP), and are available at http://esa.un.org/unpd/ppp.
This paper focuses on the methods used to convert probabilistic projections of and TFR to probabilistic projections of agespecific mortality and fertility rates. Some limitations of the methods used for the 2014 PPP are identified, and several improvements are proposed to overcome them. The methods presented in this paper are implemented in an open source R package called bayesPop (Ševčíková & Raftery, 2014; Ševčíková et al., 2014).
The paper is organized as follows. In Section 2 we describe the current method in PPP for projecting agespecific mortality rates and our proposed improvements. In Section 2.1 we outline the Probabilistic LeeCarter method used in the 2014 PPP. In the rest of Section 2 we propose several improvements to overcome limitations of this method. These include a new Coherent Kannisto Method for joint projection of future agespecific mortality rates at very high ages that avoids unrealistic crossovers between the sexes (Section 2.2), application of the Coherent LeeCarter method to avoid crossovers at lower ages (Section 2.3), new methods for avoiding jumpoff bias (Section 2.4), and application of the Rotated LeeCarter method to reflect the fact that when mortality rates are low, they tend to decline faster at older than at younger ages (Section 2.5). In Section 3, we describe the current method in PPP for projecting agespecific fertility rates and our proposed improvements. We conclude with a discussion in Section 4.
2 AgeSpecific Mortality Rates for Probabilistic Population Projections
2.1 Probabilistic LeeCarter Method
Our methodology is based on the LeeCarter model (Lee & Carter, 1992):
where is the mortality rate for age and time period . The quantity represents the baseline pattern of mortality by age over time, and is the average rate of change in mortality rate by age group for a unit change in the mortality index . The parameter is a timevarying index of the overall level of mortality, and is the residual at age and time . Throughout this paper, denotes the natural logarithm.
For a given matrix of rates the model is estimated by a least squares method. The baseline mortality pattern is estimated as the average of over the past time periods with observed data. Since the model is underdetermined, is identified by setting , where the sum is over all ages or age groups . Also, is identified by setting , where is the number of past time periods for which data are available. The estimates are then
(1)  
(2)  
(3) 
To forecast , one needs to project into the future. To project , the LeeCarter method uses a random walk with drift:
Lee & Miller (2001) proposed replacing the step of projecting by itself by matching future to future projected .
Current calculations are done using a highest age or open interval of 85. For projections one needs to extend mortality rates to higher ages , usually beyond , because mortality rates are expected broadly to decline over time in the future, so there will be larger numbers of people at higher ages. For extending the force of mortality at older age groups, the Kannisto model provides a robust way to fit available mortality rates from age 80 to 100, and to extrapolate mortality rates up to age 130 in a way that is consistent with empirical observations on oldestold mortality (Thatcher et al., 1998).
The Bayesian probabilistic projections of life expectancy (Raftery et al., 2013; Raftery et al., 2014) provide us with a set of future trajectories of female and male , representing a sample from the joint predictive distribution of future female and male for all countries and all future time periods. The 2014 PPP used methods for turning a trajectory of future values into a set of future agespecific mortality based on the ideas of Lee & Miller (2001) and Li & Gerland (2011); see Raftery et al. (2012). They were based on the following algorithm:
Algorithm 1
Let and denote the observed and projected time periods, respectively.

Using the Kannisto method extend to higher age groups so that for all .

For a given in each trajectory and given and , solve for future numerically using life tables. This yields a nonlinear equation which can be solved using the bisection method. More details are given in Section 2.6.

Compute mortality rates by for each trajectory and future time .
Applying these steps to all trajectories of yields a posterior predictive distribution of .
However, this procedure has a number of drawbacks. There is no assurance that the extension of to higher ages yields mortality rates that are coherent between males and females. Similarly, the predicted can lead to unwanted crossovers between female and male mortality rates, since they are obtained independently for each sex. In the following sections, we present solutions to these and other limitations of the simple algorithm above, and give more details about Step 3.
2.2 Coherent Kannisto Method
A sexindependent extension of the observed mortality rates to higher age categories can lead to unrealistic crossovers at higher ages. We propose a modification of the Kannisto method that treats male and female mortality rates jointly. In this section, for simplicity we omit the time index .
The original Kannisto model has the form
where is a random perturbation with mean zero. The model is usually estimated independently for each sex, assuming independence across ages and normality of the , using a maximum likelihood method (Thatcher et al., 1998; Wilmoth et al., 2007). This yields sexspecific parameter estimates , , , .
We suggest modifying this by forcing the sexspecific parameters and to be equal (i.e. ), but still allowing the parameters and to differ between the sexes:
This leads to the following model:
where if and 0 otherwise.
To estimate the parameters, we fit the model to the observed for ages 80–99 by ordinary leastsquares regression, which corresponds to maximum likelihood under the assumptions of independence and normality of the . There are four age groups in the data used for fitting the model, and thus eight points in total for both sexes. Then,
Figure 1 shows the resulting for old ages for Brazil and Lithuania in the last observed time period. From the left panels we see that there are crossovers using the classic Kannisto method, which is unrealistic. However, male mortality stays above female mortality in the coherent version, as can be seen in the right panels; this is more realistic.
2.3 Coherent LeeCarter Method
We will adopt an extension of the LeeCarter method suggested by Li & Lee (2005), the socalled coherent LeeCarter method. It takes into account the fact that mortality patterns for closely related populations are expected to be similar. In our application, these related populations will be males and females in the same country, since there is no expectation that the life expectancy will diverge between such groups. Thus, the LeeCarter method is extended by two requirements:
(4) 
where and denotes male and female sex, respectively. This ensures that the rates of change of the future mortality rates are the same for the two sexes, and thus avoids crossovers.
2.4 Avoiding Jumpoff Bias
Mortality rates in the last period of the historical data used for estimation (or jumpoff period) are commonly referred to as jumpoff rates (Booth et al., 2006). Often there is a mismatch between fitted rates for the last period and the actual rates (jumpoff bias). As a result, a discontinuity between the actual rates in the jumpoff period and the rates projected in the first projection period may occur.
A possible solution to avoid jumpoff bias is to constrain the model in such a way that passes through zero in the jumpoff period , and to use only from the last fitting period to obtain (Lee & Miller, 2001):
(5) 
A disadvantage of this solution is that in cases where the mortality rates are bumpy in the jumpoff period (i.e. not smooth across ages), this “bumpiness” propagates into the future. In general for projections, we suggest using the agespecific mortality rates from the last fitting period and smoothing them over age if necessary (e.g. for small populations with few deaths in some age groups) while preserving the value for the youngest age group:
(6) 
Figure 2 shows the resulting difference in projected to for two countries using the three different methods of computing , namely equations (1), (5) and (6). As can be seen in the case of Bangladesh, the smoothing step removes bumps whereas the averaging method does not.
Figure 3 shows the impact of the methods on as time series for Bangladesh for three different age groups. Using the average results in jumpoffs for the 5–9 and 95–99 age groups (blue curve). If the latest raw are used, the jumpoffs are eliminated (grey curve). A smoothed version creates a new jumpoff for the age group 75–79 (red curve).
This shows that there is a tradeoff between bumpy mortality rates over ages in later projection years and no jumpoffs, and smooth mortality rates with no jumpoffs. Our solution is to decide on a countryspecific basis which method is more appropriate.
2.5 Rotated LeeCarter Method
Li et al. (2013) focused on the fact that in more developed regions, once countries have already reached a high level of life expectancy at birth, the mortality decline decelerates at younger ages and accelerates at old ages. This change in the pace of mortality decline by age cannot be captured by the original LeeCarter method, since this constrains the rate of change to be constant over time. They proposed instead rotating the over time to a socalled ultimate , denoted by , which is computed as follows.
Let
Then
(7) 
with scaled to sum to unity over all ages.
The rotation is dependent on , and so the resulting also becomes timedependent. The rotation finishes at a certain level of live expectancy, denoted by . Li et al. (2013) recommend using . Using the smooth weight function
the rotated at time , denoted by , is derived as:
(8) 
Figure 4 shows the results for Japan as an example. The original is shown by the black curve. The ultimate , to be reached at life expectancy of 102, is in red. The remaining curves show the change over time starting with yellow and continuing towards the red curve.
2.6 Computing Life Tables
Step 3 in Algorithm 1 calls for matching future to projected . This is a nonlinear equation in . It is solved by an iterative nonlinear procedure in which for given values of , and a life table is produced, and the resulting life expectancy is computed and compared with the projected . We used a bisection method to solve the nonlinear equation. This is simple and robust and involves relatively few iterations. It would be possible to use a nonlinear solution method that is more efficient computationally, but the computational gains would be modest and this could make the method much more complex.
In the process of computing life tables, the conversion of mortality rates to probabilities of dying follows the approach used by the United Nations to compute abridged life tables. This is computed by the LIFTB function in Mortpak (United Nations, 1988, 2013a), where at a given time point the probability of dying for an individual between age and is:
(9) 
with being the length of the age interval and being the average number of years lived between ages and by those dying in the interval. With being the number of survivors at age , we have
(10)  
(11)  
(12) 
where denotes the number of deaths between ages and and denotes the number of personyears lived between ages and . The expectation of life at age (in years) is given by
where is the number of personyears lived at age and older.
For ages 15 and over, the expression for is derived from the Greville (1943) approach to calculating agespecific separation factors based on the age pattern of the mortality rates themselves with:
For ages 5 and 10, and for ages under 5, values from the Coale and Demeny West region relationships are used for (Coale & Demeny, 1966).^{2}^{2}2The Coale and Demeny West region formulae are used as follows. When , then for males and for females; for males and for females. When , ) for males and for females; for males and for females.
2.7 Summary of Improved Algorithm
We now summarize the modifications described in the previous sections by proposing an improved algorithm for deriving the agespecific mortality rates for potential use in future probabilistic population projections.
Algorithm 2
As before, let and denote the observed and projected time periods, respectively. Also, let be an index to distinguish sexspecific measures.

Using the Coherent Kannisto Method from Section 2.2, extend to higher age categories with for all .

Choose a method to estimate , i.e. one of equations (1), (5) or (6), depending on country specifics.^{3}^{3}3In the bayesPop package this countryspecific set of options is controlled through two dummy variables in the vwBaseYear2012 dataset: (1) whether the most recent estimate of age mortality pattern should be used (LatestAgeMortalityPattern) and (2) whether it should be smoothed (SmoothLatestAgeMortalityPattern). See vwBaseYear2012 in R. Do the estimation for each sex , obtaining .

Given and from Step 3, set .

Compute the ultimate as in equation (7).

For a combined in each trajectory, compute as in equation (8).

For a given sexspecific in each trajectory and given and , solve for future numerically using life tables. This yields a nonlinear equation which is solved using the bisection method, as described in Section 2.6.

For each trajectory, time and sex , compute mortality rates by .

Since the previous step does not comply with equation (4) and thus can lead to crossovers in high ages, an additional constraint is added:
If then
Figure 5 shows the resulting probabilistic projection of for the period 20952100 for both sexes in two selected countries. In addition to the marginal distribution for Kazakhstan in the right panel of Figure 5, its joint distribution for males and females is shown in Figure 6 on a logarithmic scale. Points below the solid line indicate crossovers in the individual trajectories. It can be seen that only a few trajectories experience crossovers when mortality is low, i.e. in young ages, suggesting a low (but nonzero) probability for such an event, while there are no crossovers for high mortality, i.e. in old ages. We observed similar results for most countries. Figure 7 shows the same joint distribution for selected age groups on a normal scale.
Exceptions
For about 50 countries, insufficient detailed data about mortality by age and sex are available between 1950 and 2010 (United Nations, 2013c). Therefore, the age patterns of mortality are based on model life tables (e.g., CoaleDemeny). For these countries a model associated with one of the regional model life tables is used (see Table 2 page 18 in Li & Gerland (2011)).^{4}^{4}4In the bayesPop package this countryspecific set of options is controlled through two variables in the vwBaseYear2012 dataset: (1) the type of age mortality pattern used for the estimation period (AgeMortalityType with the option ”Model life tables”) and (2) the specific mortality pattern used (AgeMortalityPattern with options like ”CD West”).
In addition, for about 40 countries with a generalized HIV/AIDS epidemic, age patterns of mortality since the 1980s have been affected by the impact of AIDS mortality (especially before the scaling up of antiretroviral treatment starting in 2005). For these countries the application of the conventional LeeCarter approach is inappropriate.^{5}^{5}5In the BayesPop package this specificset of countries are identified through a dummy variable (WPPAIDS) in the vwBaseYear2012 dataset. Instead, we introduce a modification where steps 2–6 in Algorithm 2 are replaced by the following steps:

Start with the most recent (affected by impact of HIV/AIDS on mortality) and smooth it as in equation (6), obtaining .

Compute an ultimate (or “AIDSfree” target) , denoted by , which is a smoothed average of historical up to 1985 (i.e., prior to the start of the impact caused by HIV/AIDS on mortality), denoted by :
(13) 
For each interpolate from to assuming that in the long run the excess mortality due to the HIV/AIDS epidemic disappears (or reaches a very low endemic level with negligible mortality impact) both as a result of decreased HIV prevalence, improved access to treatment and survival with treatment.

During the projections, pick an by moving along the interpolated line of the corresponding , so that is reached by 2100.

As above, is associated with one of the regional model life tables.
An example of the resulting projected median agespecific mortality rates for Botswana, a country with a generalized HIV/AIDS epidemic, is shown in Figure 8.
There has been recent progress in the modelling of age patterns of mortality for countries with generalized HIV/AIDS (Sharrow et al., 2014). This could provide additional options to better incorporate the uncertainty about future HIV prevalence, expanded access to treatment, underlying age mortality patterns, and their interaction on overall mortality by age into probabilistic population projections. Further calibration and validation of these models using empirical estimates from cohort studies (Zaba et al., 2007; Reniers et al., 2014) will be important in this context.
3 AgeSpecific Fertility Rates for Probabilistic Population Projections
3.1 WPP 2012 Method of Projecting AgeSpecific Fertility Rates
The United Nations probabilistic population projections released in 2014 (Gerland et al., 2014) used a set of projected agespecific fertility rates for each country obtained by combining probabilistic projections of the total fertility rate with deterministic projections of age patterns of fertility as used in the 2012 revision of the World Population Prospects (United Nations, 2014).
For highfertility and mediumfertility countries, future age patterns of fertility were obtained by interpolating linearly between a starting proportionate age pattern of fertility and a target model pattern. The target model pattern was chosen from among 15 proportionate age patterns of fertility, with mean age at childbearing varying between 24 and 28.5 years. The target pattern was held constant once the country reaches its lowest fertility level, or by 20452050 onward.
For low fertility countries, a similar approach was used. This projected future agespecific fertility patterns by assuming that they would reach a target model pattern by 20252030. This target was chosen from among five target age patterns of fertility either for the market economies of Europe (with mean age of childbearing varying between 28 and 32 years) or for countries with economies in transition (with mean age of childbearing varying between 26 and 30 years). Once the model pattern was reached, it was assumed to remain constant until the end of the projection period. In some instances, a modified LeeCarter approach (Li & Gerland, 2009) was used to extrapolate the most recent set of proportionate agespecific fertility rates using the rates of change from countryspecific historical trends.
All the trajectories making up the probabilistic projection of fertility for a given country used the same age pattern of fertility. The choice of target pattern of fertility for a given country, from among the set of model patterns considered, was driven by countryspecific expert opinion about future trends and normative assumptions. No global or regional convergence in age patterns of fertility was imposed.
Figure 9 shows the results of the projections for the Mean Age of Childbearing (MAC) for countries in Eastern Asia from the 2012 Revision of the World Population Prospects.
Overall, the method for projecting agespecific patterns of fertility in the 2012 Revision (as well as in previous revisions) has several limitations. First, no global or regional convergence has been imposed despite the overall convergence in total fertility rates observed in the projection period up to 2100. Second, the time point when the target agespecific pattern is reached is not related to the projected total fertility rates. Third, expert assumptions on the target age pattern and method used for individual countries introduce diversity in the agespecific trends are difficult to explain (see Figure 9 — Mongolia and Democratic People’s Republic of Korea were done by Analyst 1, all other countries by Analyst 2). Finally, since the analysts have used at least two different methods and 25 target age patterns of fertility, the documentation of the decisions made for individual countries have been challenging.Â
3.2 Convergence Method for Projecting AgeSpecific Fertility Rates
We now propose a new method for projecting agespecific fertility rates, to overcome some of the limitations of the existing method used in WPP 2012. This builds on the approach adopted in sets of projections of married or inunion women of reproductive age (MWRA) (United Nations, 2013b). Beginning from the most recent observation of the age pattern of fertility in the base period of projection, the projected age patterns of fertility are based on the past national trend combined with the trend towards the global model age pattern of fertility. The projection method is implemented on the proportionate agespecific fertility rates (PASFR) covering seven age groups from 1519 to 4549. The final projection of PASFRs for each age group is a weighted average of two preliminary projections:

the first preliminary projection, assuming that the PASFRs converge to the global model pattern, see Section 3.2.1; and

the second preliminary projection, assuming the observed national trend in PASFRs continues into the indefinite future, see Section 3.2.2.
The method is applied to all the trajectories that make up of the probabilistic projection of total fertility rate for all countries, based on the historical data in the WPP 2012 Revision (Gerland et al., 2014).
We now define the preliminary projections that constitute our overall projection. We use different notation than in Section 2, so the same symbol may be used to denote different quantities in the two sections.
3.2.1 Trend towards the global model pattern
Let denote the base period of a projection and the year when the global model pattern is reached. For , the proportion of the interval that has elapsed at time is
Section 3.2.4 below gives details about how to estimate .
Let denote PASFR at the base period , and let denote PASFR of the global model pattern.^{6}^{6}6In the bayesPop package the global model pattern is created as an average of most recent PASFRs for a set of countries (selected through a dummy variable in the vwBaseYear2012 dataset). For the purpose of the current analysis, the low fertility countries selected have already reached their Phase III and represent later childbearing patterns with mean age at childbearing close to or above 30 years in 20102015: Austria, the Czech Republic, Denmark, France, Germany, Japan, the Netherlands, Norway and the Republic of Korea. The specification of the countries used for the global model pattern can be changed in input file. The projections at time of PASFR towards the global model pattern, denoted by , is obtained by:
(14) 
Then is renormalized so that it sums to unity for all time periods .
3.2.2 Continuing of observed national trend
Let denote the number of 5year periods over which the model is fitted. Then is the starting time period of the estimation and is PASFR at . is the projected PASFR at time , assuming the past trend was to continue into the future under the following rule:
(15) 
As above, should be scaled to sum to unity for all . Note that in our implementation we use .
3.2.3 Resulting projection
Projected PASFR at time , , is calculated as:
(16) 
Resulting is renormalized to sums to unity for all time periods .
3.2.4 Estimating the time period of reaching global pattern
We assume that the transition from the most recent age pattern of fertility to the global model age pattern of fertility is dependent on the timing when the total fertility rate (TFR) enters Phase III, i.e. when the fertility transition is completed and the country reaches low fertility. For the countries in Phase III, a time series model to project TFR was used that assumed that in the long run fertility would approach and fluctuate around countryspecific ultimate fertility levels based on a Bayesian hierarchical model (Raftery et al., 2014). The time series model uses the empirical evidence from lowfertility countries that have experienced fertility increases from a subreplacement level after a completed fertility transition. At the same time, based on the empirical evidence on the postponement of childbearing in lowfertility countries, profound shifts to later start of childbearing and an increase in the mean age of childbearing are still taking place several periods after the start of Phase III (see Figure 10). The timing and speed of the postponement of childbearing in Phase III is countryspecific and in this paper we implement the assumption that the transition to later childbearing pattern is completed when total fertility approaches countryspecific ultimate fertility levels.
To be more specific, we assume that the time of a completion of the transition to a global model pattern corresponds to the time point , when TFR reaches the ultimate fertility level of that country. In probabilistic projections of TFR, we approximate the ultimate fertility level, denoted by , by the median TFR in the last projection period , e.g. 20952100, if TFR is in Phase III:
(17) 
Then for each TFR trajectory, is the earliest time period, at which the TFR is larger or equal to :
(18) 
where denotes the start of Phase III. For the estimation of , we will now distinguish two cases, depending if is smaller or larger than the end period .
Case 1:
In this case, is either observed () or projected within the projecting period (). In both cases, if exists,
(19) 
This includes a situation where for . In such a case, the global pattern is reached quickly, namely in two 5year periods.
If for all , then does not exist. In such a case, is set to the end of the projection period, but at least five 5year periods after :
(20) 
Case 2:
In this case, is unknown, i.e. the TFR trajectory has not reached Phase III at . Thus, we will make an estimate of , denoted by , and then simply apply
(21) 
If the TFR at is low, namely , we assume that . Otherwise, we approximate by a linear extrapolation of TFR from the last four time periods and determine when such line reaches 1.8, with an upper limit of .
3.2.5 Exception for late childbearing pattern
Since trajectories for some countries have already observed or – as projected by the algorithm described above – will in near future reach higher MAC than the MAC associated with the global model pattern, we assume that for a given country’s trajectory once the maximum MAC is reached in the convergence period the associated PASFR pattern is kept constant for the remaining projection periods. This assumption enables to keep trajectoryspecific patterns of late childbearing for trajectories after the Phase III, thus already with low total fertility (see Figure 11 for example of the Czech Republic). Note that this rule is applied only in Case 1 above.
3.3 Results of the convergence method applied to probabilistic projections
For the 2012 Revision, agespecific fertility estimates are based on empirical data for all countries of the world for the period up to 2010 (or up to 20102015 for 37 countries with empirical data up to 2011 or 2012; Gerland et al., 2014). Using the probabilistic projections of TFR, each TFR trajectory has a specific start of Phase III and therefore the timing of convergence to the global model pattern is trajectoryspecific. This yields a set of trajectories of PASFR (although not probabilistic) which in turn, when combined with the probabilistic TFR, yield probabilistic projection of agespecific fertility rates.
Figure 11 shows an example of the results for PASFR in Niger, Bangladesh and Czech Republic for selected age groups over time. Figure 12 shows an example of the probabilistic results of agespecific fertility rates for Ethiopia, Nepal and Japan at the end of projection period in 20952100.
Figure 13 shows the development of PASFR for Uganda, India and Germany over time from 20052010 to 20952100. Here, the methodology was applied to the deterministic projection of TFR from WPP 2012.
In Figure 9 we showed projections of MAC from WPP 2012. This can be compared to Figure 14 where the same measure is shown after applying the new methodology to the TFR of WPP 2012.
Overall, the new method we propose improves on the current methodology in several ways. First, in the very long term (after 2100) the age patterns of fertility are converging to one global pattern, while retaining specific late childbearing patterns for several countries that reach such patterns in the current period or in the near future. Second, the projections of the age pattern of fertility are now linked to projections of the total fertility rate. Finally, for each probabilistic trajectory, the time when the target age pattern is reached depends on the trajectoryspecific total fertility rate.
4 Discussion
We have described the methods used for converting projected life expectancies at birth and total fertility rates to agespecific mortality and fertility rates in the UN’s 2014 probabilistic population projections. We have identified some limitations of these methods and have proposed several improvements to overcome them. These include a new coherent Kannisto method to avoid crossovers in mortality rates between the sexes at very high ages. They also include the application of a coherent LeeCarter method, methods for avoiding jumpoff bias, and a rotated LeeCarter method to reflect the fact that at high life expectancies, mortality rates tend to decline faster at higher than at lower ages.
It should be noted that the 2014 PPP takes account of uncertainty about the overall level of fertility as measured by the TFR, and also about the overall level of mortality as measured by . Conditional on TFR and , however, the projected vital rates are deterministic. There is thus a missing component of uncertainty, and it would be desirable to extend the methods used to take account of this, particularly of uncertainty about the future mean age at childbearing (Ediev, 2013).
5 Acknowledgements
This research was supported by NIH grants R01 HD054511 and R01 HD070936. The views expressed in this article are those of the authors and do not necessarily reflect those of NIH or the United Nations.
References
 Alkema et al. (2011) Alkema, L., Raftery, A. E., Gerland, P., Clark, S. J., Pelletier, F., Buettner, T., & Heilig, G. K. (2011). Probabilistic projections of the total fertility rate for all countries. Demography, 48, 815–839.
 Booth et al. (2006) Booth, H., Hyndman, R. J., Tickle, L., & de Jong, P. (2006). LeeCarter mortality forecasting: A multicountry comparison of variants and extensions. Demographic Research, 15.
 Coale & Demeny (1966) Coale, A. J. & Demeny, P. G. (1966). Regional Model Life Tables and Stable Populations. Princeton, N.J.,: Princeton University Press.
 Ediev (2013) Ediev, D. M. (2013). Comparative importance of the fertility model, the total fertility, the mean age and the standard deviation of age at childbearing in population projections. Presented at the Meeting of the International Union for the Scientific Study of Population, Busan, Korea. http://iussp.org/sites/default/files/event_call_for_papers/TF%20MS%20SD_what%20matters_StWr.pdf.
 Gerland et al. (2014) Gerland, P., Raftery, A. E., Ševčíková, H., Li, N., Gu, D., Spoorenberg, T., Alkema, L., Fosdick, B. K., Chunn, J. L., Lalic, N., Bay, G., Buettner, T., Heilig, G. K., & Wilmoth, J. (2014). World population stabilization unlikely this century. Science, 346, 234–237.
 Greville (1943) Greville, T. N. (1943). Short methods of constructing abridged life tables. The Record of the American Institute of Actuaries, XXXII, part 1, 29–42.
 Lee & Carter (1992) Lee, R. D. & Carter, L. (1992). Modeling and forecasting the time series of US mortality. Journal of the American Statistical Association, 87, 659–671.
 Lee & Miller (2001) Lee, R. D. & Miller, T. (2001). Evaluating the performance of the LeeCarter method for forecasting mortality. Demography, 38, 537–549.
 Li & Gerland (2009) Li, N. & Gerland, P. (2009). Modelling and projecting the postponement of childbearing in lowfertility countries. Presented at theXXVI IUSSP International Population Conference. iussp2009.princeton.edu/papers/90315.
 Li & Gerland (2011) Li, N. & Gerland, P. (2011). Modifying the LeeCarter method to project mortality changes up to 2100. Presented at the Annual Meeting of Population Association of America. http://paa2011.princeton.edu/abstracts/110555.
 Li & Lee (2005) Li, N. & Lee, R. D. (2005). Coherent mortality forecasts for a group of populations: An extension of the LeeCarter method. Demography, 42, 575–594.
 Li et al. (2013) Li, N., Lee, R. D., & Gerland, P. (2013). Extending the LeeCarter method to model the rotation of age patterns of mortality decline for longterm projections. Demography, 50, 2037–2051.
 Raftery et al. (2014) Raftery, A. E., Alkema, L., & Gerland, P. (2014). Bayesian population projections for the United Nations. Statistical Science, 29, 58–68.
 Raftery et al. (2013) Raftery, A. E., Chunn, J. L., Gerland, P., & Ševčíková, H. (2013). Bayesian probabilistic projections of life expectancy for all countries. Demography, 50, 777–801.
 Raftery et al. (2014) Raftery, A. E., Lalic, N., & Gerland, P. (2014). Joint probabilistic projection of female and male life expectancy. Demographic Research, 30, 795–822.
 Raftery et al. (2012) Raftery, A. E., Li, N., Ševčíková, H., Gerland, P., & Heilig, G. K. (2012). Bayesian probabilistic population projections for all countries. Proceedings of the National Academy of Sciences, 109, 13915–13921.
 Reniers et al. (2014) Reniers, G., Slaymaker, E., NakiyingiMiiro, J., Nyamukapa, C., Crampin, A. C., Herbst, K., Urassa, M., Otieno, F., Gregson, S., Sewe, M., Michael, D., Lutalo, T., Hosegood, V., Kasamba, I., Price, A., Nabukalu, D., Mclean, E., Zaba, B., & Network”, A. (2014). Mortality trends in the era of antiretroviral therapy: evidence from the Network for Analysing Longitudinal Population based HIV/AIDS data on Africa (ALPHA). AIDS, 28, S533–S542.
 Ševčíková & Raftery (2014) Ševčíková, H. & Raftery, A. (2014). bayesPop: Probabilistic Population Projection. R package version 5.22.
 Ševčíková et al. (2014) Ševčíková, H., Raftery, A., & Gerland, P. (2014). Bayesian probabilistic population projections: Do it yourself. Presented at the annual meeting of Population Association of America. http://paa2014.princeton.edu/abstracts/141301.
 Sharrow et al. (2014) Sharrow, D., Clark, S., & Raftery, A. (2014). Modeling agespecific mortality for countries with generalized HIV epidemics. PLoS One, 9.
 Thatcher et al. (1998) Thatcher, A. R., Kannisto, V., & Vaupel, J. W. (1998). The Force of Mortality at Ages 80 to 120, volume 5 of Odense Monographs on Population Aging Series. Odense, Denmark: Odense University Press.
 United Nations (1988) United Nations (1988). MortPak – The United Nations Software Package for Mortality Measurement. Bachoriented Software for the Mainframe Computer, volume ST/ESA/SER.R/78. New York: United Nations.
 United Nations (2013a) United Nations (2013a). MortPak for Windows Version 4.3 – The United Nations Software Package for Demographic Mortality Measurement.
 United Nations (2013b) United Nations (2013b). National, Regional and Global Estimates and Projections of the Number of Women Aged 15 to 49 Who Are Married or in a Union, 19702030, Technical Paper 2013/2. New York, NY: Population Division, Dept. of Economic and Social Affairs, United Nations.
 United Nations (2013c) United Nations (2013c). World Population Prospects: The 2012 Revision  Online and DVD Edition  Data Sources and Meta Information (POP/DB/WPP/Rev.2012/F02). New York: Population Division, Dept. of Economic and Social Affairs.
 United Nations (2014) United Nations (2014). World Population Prospects: The 2012 Revision, Methodology of the United Nations Population Estimates and Projections. ESA/P/WP.235. New York, NY: Population Division, Dept. of Economic and Social Affairs, United Nations.
 Wilmoth et al. (2007) Wilmoth, J., Andreev, K., Jdanov, D., & Glei, D. (2007). Methods protocol for the Human Mortality Database. Online publication of the Human Mortality Database. http://www.mortality.org/Public/Docs/MethodsProtocol.pdf.
 Zaba et al. (2007) Zaba, B., Marston, M., Crampin, A. C., Isingo, R., Biraro, S., Barnighausen, T., Lopman, B., Lutalo, T., Glynn, J. R., & Todd, J. (2007). Agespecific mortality patterns in HIVinfected individuals: A comparative analysis of African community study data. Aids, 21 Suppl 6, S87–S96.