*Musterlösung Hausaufgabe 6
	

clear
set mem 100m

use q31 q32 hhinc4 hh2b hh2a gdppcap1 hdivalue homicide s_cntry s_respnr ///
  using data/eqls_4.dta

tab q32, gen(mar)

gen men= hh2a==1 if !mi(hh2a)
rename hh2b alter

iis s_cntry

// Note the difference in s_u!
xtreg q31 alter men mar2-mar4, re
eststo
xtreg q31 alter men mar2-mar4 gdppcap1, re
eststo

esttab, stats(sigma_u sigma_e rho r2_o N )
	
xtreg q31 alter men mar2-mar4 hhinc4, re
eststo
xtreg q31 alter men mar2-mar4 hhinc4 gdppcap1, re
eststo

esttab, stats(sigma_u sigma_e rho r2_o N )

xtreg q31 alter men mar2-mar4 hhinc4 hdivalue1, re
xtreg q31 alter men mar2-mar4 hhinc4 gdppcap1 hdivalue , re

xtreg q31 alter men mar2-mar4 hhinc4 homicide, re
xtreg q31 alter men mar2-mar4 hhinc4 gdppcap1 hdivalue homicide, re


// Test 1
xtreg q31 alter men mar2-mar4, re
predict u1, u
xtreg q31 alter men mar2-mar4 gdppcap1, re
predict u2, u

tw ///
  || kdensity u1 || kdensity u2

// Test 2
xtreg q31 alter men mar2-mar4 hhinc4, re
predict u3, u
xtreg q31 alter men mar2-mar4 hhinc4 gdppcap1, re
predict u4, u

tw ///
  || kdensity u3 || kdensity u4


	exit