# delimit ;

clear ;

set mem 500m;

use "c:\ErikMain\race_signaling\race_merge_into_nber_files.dta", clear;

drop year;

gen year = year_first;

replace year = 1986 if year == 86 ;
replace year = 1987 if year == 87 ;
replace year = 1988 if year == 88 ;
replace year = 1989 if year == 89 ;
replace year = 1990 if year == 90 ;
replace year = 1991 if year == 91 ;
replace year = 1992 if year == 92 ;
replace year = 1993 if year == 93 ;
replace year = 1994 if year == 94 ;
replace year = 1995 if year == 95 ;

gen id_1 = id;

#delimit ;

keep id_1 as_comp1 as_comp2 state popsize smsastat hispanic black white asian count;

sort id_1;

save cex_merge.dta, replace; 


#delimit ;


cd c:\ErikMain\lifecycle_consumption;

use merged_file, clear;

/* add a date variable so as to merge with the CPI files */; 

gen year=floor(file/10);
gen quarter=file-year*10;
gen date=yq(year,quarter);
format date %tq;
sort date;

#delimit ; 

/*create a unique identifier*/

egen id=group(newid date);

/*define head*/;
gen head=relation==1;
gen temp=relation==1 if marital==1 & sex==1;
egen spouse=max(temp) if marital==1 , by(id);
drop temp;
replace head=0 if sex==2 & marital ==1 & spouse==1;
egen temp=max(age) if head==1, by(id);
replace head=0 if age<temp;
drop temp;
egen temp=min(sex) if head==1, by(id);
replace head=0 if sex>temp;
drop temp;
egen tag=tag(head age id sex);
replace head=0 if tag==0 & head==1;


/* merge with the CPI file */;

#delimit ;

sort date ;
merge date using cpi_u, nokeep;
tab _merge;
drop _merge;



/* keep only heads */ ; 

#delimit ;

keep if head==1 & year >= 1986;



/* merge in additional CEX data */ ; 


gen id_1 = newid;
sort id_1;

merge id_1 using cex_merge; 
tab _merge;

keep if _merge == 3;
drop _merge;


/*generate real expenditure (2005 dollars)*/;

#delimit ;

foreach var of varlist 
wages		bus	farm	rents	div	interest	pension	socsec	ssi	unemp	workcomp	welfare 
scholar	foodstmp	gvpremia	rrpremia	sspremia	fedtax	statax	pproptax	othtax	
nontax	foodhome	foodout	foodwork	tobacco	alcohol	niteclub	clothes	tailors	jewelry	
toiletry	hlthbeau	renthome	rentothr	furnish	housuppl	elect	gas	water	homefuel	telephon	
servants	drugs	orthopd	doctors	hospital	nurshome	helthins	busiserv	lifeins	autos	parts	
carservs	gasoline	tolls	autoins	masstran	othtrans	airfare	books	pubs	recsport	othrec	
gambling	highedu	lowedu	othedu	charity	insrefnd	intauto	intoth	rentnpay	homevalu	
homeval2	ohint	ohtax	ohmaint	ohprinc	ohlump	ohsold	ohbuy	ohmort1	ohmort2	ohtrans	
housadd	checkng	dcheckng	saving	dsaving	securi	dsecuri	carloan	tradein	carsold	
carprinc	investb	pnpremia	sepremia	deltaiou	owe1	owe1q1	owe5	owe5q1	give	receive	
give2	lumpsums	othernet {;
gen `var'_real=`var'/cpi_u_4q*194.1;
};



/* convert consumption expenditures from annual to quarterly */;

#delimit ;

foreach var of varlist  foodhome foodout foodwork	tobacco	alcohol	niteclub	clothes	tailors	jewelry	
toiletry	hlthbeau	renthome	rentothr	furnish	housuppl	elect	gas	water	homefuel	telephon	
servants	drugs	orthopd	doctors	hospital	nurshome	helthins	busiserv	lifeins	autos	parts	
carservs	gasoline	tolls	autoins	masstran	othtrans	airfare	books	pubs	recsport	othrec	
gambling	highedu	lowedu	othedu	charity	insrefnd	intauto	intoth	
carprinc	 {;
gen `var'_adj =`var'_real/count;
};



#delimit ;

gen age_4 = .;
replace age_4 = 26 if age >=26 & age <= 29;
replace age_4 = 30 if age >=30 & age <= 33;
replace age_4 = 34 if age >=34 & age <= 37;
replace age_4 = 38 if age >=38 & age <= 41;
replace age_4 = 42 if age >=42 & age <= 45;
replace age_4 = 46 if age >=46 & age <= 49;
replace age_4 = 50 if age >=50 & age <= 53;
replace age_4 = 54 if age >=54 & age <= 57;
replace age_4 = 58 if age >=58 & age <= 61;
replace age_4 = 62 if age >=62 & age <= 65;
replace age_4 = 66 if age >=66 & age <= 69;
replace age_4 = 70 if age >=70 & age <= 73;

#delimit ;

gen hhsize = . ;
replace hhsize = 1 if famsize >= 0.5 & famsize < 1.5; 
replace hhsize = 2 if famsize >= 1.5 & famsize < 2.5; 
replace hhsize = 3 if famsize >= 2.5 & famsize < 3.5; 
replace hhsize = 4 if famsize >= 3.5 & famsize < 4.5; 
replace hhsize = 5 if famsize >= 4.5 & famsize < 5.5; 
replace hhsize = 6 if famsize >= 5.5 & famsize < 6.5; 
replace hhsize = 7 if famsize >= 6.5 & famsize < 7.5; 
replace hhsize = 8 if famsize >= 7.5 & famsize < 8.5; 
replace hhsize = 9 if famsize >= 8.5 & famsize < 9.5; 
replace hhsize = 10 if famsize >= 9.5 & famsize < 10.5; 
replace hhsize = 11 if famsize >= 10.5 & famsize < .; 

#delimit ;

gen famsize_1 = hhsize == 1;
gen famsize_2 = hhsize == 2;
gen famsize_3 = hhsize == 3;
gen famsize_4 = hhsize == 4;
gen famsize_5 = hhsize == 5;
gen famsize_6 = hhsize >= 6;

gen adults = as_comp1 + as_comp2;

gen adults_1 = adults >= 1   & adults < 1.5;
gen adults_2 = adults >= 1.5 & adults < 2.5;
gen adults_3 = adults >= 2.5 & adults < 3.5;
gen adults_4 = adults >= 3.5 & adults < 4.5;
gen adults_5 = adults >= 4.5 & adults < 5.5;
gen adults_6 = adults >= 5.5;

replace adults_1 = . if adults == .;
replace adults_2 = . if adults == .;
replace adults_3 = . if adults == .;
replace adults_4 = . if adults == .;
replace adults_5 = . if adults == .;
replace adults_6 = . if adults == .;



#delimit ;

gen married = marital == 1;
replace married = . if marital == .;

gen weight = totwt;

#delimit ;

gen ed_a = 0;
gen ed_b = 0;
gen ed_c = 0;
gen ed_d = 0;


replace ed_a = 1 if educatio > 0 & educatio <= 11 ;

replace ed_b = 1 if educatio == 12 & year < 1996;
replace ed_b = 1 if (educatio == 38 | educatio == 39) & year >= 1996;

replace ed_c = 1 if (educatio == 21 | educatio == 22 | educatio == 23) &  year < 1996 ;
replace ed_c = 1 if (educatio == 40 | educatio == 41 | educatio == 42) &  year >= 1996;

replace ed_d = 1 if (educatio == 24 | educatio == 31 | educatio == 32) &  year < 1996 ;
replace ed_d = 1 if (educatio == 43 | educatio == 44 | educatio == 45 | educatio == 46) & year >= 1996 ;

replace ed_a = . if educatio == . | educatio == 0 ;
replace ed_b = . if educatio == . | educatio == 0 ;
replace ed_c = . if educatio == . | educatio == 0 ;
replace ed_d = . if educatio == . | educatio == 0 ;


gen age_sq = age^2;

gen male = .;
replace male = 1 if sex == 1;
replace male = 0 if sex == 2;

gen northeast = region >= 1 & region < 1.5;
gen midwest = region >= 1.5 & region < 2.5;
gen south = region >= 2.5 & region < 3.5;
gen west = region >= 3.5 & region <= 4;

replace northeast = . if region == .;
replace midwest = . if region == .;
replace south = . if region == .;
replace west = . if region == .;

#delimit ;

egen income = rsum(wages_real	bus_real farm_real rents_real div_real interest_real pension_real socsec_real ssi_real unemp_real 
workcomp_real welfare_real) ;

gen ln_income = ln(income);
replace ln_income = 0 if income <= 0;
gen income_zero = income <= 0;

#delimit ;

gen income_sq = income^2;
gen income_cube = income^3;


gen year_84 = year == 1984 ;
gen year_85 = year == 1985 ;
gen year_86 = year == 1986 ;
gen year_87 = year == 1987 ;
gen year_88 = year == 1988 ;
gen year_89 = year == 1989 ;
gen year_90 = year == 1990 ;
gen year_91 = year == 1991 ;
gen year_92 = year == 1992 ;
gen year_93 = year == 1993 ;
gen year_94 = year == 1994 ;
gen year_95 = year == 1995 ;
gen year_96 = year == 1996 ;
gen year_97 = year == 1997 ;
gen year_98 = year == 1998 ;
gen year_99 = year == 1999 ;
gen year_00 = year == 2000 ;
gen year_01 = year == 2001 ;
gen year_02 = year == 2002 ;
gen year_03 = year == 2003 ;


gen urban = blsurbn == 1;

#delimit ; 


/*  visible expenditure measure */  ; 


gen clothes_pcare_adj = clothes_adj + toiletry_adj + hlthbeau_adj + jewelry_adj + tailors_adj; 
gen ln_clothes_pcare = ln(clothes_pcare_adj);
replace ln_clothes_pcare = 0 if clothes_pcare_adj == 0;


gen pcare_adj = toiletry_adj + hlthbeau_adj ; 
gen ln_pcare = ln(pcare_adj);
replace ln_pcare = 0 if pcare_adj == 0;


gen clothing_adj = clothes_adj + jewelry_adj + tailors_adj; 
gen ln_clothing = ln(clothing_adj);
replace ln_clothing = 0 if clothing_adj == 0;


gen totcar_adj = autos_adj + carservs_adj + carprinc_adj + parts_adj;
gen ln_totcar = ln(totcar_adj);
replace ln_totcar = 0 if totcar_adj == 0;



gen car_adj = autos_adj;
gen ln_car = ln(car_adj);
replace ln_car = 0 if car_adj == 0;


