cap program drop Smoking_Reduction_intervene0
program define Smoking_Reduction_intervene0
	syntax [varlist] [if] [in], is_startyr(integer)


	marksample touse, zeroweight strok novarlist 
* Smoking Reduction Initial Intervention
* Parameter: init_intervene_smk_rd (scalar)
*	Treats smoking in a fraction of the population given by the init_intervene_smk_rd parameter



	
		replace treat_now = smoken == 1 & died == 0 if `touse'
		replace smoken = 0 if treat_now & normal(x_smkstat2004) <= init_intervene_smk_rd & `touse'
		replace treat_effective = old_smoken == 1 & smoken == 0 if `touse'
		replace lsmoken = smoken
		
	end