cap program drop nsmoker_intervene0 program define nsmoker_intervene0 syntax [varlist] [if] [in], is_startyr(integer) * nsmoker Initial Intervention * Only selects persons that have never smoked and are from households where there are 1 or 2 non smokers bys hhid: egen inc = total(smokev == 0 & smoken == 0) keep if inc == 1 | inc == 2 replace weight = 0 if !(smokev == 0 & smoken == 0) end