gen totvis_full = clothes_pcare_adj + totcar_adj; 
gen ln_totvis_full = ln(totvis_full);

#delimit ;

replace ln_totvis_full = 0 if totvis_full == 0;


gen totvis = clothes_pcare_adj + car_adj; 
gen ln_totvis = ln(totvis);
replace ln_totvis = 0 if totvis == 0;


/*  All other expenditure measure */; 


gen food_adj = foodhome_adj + foodout_adj + foodwork_adj;  
gen ln_food = ln(food_adj); 
replace ln_food = 0 if food_adj == 0;

#delimit ;

gen ln_foodaway = ln(foodout_adj) ;
replace ln_foodaway = 0 if foodout_adj == 0 ;


gen rent_adj = renthome_adj + (homeval2_real) * 3;  
gen ln_rent = ln(rent_adj);
replace ln_rent = 0 if rent_adj == 0;


gen health_adj = drugs_adj + orthopd_adj + doctors_adj + hospital_adj + nurshome_adj + helthins_adj;  
gen ln_health = ln(health_adj);
replace ln_health = 0 if health_adj == 0;


gen utility_adj = elect_adj + gas_adj + water_adj + homefuel_adj + telephon_adj;
gen ln_utility = ln(utility_adj);
replace ln_utility = 0 if utility_adj == 0;


gen education_adj = highedu_adj + lowedu_adj + othedu_adj; 
gen ln_education = ln(education_adj); 
replace ln_education = 0 if education_adj == 0;



gen car_serv_adj = carservs_adj + carprinc_adj + parts_adj;
gen ln_car_serv = ln(car_serv_adj);
replace ln_car_serv = 0 if car_serv_adj == 0;

gen transport_other_adj = gasoline_adj + tolls_adj + autoins_adj + masstran_adj + othtrans_adj ;
gen ln_transport_other = ln(transport_other_adj); 
replace ln_transport_other = 0 if transport_other_adj == 0;

gen alc_tob_adj = tobacco_adj + alcohol_adj + niteclub_adj ;
gen ln_alc_tob = ln(alc_tob_adj);
replace ln_alc_tob = 0 if alc_tob_adj == 0;

gen ent_serv_adj = othrec_adj + airfare_adj + books_adj + pubs_adj ; 
gen ln_ent_serv = ln(ent_serv_adj); 
replace ln_ent_serv = 0 if ent_serv_adj == 0; 
 

#delimit ;

gen other_adj =  busiserv_adj + lifeins_adj + gambling_adj + charity_adj + housuppl_adj + servants_adj + rentothr_adj; 
gen ln_other = ln(other_adj); 
replace ln_other = 0 if other_adj == 0; 

gen ent_dur_adj = recsport_adj;
gen ln_ent_dur = ln(ent_dur_adj); 
replace ln_ent_dur = 0 if ent_dur_adj == 0; 

gen ln_furnish = ln(furnish_adj); 
replace ln_furnish = 0 if furnish_adj == 0; 



/* gen total expenditures */ ; 

#delimit ;

gen totexp = totvis + rent_adj +  food_adj + utility_adj + car_serv_adj + transport_other_adj +
ent_serv_adj + health_adj +  furnish_adj + education_adj +  ent_dur_adj + alc_tob_adj  + other_adj ;

gen ln_totexp = ln(totexp);
replace ln_totexp = 0 if totexp == 0;

gen totexp_sq = totexp^2;
gen totexp_cube = totexp^3;


/* gen expenditure share and fraction positive shares */ ; 

#delimit ;

gen vis_share = totvis/totexp;
gen rent_share = rent_adj/totexp;
gen food_share = food_adj/totexp;
gen utility_share = utility_adj/totexp;
gen transport_other_share = transport_other_adj/totexp;
gen car_serv_share = car_serv_adj/totexp;
gen ent_serv_share = ent_serv_adj/totexp;
gen health_share = health_adj/totexp;
gen furnish_share = furnish_adj/totexp;
gen education_share = education_adj/totexp;
gen ent_dur_share = ent_dur_adj/totexp;
gen alc_tob_share = alc_tob_adj/totexp;
gen other_share = other_adj/totexp;

gen totvis_pos = totvis > 0;
gen rent_pos = rent_adj > 0;
gen food_pos = food_adj > 0;
gen utility_pos = utility_adj > 0;
gen transport_other_pos = transport_other_adj > 0;
gen car_serv_pos = car_serv_adj > 0;
gen ent_serv_pos = ent_serv_adj > 0;
gen health_pos = health_adj > 0;
gen furnish_pos = furnish_adj > 0;
gen education_pos = education_adj > 0;
gen ent_dur_pos = ent_dur_adj > 0;
gen alc_tob_pos = alc_tob_adj > 0;
gen other_pos = other_adj > 0;


/*  Other controls */ ; 



gen liquid_wealth = checkng + saving + securi ;

gen liquid_pos = liquid_wealth > 0 & liquid_wealth < . ;
gen ln_liquid = ln(liquid_wealth) ;
replace ln_liquid = 0 if liquid_wealth == 0 ;

gen debt = owe1;
gen debt_pos = debt > 0 & debt < .;
gen ln_debt = ln(debt);
replace ln_debt = 0 if debt == 0 ;


/*  Restricting the sample to people with a clear race identifier */ ; 


#delimit ;

keep if  (black == 1 | black == 0)  & ( white == 1 | white ==0 ) & (hispanic == 1 | hispanic == 0) & (asian == 1 | asian == 0); 

             /*   -------------------- state incomes ------------------------------ */ ;

#delimit ;

/* generate income of black men, white men, hispanic men and all men by state */;

gen 	st_inc_bm =	.			;
replace	st_inc_bm =	21647	if  state ==	1	;
replace	st_inc_bm =	37209	if  state ==	2	;
replace	st_inc_bm =	30806	if  state ==	4	;
replace	st_inc_bm =	19141	if  state ==	5	;
replace	st_inc_bm =	31546	if  state ==	6	;
replace	st_inc_bm =	28334	if  state ==	8	;
replace	st_inc_bm =	29852	if  state ==	9	;
replace	st_inc_bm =	28838	if  state ==	10	;
replace	st_inc_bm =	26149	if  state ==	11	;
replace	st_inc_bm =	24122	if  state ==	12	;
replace	st_inc_bm =	27991	if  state ==	13	;
replace	st_inc_bm =	33795	if  state ==	15	;
replace	st_inc_bm =	.	if  state ==	16	;
replace	st_inc_bm =	26913	if  state ==	17	;
replace	st_inc_bm =	26248	if  state ==	18	;
replace	st_inc_bm =	24162	if  state ==	19	;
replace	st_inc_bm =	30879	if  state ==	20	;
replace	st_inc_bm =	25573	if  state ==	21	;
replace	st_inc_bm =	20925	if  state ==	22	;
replace	st_inc_bm =	.	if  state ==	23	;
replace	st_inc_bm =	31610	if  state ==	24	;
replace	st_inc_bm =	28860	if  state ==	25	;
replace	st_inc_bm =	27083	if  state ==	26	;
replace	st_inc_bm =	26506	if  state ==	27	;
replace	st_inc_bm =	21129	if  state ==	28	;
replace	st_inc_bm =	25865	if  state ==	29	;
replace	st_inc_bm =	25894	if  state ==	31	;
replace	st_inc_bm =	29488	if  state ==	32	;
replace	st_inc_bm =	.	if  state ==	33	;
replace	st_inc_bm =	31612	if  state ==	34	;
replace	st_inc_bm =	36467	if  state ==	35	;
replace	st_inc_bm =	26195	if  state ==	36	;
replace	st_inc_bm =	25924	if  state ==	37	;
replace	st_inc_bm =	26452	if  state ==	39	;
replace	st_inc_bm =	25219	if  state ==	40	;
replace	st_inc_bm =	33896	if  state ==	41	;
replace	st_inc_bm =	23615	if  state ==	42	;
replace	st_inc_bm =	24935	if  state ==	45	;
replace	st_inc_bm =	.	if  state ==	46	;
replace	st_inc_bm =	25251	if  state ==	47	;
replace	st_inc_bm =	27454	if  state ==	48	;
replace	st_inc_bm =	.	if  state ==	50	;
replace	st_inc_bm =	28832	if  state ==	51	;
replace	st_inc_bm =	39497	if  state ==	53	;
replace	st_inc_bm =	.	if  state ==	54	;
replace	st_inc_bm =	24039	if  state ==	55	;
replace	st_inc_bm =	29163	if  state ==	44	;
replace	st_inc_bm =	.	if  state ==	38	;
replace	st_inc_bm =	.	if  state ==	30	;
replace	st_inc_bm =	.	if  state ==	56	;
replace	st_inc_bm =	24906	if  state ==	49	;

gen 	st_inc_wm =	.			;
replace	st_inc_wm =	39888	if  state ==	1	;
replace	st_inc_wm =	50685	if  state ==	2	;
replace	st_inc_wm =	44097	if  state ==	4	;
replace	st_inc_wm =	33617	if  state ==	5	;
replace	st_inc_wm =	51700	if  state ==	6	;
replace	st_inc_wm =	46131	if  state ==	8	;
replace	st_inc_wm =	54255	if  state ==	9	;
replace	st_inc_wm =	43707	if  state ==	10	;
replace	st_inc_wm =	69075	if  state ==	11	;
replace	st_inc_wm =	42081	if  state ==	12	;
replace	st_inc_wm =	46054	if  state ==	13	;
replace	st_inc_wm =	41320	if  state ==	15	;
replace	st_inc_wm =	38406	if  state ==	16	;
replace	st_inc_wm =	49056	if  state ==	17	;
replace	st_inc_wm =	39944	if  state ==	18	;
replace	st_inc_wm =	38424	if  state ==	19	;
replace	st_inc_wm =	41807	if  state ==	20	;
replace	st_inc_wm =	37729	if  state ==	21	;
replace	st_inc_wm =	41325	if  state ==	22	;
replace	st_inc_wm =	36878	if  state ==	23	;
replace	st_inc_wm =	52261	if  state ==	24	;
replace	st_inc_wm =	48588	if  state ==	25	;
replace	st_inc_wm =	46220	if  state ==	26	;
replace	st_inc_wm =	45615	if  state ==	27	;
replace	st_inc_wm =	36881	if  state ==	28	;
replace	st_inc_wm =	40209	if  state ==	29	;
replace	st_inc_wm =	39923	if  state ==	31	;
replace	st_inc_wm =	45274	if  state ==	32	;
replace	st_inc_wm =	45499	if  state ==	33	;
replace	st_inc_wm =	55801	if  state ==	34	;
replace	st_inc_wm =	39917	if  state ==	35	;
replace	st_inc_wm =	49755	if  state ==	36	;
replace	st_inc_wm =	41689	if  state ==	37	;
replace	st_inc_wm =	43447	if  state ==	39	;
replace	st_inc_wm =	38460	if  state ==	40	;
replace	st_inc_wm =	42102	if  state ==	41	;
replace	st_inc_wm =	43409	if  state ==	42	;
replace	st_inc_wm =	41027	if  state ==	45	;
replace	st_inc_wm =	35620	if  state ==	46	;
replace	st_inc_wm =	38024	if  state ==	47	;
replace	st_inc_wm =	48618	if  state ==	48	;
replace	st_inc_wm =	39136	if  state ==	50	;
replace	st_inc_wm =	49495	if  state ==	51	;
replace	st_inc_wm =	45622	if  state ==	53	;
replace	st_inc_wm =	31535	if  state ==	54	;
replace	st_inc_wm =	42977	if  state ==	55	;
replace	st_inc_wm =	43832	if  state ==	44	;
replace	st_inc_wm =	.	if  state ==	38	;
replace	st_inc_wm =	.	if  state ==	30	;
replace	st_inc_wm =	.	if  state ==	56	;
replace	st_inc_wm =	40978	if  state ==	49	;

