/*****************************************************************************
 File Name:             figures1-3.do         
 Project:				SPOON Guatemala
 Created by:            Solis Winters
 Created on:			12/07/2022

* Description: This do file replicates Figures S1-S3 in the Appendix
*****************************************************************************/

*---------------------------------------------------*
*		 	FIGURE S1. CHILD NUTRITION				*
*---------------------------------------------------*


use "${processed}/child_target_analysis.dta", clear
drop if comunidad_elegible==0

	gen t1_t2 =tratamiento_hogar if tratamiento_hogar!=3
	gen t1_control=tratamiento_hogar if tratamiento_hogar!=2
	gen t2_control=tratamiento_hogar if tratamiento_hogar!=1

	* HAZ
	ksmirnov _zlen, by(t1_control)
	ksmirnov _zlen, by(t2_control)
	ksmirnov _zlen, by(t1_t2)

		twoway (kdensity _zlen if tratamiento_hogar==1) (kdensity _zlen if tratamiento_hogar==2) (kdensity _zlen if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Height-for-age z-score (HAZ)") ytitle("Density") text(0.35 1 "{bf:K-S test p-values}" "T1 v. C: 0.059" "T2 v. C: 0.670" "T1 v. T2: 0.809", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/_zlen.gph", replace) name(_zlen, replace)

	* WAZ
	
	ksmirnov _zwei, by(t1_control)
	ksmirnov _zwei, by(t2_control)
	ksmirnov _zwei, by(t1_t2)
	
	twoway (kdensity _zwei if tratamiento_hogar==1) (kdensity _zwei if tratamiento_hogar==2) (kdensity _zwei if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Weight-for-age z-score (WAZ)") ytitle("Density") text(0.45 3 "{bf:K-S test p-values}" "T1 v. C: 0.194" "T2 v. C: 0.421" "T1 v. T2: 0.639", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/_zwei.gph", replace) name(_zwei, replace)
	
	* BAZ
	
	ksmirnov _zbmi, by(t1_control)
	ksmirnov _zbmi, by(t2_control)
	ksmirnov _zbmi, by(t1_t2)
	
	twoway (kdensity _zbmi if tratamiento_hogar==1) (kdensity _zbmi if tratamiento_hogar==2) (kdensity _zbmi if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("BMI-for-age z-score (BAZ)") ytitle("Density") text(0.45 5 "{bf:K-S test p-values}" "T1 v. C: 0.598" "T2 v. C: 0.636" "T1 v. T2: 0.240", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/_zbmi.gph", replace) name(_zbmi, replace)
	
	* Hemoglobin
	
	ksmirnov adj_haemoglobin, by(t1_control)
	ksmirnov adj_haemoglobin, by(t2_control)
	ksmirnov adj_haemoglobin, by(t1_t2)
	
	twoway (kdensity adj_haemoglobin if tratamiento_hogar==1) (kdensity adj_haemoglobin if tratamiento_hogar==2) (kdensity adj_haemoglobin if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Adjusted hemoglobin (g/dL)") ytitle("Density") text(0.35 17 "{bf:K-S test p-values}" "T1 v. C: 0.183" "T2 v. C: 0.348" "T1 v. T2: 0.755", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/adj_haemoglobin.gph", replace) name(adj_haemoglobin, replace)
	
	* Combined graph
	cap net install grc1leg,from( http://www.stata.com/users/vwiggins/)
	grc1leg _zlen _zwei _zbmi adj_haemoglobin, rows(2) col(2) scheme(538w)
	graph export "$tables_sup/fig_S1.nutrition.jpg", as(jpg) quality(90) replace

*---------------------------------------------------*
*		 	FIGURE S2. IYCF PRACTICES				*
*---------------------------------------------------*

	* Practice index
	ksmirnov practice_index, by(t1_control)
	ksmirnov practice_index, by(t2_control)
	ksmirnov practice_index, by(t1_t2)
	
	twoway (kdensity practice_index if tratamiento_hogar==1) (kdensity practice_index if tratamiento_hogar==2) (kdensity practice_index if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Practice index") ytitle("Density") text(3.5 0.9 "{bf:K-S test p-values}" "T1 v. C: 0.002" "T2 v. C: 0.489" "T1 v. T2: 0.267", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/practices_index_sub.gph", replace) name(practices_index_sub, replace)

	* Practices subscale: breastfeeding
	
	ksmirnov practice_index_bf, by(t1_control)
	ksmirnov practice_index_bf, by(t2_control)
	ksmirnov practice_index_bf, by(t1_t2)
	
	twoway (kdensity practice_index_bf if tratamiento_hogar==1) (kdensity practice_index_bf if tratamiento_hogar==2) (kdensity practice_index_bf if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Practices subscale: breastfeeding") ytitle("Density") text(3.5 0.2 "{bf:K-S test p-values}" "T1 v. C: 0.000" "T2 v. C: 0.282" "T1 v. T2: 0.163", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/practice_index_bf.gph", replace) name(practice_index_bf, replace)

	* Practices subscale: solid foods
	
	ksmirnov practice_index_food, by(t1_control)
	ksmirnov practice_index_food, by(t2_control)
	ksmirnov practice_index_food, by(t1_t2)
	
	twoway (kdensity practice_index_food if tratamiento_hogar==1) (kdensity practice_index_food if tratamiento_hogar==2) (kdensity practice_index_food if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Practices subscale: solid foods") ytitle("Density") text(1.75 0.2 "{bf:K-S test p-values}" "T1 v. C: 0.780" "T2 v. C: 0.845" "T1 v. T2: 1.000", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/practice_index_food.gph", replace) name(practice_index_food, replace)

	* Practices subscale: supplementation
	ksmirnov practice_index_sup, by(t1_control)
	ksmirnov practice_index_sup, by(t2_control)
	ksmirnov practice_index_sup, by(t1_t2)
	
	twoway (kdensity practice_index_sup if tratamiento_hogar==1) (kdensity practice_index_sup if tratamiento_hogar==2) (kdensity practice_index_sup if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Practices subscale: supplementation") ytitle("Density") text(2.5 0.9 "{bf:K-S test p-values}" "T1 v. C: 0.798" "T2 v. C: 1.000" "T1 v. T2: 0.557", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/practice_index_sup.gph", replace) name(practice_index_sup, replace)

	* Combined graph
	grc1leg practices_index_sub practice_index_bf practice_index_food practice_index_sup, rows(2) col(2) scheme(538w)
	graph export "$tables_sup/fig_S2.practices.jpg", as(jpg) quality(90) replace
		

*---------------------------------------------------*
*			FIGURE S3. CAREGIVER KNOWLEDGE			*
*---------------------------------------------------*
use "${processed}/caregiver_target_analysis.dta", clear
drop if comunidad_elegible==0

	gen t1_t2 =tratamiento_hogar if tratamiento_hogar!=3
	gen t1_control=tratamiento_hogar if tratamiento_hogar!=2
	gen t2_control=tratamiento_hogar if tratamiento_hogar!=1

	* Knowledge index
	ksmirnov knowledge_indexV1, by(t1_control)
	ksmirnov knowledge_indexV1, by(t2_control)
	ksmirnov knowledge_indexV1, by(t1_t2)
	
	twoway (kdensity knowledge_indexV1 if tratamiento_hogar==1) (kdensity knowledge_indexV1 if tratamiento_hogar==2) (kdensity knowledge_indexV1 if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Knowledge index") ytitle("Density") text(4.5 0.7 "{bf:K-S test p-values}" "T1 v. C: 0.000" "T2 v. C: 0.000" "T1 v. T2: 0.220", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/knowledge_indexV1.gph", replace) name(knowledge_indexV1, replace)

	* Knowledge subscale: breastfeeding 
	ksmirnov knowledge_indexV1_bf, by(t1_control)
	ksmirnov knowledge_indexV1_bf, by(t2_control)
	ksmirnov knowledge_indexV1_bf, by(t1_t2)
	
	twoway (kdensity knowledge_indexV1_bf if tratamiento_hogar==1) (kdensity knowledge_indexV1_bf if tratamiento_hogar==2) (kdensity knowledge_indexV1_bf if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Knowledge subscale: breastfeeding") ytitle("Density") text(3.5 0.9 "{bf:K-S test p-values}" "T1 v. C: 0.000" "T2 v. C: 0.000" "T1 v. T2:  0.635", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/knowledge_indexV1_bf.gph", replace) name(knowledge_indexV1_bf, replace)
	

	* Knowledge subscale: solid foods
	ksmirnov knowledge_indexV1_food, by(t1_control)
	ksmirnov knowledge_indexV1_food, by(t2_control)
	ksmirnov knowledge_indexV1_food, by(t1_t2)
	
	twoway (kdensity knowledge_indexV1_food if tratamiento_hogar==1) (kdensity knowledge_indexV1_food if tratamiento_hogar==2) (kdensity knowledge_indexV1_food if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Knowledge subscale: solid foods") ytitle("Density") text(3.5 0.7 "{bf:K-S test p-values}" "T1 v. C: 0.000" "T2 v. C: 0.000" "T1 v. T2:  0.044", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/knowledge_indexV1_food.gph", replace) name(knowledge_indexV1_food, replace)

	* Knowledge subscale: health/health care
	ksmirnov knowledge_indexV1_health, by(t1_control)
	ksmirnov knowledge_indexV1_health, by(t2_control)
	ksmirnov knowledge_indexV1_health, by(t1_t2)
	
	twoway (kdensity knowledge_indexV1_health if tratamiento_hogar==1) (kdensity knowledge_indexV1_health if tratamiento_hogar==2) (kdensity knowledge_indexV1_health if tratamiento_hogar==3), ///
	note(" ") legend(label( 1 "T1") label(2 "T2") label(3 "Control") rows(1)) xtitle("Knowledge subscale: health/health care") ytitle("Density") text(3.5 0.7 "{bf:K-S test p-values}" "T1 v. C: 0.000" "T2 v. C: 0.000" "T1 v. T2:  0.689", placement(c) justification(left) size(vsmall)) scheme(538w) saving("$tables_sup/knowledge_indexV1_health.gph", replace) name(knowledge_indexV1_health, replace)
	
	* Combined graph
	grc1leg knowledge_indexV1 knowledge_indexV1_bf knowledge_indexV1_food knowledge_indexV1_health, rows(2) col(2) scheme(538w)
	graph export "$tables_sup/fig_S3.knowledge.jpg", as(jpg) quality(90) replace