*********************************************************
** 7/22/2016 Revised to person-wave
** Create variable indicate people survive for four waves
*********************************************************

set more off

qui include ../fem_env.do

use "$dua_rand_hrs/bio_imp_all.dta", clear

gen alive = (iwstat == 1 | iwstat == 4)

bys hhidpn: egen alive_total = total(alive)

recode age (51/64=1) (65/80=2) (nonmiss=4), gen(agecat3)
replace agecat3 = 3 if age > 80
recode agecat3 (4 = .)

preserve

duplicates drop hhidpn, force

gen indicator = (wave == 8 & alive_total >=3)
tab indicator

gen indicator2 = (wave == 10 & alive_total >=2)
tab indicator2

replace hibperx = 1 if hibperx == 1 & hibpe == 1
replace hibperx = 2 if hibperx == 0 & hibpe == 1

replace diaberx = 1 if diaberx == 1 & diabe == 1
replace diaberx = 2 if diaberx == 0 & diabe == 1

local bioimp hdl a1c sysbp crp

keep if indicator == 1

foreach v of local bioimp {
	gen `v'_mp = 0 if `v'_nowave == 1
	replace `v'_mp = 1 if `v'_onewave == 1
	replace `v'_mp = 2 if `v'_twowaves == 1
	replace `v'_mp = 3 if `v'_threewaves == 1

tab `v'_mp, m


  tab agecat3 `v'_mp , col
  tab agecat3 `v'_mp [aw=wtresp] , col
  tab agecat3 `v'_mp [aw=wtresp] , col
  
  tab male `v'_mp , col 
  tab male `v'_mp [aw=wtresp] , col
  tab male `v'_mp [aw=wtresp] , col
   
  tab race `v'_mp  , col
  tab race `v'_mp [aw=wtresp] , col
  tab race `v'_mp [aw=wtresp] , col
  
  tab educ `v'_mp , col
  tab educ `v'_mp [aw=wtresp], col
  tab educ `v'_mp [aw=wtresp], col
  
  tab hearte `v'_mp , col
  tab hearte `v'_mp [aw=wtresp], col
  tab hearte `v'_mp [aw=wtresp], col
  
  tab stroke `v'_mp , col
  tab stroke `v'_mp [aw=wtresp], col
  tab stroke `v'_mp [aw=wtresp], col
  
  tab rxchol `v'_mp , col
  tab rxchol `v'_mp [aw=wtresp], col
  tab rxchol `v'_mp [aw=wtresp], col
  
  tab diabe `v'_mp , col
  tab diabe `v'_mp [aw=wtresp], col
  tab diabe `v'_mp [aw=wtresp], col
  
  tab hibpe `v'_mp , col
  tab hibpe `v'_mp [aw=wtresp], col
  tab hibpe `v'_mp [aw=wtresp], col
  
  tab diaberx `v'_mp, col
  tab diaberx `v'_mp [aw=wtresp], col
  tab diaberx `v'_mp [aw=wtresp], col
  
  tab hibperx `v'_mp  , col
  tab hibperx `v'_mp [aw=wtresp], col
  tab hibperx `v'_mp [aw=wtresp], col
  
  tab smkstat `v'_mp , col
  tab smkstat `v'_mp [aw=wtresp], col
  tab smkstat `v'_mp [aw=wtresp], col
  
  tab overwt `v'_mp, col 
  tab overwt `v'_mp [aw=wtresp], col
  tab overwt `v'_mp [aw=wtresp], col
   
  tab obese `v'_mp , col
  tab obese `v'_mp [aw=wtresp], col
  tab obese `v'_mp [aw=wtresp], col

}

restore