capt log close

log using logs/population_analysis_out.log, replace


/*
This program pulls in the cohort analysis output by year and puts it out into csv files.
September 25, 2008
*/

global tabdata "/zeno/a/DOL/Tables_ahg"
global tabdata "/zeno/a/DOL/Output_prev2"
global tabdata "/zeno/a/DOL/Output_prev2"

#d ; 
	global sumvar0  dage5164 dage6574 dage7584 dage85p male hispan black hsless college;
	global sumvar1  pop qaly meanage pop_healthy pop_disability 
			pop_iadl1 pop_adl12 pop_adl3 pop_nhm pop_work pop_ssclaim pop_diclaim pop_ssiclaim  pop_lfhibpe pop_lfdiabe pop_lfsmoken pop_lfobese;
	global sumvar2  died cancre diabe hearte hibpe lunge stroke anycond cond3p iadl1 
			adl12 adl3 nhmliv obese overwt smokev smoken $addvars wlth_nonzero; 
	global sumvar3	 ttl_earn ttl_ssben ttl_diben ttl_ssiben;
	global sumvar4   work ry_earn ssclaim ssben diclaim diben ssiclaim ssiben anyhi;
	global sumvar5	 hhwealth ftax stax ctax hoasi hmed gross net;
	global sumvar6   raime totmd caremd oopmd caidmd hsptim hspnit doctim;
	global sumvar7   ttl_totmd ttl_oopmd ttl_caremd ttl_caidmd ttl_hsptim ttl_hspnit ttl_doctim;
	global sumvar8	 ttl_raime ttl_hhwealth ttl_ftax ttl_stax ttl_ctax ttl_hoasi ttl_hmed ttl_gross ttl_net;
	global sumvar9   mhhwealth mraime mry_earn mry_earnwork;
	global sumvar10  tot65 oop65 fracss62 fracss65;
	global sumvar11  treat_now treat_new diabe_new hibpe_new;
	
	global sumvar12 totmd65l caremd65l oopmd65l caidmd65l hsptim65l hspnit65l doctim65l ;
	global sumvar13 totmd65p caremd65p oopmd65p caidmd65p hsptim65p hspnit65p doctim65p ;
	global sumvar14 ttl_totmdl ttl_oopmdl ttl_caremdl ttl_caidmdl ttl_hsptiml ttl_hspnitl ttl_doctiml;
	global sumvar15 ttl_totmdp ttl_oopmdp ttl_caremdp ttl_caidmdp ttl_hsptimp ttl_hspnitp ttl_doctimp;
	
	global sumvars year $sumvar0 $sumvar1 $sumvar2 $sumvar3 $sumvar4 $sumvar5 $sumvar6 
	$sumvar7 $sumvar8 $sumvar9 $sumvar10 $sumvar11 ;

    global sumvars year pop cancre diabe hearte hibpe 
    lunge stroke smokev smoken ttl_totmd ttl_caremd ttl_caidmd caremd;
#d cr

clear

local status_quo_file "/zeno/a/DOL/Output_ahg/51p_sstatus_quo_c.dta"  
local obesity_r_file = "/zeno/a/DOL/Output_ahg/51p_sobese_r_c.dta"  
local nosmoke_r_file = "/zeno/a/DOL/Output_ahg/51p_snosmoke_r_c.dta"
local shareprev_file = "/zeno/a/DOL/Output_ahg/51p_sshareprev_c.dta"
local invobese_r_file = "/zeno/a/DOL/Output_ahg/51p_sinvobese_r_c.dta"

local status_quo_file "/zeno/a/DOL/Output_prev2/51p_sstatus_quo_c.dta" 
local ebmall_file "/zeno/a/DOL/Output_prev2/51p_sebmall_c.dta" 
local Hbp_R_file "/zeno/a/DOL/Output_prev2/51p_sHbp_R_c.dta" 
local Diabetes_R_file "/zeno/a/DOL/Output_prev2/51p_sDiabetes_R_c.dta" 
local Smoking_R_file "/zeno/a/DOL/Output_prev2/51p_sSmoking_R_c.dta" 

local Diab_Hibpe_file "/zeno/a/DOL/Output_prev2/51p_sDiab_Hibpe_c.dta" 

local adj_totmd65l = 1.31
local adj_caremd65l = 1.59
local adj_caidmd65l = 2.53
local adj_totmd65p = 1.14
local adj_caremd65p = 1.11
local adj_caidmd65p = 1.64

* foreach i in status_quo obesity_r nosmoke_r shareprev invobese_r{
* foreach i in status_quo ebmall {
* foreach i in Hbp_R Diabetes_R Smoking_R {
    foreach i in Diab_Hibpe {
	clear
	use "``i'_file'"
	
	foreach v in totmd caremd caidmd {
		replace ttl_`v'l = ttl_`v'l * `adj_`v'65l'
		replace ttl_`v'p = ttl_`v'p * `adj_`v'65p'
		replace ttl_`v' = (ttl_`v'l + ttl_`v'p)/10^9
	}
	
	gen tot65 = .
	gen oop65 = .
		
	outsheet $sumvars using "$tabdata//`i'.csv", c names replace
}

log close

capt log close