This vignette explains how the value of lives lost due to infection in a pandemic scenario is calculated in daedalus.
The main focus is on explaining how the value of a statistical life (VSL, also referred to hereafter as ‘life-value’) is calculated in daedalus. The VSL is multiplied by the number of deaths to get the total value of life-years lost due to the pandemic.
The methodology used here has been chosen so as to be generalisable for the widest possible range of countries and territories, while using publicly available data.
This approach balances the ability to run daedalus scenarios for a large number of countries, and the effort required to obtain and update official life-values from national governments.
Note that the VSL values used in a call to
daedalus()
for any country can always be substituted with
values considered to be more appropriate (e.g., from a national health
ministry).
Value of a life year
The value of a life year (VLY) in a country is assumed to be the gross national income (GNI) per capita scaled for purchasing power parity in international dollars.
Country GNI values are for the most recent year available, and are taken from the World Bank’s World Development Indicators.
Note that there is no discounting applied to older age groups, that is, all remaining years of life are valued the same, and are not adjusted for disabilities or quality.
Life expectancy
Life-expectancy data in terms of years for each country and territory is initially accessed for age-bins and separated by sex.
The life-expectancy data is then aggregated by the DAEDALUS age bins: 0 – 4, 5 – 19 (school age), 20 – 65 (working age), 65+ (retired or not working).
Note that the aggregation of is a simple mean across all age-bins in the raw data; a more appropriate approach which we aim to add in future would be a weighted mean that takes into account the size of each age-bin in the raw data.
This results in four life-expectancy values for each country, each corresponding to one of the four DAEDALUS age groups.
Value of a statistical life
Many policies aim to improve longevity, decreasing the risk of death in each year. The value of these risk reductions is often expressed as the value per statistical life (VSL); at times a value per statistical life year (VSLY) may be used.
The VSL concept is widely misunderstood. It is not the value that the analyst, the government, or the individual places on saving an identified life with certainty. Instead, it reflects individuals’ willingness to exchange money for a small change in their own risk, such as a 1 in 10,000 decrease in the chance of dying in a specific year (Robinson et al. 2019).
The value of a statistical life for an individual of country in age group denoted as is the product of the country-specific VLY () and the country- and age-specific life expectancy (), and depends on the country and age-group of the individual; this can be expressed as:
In daedalus, the country- and age-specifc VSL values can be
accessed from the package data object daedalus::life_value
.
These values are also automatically accessed and included when a
<daedalus_country>
object is created, under the name
vsl
.
# load library for examples
library(daedalus)
# create a <daedalus_country>
x <- daedalus_country("China")
# access VSL values using either a helper function or
# the subsetting operator
get_data(x, "vsl")
#> [1] 1322296.0 1166990.6 674010.6 177462.6
x$vsl
#> [1] 1322296.0 1166990.6 674010.6 177462.6
Value of lives lost
In a daedalus()
scenario, the total value of lives lost
is the sum of the products of the country- and age-specific VSL and the
number of deaths in each age-group in that country:
where is the total number of deaths in age group , with (as detailed above).
For users’ convenience, the total value of lives lost is reported in million dollars, to bring it in line with the scaling used for other pandemic costs.