cd "C:\Users\lorenzo\Desktop\Development Economics 1112\Nigeria Data"
set more off
use sect1.dta, clear

tab s1q12
tab s1q12 zone
tab s1q12 zone, mis

histogram s1q12, discrete percent addlabel xlabel(, labels valuelabel) legend(on)
histogram s1q12, discrete percent addlabel xlabel(, labels valuelabel) by(, legend(on)) by(zone)

tab s1q2, nol
gen female = s1q2
recode female 1=0 2=1
table s1q12 zone, c(m female) format(%3.2f)

gen religion = s1q12 if s1q12<=3
tab religion

label var religion "individual religion"
label define indrel 1 "chris" 2 "islam" 3 "trad"
label values religion indrel

gen age = s1q4
replace age = . if age>108
sort hhid
by hhid : egen Mage = mean(age) 

sort hhid s1q3
by hhid: gen temp = age if s1q3==2
by hhid: egen spouseage = max(temp)
scatter spouseage age if s1q3==1
twoway (scatter spouseage age if s1q3==1) (lfit age age if s1q3==1) 

sum age spouseage if s1q3==1
sum age spouseage if s1q3==1, det


sort hhid indiv
save temp0.dta, replace


use sect11a1.dta, clear
sort hhid plotid
by hhid: egen totland = total(s11aq4c)
sum totland, det
kdensity totland
drop if totland==.
drop if totland==0
drop if totland>120000
gen share = s11aq4c/totland

collapse totland, by(hhid)
sort hhid
save temp1.dta, replace

use temp0.dta, clear
merge hhid using temp1.dta
tab _merge
drop _merge

erase temp0.dta
erase temp1.dta
sort hhid indiv
save temp2.dta, replace


use sect2.dta, clear
sort hhid indiv
save sect2.dta, replace

use temp2, clear
merge hhid indiv using sect2.dta
tab _merge
drop _merge
erase temp2.dta

gen schooling = s2q4
recode schooling 1=1 2=0
save mydata.dta, replace
 
preserve
collapse schooling, by(hhid)
histogram schooling, percent
save temp.dta, replace
restore

sort hhid female
by hhid : egen malesch = mean(schooling) if female==0
by hhid : egen femalesch = mean(schooling) if female==1
by hhid female : gen Nmale = _N if female==0
by hhid female : gen Nfemale = _N if female==1 

collapse totland malesch femalesch Nmale Nfemale, by(zone hhid)  

tab zone, gen(z)

reg totland Nmale Nfemale malesch femalesch z2-z6, robust
reg totland Nmale Nfemale malesch femalesch, robust