gen 	st_inc_hm =	.			;
replace	st_inc_hm =	32955	if  state ==	1	;
replace	st_inc_hm =	46671	if  state ==	2	;
replace	st_inc_hm =	26377	if  state ==	4	;
replace	st_inc_hm =	25874	if  state ==	5	;
replace	st_inc_hm =	27770	if  state ==	6	;
replace	st_inc_hm =	40037	if  state ==	8	;
replace	st_inc_hm =	33300	if  state ==	9	;
replace	st_inc_hm =	31657	if  state ==	10	;
replace	st_inc_hm =	27022	if  state ==	11	;
replace	st_inc_hm =	30201	if  state ==	12	;
replace	st_inc_hm =	30558	if  state ==	13	;
replace	st_inc_hm =	32722	if  state ==	15	;
replace	st_inc_hm =	24403	if  state ==	16	;
replace	st_inc_hm =	30387	if  state ==	17	;
replace	st_inc_hm =	30961	if  state ==	18	;
replace	st_inc_hm =	23749	if  state ==	19	;
replace	st_inc_hm =	31289	if  state ==	20	;
replace	st_inc_hm =	29029	if  state ==	21	;
replace	st_inc_hm =	36517	if  state ==	22	;
replace	st_inc_hm =	.	if  state ==	23	;
replace	st_inc_hm =	35464	if  state ==	24	;
replace	st_inc_hm =	25425	if  state ==	25	;
replace	st_inc_hm =	35356	if  state ==	26	;
replace	st_inc_hm =	27090	if  state ==	27	;
replace	st_inc_hm =	35247	if  state ==	28	;
replace	st_inc_hm =	31191	if  state ==	29	;
replace	st_inc_hm =	26480	if  state ==	31	;
replace	st_inc_hm =	29606	if  state ==	32	;
replace	st_inc_hm =	.	if  state ==	33	;
replace	st_inc_hm =	32158	if  state ==	34	;
replace	st_inc_hm =	26810	if  state ==	35	;
replace	st_inc_hm =	27617	if  state ==	36	;
replace	st_inc_hm =	25643	if  state ==	37	;
replace	st_inc_hm =	35228	if  state ==	39	;
replace	st_inc_hm =	25768	if  state ==	40	;
replace	st_inc_hm =	23143	if  state ==	41	;
replace	st_inc_hm =	31342	if  state ==	42	;
replace	st_inc_hm =	34215	if  state ==	45	;
replace	st_inc_hm =	.	if  state ==	46	;
replace	st_inc_hm =	28284	if  state ==	47	;
replace	st_inc_hm =	26862	if  state ==	48	;
replace	st_inc_hm =	.	if  state ==	50	;
replace	st_inc_hm =	39968	if  state ==	51	;
replace	st_inc_hm =	27243	if  state ==	53	;
replace	st_inc_hm =	.	if  state ==	54	;
replace	st_inc_hm =	27673	if  state ==	55	;
replace	st_inc_hm =	23190	if  state ==	44	;
replace	st_inc_hm =	.	if  state ==	38	;
replace	st_inc_hm =	.	if  state ==	30	;
replace	st_inc_hm =	.	if  state ==	56	;
replace	st_inc_hm =	27639	if  state ==	49	;

gen 	st_inc_m =	.			;
replace	st_inc_m =	35140	if  state ==	1	;
replace	st_inc_m =	47267	if  state ==	2	;
replace	st_inc_m =	38277	if  state ==	4	;
replace	st_inc_m =	31239	if  state ==	5	;
replace	st_inc_m =	41573	if  state ==	6	;
replace	st_inc_m =	44559	if  state ==	8	;
replace	st_inc_m =	50882	if  state ==	9	;
replace	st_inc_m =	40314	if  state ==	10	;
replace	st_inc_m =	41061	if  state ==	11	;
replace	st_inc_m =	37308	if  state ==	12	;
replace	st_inc_m =	39685	if  state ==	13	;
replace	st_inc_m =	38313	if  state ==	15	;
replace	st_inc_m =	36674	if  state ==	16	;
replace	st_inc_m =	43694	if  state ==	17	;
replace	st_inc_m =	38824	if  state ==	18	;
replace	st_inc_m =	37602	if  state ==	19	;
replace	st_inc_m =	40172	if  state ==	20	;
replace	st_inc_m =	36957	if  state ==	21	;
replace	st_inc_m =	36092	if  state ==	22	;
replace	st_inc_m =	36902	if  state ==	23	;
replace	st_inc_m =	45683	if  state ==	24	;
replace	st_inc_m =	46221	if  state ==	25	;
replace	st_inc_m =	43783	if  state ==	26	;
replace	st_inc_m =	44124	if  state ==	27	;
replace	st_inc_m =	31323	if  state ==	28	;
replace	st_inc_m =	38659	if  state ==	29	;
replace	st_inc_m =	38631	if  state ==	31	;
replace	st_inc_m =	41018	if  state ==	32	;
replace	st_inc_m =	45466	if  state ==	33	;
replace	st_inc_m =	49793	if  state ==	34	;
replace	st_inc_m =	32781	if  state ==	35	;
replace	st_inc_m =	42501	if  state ==	36	;
replace	st_inc_m =	37609	if  state ==	37	;
replace	st_inc_m =	41607	if  state ==	39	;
replace	st_inc_m =	35717	if  state ==	40	;
replace	st_inc_m =	39934	if  state ==	41	;
replace	st_inc_m =	41302	if  state ==	42	;
replace	st_inc_m =	36548	if  state ==	45	;
replace	st_inc_m =	34286	if  state ==	46	;
replace	st_inc_m =	35869	if  state ==	47	;
replace	st_inc_m =	38940	if  state ==	48	;
replace	st_inc_m =	38940	if  state ==	50	;
replace	st_inc_m =	44820	if  state ==	51	;
replace	st_inc_m =	44090	if  state ==	53	;
replace	st_inc_m =	31421	if  state ==	54	;
replace	st_inc_m =	41468	if  state ==	55	;
replace	st_inc_m =	41577	if  state ==	44	;
replace	st_inc_m =	.	if  state ==	38	;
replace	st_inc_m =	.	if  state ==	30	;
replace	st_inc_m =	.	if  state ==	56	;
replace	st_inc_m =	39492	if  state ==	49	;

#delimit ;

/* generate standard deviation of incomes for black men, white men, hispanic men and all men by state */;

gen 	st_sdinc_bm =	.			;
replace	st_sdinc_bm =	21051	if  state ==	1	;
replace	st_sdinc_bm =	45697	if  state ==	2	;
replace	st_sdinc_bm =	23711	if  state ==	4	;
replace	st_sdinc_bm =	16936	if  state ==	5	;
replace	st_sdinc_bm =	35509	if  state ==	6	;
replace	st_sdinc_bm =	28334	if  state ==	8	;
replace	st_sdinc_bm =	26572	if  state ==	9	;
replace	st_sdinc_bm =	26394	if  state ==	10	;
replace	st_sdinc_bm =	29001	if  state ==	11	;
replace	st_sdinc_bm =	26049	if  state ==	12	;
replace	st_sdinc_bm =	31316	if  state ==	13	;
replace	st_sdinc_bm =	26252	if  state ==	15	;
replace	st_sdinc_bm =	.	if  state ==	16	;
replace	st_sdinc_bm =	28571	if  state ==	17	;
replace	st_sdinc_bm =	23618	if  state ==	18	;
replace	st_sdinc_bm =	17721	if  state ==	19	;
replace	st_sdinc_bm =	22726	if  state ==	20	;
replace	st_sdinc_bm =	19990	if  state ==	21	;
replace	st_sdinc_bm =	19873	if  state ==	22	;
replace	st_sdinc_bm =	.	if  state ==	23	;
replace	st_sdinc_bm =	34881	if  state ==	24	;
replace	st_sdinc_bm =	36597	if  state ==	25	;
replace	st_sdinc_bm =	27057	if  state ==	26	;
replace	st_sdinc_bm =	27883	if  state ==	27	;
replace	st_sdinc_bm =	23987	if  state ==	28	;
replace	st_sdinc_bm =	21715	if  state ==	29	;
replace	st_sdinc_bm =	18204	if  state ==	31	;
replace	st_sdinc_bm =	31453	if  state ==	32	;
replace	st_sdinc_bm =	.	if  state ==	33	;
replace	st_sdinc_bm =	38829	if  state ==	34	;
replace	st_sdinc_bm =	29289	if  state ==	35	;
replace	st_sdinc_bm =	28516	if  state ==	36	;
replace	st_sdinc_bm =	25063	if  state ==	37	;
replace	st_sdinc_bm =	27730	if  state ==	39	;
replace	st_sdinc_bm =	21761	if  state ==	40	;
replace	st_sdinc_bm =	52965	if  state ==	41	;
replace	st_sdinc_bm =	26440	if  state ==	42	;
replace	st_sdinc_bm =	23195	if  state ==	45	;
replace	st_sdinc_bm =	.	if  state ==	46	;
replace	st_sdinc_bm =	27638	if  state ==	47	;
replace	st_sdinc_bm =	31017	if  state ==	48	;
replace	st_sdinc_bm =	.	if  state ==	50	;
replace	st_sdinc_bm =	26812	if  state ==	51	;
replace	st_sdinc_bm =	72942	if  state ==	53	;
replace	st_sdinc_bm =	.	if  state ==	54	;
replace	st_sdinc_bm =	23269	if  state ==	55	;
replace	st_sdinc_bm =	26559	if  state ==	44	;
replace	st_sdinc_bm =	.	if  state ==	38	;
replace	st_sdinc_bm =	.	if  state ==	30	;
replace	st_sdinc_bm =	.	if  state ==	56	;
replace	st_sdinc_bm =	18628	if  state ==	49	;

