* compute age-specific moratlity rates from FEM output for each demographic group

global scen : env SCEN

log using logs/${scen}_mortality_rates.log, replace

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

foreach gender in m f {
	       forvalues agegrp=55/99 {
	       local pop `gender'_`agegrp'
	       gen exposure_`pop' = 2*startpop_`pop' - n_died_`pop'
	}
}

collapse (mean) exposure* n_died* startpop*, by(year)

keep year n_died_* exposure_* startpop_* 

reshape long n_died exposure startpop, i(year) j(grp) string

drop if strlen(grp)!=5 | strmatch("_",substr(grp,3,1))==0
tab grp, m
gen male = strmatch("_m_",substr(grp,1,3))
gen agegrp = real(substr(grp,4,2))
desc
sort male agegrp year

tab male, m
tab agegrp year, m
sum exposure

save data/${scen}_mortality_rates.dta, replace

exit, clear