Vaccinating the oldest against COVID-19 saves both the most lives and most years of life.

Many competing criteria are under consideration for prioritizing COVID-19 vaccination. Two criteria based on age are demographic: lives saved and years of future life saved. Vaccinating the very old against COVID-19 saves the most lives, but, since older age is accompanied by falling life expectancy, it is widely supposed that these two goals are in conflict. We show this to be mistaken. The age patterns of COVID-19 mortality are such that vaccinating the oldest first saves the most lives and, surprisingly, also maximizes years of remaining life expectancy. We demonstrate this relationship empirically in the United States, Germany, and South Korea and with mathematical analysis of life tables. Our age-risk results, under usual conditions, also apply to health risks.

Margaret Keenan received the first vaccination against COVID-19 given in the United Kingdom just a week before her 91st birthday (1). In contexts where vaccine allocation is based on age, prioritizing the very old, who have the highest risk of dying, averts the most deaths. But, from the point of view of saving years of remaining life, it seems that it might be better to target younger people, who are less likely to die quickly of another cause if they can be protected from the coronavirus.

Common intuition suggests a trade-off between saving the maximum number of lives and saving the most future life. Such a trade-off is central to philosophical and ethical discussions (2–4), and epidemiological and public health studies often report both numbers of deaths and numbers of years of life lost (YLL) (5–7) in quests for balanced health intervention priorities.

Here we show that this intuition about demographic trade-offs is wrong in the case of COVID-19. While it is true that remaining life expectancy declines with age, this decline is overwhelmed by the exponentially increasing risk of death. For COVID-19, it turns out that vaccinating first the oldest old saves the most lives and also the most life left. Early guidelines from the Centers for Disease Control and Prevention list those aged “65+” as one priority group among several. But, even in terms of YLL, the oldest old in this group should be first in line.

Naturally, many complex considerations besides age impinge on vaccination strategies, and we defer to epidemiologists, ethicists, and experts in health-related sciences for comprehensive studies that include transmission dynamics and other priorities. Our aim in this report is strictly demographic. We take one—narrow but central—aspect of the picture, and show that formal demography can avert one common misunderstanding.

Concepts

We write for the expected lives saved per vaccination of a person aged . Each death averted at age prevents the loss of years of potential life. The product gives the expected years of life saved per vaccination.

The age pattern of mortality reduction per vaccination depends on the age pattern of vaccine efficacy and on the age pattern of mortality risk. (The mortality risk reflects the combined effects of infection risk and the infection fatality rate [IFR].) Early trial results suggest that there are not yet detectable differences in vaccine efficacy by age (8). COVID-19 mortality risk rises exponentially with age at nearly the same rate as all-cause mortality (9, 10).

These two stylized facts together allow expected person-years saved per vaccination to be expressed aswith vaccine efficacy being a constant over age, and COVID-19 mortality being a constant fraction of all-cause mortality, .

Observations

Fig. 1 shows age profiles of COVID-19 mortality rates, remaining life expectancy, and person-years saved per vaccination for the United States, Germany, and South Korea. COVID-19 mortality rates (Fig. 1) grow at a nearly constant exponential rate in all three countries (9). Within an age-based framework, vaccination of the oldest group will clearly maximize , lives saved.

Fig. 1.

Age patterns of (A) COVID-19 death rates, (B) remaining life expectancy, and (C) person-years saved per effective vaccination in the United States, Germany, and South Korea. In all three countries, within an age-based framework, vaccination of the oldest group will maximize person-years saved per effective vaccination. (Details are provided in ).

Remaining life expectancies (Fig. 1) decline fairly linearly in the same manner in all three countries, flattening out only at very old ages.

Years saved per vaccination (Fig. 1) rise with age, with the oldest age saving the most remaining life. Some decline in the rate of increase occurs at older ages in Germany and South Korea, but the maximum is still reached at the oldest ages.

The rate of exponential increase with age of COVID-19 mortality is similar in all three countries, about 11% per year. The differences between age groups are enormous. Assuming equal effectiveness over age, vaccines given to 90-year-olds would save 3 times more lives than would the same doses given to 80-year-olds—and more than 80 times as many lives were the doses given to 50-year-olds.

In terms of maximizing person-years of remaining life, vaccinating a 90-year-old in the United States would be expected to save twice as many person-years as vaccinating a 75-year-old, and 6 times as many as vaccinating a 50-year-old.

Mathematical Generalization

The pattern of increasing years of life saved is not an accidental feature of the particular mortality schedules presented in Fig. 1. It is a general mathematical property holding for a class of mortality models including all of the main models favored by demographers for describing rising age-specific rates from midadulthood onward. In , we establish inequality (Eq.) which guarantees that years of life saved increases with age. A prime example satisfying the inequality is the Gompertz model which posits exponentially rising mortality with age and gives good fits to the empirical schedules in Fig. 1. Other models satisfying the inequality are the Weibull models with rates increasing like a power of age, and logistic models with rates that begin to diverge from Gompertz patterns at extreme ages and level out into mortality plateaus.

