use "C:\Users\smendoli\Dropbox\Torino new\Panel data course\Understanding society teaching dataset\stata\stata13\longitudinal_td.dta"
ta wave
lookfor qualif
lookfor satisf
ta sclfsato, m
xtset pidp wave
xttab  sclfsato
tab jbstat, gen(job)


g outlforce=0
replace outlforce=1 if jbstat>=4
replace outlforce=. if jbstat==.

g employed=0
replace employed=1 if jbstat==2
replace employed=. if jbstat==.

g self_employed=0
replace self_employed=1 if jbstat==1
replace self_employed=. if jbstat==.

g unemployed=0
replace unemployed=1 if jbstat==3
replace unemployed=. if jbstat==.

ta mstat_dv, gen (marit)

g married=0
replace married=1 if marit2==1| marit5==1
replace married=. if mstat==.
ren marit3 separated
ren marit4 widow

xtset pidp wave

reg sclfsato age_dv i.sex_dv married separated widow i.hiqual_dv self_employed unemployed outlforce
est store OLS

xtreg sclfsato age_dv i.sex_dv married separated widow i.hiqual_dv self_employed unemployed outlforce, fe
est store FE

coefplot OLS FE, keep (married separated widow self_employed unemployed outlforce) xline (0)


xtreg sclfsato age_dv i.sex_dv married separated widow i.hiqual_dv self_employed unemployed outlforce, re
est store RE

estimates table FE RE, b se 

hausman FE RE, sigmamore