gen 	st_sdinc_wm =	.			;
replace	st_sdinc_wm =	39466	if  state ==	1	;
replace	st_sdinc_wm =	41746	if  state ==	2	;
replace	st_sdinc_wm =	45349	if  state ==	4	;
replace	st_sdinc_wm =	32362	if  state ==	5	;
replace	st_sdinc_wm =	51258	if  state ==	6	;
replace	st_sdinc_wm =	46131	if  state ==	8	;
replace	st_sdinc_wm =	50977	if  state ==	9	;
replace	st_sdinc_wm =	38658	if  state ==	10	;
replace	st_sdinc_wm =	70168	if  state ==	11	;
replace	st_sdinc_wm =	43983	if  state ==	12	;
replace	st_sdinc_wm =	44638	if  state ==	13	;
replace	st_sdinc_wm =	37909	if  state ==	15	;
replace	st_sdinc_wm =	38406	if  state ==	16	;
replace	st_sdinc_wm =	45530	if  state ==	17	;
replace	st_sdinc_wm =	36136	if  state ==	18	;
replace	st_sdinc_wm =	34654	if  state ==	19	;
replace	st_sdinc_wm =	41554	if  state ==	20	;
replace	st_sdinc_wm =	39274	if  state ==	21	;
replace	st_sdinc_wm =	40958	if  state ==	22	;
replace	st_sdinc_wm =	36294	if  state ==	23	;
replace	st_sdinc_wm =	46933	if  state ==	24	;
replace	st_sdinc_wm =	46837	if  state ==	25	;
replace	st_sdinc_wm =	42146	if  state ==	26	;
replace	st_sdinc_wm =	45526	if  state ==	27	;
replace	st_sdinc_wm =	33835	if  state ==	28	;
replace	st_sdinc_wm =	39684	if  state ==	29	;
replace	st_sdinc_wm =	35667	if  state ==	31	;
replace	st_sdinc_wm =	42210	if  state ==	32	;
replace	st_sdinc_wm =	43213	if  state ==	33	;
replace	st_sdinc_wm =	52294	if  state ==	34	;
replace	st_sdinc_wm =	34715	if  state ==	35	;
replace	st_sdinc_wm =	50723	if  state ==	36	;
replace	st_sdinc_wm =	40379	if  state ==	37	;
replace	st_sdinc_wm =	40354	if  state ==	39	;
replace	st_sdinc_wm =	39020	if  state ==	40	;
replace	st_sdinc_wm =	40650	if  state ==	41	;
replace	st_sdinc_wm =	42314	if  state ==	42	;
replace	st_sdinc_wm =	35807	if  state ==	45	;
replace	st_sdinc_wm =	33444	if  state ==	46	;
replace	st_sdinc_wm =	41936	if  state ==	47	;
replace	st_sdinc_wm =	46406	if  state ==	48	;
replace	st_sdinc_wm =	36128	if  state ==	50	;
replace	st_sdinc_wm =	47294	if  state ==	51	;
replace	st_sdinc_wm =	43212	if  state ==	53	;
replace	st_sdinc_wm =	30120	if  state ==	54	;
replace	st_sdinc_wm =	37463	if  state ==	55	;
replace	st_sdinc_wm =	43378	if  state ==	44	;
replace	st_sdinc_wm =	.	if  state ==	38	;
replace	st_sdinc_wm =	.	if  state ==	30	;
replace	st_sdinc_wm =	.	if  state ==	56	;
replace	st_sdinc_wm =	39610	if  state ==	49	;

gen 	st_sdinc_hm =	.			;
replace	st_sdinc_hm =	44960	if  state ==	1	;
replace	st_sdinc_hm =	47670	if  state ==	2	;
replace	st_sdinc_hm =	28239	if  state ==	4	;
replace	st_sdinc_hm =	32383	if  state ==	5	;
replace	st_sdinc_hm =	29900	if  state ==	6	;
replace	st_sdinc_hm =	40037	if  state ==	8	;
replace	st_sdinc_hm =	44487	if  state ==	9	;
replace	st_sdinc_hm =	28640	if  state ==	10	;
replace	st_sdinc_hm =	35565	if  state ==	11	;
replace	st_sdinc_hm =	33436	if  state ==	12	;
replace	st_sdinc_hm =	27634	if  state ==	13	;
replace	st_sdinc_hm =	24441	if  state ==	15	;
replace	st_sdinc_hm =	24403	if  state ==	16	;
replace	st_sdinc_hm =	33344	if  state ==	17	;
replace	st_sdinc_hm =	33757	if  state ==	18	;
replace	st_sdinc_hm =	17827	if  state ==	19	;
replace	st_sdinc_hm =	30057	if  state ==	20	;
replace	st_sdinc_hm =	29144	if  state ==	21	;
replace	st_sdinc_hm =	38038	if  state ==	22	;
replace	st_sdinc_hm =	.	if  state ==	23	;
replace	st_sdinc_hm =	40187	if  state ==	24	;
replace	st_sdinc_hm =	24417	if  state ==	25	;
replace	st_sdinc_hm =	38353	if  state ==	26	;
replace	st_sdinc_hm =	24320	if  state ==	27	;
replace	st_sdinc_hm =	34185	if  state ==	28	;
replace	st_sdinc_hm =	39450	if  state ==	29	;
replace	st_sdinc_hm =	19648	if  state ==	31	;
replace	st_sdinc_hm =	32444	if  state ==	32	;
replace	st_sdinc_hm =	.	if  state ==	33	;
replace	st_sdinc_hm =	33924	if  state ==	34	;
replace	st_sdinc_hm =	25113	if  state ==	35	;
replace	st_sdinc_hm =	32522	if  state ==	36	;
replace	st_sdinc_hm =	24808	if  state ==	37	;
replace	st_sdinc_hm =	34044	if  state ==	39	;
replace	st_sdinc_hm =	23162	if  state ==	40	;
replace	st_sdinc_hm =	22362	if  state ==	41	;
replace	st_sdinc_hm =	36267	if  state ==	42	;
replace	st_sdinc_hm =	38009	if  state ==	45	;
replace	st_sdinc_hm =	.	if  state ==	46	;
replace	st_sdinc_hm =	27386	if  state ==	47	;
replace	st_sdinc_hm =	30777	if  state ==	48	;
replace	st_sdinc_hm =	.	if  state ==	50	;
replace	st_sdinc_hm =	44431	if  state ==	51	;
replace	st_sdinc_hm =	25915	if  state ==	53	;
replace	st_sdinc_hm =	.	if  state ==	54	;
replace	st_sdinc_hm =	29724	if  state ==	55	;
replace	st_sdinc_hm =	29766	if  state ==	44	;
replace	st_sdinc_hm =	.	if  state ==	38	;
replace	st_sdinc_hm =	.	if  state ==	30	;
replace	st_sdinc_hm =	.	if  state ==	56	;
replace	st_sdinc_hm =	22870	if  state ==	49	;

gen 	st_sdinc_m =	.			;
replace	st_sdinc_m =	36691	if  state ==	1	;
replace	st_sdinc_m =	41543	if  state ==	2	;
replace	st_sdinc_m =	41034	if  state ==	4	;
replace	st_sdinc_m =	31220	if  state ==	5	;
replace	st_sdinc_m =	45295	if  state ==	6	;
replace	st_sdinc_m =	44964	if  state ==	8	;
replace	st_sdinc_m =	49762	if  state ==	9	;
replace	st_sdinc_m =	36650	if  state ==	10	;
replace	st_sdinc_m =	52077	if  state ==	11	;
replace	st_sdinc_m =	40565	if  state ==	12	;
replace	st_sdinc_m =	41073	if  state ==	13	;
replace	st_sdinc_m =	34604	if  state ==	15	;
replace	st_sdinc_m =	36674	if  state ==	16	;
replace	st_sdinc_m =	43229	if  state ==	17	;
replace	st_sdinc_m =	35647	if  state ==	18	;
replace	st_sdinc_m =	34087	if  state ==	19	;
replace	st_sdinc_m =	39854	if  state ==	20	;
replace	st_sdinc_m =	38686	if  state ==	21	;
replace	st_sdinc_m =	38486	if  state ==	22	;
replace	st_sdinc_m =	36408	if  state ==	23	;
replace	st_sdinc_m =	44677	if  state ==	24	;
replace	st_sdinc_m =	45611	if  state ==	25	;
replace	st_sdinc_m =	41209	if  state ==	26	;
replace	st_sdinc_m =	44637	if  state ==	27	;
replace	st_sdinc_m =	31629	if  state ==	28	;
replace	st_sdinc_m =	38874	if  state ==	29	;
replace	st_sdinc_m =	34712	if  state ==	31	;
replace	st_sdinc_m =	40454	if  state ==	32	;
replace	st_sdinc_m =	43096	if  state ==	33	;
replace	st_sdinc_m =	50433	if  state ==	34	;
replace	st_sdinc_m =	30887	if  state ==	35	;
replace	st_sdinc_m =	46969	if  state ==	36	;
replace	st_sdinc_m =	37742	if  state ==	37	;
replace	st_sdinc_m =	39563	if  state ==	39	;
replace	st_sdinc_m =	36905	if  state ==	40	;
replace	st_sdinc_m =	39984	if  state ==	41	;
replace	st_sdinc_m =	41497	if  state ==	42	;
replace	st_sdinc_m =	33934	if  state ==	45	;
replace	st_sdinc_m =	32777	if  state ==	46	;
replace	st_sdinc_m =	40252	if  state ==	47	;
replace	st_sdinc_m =	41633	if  state ==	48	;
replace	st_sdinc_m =	38940	if  state ==	50	;
replace	st_sdinc_m =	44431	if  state ==	51	;
replace	st_sdinc_m =	43576	if  state ==	53	;
replace	st_sdinc_m =	30846	if  state ==	54	;
replace	st_sdinc_m =	36912	if  state ==	55	;
replace	st_sdinc_m =	42201	if  state ==	44	;
replace	st_sdinc_m =	.	if  state ==	38	;
replace	st_sdinc_m =	.	if  state ==	30	;
replace	st_sdinc_m =	.	if  state ==	56	;
replace	st_sdinc_m =	38348	if  state ==	49	;

#delimit ;