A corollary applies to health status, assuming the popular “proportional hazards” approach. Vaccinating those at age who are in poorer health will save more years of life than vaccinating healthy people of the same age.

Discussion

Allocating scarce COVID-19 vaccine doses involves many complex trade-offs. However, a conflict between minimizing the count of deaths and maximizing remaining life is not one of them.

Contrary views are widespread. The World Health Organization, for example, rejected the use of YLL as an allocation criterion in part because “a priority approach relying on YLL could be viewed as disrespectful to older people by failing to address their disproportionately higher risk of death” (11). Our empirical analysis shows it is easier than thought to set such fears aside and to give vaccine priority to the oldest old and those in the most vulnerable states of health.

The assumptions of our formal investigation can be relaxed somewhat. COVID-19 mortality rates rise slightly faster with age than all-cause mortality rates [about 1% faster (9)], so our results would continue to hold if it is found that vaccination effectiveness declines slightly with age. Likewise, the period mortality rates in Fig. 1 presumably rise a little more steeply than the cohort rates driving actual future years lived. In situations where those in long-term care have already been prioritized, a relevant question for further prioritization is the steepness of the mortality age gradient for the remaining population. This is an area of active research (10).

The age-based mathematical analysis can be extended, using common assumptions, to other risk factors such as health status. Under the proportional hazards Gompertz model, giving priority to people most at risk due to combinations of age and health status will save both the most lives and life left.

Our findings give further support to announced priority rankings that vaccinate first those most at risk for dying from COVID-19. Our analysis is not a replacement for more comprehensive modeling accounting for both the direct and indirect effects of vaccines. Rather, in the case of COVID-19, we show that a general property of life tables can be invoked to wipe away the potential conflict between saving lives and maximizing future years of life. Giving COVID-19 vaccines to those facing the highest risk of death from the disease accomplishes both of these goals.

Materials and Methods

Data and Methods for Fig. 1.

COVID-19 death rates by age are calculated per capita, using cumulative data from late fall 2020 based on data from the Institut National d’Etudes Demographiques (INED) (12). The last week of reporting available was November 24 in Germany, November 28 in the United States, and December 2 in South Korea.

Life expectancy profiles are from the Human Mortality Database (HMD) (13), using recent period life tables (2017 for the United States and Germany, and 2018 for South Korea). Life expectancy was interpolated to the midpoint of all COVID-19 mortality reporting age groups. The average age of death from COVID-19 for the open interval was approximated using all-cause HMD life tables as follows: (92.0) for the United States, (101.9) for Germany, and (89.5) for South Korea.

Mathematical Proofs.

We use continuous notation for the hazard and cumulative hazard , with , , and . As in Eq. , years of life saved are assumed to be proportional to the product given bywith the integral and its denominator defining , life expectancy at age .

Proposition: Suppose that the hazard is twice differentiable and for all ages ,Then is monotone increasing in age.

Proof: Consider the inverse function for the nondecreasing function and let be the composition of the hazard with this inverse function. That is, when , and when . Using to change the variable of integration from to (and absorbing the factor that appears outside the integral), we find that the value of corresponding to is given byFor ,Write for the fraction in brackets inside the integral. Define ages by , , , and , and put . From the definition of in terms of the inverse function of , we can write as .

For small , we haveand similarly for the relationship between and . ThusBy assumption, the last integrand is greater than zero. That assures us that is positive for all and all , making greater than and thus making a monotone increasing function of . Q.E.D.

As stated in Mathematical Generalization, most demographic models in common use for later adult age-specific mortality satisfy the conditions of the proposition. For exponentially increasing Gompertz hazards given by for and , the right-hand side of Eq. vanishes. For logistic hazards with , the right-hand side is negative. For Weibull hazards with with , it is also negative. Thus the proposition holds for all these models. Artificial hazard functions like for some values of fail to satisfy the condition for some , but these have little empirical relevance.

Remark: In the presence of a proportional hazard multiplier , is replaced by a product , and is replaced by a product . Using the inverse of the original, baseline cumulative hazard to implement the change of variables, under the rescaled hazard rates, in the counterpart of Eq. , in the numerator becomes , and in the denominator becomes , but in the denominator remains unchanged, since it derives from the change of variables. Under a further change of variables from to , we haveAs long as and therefore are monotone increasing functions and , is smaller than for all , and therefore is bigger than . At any fixed age, those with a higher risk multiplier have higher values for years of life saved.

For Gompertz hazards , a stronger property holds. A person with proportional hazards multiplier faces hazardsIn other words, such a person experiences the same hazards as would a person years older who has no multiplier, implying that orderings of risk by and by come out to be the same.