*****************************************
* Examine how relationship between smoking and mortality change by age
*****************************************

clear
set more off
set mem 500m

global outdata "/zeno/d/zeno_a/DOL/Output_Dec/51p"

global slist status_quo smk80

foreach s in $slist{
	use "$outdata//`s'/51p_s`s'_y2004_rep1.dta
	keep entry hhidpn age fsmokev fsmoken weight died year

	save tmp, replace
	forvalues i = 2006(2)2050 {
		use "$outdata//`s'/51p_s`s'_y`i'_rep1.dta", clear
		keep entry hhidpn age fsmokev fsmoken weight died pdied year
		append using tmp
		save tmp,replace
	}
	save "$outdata/`s'_sub.dta", replace
	erase tmp.dta
}
exit

* Check mortality by smoking status
use "$outdata/status_quo_sub.dta"
	foreach x in age fsmokev fsmoken weight died pdied {
		ren `x' `x'_o
	}
	
foreach x in pop died pdied {
	ren `x' `x'_s
}
sort entry fsmokev year, stable
merge entry fsmokev year using "$outdata/status_quo_collapse.dta", sort