gen	state_1 =	0		;
gen	state_2 =	0		;
gen	state_4 =	0		;
gen	state_5 =	0		;
gen	state_6 =	0		;
gen	state_8 =	0		;
gen	state_9 =	0		;
gen	state_10 =	0		;
gen	state_11 =	0		;
gen	state_12 =	0		;
gen	state_13 =	0		;
gen	state_15 =	0		;
gen	state_16 =	0		;
gen	state_17 =	0		;
gen	state_18 =	0		;
gen	state_19 =	0		;
gen	state_20 =	0		;
gen	state_21 =	0		;
gen	state_22 =	0		;
gen	state_23 =	0		;
gen	state_24 =	0		;
gen	state_25 =	0		;
gen	state_26 =	0		;
gen	state_27 =	0		;
gen	state_28 =	0		;
gen	state_29 =	0		;
gen	state_31 =	0		;
gen	state_32 =	0		;
gen	state_33 =	0		;
gen	state_34 =	0		;
gen	state_35 =	0		;
gen	state_36 =	0		;
gen	state_37 =	0		;
gen	state_39 =	0		;
gen	state_40 =	0		;
gen	state_41 =	0		;
gen	state_42 =	0		;
gen	state_45 =	0		;
gen	state_46 =	0		;
gen	state_47 =	0		;
gen	state_48 =	0		;
gen	state_50 =	0		;
gen	state_51 =	0		;
gen	state_53 =	0		;
gen	state_54 =	0		;
gen	state_55 =	0		;
gen	state_44 =	0		;
gen	state_38 =	0		;
gen	state_30 =	0		;
gen	state_56 =	0		;
gen	state_49 =	0		;

replace	state_1 =	1	if  state ==	1	;
replace	state_2 =	1	if  state ==	2	;
replace	state_4 =	1	if  state ==	4	;
replace	state_5 =	1	if  state ==	5	;
replace	state_6 =	1	if  state ==	6	;
replace	state_8 =	1	if  state ==	8	;
replace	state_9 =	1	if  state ==	9	;
replace	state_10 =	1	if  state ==	10	;
replace	state_11 =	1	if  state ==	11	;
replace	state_12 =	1	if  state ==	12	;
replace	state_13 =	1	if  state ==	13	;
replace	state_15 =	1	if  state ==	15	;
replace	state_16 =	1	if  state ==	16	;
replace	state_17 =	1	if  state ==	17	;
replace	state_18 =	1	if  state ==	18	;
replace	state_19 =	1	if  state ==	19	;
replace	state_20 =	1	if  state ==	20	;
replace	state_21 =	1	if  state ==	21	;
replace	state_22 =	1	if  state ==	22	;
replace	state_23 =	1	if  state ==	23	;
replace	state_24 =	1	if  state ==	24	;
replace	state_25 =	1	if  state ==	25	;
replace	state_26 =	1	if  state ==	26	;
replace	state_27 =	1	if  state ==	27	;
replace	state_28 =	1	if  state ==	28	;
replace	state_29 =	1	if  state ==	29	;
replace	state_31 =	1	if  state ==	31	;
replace	state_32 =	1	if  state ==	32	;
replace	state_33 =	1	if  state ==	33	;
replace	state_34 =	1	if  state ==	34	;
replace	state_35 =	1	if  state ==	35	;
replace	state_36 =	1	if  state ==	36	;
replace	state_37 =	1	if  state ==	37	;
replace	state_39 =	1	if  state ==	39	;
replace	state_40 =	1	if  state ==	40	;
replace	state_41 =	1	if  state ==	41	;
replace	state_42 =	1	if  state ==	42	;
replace	state_45 =	1	if  state ==	45	;
replace	state_46 =	1	if  state ==	46	;
replace	state_47 =	1	if  state ==	47	;
replace	state_48 =	1	if  state ==	48	;
replace	state_50 =	1	if  state ==	50	;
replace	state_51 =	1	if  state ==	51	;
replace	state_53 =	1	if  state ==	53	;
replace	state_54 =	1	if  state ==	54	;
replace	state_55 =	1	if  state ==	55	;
replace	state_44 =	1	if  state ==	44	;
replace	state_38 =	1	if  state ==	38	;
replace	state_30 =	1	if  state ==	30	;
replace	state_56 =	1	if  state ==	56	;
replace	state_49 =	1	if  state ==	49	;



#delimit ; 



gen ln_st_inc_mean_wm = ln(st_inc_wm);
gen ln_st_inc_mean_bm = ln(st_inc_bm);
gen ln_st_inc_mean_hm = ln(st_inc_hm);
gen ln_st_inc_mean_m = ln(st_inc_m);

gen ln_st_inc_meanbyrace = .;

replace ln_st_inc_meanbyrace = ln_st_inc_mean_wm if white == 1;
replace ln_st_inc_meanbyrace = ln_st_inc_mean_bm if black == 1;
replace ln_st_inc_meanbyrace = ln_st_inc_mean_hm if hispanic == 1; 


#delimit ;


gen ln_st_sdinc_mean_wm = ln(st_sdinc_wm);
gen ln_st_sdinc_mean_bm = ln(st_sdinc_bm);
gen ln_st_sdinc_mean_hm = ln(st_sdinc_hm);
gen ln_st_sdinc_mean_m = ln(st_sdinc_m);

#delimit ;

gen sharp_wm =  ln_st_sdinc_mean_wm - ln_st_inc_mean_wm  ;
gen sharp_bm = ln_st_sdinc_mean_bm - ln_st_inc_mean_bm  ;
gen sharp_hm = ln_st_sdinc_mean_hm - ln_st_inc_mean_hm ;
gen sharp_m = ln_st_sdinc_mean_m  - ln_st_inc_mean_m ;

gen sharp_byrace = .;

replace sharp_byrace = sharp_wm if white == 1;
replace sharp_byrace = sharp_bm if black == 1;
replace sharp_byrace = sharp_hm if hispanic == 1; 


#delimit ;

/* enter in race share in state */; 

gen share_wm =	0.722	if state == 	1	;
replace share_wm = 	0.76	if state == 	2	;
replace share_wm = 	0.661	if state == 	4	;
replace share_wm = 	0.807	if state == 	5	;
replace share_wm = 	0.496	if state == 	6	;
replace share_wm = 	0.78	if state == 	8	;
replace share_wm = 	0.828	if state == 	9	;
replace share_wm = 	0.74	if state == 	10	;
replace share_wm = 	0.341	if state == 	11	;
replace share_wm = 	0.65	if state == 	12	;
replace share_wm = 	0.634	if state == 	13	;
replace share_wm = 	0.32	if state == 	15	;
replace share_wm = 	.	if state == 	16	;
replace share_wm = 	0.715	if state == 	17	;
replace share_wm = 	0.894	if state == 	18	;
replace share_wm = 	0.936	if state == 	19	;
replace share_wm = 	0.855	if state == 	20	;
replace share_wm = 	0.913	if state == 	21	;
replace share_wm = 	0.704	if state == 	22	;
replace share_wm = 	0.979	if state == 	23	;
replace share_wm = 	0.646	if state == 	24	;
replace share_wm = 	0.868	if state == 	25	;
replace share_wm = 	0.828	if state == 	26	;
replace share_wm = 	0.905	if state == 	27	;
replace share_wm = 	0.63	if state == 	28	;
replace share_wm = 	0.869	if state == 	29	;
replace share_wm = 	0.904	if state == 	31	;
replace share_wm = 	0.703	if state == 	32	;
replace share_wm = 	.	if state == 	33	;
replace share_wm = 	0.695	if state == 	34	;
replace share_wm = 	0.473	if state == 	35	;
replace share_wm = 	0.669	if state == 	36	;
replace share_wm = 	0.726	if state == 	37	;
replace share_wm = 	0.87	if state == 	39	;
replace share_wm = 	0.795	if state == 	40	;
replace share_wm = 	0.852	if state == 	41	;
replace share_wm = 	0.87	if state == 	42	;
replace share_wm = 	0.705	if state == 	45	;
replace share_wm = 	.	if state == 	46	;
replace share_wm = 	0.806	if state == 	47	;
replace share_wm = 	0.534	if state == 	48	;
replace share_wm = 	.	if state == 	50	;
replace share_wm = 	0.723	if state == 	51	;
replace share_wm = 	0.855	if state == 	53	;
replace share_wm = 	.	if state == 	54	;
replace share_wm = 	0.908	if state == 	55	;

gen share_bm =	0.257	if state == 	1	;
replace share_bm = 	0.05	if state == 	2	;
replace share_bm = 	0.027	if state == 	4	;
replace share_bm = 	0.151	if state == 	5	;
replace share_bm = 	0.062	if state == 	6	;
replace share_bm = 	0.037	if state == 	8	;
replace share_bm = 	0.081	if state == 	9	;
replace share_bm = 	0.197	if state == 	10	;
replace share_bm = 	0.554	if state == 	11	;
replace share_bm = 	0.133	if state == 	12	;
replace share_bm = 	0.311	if state == 	13	;
replace share_bm = 	0.03	if state == 	15	;
replace share_bm = 	.	if state == 	16	;
replace share_bm = 	0.136	if state == 	17	;
replace share_bm = 	0.062	if state == 	18	;
replace share_bm = 	0.017	if state == 	19	;
replace share_bm = 	0.067	if state == 	20	;
replace share_bm = 	0.063	if state == 	21	;
replace share_bm = 	0.253	if state == 	22	;
replace share_bm = 	0.002	if state == 	23	;
replace share_bm = 	0.268	if state == 	24	;
replace share_bm = 	0.045	if state == 	25	;
replace share_bm = 	0.113	if state == 	26	;
replace share_bm = 	0.036	if state == 	27	;
replace share_bm = 	0.349	if state == 	28	;
replace share_bm = 	0.091	if state == 	29	;
replace share_bm = 	0.032	if state == 	31	;
replace share_bm = 	0.054	if state == 	32	;
replace share_bm = 	.	if state == 	33	;
replace share_bm = 	0.122	if state == 	34	;
replace share_bm = 	0.016	if state == 	35	;
replace share_bm = 	0.13	if state == 	36	;
replace share_bm = 	0.207	if state == 	37	;
replace share_bm = 	0.095	if state == 	39	;
replace share_bm = 	0.078	if state == 	40	;
replace share_bm = 	0.017	if state == 	41	;
replace share_bm = 	0.084	if state == 	42	;
replace share_bm = 	0.273	if state == 	45	;
replace share_bm = 	.	if state == 	46	;
replace share_bm = 	0.164	if state == 	47	;
replace share_bm = 	0.111	if state == 	48	;
replace share_bm = 	.	if state == 	50	;
replace share_bm = 	0.191	if state == 	51	;
replace share_bm = 	0.031	if state == 	53	;
replace share_bm = 	.	if state == 	54	;
replace share_bm = 	0.045	if state == 	55	;

