set more off

sum dmage dmeduc married black nlbnl smoke dbirwt gestat adeqcode2 adeqcode3 novisit pretri2 pretri3 wtgainnet bigcity bigcnty

sum cigar if smoke & cigar<99

// drop if statenat==5 | statenat==15 | statenat==33 | statenat==42
// drop california, indiana, new york, south dakota observations

regr dbirwt smoke male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3, robust

regr dbirwt smoke cigar male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3 if cigar<99, robust

vce

xi: regr dbirwt smoke male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3 i.nlbnl, robust

xi: regr dbirwt smoke cigar male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3 i.nlbnl if cigar<99, robust

vce

xi: regr dbirwt smoke male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3 i.nlbnl i.stateres, robust

xi: regr dbirwt smoke cigar male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3 i.nlbnl i.stateres if cigar<99, robust

vce

xi: regr dbirwt smoke male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3 wtgainnet wtgainsq i.nlbnl i.stateres, robust

xi: regr dbirwt smoke cigar male dmage agesq hsgrad somecoll collgrad married black year* adeqcode2 adeqcode3 novisit pretri2 pretri3 wtgainnet wtgainsq i.nlbnl i.stateres if cigar<99, robust

vce