global scen : env SCEN
log using logs/${scen}_cohort_e65.log, replace

use ../../output/$scen/${scen}_by_rep.dta

keep agedied_*_65p_rbyr* n_died_*_65p_rbyr* year rep

reshape long agedied n_died, i(year rep) j(lbl) string
gen male= substr(lbl,1,3)=="_m_"
gen rbyr = substr(lbl,12,4)
destring(rbyr), replace

* sum years lived for those who died after 65 in each year and sum total number of deaths in each year
collapse (sum) agedied n_died, by(male rbyr)

* number of years lived after age 65; assuming death 1 year into the last time step, so subtract 1 year from the final age 
gen e65 = agedied / n_died - 1.0 - 65.0

save data/${scen}_e65.dta