gen share_hm =	0.014	if state == 	1	;
replace share_hm = 	0.04	if state == 	2	;
replace share_hm = 	0.27	if state == 	4	;
replace share_hm = 	0.026	if state == 	5	;
replace share_hm = 	0.331	if state == 	6	;
replace share_hm = 	0.151	if state == 	8	;
replace share_hm = 	0.068	if state == 	9	;
replace share_hm = 	0.043	if state == 	10	;
replace share_hm = 	0.078	if state == 	11	;
replace share_hm = 	0.197	if state == 	12	;
replace share_hm = 	0.04	if state == 	13	;
replace share_hm = 	0.05	if state == 	15	;
replace share_hm = 	.	if state == 	16	;
replace share_hm = 	0.117	if state == 	17	;
replace share_hm = 	0.031	if state == 	18	;
replace share_hm = 	0.027	if state == 	19	;
replace share_hm = 	0.052	if state == 	20	;
replace share_hm = 	0.016	if state == 	21	;
replace share_hm = 	0.027	if state == 	22	;
replace share_hm = 	0.006	if state == 	23	;
replace share_hm = 	0.049	if state == 	24	;
replace share_hm = 	0.054	if state == 	25	;
replace share_hm = 	0.03	if state == 	26	;
replace share_hm = 	0.03	if state == 	27	;
replace share_hm = 	0.012	if state == 	28	;
replace share_hm = 	0.023	if state == 	29	;
replace share_hm = 	0.042	if state == 	31	;
replace share_hm = 	0.189	if state == 	32	;
replace share_hm = 	.	if state == 	33	;
replace share_hm = 	0.129	if state == 	34	;
replace share_hm = 	0.423	if state == 	35	;
replace share_hm = 	0.144	if state == 	36	;
replace share_hm = 	0.039	if state == 	37	;
replace share_hm = 	0.017	if state == 	39	;
replace share_hm = 	0.044	if state == 	40	;
replace share_hm = 	0.087	if state == 	41	;
replace share_hm = 	0.027	if state == 	42	;
replace share_hm = 	0.015	if state == 	45	;
replace share_hm = 	.	if state == 	46	;
replace share_hm = 	0.016	if state == 	47	;
replace share_hm = 	0.32	if state == 	48	;
replace share_hm = 	.	if state == 	50	;
replace share_hm = 	0.044	if state == 	51	;
replace share_hm = 	0.051	if state == 	53	;
replace share_hm = 	.	if state == 	54	;
replace share_hm = 	0.029	if state == 	55	;

gen 	st_inc_b_all =	.			;
replace	st_inc_b_all =	17231	if state == 	1	;
replace	st_inc_b_all =	29202	if state == 	2	;
replace	st_inc_b_all =	24383	if state == 	4	;
replace	st_inc_b_all =	15288	if state == 	5	;
replace	st_inc_b_all =	26260	if state == 	6	;
replace	st_inc_b_all =	22713	if state == 	8	;
replace	st_inc_b_all =	25127	if state == 	9	;
replace	st_inc_b_all =	24055	if state == 	10	;
replace	st_inc_b_all =	23707	if state == 	11	;
replace	st_inc_b_all =	19414	if state == 	12	;
replace	st_inc_b_all =	21973	if state == 	13	;
replace	st_inc_b_all =	26871	if state == 	15	;
replace	st_inc_b_all =	.	if state == 	16	;
replace	st_inc_b_all =	21598	if state == 	17	;
replace	st_inc_b_all =	22165	if state == 	18	;
replace	st_inc_b_all =	19808	if state == 	19	;
replace	st_inc_b_all =	22505	if state == 	20	;
replace	st_inc_b_all =	19127	if state == 	21	;
replace	st_inc_b_all =	15599	if state == 	22	;
replace	st_inc_b_all =	.	if state == 	23	;
replace	st_inc_b_all =	27076	if state == 	24	;
replace	st_inc_b_all =	22053	if state == 	25	;
replace	st_inc_b_all =	21698	if state == 	26	;
replace	st_inc_b_all =	21576	if state == 	27	;
replace	st_inc_b_all =	16198	if state == 	28	;
replace	st_inc_b_all =	20488	if state == 	29	;
replace	st_inc_b_all =	19989	if state == 	31	;
replace	st_inc_b_all =	22913	if state == 	32	;
replace	st_inc_b_all =	.	if state == 	33	;
replace	st_inc_b_all =	25132	if state == 	34	;
replace	st_inc_b_all =	27775	if state == 	35	;
replace	st_inc_b_all =	21290	if state == 	36	;
replace	st_inc_b_all =	20636	if state == 	37	;
replace	st_inc_b_all =	21281	if state == 	39	;
replace	st_inc_b_all =	18514	if state == 	40	;
replace	st_inc_b_all =	24407	if state == 	41	;
replace	st_inc_b_all =	20987	if state == 	42	;
replace	st_inc_b_all =	19330	if state == 	45	;
replace	st_inc_b_all =	.	if state == 	46	;
replace	st_inc_b_all =	19803	if state == 	47	;
replace	st_inc_b_all =	22138	if state == 	48	;
replace	st_inc_b_all =	.	if state == 	50	;
replace	st_inc_b_all =	23325	if state == 	51	;
replace	st_inc_b_all =	31399	if state == 	53	;
replace	st_inc_b_all =	.	if state == 	54	;
replace	st_inc_b_all =	19179	if state == 	55	;
replace	st_inc_b_all =	22544	if state == 	44	;
replace	st_inc_b_all =	.	if state == 	38	;
replace	st_inc_b_all =	.	if state == 	30	;
replace	st_inc_b_all =	.	if state == 	56	;
replace	st_inc_b_all =	20899	if state == 	49	;

gen 	st_inc_w_all =	.			;
replace	st_inc_w_all =	27785	if state == 	1	;
replace	st_inc_w_all =	37397	if state == 	2	;
replace	st_inc_w_all =	31927	if state == 	4	;
replace	st_inc_w_all =	24798	if state == 	5	;
replace	st_inc_w_all =	37835	if state == 	6	;
replace	st_inc_w_all =	35449	if state == 	8	;
replace	st_inc_w_all =	39481	if state == 	9	;
replace	st_inc_w_all =	32925	if state == 	10	;
replace	st_inc_w_all =	56905	if state == 	11	;
replace	st_inc_w_all =	31228	if state == 	12	;
replace	st_inc_w_all =	32771	if state == 	13	;
replace	st_inc_w_all =	31996	if state == 	15	;
replace	st_inc_w_all =	26718	if state == 	16	;
replace	st_inc_w_all =	34871	if state == 	17	;
replace	st_inc_w_all =	28469	if state == 	18	;
replace	st_inc_w_all =	28074	if state == 	19	;
replace	st_inc_w_all =	30061	if state == 	20	;
replace	st_inc_w_all =	26207	if state == 	21	;
replace	st_inc_w_all =	27806	if state == 	22	;
replace	st_inc_w_all =	26954	if state == 	23	;
replace	st_inc_w_all =	38621	if state == 	24	;
replace	st_inc_w_all =	36375	if state == 	25	;
replace	st_inc_w_all =	32362	if state == 	26	;
replace	st_inc_w_all =	33625	if state == 	27	;
replace	st_inc_w_all =	26052	if state == 	28	;
replace	st_inc_w_all =	29440	if state == 	29	;
replace	st_inc_w_all =	28538	if state == 	31	;
replace	st_inc_w_all =	33745	if state == 	32	;
replace	st_inc_w_all =	33245	if state == 	33	;
replace	st_inc_w_all =	39187	if state == 	34	;
replace	st_inc_w_all =	28651	if state == 	35	;
replace	st_inc_w_all =	35119	if state == 	36	;
replace	st_inc_w_all =	30461	if state == 	37	;
replace	st_inc_w_all =	30638	if state == 	39	;
replace	st_inc_w_all =	26707	if state == 	40	;
replace	st_inc_w_all =	30165	if state == 	41	;
replace	st_inc_w_all =	30764	if state == 	42	;
replace	st_inc_w_all =	29260	if state == 	45	;
replace	st_inc_w_all =	26235	if state == 	46	;
replace	st_inc_w_all =	27027	if state == 	47	;
replace	st_inc_w_all =	34466	if state == 	48	;
replace	st_inc_w_all =	28785	if state == 	50	;
replace	st_inc_w_all =	35202	if state == 	51	;
replace	st_inc_w_all =	32738	if state == 	53	;
replace	st_inc_w_all =	21684	if state == 	54	;
replace	st_inc_w_all =	31988	if state == 	55	;
replace	st_inc_w_all =	33043	if state == 	44	;
replace	st_inc_w_all =	25844	if state == 	38	;
replace	st_inc_w_all =	24222	if state == 	30	;
replace	st_inc_w_all =	27401	if state == 	56	;
replace	st_inc_w_all =	27437	if state == 	49	;

gen 	st_inc_h_all =	.			;
replace	st_inc_h_all =	17231	if state == 	1	;
replace	st_inc_h_all =	29202	if state == 	2	;
replace	st_inc_h_all =	24383	if state == 	4	;
replace	st_inc_h_all =	15288	if state == 	5	;
replace	st_inc_h_all =	26260	if state == 	6	;
replace	st_inc_h_all =	22713	if state == 	8	;
replace	st_inc_h_all =	25127	if state == 	9	;
replace	st_inc_h_all =	24055	if state == 	10	;
replace	st_inc_h_all =	23707	if state == 	11	;
replace	st_inc_h_all =	19414	if state == 	12	;
replace	st_inc_h_all =	21973	if state == 	13	;
replace	st_inc_h_all =	26871	if state == 	15	;
replace	st_inc_h_all =	.	if state == 	16	;
replace	st_inc_h_all =	21598	if state == 	17	;
replace	st_inc_h_all =	22165	if state == 	18	;
replace	st_inc_h_all =	19808	if state == 	19	;
replace	st_inc_h_all =	22505	if state == 	20	;
replace	st_inc_h_all =	19127	if state == 	21	;
replace	st_inc_h_all =	15599	if state == 	22	;
replace	st_inc_h_all =	.	if state == 	23	;
replace	st_inc_h_all =	27076	if state == 	24	;
replace	st_inc_h_all =	22053	if state == 	25	;
replace	st_inc_h_all =	21698	if state == 	26	;
replace	st_inc_h_all =	21576	if state == 	27	;
replace	st_inc_h_all =	16198	if state == 	28	;
replace	st_inc_h_all =	20488	if state == 	29	;
replace	st_inc_h_all =	19989	if state == 	31	;
replace	st_inc_h_all =	22913	if state == 	32	;
replace	st_inc_h_all =	.	if state == 	33	;
replace	st_inc_h_all =	25132	if state == 	34	;
replace	st_inc_h_all =	27775	if state == 	35	;
replace	st_inc_h_all =	21290	if state == 	36	;
replace	st_inc_h_all =	20636	if state == 	37	;
replace	st_inc_h_all =	21281	if state == 	39	;
replace	st_inc_h_all =	18514	if state == 	40	;
replace	st_inc_h_all =	24407	if state == 	41	;
replace	st_inc_h_all =	20987	if state == 	42	;
replace	st_inc_h_all =	19330	if state == 	45	;
replace	st_inc_h_all =	.	if state == 	46	;
replace	st_inc_h_all =	19803	if state == 	47	;
replace	st_inc_h_all =	22138	if state == 	48	;
replace	st_inc_h_all =	.	if state == 	50	;
replace	st_inc_h_all =	23325	if state == 	51	;
replace	st_inc_h_all =	31399	if state == 	53	;
replace	st_inc_h_all =	.	if state == 	54	;
replace	st_inc_h_all =	19179	if state == 	55	;
replace	st_inc_h_all =	22544	if state == 	44	;
replace	st_inc_h_all =	.	if state == 	38	;
replace	st_inc_h_all =	.	if state == 	30	;
replace	st_inc_h_all =	.	if state == 	56	;
replace	st_inc_h_all =	20899	if state == 	49	;

