** three-level logistic regression lab       **
** use the guatemala data set on the website **

*EDA*
codebook cluster
sort cluster
tab cluster

summ immun

* analyses *
gllamm immun kid2p indNoSpa indSpa  momEdPri momEdSec husEdPri husEdSec husEdDK rural pcInd81, family(binomial) link(logit) i(mom cluster) nip(5)

matrix a=e(b) 
gllamm immun kid2p indNoSpa indSpa  momEdPri momEdSec husEdPri husEdSec husEdDK rural pcInd81, family(binomial) link(logit) i(mom cluster) from(a) adapt
gllamm, eform

estimates store model1 

*with fewer covariates 
matrix a=e(b)
gllamm immun kid2p rural pcInd81, family(binomial) link(logit) i(mom cluster) from(a) skip adapt
gllamm, eform
   
estimates store model2

*random intercept and random slope 
gen cons=1
eq inter: cons
eq slope: kid2p
matrix a = e(b)
matrix a = (a,.2,0)
gllamm immun kid2p rural pcInd81, family(binomial) link(logit) i(mom cluster) nrf(1 2) eqs(inter inter slope) nip(8 4 4) from(a) copy adapt eform
estimates store model3


* cross-level interaction
* model 1
gen cons=1
eq inter: cons
eq slope: kid2p
gllamm immun kid2p, family(binomial) link(logit) i(cluster) nrf(2) eqs(inter slope) nip(4 4) adapt eform


* model 2
eq inter: cons
eq slope: kid2p
gllamm immun kid2p rural pcInd81, family(binomial) link(logit) i(cluster) nrf(2) eqs(inter slope) nip(4 4) adapt eform

* model 3
gen int_2p_ru = kid2p * rural
gen int_2p_pc = kid2p * pcInd81

eq inter: cons
eq slope: kid2p
gllamm immun kid2p rural pcInd81 int_2p_ru int_2p_pc, family(binomial) link(logit) i(cluster) nrf(2) eqs(inter slope) nip(4 4) adapt eform