/*
cd "C:\Users\chkim\OneDrive - The University of Kansas\Lectures\Soc910_AdvStat\examples"

use soc910_1_data1.dta, clear

* sample limitation

keep if relate==1 | relate==2 
keep if age>=25 & age<=54
gen married=marst==1
drop if (age_sp<25 | age_sp>54) & married==1

sort serial
by serial: egen nadult = count(serial)
keep if nadult==1 | nadult==2

keep if age>=35 & age<=44

save, replace
*/

cd "C:\Users\chkim\OneDrive - The University of Kansas\Lectures\Soc910_AdvStat\examples"

use soc910_1_data1.dta, clear

* ===========
* Recode
* ===========
numlabel, add

gen eqincome  = ftotinc/sqrt(famsize)
gen lneqincome= ln(eqincome)
gen lnpincome = ln(inctot) if incearn>0 
gen lnspincome= ln(inctot_sp)

gen female=sex==2
gen single=marst==6

gen     races=5
replace races=1 if raced==100 & hispan==0
replace races=2 if raced==200 & hispan==0 
replace races=4 if raced>=400 & raced<=679 & hispan==0
replace races=3 if hispan != 0 

* education
gen     edu=1
replace edu=2 if educd==63 | educd==64
replace edu=3 if educd>=65 & educd<=100
replace edu=4 if educd==101
replace edu=5 if educd==114
replace edu=6 if educd==115
replace edu=7 if educd==116

* spouse education
gen     edu_sp=1
replace edu_sp=2 if educd_sp==63 | educd_sp==64
replace edu_sp=3 if educd_sp>=65 & educd_sp<=100
replace edu_sp=4 if educd_sp==101
replace edu_sp=5 if educd_sp==114
replace edu_sp=6 if educd_sp==115
replace edu_sp=7 if educd_sp==116
replace edu_sp=0 if married==0

* assortative mating
gen     assedu=-1 if edu_sp<edu  & marst==1
replace assedu=0  if edu_sp==edu & marst==1
replace assedu=1  if edu_sp>edu  & marst==1
replace assedu=2  if marst==6

gen marryup=assedu==1
gen marrydown=assedu==-1

gen age2=age*age
gen immig=yrimmig!=0

set more off
tab edu, gen(edu)


* ===========
* regression models 
* ===========	
sum inctot lnpincome, d

reg inctot    female single i.edu i.races age age2, if inctot > 0 
reg lnpincome female single i.edu i.races age age2

reg lneqincome lnpincome single i.edu i.races age age2 if female==0
reg lneqincome lnpincome single i.edu i.races age age2 if female==1

reg lneqincome lnspincome single i.edu i.races age age2 if female==0
reg lneqincome lnspincome single i.edu i.races age age2 if female==1

reg lnpincome  female i.edu i.races age age2 
reg lnpincome  i.edu i.races i.female##c.age i.female##age2 


* ===========
* weighting
* ===========	

reg lnpincome  female i.edu i.races age age2 
reg lnpincome  female i.edu i.races age age2 [aw = perwt]
reg lnpincome  female i.edu i.races age age2 [pw = perwt] 
reg lnpincome  female i.edu i.races age age2 [fw = perwt]