/*****************************************************************************
This do file replicates Tables S4-S10 in the Supplementary Material of the 
impact evaluation paper of SPOON Guatemala

Updated on 6/17/2024

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

*---------------------------------------------------------------------------*
*	TABLE S4. Marginal effects on overweight/obesity, logistic regression 	*
*---------------------------------------------------------------------------*

use "${processed}/child_target_analysis.dta", clear

	* Generate and label interaction variables 
	xtile age_q2=cg_age, n(2) // median cutoff
	recode age_q2 (1=0) (2=1)
	xtile decisiones_q2=decisiones_bl, n(2) // median cutoff
	recode decisiones_q2 (1=0) (2=1)
	gen grit_q2 = 0 if grit_bl<=3 & grit_bl!=. // cutoff at 3
	replace grit_q2 = 1 if grit_bl>3 & grit_bl!=. 
	gen rosenberg_q2 = 0 if rosenberg_bl<15 & rosenberg_bl!=. // cutoff at 15
	replace rosenberg_q2 = 1 if rosenberg_bl>=15 & rosenberg_bl!=.
	xtile wealth_q2=wealth_bl, n(2) // median cutoff
	recode wealth_q2 (1=0) (2=1)
	
	* Generate age-sex dummies
	tab child_age_cat child_sex
	tab child_age_cat, gen(child_age)
	gen child_age1_m = (child_age1==1 & child_sex==1)
	gen child_age1_f = (child_age1==1 & child_sex==0)
	gen child_age2_m = (child_age2==1 & child_sex==1)
	gen child_age2_f = (child_age2==1 & child_sex==0)
	gen child_age3_m = (child_age3==1 & child_sex==1)
	gen child_age3_f = (child_age3==1 & child_sex==0)
	gen child_age4_m = (child_age4==1 & child_sex==1)
	gen child_age4_f = (child_age4==1 & child_sex==0)
	gen child_age5_m = (child_age5==1 & child_sex==1)
	gen child_age5_f = (child_age5==1 & child_sex==0)
	gen child_age6_m = (child_age6==1 & child_sex==1)
	gen child_age6_f = (child_age6==1 & child_sex==0)
	gen child_age7_m = (child_age7==1 & child_sex==1)
	gen child_age7_f = (child_age7==1 & child_sex==0)
	gen child_age8_m = (child_age8==1 & child_sex==1)
	gen child_age8_f = (child_age8==1 & child_sex==0)
	gen child_age9_m = (child_age9==1 & child_sex==1)
	gen child_age9_f = (child_age9==1 & child_sex==0)
	gen child_age10_m = (child_age10==1 & child_sex==1)
	gen child_age10_f = (child_age10==1 & child_sex==0)
	gen child_age11_m = (child_age11==1 & child_sex==1)
	gen child_age11_f = (child_age11==1 & child_sex==0)
	gen child_age12_m = (child_age12==1 & child_sex==1)
	gen child_age12_f = (child_age12==1 & child_sex==0)

	lab var child_age1_m "Child age: <30 months, male"
	lab var child_age2_m "Child age: 30-31 months, male"
	lab var child_age3_m "Child age: 32 months, male"
	lab var child_age4_m "Child age: 33 months, male"
	lab var child_age5_m "Child age: 34 months, male"
	lab var child_age6_m "Child age: 35 months, male"
	lab var child_age7_m "Child age: 36 months, male"
	lab var child_age8_m "Child age: 37 months, male"
	lab var child_age9_m "Child age: 38 months, male"
	lab var child_age10_m "Child age: 39 months, male"
	lab var child_age11_m "Child age: 40-41 months, male"
	lab var child_age12_m "Child age: >41 months, male"

	lab var child_age1_f "Child age: <30 months, female"
	lab var child_age2_f "Child age: 30-31 months, female"
	lab var child_age3_f "Child age: 32 months, female"
	lab var child_age4_f "Child age: 33 months, female"
	lab var child_age5_f "Child age: 34 months, female"
	lab var child_age6_f "Child age: 35 months, female"
	lab var child_age7_f "Child age: 36 months, female"
	lab var child_age8_f "Child age: 37 months, female"
	lab var child_age9_f "Child age: 38 months, female"
	lab var child_age10_f "Child age: 39 months, female"
	lab var child_age11_f "Child age: 40-41 months, female"
	lab var child_age12_f "Child age: >41 months, female"

	lab var wealth_bl "Wealth"

** model 1) base model with no controls (unadjusted difference in means)
	logit bmi_2sd T1 T2, cluster(communityid)
	eststo margin: margins, dydx(*) post
	test T1=T2
	outreg2 [margin] using "$tables_sup/tableS4.xls", label dec(3) replace excel addstat("Significance test T1=T2:",r(p)) 

** model 2) + endline child age-sex dummies, leaving out youngest group for each sex, and child order
	* combine omitted age categories for logit: 
	gen child_age2_3_f = child_age2_f==1 | child_age3_f==1 
	gen child_age5_6_f = child_age5_f==1 | child_age6_f==1 
	gen child_age9_10_m = child_age9_m==1 | child_age10_m==1 
	gen child_age11_12_m = child_age11_m==1 | child_age12_m==1 

	logit bmi_2sd T1 T2 child_age2_m child_age2_3_f child_age3_m child_age4_m child_age4_f child_age5_m child_age5_6_f child_age6_m child_age7_m child_age7_f child_age8_f child_age9_10_m child_age9_f child_age10_f child_age11_12_m child_age11_f child_age12_f child_order, cluster(communityid) 
	eststo margin: margins, dydx(*) post // report marginal effects
	test T1=T2
	outreg2 [margin] using "$tables_sup/tableS4.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 


** model 3) + primary caregiver's education in years 
*			+ primary caregiver's age in years
*			+ primary caregiver's baseline decision-making power 
*			+ primary caregiver's baseline  grit score 
*			+ primary caregiver's baseline Rosenberg (self-esteem) score
*			+ baseline household wealth 
*			+ baseline household demographic composition (proportion of household members by age/sex categories)
		
	logit bmi_2sd T1 T2 child_age2_m child_age2_3_f child_age3_m child_age4_m child_age4_f child_age5_m child_age5_6_f child_age6_m child_age7_m child_age7_f child_age8_f child_age9_10_m child_age9_f child_age10_f child_age11_12_m child_age11_f child_age12_f child_order cg_age cg_educ decisiones_bl decisiones_bl_imputed grit_bl grit_bl_imputed rosenberg_bl rosenberg_bl_imputed wealth_bl hh_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 
	eststo margin: margins, dydx(*) post
	test T1=T2
	outreg2 [margin] using "$tables_sup/tableS4.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 	

	
*-----------------------------------------------------------------------*
*	TABLES S5-S10. Hetereogeneity analysis								*
*-----------------------------------------------------------------------*

	* Run regressions and output results to Excel files for each outcome-interaction combination 
	global outcomes _zlen hfa_2sd hfa_3sd _zbmi bmi_2sd

	foreach var of global outcomes {

	** model 1) base model with no controls (unadjusted difference in means)
	reg `var' T1 T2 1.T1#1.cg_educ 1.T2#1.cg_educ cg_educ, cluster(communityid)
	test T1=T2
	outreg2 using "$tables_sup/`var' - Educ interaction.xls", label dec(3) replace excel addstat("Significance test T1=T2:",r(p))  

	reg `var' T1 T2 1.T1#1.age_q2 1.T2#1.age_q2 age_q2, cluster(communityid)
	test T1=T2
	outreg2 using "$tables_sup/`var' - Age interaction.xls", label dec(3) replace excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.decisiones_q2 1.T2#1.decisiones_q2 decisiones_q2 decisiones_bl_imputed, cluster(communityid)
	test T1=T2
	outreg2 using "$tables_sup/`var' - Decision-making power interaction.xls", label dec(3) replace excel addstat("Significance test T1=T2:",r(p))  

	reg `var' T1 T2 1.T1#1.grit_q2 1.T2#1.grit_q2 grit_q2 grit_bl_imputed, cluster(communityid)
	test T1=T2
	outreg2 using "$tables_sup/`var' - Grit interaction.xls", label dec(3) replace excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.rosenberg_q2 1.T2#1.rosenberg_q2 rosenberg_q2 rosenberg_bl_imputed, cluster(communityid)
	test T1=T2
	outreg2 using "$tables_sup/`var' - Rosenberg interaction.xls", label dec(3) replace excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.wealth_q2 1.T2#1.wealth_q2 wealth_q2 hh_imputed, cluster(communityid)
	test T1=T2
	outreg2 using "$tables_sup/`var' - Wealth interaction.xls", label dec(3) replace excel addstat("Significance test T1=T2:",r(p)) 


	** model 2) + endline child age-sex dummies, leaving out youngest group for each sex, and child order
	reg `var' T1 T2 1.T1#1.cg_educ 1.T2#1.cg_educ cg_educ child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Educ interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.age_q2 1.T2#1.age_q2 age_q2 child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Age interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.decisiones_q2 1.T2#1.decisiones_q2 decisiones_q2 decisiones_bl_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Decision-making power interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.grit_q2 1.T2#1.grit_q2 grit_q2 grit_bl_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Grit interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.rosenberg_q2 1.T2#1.rosenberg_q2 rosenberg_q2 rosenberg_bl_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Rosenberg interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.wealth_q2 1.T2#1.wealth_q2 wealth_q2 hh_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Wealth interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 

	** model 3) + primary caregiver's education in years 
	*			+ primary caregiver's age in years
	*			+ primary caregiver's baseline decision-making power 
	*			+ primary caregiver's baseline  grit score 
	*			+ primary caregiver's baseline Rosenberg (self-esteem) score
	*			+ baseline household wealth 
	*			+ baseline household demographic composition (proportion of household members by age/sex categories)

	reg `var' T1 T2 1.T1#1.cg_educ 1.T2#1.cg_educ cg_educ child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order cg_age decisiones_bl decisiones_bl_imputed grit_bl grit_bl_imputed rosenberg_bl rosenberg_bl_imputed wealth_bl hh_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Educ interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 

	reg `var' T1 T2 1.T1#1.age_q2 1.T2#1.age_q2 age_q2 child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order cg_educ decisiones_bl decisiones_bl_imputed grit_bl grit_bl_imputed rosenberg_bl rosenberg_bl_imputed wealth_bl hh_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Age interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 
	 
	reg `var' T1 T2 1.T1#1.decisiones_q2 1.T2#1.decisiones_q2 decisiones_q2 decisiones_bl_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order cg_age cg_educ grit_bl grit_bl_imputed rosenberg_bl rosenberg_bl_imputed wealth_bl hh_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Decision-making power interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 
	 
	reg `var' T1 T2 1.T1#1.grit_q2 1.T2#1.grit_q2 grit_q2 grit_bl_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order cg_age cg_educ decisiones_bl decisiones_bl_imputed rosenberg_bl rosenberg_bl_imputed wealth_bl hh_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Grit interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 
	 
	reg `var' T1 T2 1.T1#1.rosenberg_q2 1.T2#1.rosenberg_q2 rosenberg_q2 rosenberg_bl_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order cg_age cg_educ decisiones_bl decisiones_bl_imputed grit_bl grit_bl_imputed wealth_bl hh_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Rosenberg interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 
	 
	reg `var' T1 T2 1.T1#1.wealth_q2 1.T2#1.wealth_q2 wealth_q2 hh_imputed child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order cg_age cg_educ decisiones_bl decisiones_bl_imputed grit_bl grit_bl_imputed rosenberg_bl rosenberg_bl_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 
	test T1=T2
	outreg2 using "$tables_sup/`var' - Wealth interaction.xls", label dec(3) append excel addstat("Significance test T1=T2:",r(p)) 
	}

*-----------------------------------------------------------------------------------------------------*
*	TABLE S11. Association between breastfeeding practices and overweight/obesity, for control group  *
*-----------------------------------------------------------------------------------------------------*

keep if tratamiento_hogar==3
	
* breastfeeding practice subscale
	
	** model 1) base model with no controls (unadjusted difference in means)
	reg bmi_2sd practice_index_bf, cluster(communityid) 
	outreg2 using "$tables_sup/tableS11.xls", label dec(3) replace excel

	** model 2) + endline child age-sex dummies, leaving out youngest group for each sex, and child order
	reg bmi_2sd practice_index_bf child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order, cluster(communityid) 
	outreg2 using "$tables_sup/tableS11.xls", label dec(3) append excel

	** model 3) + primary caregiver's education in years 
	*			+ primary caregiver's age in years
	*			+ primary caregiver's baseline decision-making power 
	*			+ primary caregiver's baseline  grit score 
	*			+ primary caregiver's baseline Rosenberg (self-esteem) score
	*			+ baseline household wealth 
	*			+ baseline household demographic composition (proportion of household members by age/sex categories)

		reg bmi_2sd practice_index_bf child_age2_m child_age2_f child_age3_m child_age3_f child_age4_m child_age4_f child_age5_m child_age5_f child_age6_m child_age6_f child_age7_m child_age7_f child_age8_f child_age9_m child_age9_f child_age10_m child_age10_f child_age11_m child_age11_f child_age12_m child_age12_f child_order cg_age cg_educ decisiones_bl decisiones_bl_imputed grit_bl grit_bl_imputed rosenberg_bl rosenberg_bl_imputed wealth_bl hh_imputed prop_males0_5_bl prop_males6_18_bl prop_males19_49_bl prop_males50_bl prop_females0_5_bl prop_females6_18_bl prop_females19_49_bl prop_females50_bl, cluster(communityid) 

	outreg2 using "$tables_sup/tableS11.xls", label dec(3) append excel