***=======================================================*	 
* FEM scenario combination - lifetime
* Mar 2008
* May 2008, graph lifetime trajectory of diabetes and hibpe
* Export to Excel files
***=======================================================*	 

clear
cap clear mata
set more off
set mem 500m

***=============================*	 
* DEFINE DIRECTORIES 
global workdir "/zeno/a/zyuhui/DOL/Model"
global indata "/zeno/a/zyuhui/DOL/Input"
global indata2 "/zeno/a/zyuhui/DOL/Indata"
global netindata "/homer/c/Retire/yuhui/rdata2"
global outdata "/zeno/a/zyuhui/DOL/Output"
global netoutdata "/homer/c/Retire/yuhui/simu"

* For local
global workdir "//zeno/zeno_a/zyuhui/DOL/Model"
global indata "//zeno/zeno_a/zyuhui/DOL/Input"
global indata2 "//zeno/zeno_a/zyuhui/DOL/Indata"
global outdata "//zeno/zeno_a/zyuhui/DOL/Output"
global outdata "//zeno/zeno_a/zyuhui/DOL/Output_Prev2"
global netindata "//homer/homer_c/Retire/yuhui/rdata2"
global netoutdata "//homer/homer_c/Retire/yuhui/simu"

global workdir "//zeno/zeno_a/zyuhui/ahg/Model_yh"
global indata "//zeno/zeno_a/zyuhui/ahg/Input_yh"
global indata2 "//zeno/zeno_a/zyuhui/ahg/Indata_yh"
global outdata "//zeno/zeno_a/zyuhui/ahg/Output_yh"
global outdata "//zeno/zeno_a/zyuhui/ahg/Output_yh"
global netindata "//homer/homer_c/Retire/yuhui/rdata2"
global netoutdata "//homer/homer_c/Retire/yuhui/simu"

global workdir "//zeno/zeno_a/DOL/Othercodes_ahg"
global indata 	"//zeno/zeno_a/DOL/Input_yh"
global indata2 	"//zeno/zeno_a/DOL/Indata_yh"
global outdata 	"//zeno/zeno_a/DOL/Output_prev2"
global outdata 	"//zeno/zeno_a/DOL/Output_prev2"
global netindata "//homer/homer_c/Retire/yuhui/rdata2"
global netoutdata "//homer/homer_c/Retire/yuhui/simu"

**************************************************
* Diabetes versus status-quo
**************************************************
local sname Diabetes_R Diabetes_J Diabetes_E
local dname diabe
local sname2 diabetes
local f 5152
foreach scr in status_quo `sname' {
	use "$outdata//`f'_s`scr'_c.dta"
	keep year pop `dname'
	* replace `dname' = pop * `dname'
	* replace pop = int(pop * 10^6)
	* replace `dname' = int(`dname' * 10^6)
	replace `dname'= `dname' * 100
	foreach x in pop `dname'{
		ren `x' `x'_`scr'
	}
	format pop `dname' %9.1f
	sort year, stable
	save tmp_`scr'.dta, replace
}
drop _all
use tmp_status_quo
erase tmp_status_quo.dta
foreach scr in `sname'{
	merge year using  tmp_`scr'
	drop *merge
	sort year
	erase tmp_`scr'.dta
}

gen age = (year - 2004) + 51
label var age "Age"
outsheet age `dname'_* using "$outdata/lifetime_`dname'.csv", replace comma nol 

**************************************************
* Hypertension versus status-quo
**************************************************
drop _all
local sname Hbp_R Hbp_J Hbp_E
local dname hibpe
local sname2 hypertension

local f 5152
foreach scr in status_quo `sname' {
	use "$outdata//`f'_s`scr'_c.dta"
	keep year pop `dname'
	* replace `dname' = pop * `dname'
	* replace pop = int(pop * 10^6)
	* replace `dname' = int(`dname' * 10^6)
	replace `dname'= `dname' * 100
	foreach x in pop `dname'{
		ren `x' `x'_`scr'
	}
	
	format pop `dname' %9.1f
	
	sort year, stable
	save tmp_`scr'.dta, replace
}
drop _all
use tmp_status_quo
erase tmp_status_quo.dta
foreach scr in `sname'{
	merge year using  tmp_`scr'
	drop *merge
	sort year
	erase tmp_`scr'.dta
}

gen age = (year - 2004) + 51
label var age "Age"
outsheet age `dname'_* using "$outdata/lifetime_`dname'.csv", replace comma nol 

**************************************************
* Smoking now versus status-quo
**************************************************
drop _all
local sname Smoking_R Smoking_J Smoking_E
local dname smoken
local sname2 smoking-now

local f 5152
foreach scr in status_quo `sname' {
	use "$outdata//`f'_s`scr'_c.dta"
	keep year pop `dname'
	* replace `dname' = pop * `dname'
	* replace pop = int(pop * 10^6)
	* replace `dname' = int(`dname' * 10^6)
	replace `dname'= `dname' * 100
	foreach x in pop `dname'{
		ren `x' `x'_`scr'
	}
	
	format pop `dname' %9.1f
	
	sort year, stable
	save tmp_`scr'.dta, replace
}
drop _all
use tmp_status_quo
erase tmp_status_quo.dta
foreach scr in `sname'{
	merge year using  tmp_`scr'
	drop *merge
	sort year
	erase tmp_`scr'.dta
}

gen age = (year - 2004) + 51
label var age "Age"
outsheet age `dname'_* using "$outdata/lifetime_`dname'.csv", replace comma nol 

**************************************************
* Obesity prevalence under diabetes 
**************************************************
drop _all
local sname Obs_R Obs_J Obs_E
local dname obese
local sname2 obese

local f 5152
foreach scr in status_quo `sname' {
	use "$outdata//`f'_s`scr'_c.dta"
	keep year pop `dname'
	* replace `dname' = pop * `dname'
	* replace pop = int(pop * 10^6)
	* replace `dname' = int(`dname' * 10^6)
	replace `dname'= `dname' * 100
	foreach x in pop `dname'{
		ren `x' `x'_`scr'
	}
	
	format pop `dname' %9.1f
	
	sort year, stable
	save tmp_`scr'.dta, replace
}
drop _all
use tmp_status_quo
erase tmp_status_quo.dta
foreach scr in `sname'{
	merge year using  tmp_`scr'
	drop *merge
	sort year
	erase tmp_`scr'.dta
}

gen age = (year - 2004) + 51
label var age "Age"
outsheet age `dname'_* using "$outdata/lifetime_`dname'.csv", replace comma nol 
exit