############################################################################# # # # LIFE LIVED AND LEFT: # # ESTIMATING AGE-SPECIFIC SURVIVAL IN STABLE POPULATIONS WITH UNKNOWN AGES # # # # James W. VAUPEL jvaupel@sdu.dk # # Francisco VILLAVICENCIO fvillavicencio@health.sdu.dk # # # # October 2018 # ############################################################################# INTRODUCTION ------------ This file offers a brief overview of the data and R code (R Core Team 2018) available in the supplementary materials of the article Vaupel, J. W. and F. Villavicencio (2018). Life lived and left: Estimating age-specific survival in stable populations with unknown ages. Demographic Research 39(37): 991-1008. Files available: 1) JWV-FV-2018-LifeLivedLeft.R: R-script that can be used to reproduce the empirical results discussed in the article. We have commented the R-script as much as possible in order to help any external user to understand the code. Feel free to contact us if something is not clear. 2) HMD-SwePop5x5.txt: Data on population counts from Sweden 1751-2017, 5-year age groups (Source: Human Mortality Database 2018). 3) HMD-SweFemLtCoh5x5.txt: Data on cohort life tables, Swedish females, 1751-1924, 5-year age groups and time intervals (Source: Human Mortality Database 2018). R-SCRIPT -------- The file JWV-FV-2018-LifeLivedLeft.R can be used to reproduce in R (R Core Team 2018) the results presented in Sections 4.2 and 5 of the main manuscript. SECTION 4.2: Data simulation to test the validity of theoretical results from Equations (9*), (10*) and (11*) on the main manuscript. This section of the code can be used to reproduce the results in Tables 2 and 3 (except from small variations due to the randomness on the data simulation). SECTION 5: In this section we show an example of application with actual data that aims to recover the 1805 cohort life table of Swedish females by applying Equations (9*) and (10*) to data on population counts from the Human Mortality Database (2018). The R code automatically uploads into R the data files HMD-SwePop5x5.txt and HMD-SweFemLtCoh5x5.txt. Next, it estimates the 1805 cohort survival schedule of Swedish females using data on life left. The resulting survival curve is smoothed by combining the R functions predict() and loess(). The code also can be used to reproduce FIGURE 4 on the main manuscript. NOTE: To execute the file JWV-FV-2018-LifeLivedLeft.R it is necessary to have a computer with the open-source statistical software R installed, which can be downloaded from https://www.R-project.org/. Moreover, it is recommended to store all the files (data and code) in the same working directory. REFERENCES ---------- Human Mortality Database (2018). University of California, Berkeley, and Max Planck Institute for Demographic Research, Rostock. URL http://www.mortality.org (data downloaded on 4 June 2018). R Core Team (2018). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/.