gen 	st_inc_all =	.			;
replace	st_inc_all =	24718	if state == 	1	;
replace	st_inc_all =	34883	if state == 	2	;
replace	st_inc_all =	27648	if state == 	4	;
replace	st_inc_all =	23039	if state == 	5	;
replace	st_inc_all =	30490	if state == 	6	;
replace	st_inc_all =	33129	if state == 	8	;
replace	st_inc_all =	36865	if state == 	9	;
replace	st_inc_all =	30693	if state == 	10	;
replace	st_inc_all =	34613	if state == 	11	;
replace	st_inc_all =	27695	if state == 	12	;
replace	st_inc_all =	28670	if state == 	13	;
replace	st_inc_all =	29581	if state == 	15	;
replace	st_inc_all =	25642	if state == 	16	;
replace	st_inc_all =	31311	if state == 	17	;
replace	st_inc_all =	27832	if state == 	18	;
replace	st_inc_all =	27480	if state == 	19	;
replace	st_inc_all =	28857	if state == 	20	;
replace	st_inc_all =	25711	if state == 	21	;
replace	st_inc_all =	24075	if state == 	22	;
replace	st_inc_all =	26949	if state == 	23	;
replace	st_inc_all =	34746	if state == 	24	;
replace	st_inc_all =	34480	if state == 	25	;
replace	st_inc_all =	30780	if state == 	26	;
replace	st_inc_all =	32613	if state == 	27	;
replace	st_inc_all =	22200	if state == 	28	;
replace	st_inc_all =	28261	if state == 	29	;
replace	st_inc_all =	27682	if state == 	31	;
replace	st_inc_all =	30420	if state == 	32	;
replace	st_inc_all =	33219	if state == 	33	;
replace	st_inc_all =	35297	if state == 	34	;
replace	st_inc_all =	23656	if state == 	35	;
replace	st_inc_all =	30481	if state == 	36	;
replace	st_inc_all =	27670	if state == 	37	;
replace	st_inc_all =	29496	if state == 	39	;
replace	st_inc_all =	24980	if state == 	40	;
replace	st_inc_all =	28835	if state == 	41	;
replace	st_inc_all =	29521	if state == 	42	;
replace	st_inc_all =	26203	if state == 	45	;
replace	st_inc_all =	25318	if state == 	46	;
replace	st_inc_all =	25676	if state == 	47	;
replace	st_inc_all =	27747	if state == 	48	;
replace	st_inc_all =	28667	if state == 	50	;
replace	st_inc_all =	32393	if state == 	51	;
replace	st_inc_all =	31770	if state == 	53	;
replace	st_inc_all =	21510	if state == 	54	;
replace	st_inc_all =	30748	if state == 	55	;
replace	st_inc_all =	31467	if state == 	44	;
replace	st_inc_all =	25109	if state == 	38	;
replace	st_inc_all =	23561	if state == 	30	;
replace	st_inc_all =	26800	if state == 	56	;
replace	st_inc_all =	26667	if state == 	49	;

#delimit ;

sort state ; 

save c:\ErikMain\race_signaling\nber_cex_main_tables_figures.dta, replace ;

#delimit ;

clear;

/* segregation measures */; 

#delimit ; 

insheet state dissim_ind_bm  isolation_ind_bm using c:\ErikMain\race_signaling\segregation_index.txt;

sort state;

#delimit ;

save segregation.dta, replace;

clear;  

use c:\ErikMain\race_signaling\nber_cex_main_tables_figures.dta, clear ;

merge state using segregation;

tab _merge;

#delimit ;

drop _merge ;

#delimit ; 

clear;

insheet using c:\ErikMain\race_signaling\cps_income_bystate.txt;

sort state;

save cps_income_bystate, replace  ;

#delimit ; 

clear ; 

use c:\ErikMain\race_signaling\nber_cex_main_tables_figures.dta, clear ;

sort state ;

merge state using cps_income_bystate ;

tab _merge;

#delimit ;

drop if _merge == 2 ; 

drop _merge ; 


#delimit ; 

gen ln_st_inc_wm_1 = ln(lab_inc_main_w);
gen ln_st_inc_bm_1 = ln(lab_inc_main_b);
gen ln_st_inc_hm_1 = ln(lab_inc_main_h);
gen ln_st_inc_m_1 = ln(lab_inc_main);

gen ln_st_inc_meanbyrace_1 = .;

replace ln_st_inc_meanbyrace_1 = ln_st_inc_wm_1 if white == 1;
replace ln_st_inc_meanbyrace_1 = ln_st_inc_bm_1 if black == 1;
replace ln_st_inc_meanbyrace_1 = ln_st_inc_hm_1 if hispanic == 1; 


#delimit ; 

gen ln_st_inc_wm_2 = ln(inc_main_w);
gen ln_st_inc_bm_2 = ln(inc_main_b);
gen ln_st_inc_hm_2 = ln(inc_main_h);
gen ln_st_inc_m_2 = ln(inc_main);

gen ln_st_inc_meanbyrace_2 = .;

replace ln_st_inc_meanbyrace_2 = ln_st_inc_wm_2 if white == 1;
replace ln_st_inc_meanbyrace_2 = ln_st_inc_bm_2 if black == 1;
replace ln_st_inc_meanbyrace_2 = ln_st_inc_hm_2 if hispanic == 1; 


#delimit ; 

gen ln_st_inc_wm_3 = ln(inc_hh_main_w);
gen ln_st_inc_bm_3 = ln(inc_hh_main_b);
gen ln_st_inc_hm_3 = ln(inc_hh_main_h);
gen ln_st_inc_m_3 = ln(inc_hh_main);

gen ln_st_inc_meanbyrace_3 = .;

replace ln_st_inc_meanbyrace_3 = ln_st_inc_wm_3 if white == 1;
replace ln_st_inc_meanbyrace_3 = ln_st_inc_bm_3 if black == 1;
replace ln_st_inc_meanbyrace_3 = ln_st_inc_hm_3 if hispanic == 1; 


#delimit ; 

gen ln_st_inc_wm_4 = ln(lab_inc_alt1_w);
gen ln_st_inc_bm_4 = ln(lab_inc_alt1_b);
gen ln_st_inc_hm_4 = ln(lab_inc_alt1_h);
gen ln_st_inc_m_4 = ln(lab_inc_alt1);

gen ln_st_inc_meanbyrace_4 = .;

replace ln_st_inc_meanbyrace_4 = ln_st_inc_wm_4 if white == 1;
replace ln_st_inc_meanbyrace_4 = ln_st_inc_bm_4 if black == 1;
replace ln_st_inc_meanbyrace_4 = ln_st_inc_hm_4 if hispanic == 1; 


#delimit ; 

gen ln_st_inc_wm_5 = ln(lab_inc_alt2_w);
gen ln_st_inc_bm_5 = ln(lab_inc_alt2_b);
gen ln_st_inc_hm_5 = ln(lab_inc_alt2_h);
gen ln_st_inc_m_5 = ln(lab_inc_alt2);

gen ln_st_inc_meanbyrace_5 = .;

replace ln_st_inc_meanbyrace_5 = ln_st_inc_wm_5 if white == 1;
replace ln_st_inc_meanbyrace_5 = ln_st_inc_bm_5 if black == 1;
replace ln_st_inc_meanbyrace_5 = ln_st_inc_hm_5 if hispanic == 1; 




                /*  --------- Empirical Results for the Paper --------------- */ ; 


/* ------------------ Appendix Table A1 ------------------------------ */ ; 

#delimit ;

sum totexp totvis rent_adj  food_adj utility_adj car_serv_adj transport_other_adj  ent_serv_adj   
health_adj furnish_adj education_adj  ent_dur_adj alc_tob_adj  other_adj
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= ., 
;

#delimit ;

sum totexp totvis rent_adj  food_adj utility_adj car_serv_adj transport_other_adj  ent_serv_adj   
health_adj furnish_adj education_adj  ent_dur_adj alc_tob_adj  other_adj
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & white == 1, 
;

#delimit ;

sum totexp totvis rent_adj  food_adj utility_adj car_serv_adj transport_other_adj  ent_serv_adj   
health_adj furnish_adj education_adj  ent_dur_adj alc_tob_adj  other_adj
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & black == 1, 
;

#delimit ;

sum totexp totvis rent_adj  food_adj utility_adj car_serv_adj transport_other_adj  ent_serv_adj   
health_adj furnish_adj education_adj  ent_dur_adj alc_tob_adj  other_adj
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & hispanic == 1, 
;




#delimit ;

sum totvis_pos  rent_pos food_pos utility_pos car_serv_pos transport_other_pos  ent_serv_pos   
health_pos furnish_pos education_pos  ent_dur_pos alc_tob_pos  other_pos
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= ., 
;

#delimit ;

sum totvis_pos  rent_pos food_pos utility_pos car_serv_pos transport_other_pos  ent_serv_pos   
health_pos furnish_pos education_pos  ent_dur_pos alc_tob_pos  other_pos
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & white == 1, 
;

#delimit ;

sum totvis_pos  rent_pos food_pos utility_pos car_serv_pos transport_other_pos  ent_serv_pos   
health_pos furnish_pos education_pos  ent_dur_pos alc_tob_pos  other_pos
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & black == 1, 
;

