/******************************************************************************

	PROJECT: 		ERC SES
	AUTHOR: 		NICOLA FOSTER
	LAST UPDATED: 	31 JANUARY 2022

******************************************************************************/


*******************************************************************************
* INSTALLING SCHEMES
ssc install blindschemes, replace
ssc install scheme-burd, replace
ssc install scheme-mrc, replace
ssc install scheme-tfl, replace
ssc install scheme_scientific
ssc install scheme_tufte, replace
ssc install scheme-burd, replace


* SET SCHEME
* set scheme plotplain, permanently
* set scheme plotplainblind, permanently
* set scheme plottig, permanently
* set scheme s2gcolor, permanently
set scheme burd, permanently
graph drop _all


*******************************************************************************
* CONCENTRATION INDEX
* the compare(survey) 
conindex tbstat [aweight= Wtotps], rankvar(pca_quartile) truezero robust compare(survey)


*******************************************************************************
* CONCENTRATION CURVES
* lorenz curve - tb cases
lorenz estimate tbcase if survey==0, pvar(pca_quartile)
lorenz graph, aspectratio(1) xlabels(, grid) overlay legend(cols(1)) name(lo1)
lorenz estimate tbcase if survey==1, pvar(pca_quartile)
lorenz graph, aspectratio(1) xlabels(, grid) overlay legend(cols(1)) name(lo2)
graph combine lo1 lo2	

* lorenz curve - tb symptoms
lorenz estimate symptoms_any if survey==0, pvar(pca_quartile)
lorenz graph, aspectratio(1) xlabels(, grid) overlay legend(cols(1)) name(lo3)
lorenz estimate symptoms_any if survey==1, pvar(pca_quartile)
lorenz graph, aspectratio(1) xlabels(, grid) overlay legend(cols(1)) name(lo4)
graph combine lo3 lo4


*******************************************************************************	
* Test of significant differences
conindex tbstat, rankvar(pca_quartile) truezero robust compare(survey) graph box