* SanQoL evaluation and CEA paper regressions
* Ian Ross

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

*if running, edit paths to your own folder structure

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


*****************************************************************************
* Table 4: Effects on primary and secondary outcomes 
*****************************************************************************

*************** unadjusted  **************

eststo clear 

eststo: meglm sqwt i.treatment || c2a:, base vce(cluster c2a)
eststo: meglm v1 i.treatment || c2a:, base vce(cluster c2a)
eststo: meglm who5index i.treatment || c2a:, base vce(cluster c2a)
		
esttab using tab4_headline_unadj.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p) nogap label
esttab using tab4_headline_unadj_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label

*************** adjusted  **************

eststo clear 

eststo: meglm sqwt i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm v1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm who5index i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)

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

*****************************************************************************
* Table 5: Interactions by sex and age-group
*****************************************************************************

eststo clear 

eststo: meglm sqwt i.over60 i.treatment##i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm sqwt i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm v1 i.over60 i.treatment##i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm v1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)

eststo: meglm who5index i.over60 i.treatment##i.female wealth || c2a:, base vce(cluster c2a)

eststo: meglm who5index i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)
			
			
esttab using tab5_interactions.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p) nogap label
esttab using tab5_interactions_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label


*****************************************************************************
* Table 6: Differences in SanQoL attributes and interactions w gender and age
*****************************************************************************

* unadjusted

eststo clear 

eststo: meglm d1 i.treatment || c2a:, base vce(cluster c2a)
eststo: meglm h1 i.treatment || c2a:, base vce(cluster c2a)
eststo: meglm s1 i.treatment || c2a:, base vce(cluster c2a)
eststo: meglm t1 i.treatment || c2a:, base vce(cluster c2a)
eststo: meglm p1 i.treatment || c2a:, base vce(cluster c2a)

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

* adjusted

eststo clear 

eststo: meglm d1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm h1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm s1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm t1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm p1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)

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

* gender interaction

eststo clear 

eststo: meglm d1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a)
eststo: meglm h1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a)
eststo: meglm s1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a)
eststo: meglm t1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a)
eststo: meglm p1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a)

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


* over60 interaction

eststo clear 

eststo: meglm d1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm h1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm s1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm t1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)
eststo: meglm p1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a)

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

*******************************************************************************
* Online Appendix Table D-2: table 3 regression incl toilet type as covariate
*******************************************************************************

eststo clear

eststo: meglm sqwt i.treatment i.wsuptype i.over60 i.female wealth || c2a:, base vce(cluster c2a)

lincom 1.treatment + 2.wsuptype
		
esttab using tabd2_STvsCSB.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p) nogap label
esttab using tabd2_STvsCSB_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label

******************************************************************************
* Online Appendix Table E 1: Robustness checks for SanQoL and VAS
******************************************************************************

*** CHECKS ON SANQOL    ****
  
eststo clear 

eststo: meglm sqwt i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
xtset c2a
eststo: xtgee sqwt i.treatment i.over60 i.female wealth, base
eststo: reg sqwt i.treatment i.over60 i.female wealth, vce(cluster c2a) baselevels
eststo: meglm sqwt i.treatment respage_cat wealth i.secondaryorabove  || c2a:, base vce(cluster c2a)	
eststo: meglm sqwt i.treatment i.over60 i.female peoplesharingcub i.sharetoiletyesno i.renter || c2a:, base vce(cluster c2a)	
	
*** CHECKS ON VAS	    ****

eststo: meglm v1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
xtset c2a
eststo: xtgee v1 i.treatment i.over60 i.female wealth, base
eststo: reg v1 i.treatment i.over60 i.female wealth, vce(cluster c2a) baselevels
eststo: meglm v1 i.treatment respage_cat wealth i.secondaryorabove || c2a:, base vce(cluster c2a)	
eststo: meglm v1 i.treatment i.over60 i.female peoplesharingcub i.sharetoiletyesno i.renter || c2a:, base vce(cluster c2a)	
	
esttab using tabe1_checks_sanqol_vas.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p) nogap label
esttab using tabe1_checks_sanqol_vas_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label
	
*****************************************************************************
* Online Appendix Table E 2: Robustness checks for WHO-5
*****************************************************************************

eststo clear 

eststo: meglm who5index i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a)
xtset c2a
eststo: xtgee who5index i.treatment i.over60 i.female wealth, base
eststo: reg who5index i.treatment i.over60 i.female wealth, vce(cluster c2a) baselevels
eststo: meglm who5index i.treatment respage_cat wealth i.secondaryorabove || c2a:, base vce(cluster c2a)	
eststo: meglm who5index i.treatment i.over60 i.female partner eq_painscale eq_mobility || c2a:, base vce(cluster c2a)	
		
esttab using tabe2_checks_who5.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(3) b(3) cells(b se p) nogap label
esttab using tabe2_checks_who5_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label

************************************************************************************************
* Online Appendix Table E 3: Robustness checks for attribute level models (coefficients are odds ratios)
************************************************************************************************

eststo clear

* adjusted

eststo: meologit d1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit h1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit s1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit t1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit p1 i.treatment i.over60 i.female wealth || c2a:, base vce(cluster c2a) or

esttab using tabe3_attlevel_meologit.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) cells(b se p) nogap label eform
esttab using tabe3_attlevel_meologit_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label eform


************************************************************************************************
* Online Appendix Table E 4: Robustness checks for attribute level interactions (coefficients are odds ratios)
************************************************************************************************

* gender interaction

eststo clear 

eststo: meologit d1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a) or
eststo: meologit h1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a) or
eststo: meologit s1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a) or
eststo: meologit t1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a) or
eststo: meologit p1 i.treatment##i.female i.over60 wealth || c2a:, base vce(cluster c2a) or

* over60 interaction

eststo: meologit d1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit h1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit s1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit t1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a) or
eststo: meologit p1 i.treatment##i.over60 i.female wealth || c2a:, base vce(cluster c2a) or


esttab using tabe3_attlevel_inters_meologit.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) cells(b se p) nogap label eform
esttab using tabe3_attlevel_inters_meologit_formatted.csv, star(* 0.10 ** 0.05 *** 0.01) replace se(2) b(2) nogap label eform
 
************************************************************************************************
* Online Appendix Table E 5: Interactions with sharing toilets
************************************************************************************************

eststo clear

*sanqol
eststo: meglm sqwt i.treatment##i.sharetoiletyesno i.over60 i.female wealth  || c2a:, base vce(cluster c2a)

lincom 1.treatment + 1.sharetoiletyesno#1.treatment
lincom 1.treatment + 0.sharetoiletyesno#1.treatment

*vas
eststo: meglm v1 i.treatment##i.sharetoiletyesno  i.over60 i.female wealth || c2a:, base vce(cluster c2a)

lincom 1.treatment + 1.sharetoiletyesno#1.treatment
lincom 1.treatment + 0.sharetoiletyesno#1.treatment
		
*who-5
eststo: meglm who5index i.treatment##i.sharetoiletyesno i.over60 i.female wealth || c2a:, base vce(cluster c2a)

lincom 1.treatment + 1.sharetoiletyesno#1.treatment
lincom 1.treatment + 0.sharetoiletyesno#1.treatment

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


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

clear