#delimit ;

sum totvis_pos  rent_pos food_pos utility_pos car_serv_pos transport_other_pos  ent_serv_pos   
health_pos furnish_pos education_pos  ent_dur_pos alc_tob_pos  other_pos
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & hispanic == 1, 
;


#delimit ;

sum vis_share  rent_share food_share utility_share car_serv_share transport_other_share  ent_serv_share   
health_share furnish_share education_share  ent_dur_share alc_tob_share  other_share
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= ., 
;


#delimit ;

sum vis_share  rent_share food_share utility_share car_serv_share transport_other_share  ent_serv_share   
health_share furnish_share education_share  ent_dur_share alc_tob_share  other_share
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & white == 1, 
;


#delimit ;

sum vis_share  rent_share food_share utility_share car_serv_share transport_other_share  ent_serv_share   
health_share furnish_share education_share  ent_dur_share alc_tob_share  other_share
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & black == 1, 
;


#delimit ;

sum vis_share  rent_share food_share utility_share car_serv_share transport_other_share  ent_serv_share   
health_share furnish_share education_share  ent_dur_share alc_tob_share  other_share
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & hispanic == 1, 
;


                           /* ---- Table 1 (means) ------ */; 

#delimit ;

sum age ed_a ed_b ed_c ed_d married urban famsize adults income_zero totexp 
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= .
;

#delimit ;

sum income
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & income > 0 
;


#delimit ;

sum age ed_a ed_b ed_c ed_d married urban famsize adults income_zero totexp 
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & white == 1 
;

#delimit ;

sum income
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & income > 0 & white == 1 
;


#delimit ;

sum age ed_a ed_b ed_c ed_d married urban famsize adults income_zero totexp 
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & black == 1 
;

#delimit ;

sum income
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & income > 0 & black == 1 
;


#delimit ;

sum age ed_a ed_b ed_c ed_d married urban famsize adults income_zero totexp 
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & hispanic == 1 
;

#delimit ;

sum income
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= . & income > 0 & hispanic == 1 
;





/* ----   Table 2 -  ---- */ ; 

#delimit ;

reg ln_totvis black hispanic  
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= ., cluster(state)
;


#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube  
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)   
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1)  
& year < 2003 & midwest ~= ., cluster(state)
;


#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube 
[aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= ., cluster(state)
;


#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= ., cluster(state)
;


#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= ., cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= ., cluster(state)
;

                                 /*  robustness table A2 */; 

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & income > 0 , cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & totexp > 4000 , cluster(state)
;


#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 24 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . &  fullyr == 1, cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt ln_rent
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;

#delimit ;

xi:reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt i.occup
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;




#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & year >= 1996 , cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt << add in city size controls >>
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & year >= 1996 , cluster(state)
;

#delimit ;

xi:reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt i.popsize
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & year >= 1996 , cluster(state)
;

                              /*   Table 3   */;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & (male == 1 & married == 0), cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & (male == 0 & married == 0), cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & (married == 1), cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & (ed_a == 1), cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & (ed_b == 1), cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & (ed_c == 1), cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & (ed_d == 1), cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 35 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . , cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 35 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . , cluster(state)
;

#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 50 & age < 70 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . , cluster(state)
;

                     /*   Table 4  */ ; 

#delimit ;

reg ln_clothing black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . , cluster(state)
;

#delimit ;

reg ln_pcare black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . , cluster(state)
;


#delimit ;

reg ln_car black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . , cluster(state)
;




#delimit ;

reg ln_totcar black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= ., cluster(state)
;



#delimit ;

reg ln_clothing black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & numauto > 0, cluster(state)
;

#delimit ;

reg ln_pcare black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & numauto > 0, cluster(state)
;


#delimit ;

reg ln_car black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & numauto > 0, cluster(state)
;

#delimit ;

reg ln_totcar black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & numauto > 0 , cluster(state)
;


               /*              Table 5            */ ;


#delimit ;

reg ln_rent black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;


#delimit ;

reg ln_utility black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;

reg ln_food black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;

reg ln_transport_other black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;

#delimit ;

tobit ln_furnish black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , ll(0) 
;

tobit ln_education black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , ll(0) 
;

reg ln_ent_serv black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state) 
;

tobit ln_ent_dur black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , ll(0) 
;

tobit ln_health black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , ll(0) 
;

#delimit ;

tobit ln_alc_tob black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , ll(0) 
;

tobit ln_other black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , ll(0) 
;

         /*  Table 6 */ ; 

#delimit ;

reg ln_totvis ln_st_inc_meanbyrace income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & white == 1, cluster(state)
;  

reg ln_totvis  ln_st_inc_meanbyrace ln_rent income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & white == 1, cluster(state)
;  

reg ln_totvis  ln_st_inc_meanbyrace ln_st_inc_mean_m  ln_rent income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & white == 1, cluster(state)
;  

test ln_st_inc_meanbyrace = ln_st_inc_mean_m ; 

reg ln_totvis  ln_st_inc_meanbyrace ln_rent sharp_byrace income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & white == 1, cluster(state)
;  

         /*  Table 7 */ ; 

#delimit ;

reg ln_totvis ln_st_inc_meanbyrace income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & black == 1, cluster(state)
;  

reg ln_totvis ln_st_inc_meanbyrace income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & hispanic == 1, cluster(state)
;  


reg ln_totvis ln_st_inc_meanbyrace income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & (black == 1 | hispanic == 1), cluster(state)
;  

#delimit ;

reg ln_totvis ln_st_inc_meanbyrace ln_rent income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & (black == 1 | hispanic == 1), cluster(state)
; 


reg ln_totvis ln_st_inc_meanbyrace ln_st_inc_mean_m ln_rent income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & (black == 1 | hispanic == 1), cluster(state)
; 

test ln_st_inc_meanbyrace = ln_st_inc_mean_m ;

reg ln_totvis ln_st_inc_meanbyrace sharp_byrace ln_rent income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  & (black == 1 | hispanic == 1), cluster(state)
; 


          /* Table 8 */;


#delimit ;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state)
;

reg ln_totvis black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
state_2 state_5 state_4 state_6 state_8 state_9 state_10 state_11 state_12 state_13 state_15 state_16 
state_17 state_18 state_19 state_20 state_21 state_22 state_23 state_24 state_25 state_26 state_27 state_28 state_29 state_31
 state_32 state_33 state_34 state_35 state_36 state_37 state_39 state_40 state_41 state_42 state_45 state_46 state_47 state_48 
state_50 state_51 state_53 state_54 state_55 state_44 state_38 state_30 state_56 state_49
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(id)
;

#delimit ;

reg ln_totvis black hispanic ln_st_inc_meanbyrace  income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
state_2 state_5 state_4 state_6 state_8 state_9 state_10 state_11 state_12 state_13 state_15 state_16 
state_17 state_18 state_19 state_20 state_21 state_22 state_23 state_24 state_25 state_26 state_27 state_28 state_29 state_31
 state_32 state_33 state_34 state_35 state_36 state_37 state_39 state_40 state_41 state_42 state_45 state_46 state_47 state_48 
state_50 state_51 state_53 state_54 state_55 state_44 state_38 state_30 state_56 state_49
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(id)
;

reg ln_totvis black hispanic ln_st_inc_meanbyrace sharp_byrace  income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
state_2 state_5 state_4 state_6 state_8 state_9 state_10 state_11 state_12 state_13 state_15 state_16 
state_17 state_18 state_19 state_20 state_21 state_22 state_23 state_24 state_25 state_26 state_27 state_28 state_29 state_31
 state_32 state_33 state_34 state_35 state_36 state_37 state_39 state_40 state_41 state_42 state_45 state_46 state_47 state_48 
state_50 state_51 state_53 state_54 state_55 state_44 state_38 state_30 state_56 state_49
[aw=weight] if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(id)
;


/* table 9 */; 

#delimit ;

tobit ln_education black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , ll(0);

tobit ln_education black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
 rent_share  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , ll(0);

tobit ln_education black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
rent_share vis_share  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , ll(0);




#delimit ;

tobit ln_health black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
   [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , ll(0);

tobit ln_health black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
 rent_share  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , ll(0);

tobit ln_health black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
 rent_share vis_share  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , ll(0);




#delimit ;

reg ln_food black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
 [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state);

reg ln_food black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
 rent_share  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , cluster(state);

reg ln_food black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
 rent_share vis_share  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , cluster(state);



#delimit ;

reg ln_rent black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt 
   [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .  , cluster(state);

reg ln_rent black hispanic income_zero ln_income income income_sq income_cube ln_totexp totexp totexp_sq totexp_cube ed_b ed_c ed_d  
year_87 year_88 year_89 year_90 year_91 year_92 year_93 year_94 year_95 year_96  year_97 year_98 year_99 year_00 year_01 year_02 
age age_sq married male  famsize_2 famsize_3 famsize_4 famsize_5 famsize_6 adults_2 adults_3 adults_4 adults_5  adults_6  
urban midwest south west liquid_pos ln_liquid debt_pos ln_debt  vis_share  [aw=weight]  if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= .   , cluster(state);





* ----- Figure 1 - Engel Curve Estimation ---- */; 

#delimit ;


ivreg ln_totvis ( ln_income   = totexp totexp_sq totexp_cube ed_b ed_c ed_d ln_totexp)  [aw=weight]  
if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & white == 1 & income > 600 & income < 148000 , 
;

predict ln_totvis_hat_w 
if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & white == 1   , 
;

#delimit ;

ivreg ln_totvis (ln_income  = totexp totexp_sq totexp_cube ed_b ed_c ed_d ln_totexp)  [aw=weight]  
if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & black == 1  & income > 600 & income < 104000, 
;

predict ln_totvis_hat_b 
if age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . & black == 1 , 
;



#delimit ;

label variable ln_totvis_hat_b "Black";
label variable ln_totvis_hat_w "White";

#delimit ;

sort income;

twoway (line ln_totvis_hat_b income) (line ln_totvis_hat_w income, lpattern(_-)) if income > 500 & income < 100000 & age >= 18 & age < 50 & (black == 1 | white == 1 | hispanic == 1)  
& totexp < 100000 & ln_st_inc_meanbyrace ~= . & ln_st_inc_mean_bm ~= . & (ed_a == 1 | ed_b == 1 | ed_c == 1 | ed_d == 1) 
& year < 2003 & midwest ~= . , 
ti(Figure 1: Black/White Visible Expenditure Engel Curves)
xti(Household Total Income)
yti(Log Visible Expenditures) name(graph_main_engel, replace) 
;

graph save graph_main_engel, replace;
;