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

* Obese N Initial Intervention
*	Only selects persons that are obese but not overweight, and live in a household with or 1 or 2 obese persons. Sets obese to 0 for these persons.


        	bys hhid: egen inc = total(overwt == 0 & obese == 1)
        	keep if inc == 1 | inc == 2
        	replace weight = 0 if !(overwt == 0 & obese == 1)
        	replace obese = 0
        	replace fobese = 0
        	replace lobese = 0
        	
end