*SanQoL dev't paper regressions
* Ian Ross

********************************************
*******  		setup directories    *******
********************************************

clear all
cd "/Users/Ian/zDATA/SanQoL/input"
use dta_sanqol_moz_2019
cd "/Users/Ian/zDATA/SanQoL/output/devt"


******************************************************************************
*Table 6: Difference in ranks by user characteristics (mixed ordered logit)
******************************************************************************

**** main analysis - meologit

eststo clear 

eststo: meologit r_disg i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a) eform
eststo: meologit r_heal i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a) eform
eststo: meologit r_sham i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a) eform
eststo: meologit r_safe i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a) eform 
eststo: meologit r_priv i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a) eform

esttab using reg_ranks_meologit.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p ci) nogap label eform plain
esttab using reg_ranks_meologit_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label eform 

**** robustness check - GLMM

eststo clear 

eststo: meglm r_disg i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a)
eststo: meglm r_heal i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a)
eststo: meglm r_sham i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a)
eststo: meglm r_safe i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a)
eststo: meglm r_priv i.female i.over60 i.treatment wealth i.completed_primary || c2a:, base vce(cluster c2a)

esttab using reg_ranks_glmm.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p) nogap label
esttab using reg_ranks_glmm_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label


****************************************************************************
*Table 8: Associations b/w SanQoL attributes and user/toilet characteristics
****************************************************************************

**** main analysis - GLMM

eststo clear 

eststo: meglm d1 i.female i.over60 wealth i.floormat i.wallmat_bin i.lock_inside i.obs_smellfa i.completed_primary  || c2a:, base vce(cluster c2a)
eststo: meglm h1 i.female i.over60 wealth i.floormat i.wallmat_bin i.lock_inside i.obs_smellfa i.completed_primary || c2a:, base vce(cluster c2a)
eststo: meglm s1 i.female i.over60 wealth i.floormat i.wallmat_bin i.lock_inside i.obs_smellfa i.completed_primary || c2a:, base vce(cluster c2a)
eststo: meglm t1 i.female i.over60 wealth i.floormat i.wallmat_bin i.lock_inside i.obs_smellfa i.completed_primary || c2a:, base vce(cluster c2a)
eststo: meglm p1 i.female i.over60 wealth i.floormat i.wallmat_bin i.lock_inside i.obs_smellfa i.completed_primary || c2a:, base vce(cluster c2a)

esttab using reg_validity_glmm.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p) nogap label
esttab using reg_validity_glmm_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label

******************************************************************************************


clear