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

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

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