/****************************************************************************
*																			*
*					Code for Arouna, Mihcler, and Lokossou (2020)			*
*																			*
****************************************************************************/
cap log close
clear all

cd "C:\Users\jdmichler\Dropbox\Contract Farming\JDE CF Data\"

set pformat %5.4f

log using "CF_results.smcl", replace

*Load data from Aminou
use Benin_panel_20200923, clear

*Set up regression globals
quietly tabulate arrond, generate(_a)

global depvar rsize ryield psold inc_pc
global depvar15 rsize_15 ryield_15 psold_15 inc_pc_15
global arond _a2- _a33
global hhvar15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 agric_15 rice_forma_15 assoc_15
global trt T1C T2C T3C T2T1 T3T1 T3T2

/****************************************************************************
*																			*
*				Table 1: Attributes of existing contracts					*
*																			*
****************************************************************************/

tab cont_know
tab cont_adopt
tab cont_write1
tab cont_price1
tab cont_qual1
tab cont_quant1
tab cont_train1
tab cont_cred1

/****************************************************************************
*																			*
*					Figure 3: Outcomes by treatment group					*
*																			*
****************************************************************************/

by typ_contract, sort: sum rsize ryield psold inc_pc if year == 2016

twoway 	(kdensity rsize if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) ///
		(kdensity rsize if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) ///
		(kdensity rsize if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) ///
		(kdensity rsize if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) ///
		ytitle("Density") xtitle("Rice area (ha)") ///
		xline(.7721364, lpattern(solid) lwidth(thin) lcolor(black)) ///
		text(.18 `=.7721364-.01' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(1.007895, lpattern(solid) lwidth(thin) lcolor(gs10)) ///
		text(.18 `=1.007895-.01' "T1 mean", color(gs10) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(.9035317, lpattern(solid) lwidth(thin) lcolor(sky)) ///
		text(.18 `=.9035317-.01' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(1.030595, lpattern(solid) lwidth(thin) lcolor(turquoise)) ///
		text(.18 `=1.030595+.01' "T3 mean", color(turquoise) j(left) size(vsmall) place(ne) orient(vertical))), ///
		legend(label(1 "Control [C]") label(2 "Price [T1]") ///
		label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3]") ///
		position(6) col(2) row(2)) saving(rsize, replace)
		
		graph export "rsize.emf", as(emf) replace

twoway 	(kdensity ryield if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) ///
		(kdensity ryield if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) ///
		(kdensity ryield if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) ///
		(kdensity ryield if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) ///
		ytitle("Density") xtitle("Yield (kg/ha)") ///
		xline(1652.126, lpattern(solid) lwidth(thin) lcolor(black)) ///
		text(.00009 `=1652.126-10' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(2132.92, lpattern(solid) lwidth(thin) lcolor(gs10)) ///
		text(.00009 `=2132.92-10' "T1 mean", color(gs10) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(2036.185, lpattern(solid) lwidth(thin) lcolor(sky)) ///
		text(.00009 `=2036.185-10' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(2220.142, lpattern(solid) lwidth(thin) lcolor(turquoise)) ///
		text(.00009 `=2220.142-10' "T3 mean", color(turquoise) j(left) size(vsmall) place(nw) orient(vertical))), ///
		legend(label(1 "Control [C]") label(2 "Price [T1]") ///
		label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3] ") ///
		position(6) col(2) row(2)) saving(ryield, replace)
		
		graph export "ryield.emf", as(emf) replace

twoway 	(kdensity psold if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) ///
		(kdensity psold if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) ///
		(kdensity psold if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) ///
		(kdensity psold if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) ///
		ytitle("Density") xtitle("Market participation (%)") ///
		xline(24.95936, lpattern(solid) lwidth(thin) lcolor(black)) ///
		text(.0045 `=24.95936-.5' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(51.13158, lpattern(solid) lwidth(thin) lcolor(gs10)) ///
		text(.0045 `=51.13158-.5' "T1 mean", color(gs10) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(57.25, lpattern(solid) lwidth(thin) lcolor(sky)) ///
		text(.0045 `=57.25-.5' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(67.15613, lpattern(solid) lwidth(thin) lcolor(turquoise)) ///
		text(.0045 `=67.15613-.5' "T3 mean", color(turquoise) j(left) size(vsmall) place(nw) orient(vertical))), ///
		legend(label(1 "Control [C]") label(2 "Price [T1]") ///
		label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3] ") ///
		position(6) col(2) row(2)) saving(psold, replace)
		
		graph export "psold.emf", as(emf) replace		
	
twoway 	(kdensity inc_pc if typ_contract == 0 & year == 2016, lwidth(medium) lcolor(black) lpattern(solid)) ///
		(kdensity inc_pc if typ_contract == 1 & year == 2016, lwidth(medium) lcolor(gs10) lpattern(dash)) ///
		(kdensity inc_pc if typ_contract == 2 & year == 2016, lwidth(medthick) lcolor(sky) lpattern(shortdash_dot)) ///
		(kdensity inc_pc if typ_contract == 3 & year == 2016, lwidth(medium) lcolor(turquoise) lpattern(dash_dot) ///
		ytitle("Density") xtitle("Income per capita (US$)") ///
		xline(265.2537, lpattern(solid) lwidth(thin) lcolor(black)) ///
		text(-.0001 `=265.2537-5' "Control mean", color(black) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(378.8236, lpattern(solid) lwidth(thin) lcolor(gs10)) ///
		text(-.0001 `=378.8236+5' "T1 mean", color(gs10) j(left) size(vsmall) place(ne) orient(vertical)) ///
		xline(370.9411, lpattern(solid) lwidth(thin) lcolor(sky)) ///
		text(-.0001 `=370.9411-5' "T2 mean", color(sky) j(left) size(vsmall) place(nw) orient(vertical)) ///
		xline(458.6367, lpattern(solid) lwidth(thin) lcolor(turquoise)) ///
		text(-.0001 `=458.6367+5' "T3 mean", color(turquoise) j(left) size(vsmall) place(ne) orient(vertical))), ///
		legend(label(1 "Control [C]") label(2 "Price [T1]") ///
		label(3 "Extension & price [T2]") label(4 "Input loans, extension, & price [T3] ") ///
		position(6) col(2) row(2)) saving(inc_pc, replace)
		
		graph export "inc_pc.emf", as(emf) replace
	
gr combine rsize.gph ryield.gph psold.gph inc_pc.gph, col(2) iscale(.5) commonscheme

graph export "outcomes.emf", as(emf) replace

/****************************************************************************
*																			*
*		Table 2, B.1: Baseline characteristics and balance tests			*
*																			*
****************************************************************************/

by treated, sort: eststo: quietly estpost sum $depvar15 $hhvar15

esttab using sum_stat_treated.csv, replace cell(mean(fmt(3)) sd(par fmt(3))) p(3) label nodepvar

by treatment, sort: eststo: quietly estpost sum $depvar15 $hhvar15

esttab using sum_stat_treatment.csv, replace cell(mean(fmt(3)) sd(par fmt(3))) p(3) label nodepvar

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' treated, vce(cl groupement)
} // Compare T with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' T1C, vce(cl groupement)
} // Compare T1 with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' T2C, vce(cl groupement)
} // Compare T2 with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' T3C, vce(cl groupement)
} // Compare T3 with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' T2T1, vce(cl groupement)
} // Compare T2 with T1

foreach dvar of varlist $depvar15 $hhvar15{
eststo: quietly reg `dvar' T3T1, vce(cl groupement)
} // Compare T3 with T1

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' T3T2, vce(cl groupement)
} // Compare T3 with T2

esttab, se nostar drop(_cons) 

matrix C = r(coefs)
eststo clear
local rnames : rownames C
local models : coleq C
local models : list uniq models
local i 0

foreach name of local rnames {
   local ++i
   local j 0
   capture matrix drop b
   capture matrix drop se
   foreach model of local models {
     local ++j
     matrix tmp = C[`i', 2*`j'-1]
     if tmp[1,1]<. {
       matrix colnames tmp = `model'
       matrix b = nullmat(b), tmp
       matrix tmp[1,1] = C[`i', 2*`j']
       matrix se = nullmat(se), tmp
     }
   }
   ereturn post b
   quietly estadd matrix se
   eststo `name'
 }

esttab using balance_new.csv, cells(b(star fmt(3)) se(par fmt(3))) ///
	noobs label replace
	
/****************************************************************************
*																			*
*			Table 3: Baseline characteristics and balance tests				*
*																			*
****************************************************************************/

foreach dvar of varlist $trt {
eststo: reg `dvar' $depvar15 $hhvar15 $arond, vce(cl groupement)
	test $depvar $hhvar
} // Compare across all treatments


esttab, se nostar drop(_cons) 

esttab using balance_ftest.csv, cells(b(star fmt(3)) se(par fmt(3))) ///
	noobs label replace

/****************************************************************************
*																			*
*					Table 4: Tests for attrition bias						*
*																			*
****************************************************************************/

preserve

gen attrited = 1 if typ_contract == 4
replace attrited = 0 if attrited == .

by attrited, sort: eststo: quietly estpost sum treated $depvar15 $hhvar15

esttab using sum_stat_attrit.csv, replace cell(mean(fmt(3)) sd(par fmt(3))) p(3) label nodepvar

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' attrited if treated == 0 | treated == 1, vce(cl groupement)
} // Compare T with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' attrited if treatment == 0 | treatment == 1, vce(cl groupement)
} // Compare T1 with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' attrited if treatment == 0 | treatment == 2, vce(cl groupement)
} // Compare T2 with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' attrited if treatment == 0 | treatment == 3, vce(cl groupement)
} // Compare T3 with C

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' attrited if treatment == 1 | treatment == 2, vce(cl groupement)
} // Compare T2 with T1

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' attrited if treatment == 1 | treatment == 3, vce(cl groupement)
} // Compare T3 with T1

foreach dvar of varlist $depvar15 $hhvar15 {
eststo: quietly reg `dvar' attrited if treatment == 2 | treatment == 3, vce(cl groupement)
} // Compare T3 with T2

esttab, se nostar drop(_cons) 

matrix C = r(coefs)
eststo clear
local rnames : rownames C
local models : coleq C
local models : list uniq models
local i 0

foreach name of local rnames {
   local ++i
   local j 0
   capture matrix drop b
   capture matrix drop se
   foreach model of local models {
     local ++j
     matrix tmp = C[`i', 2*`j'-1]
     if tmp[1,1]<. {
       matrix colnames tmp = `model'
       matrix b = nullmat(b), tmp
       matrix tmp[1,1] = C[`i', 2*`j']
       matrix se = nullmat(se), tmp
     }
   }
   ereturn post b
   quietly estadd matrix se
   eststo `name'
 }

esttab using balance_attrit.csv, cells(b(star fmt(3)) se(par fmt(3))) ///
	noobs label replace

* attrition on treatment
reg attrited treatment, vce(cl groupement)
reg attrited treatment $arond, vce(cl groupement)
reg attrited treatment $arond $hhvar15, vce(cl groupement)
	
restore

/****************************************************************************
*																			*
*				Table B.2: Coefficient of variation in outcomes				*
*																			*
****************************************************************************/

by typ_contract, sort: cv2 $depvar

/****************************************************************************
*																			*
*		Table 5, B.3, B.4, & B.5: Treatment effect of contract				*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4

/* Panel A: Rice area */
/* OLS */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_1, replace): ///
	 reg rsize treated $arond, vce(cl groupement) 
		outreg2 using "CTa_jde.doc", ///
		title (Table 5a: Treatment effects on Rice area) ctitle (OLS) ///
		drop(_a*) se dec(3) label replace nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No)
	
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_2, replace): ///
	 reg rsize treated $hhvar15 $arond, vce(cl groupement) 
		outreg2 using "CTa_jde.doc", append ctitle (OLS) ///
		drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
/* ANCOVA */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_3, replace): ///
	 reg rsize treated rsize_15 $arond, vce(cl groupement) 
		outreg2 using "CTa_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* rsize_15) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No) 

ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_rsize_4, replace): ///
	 reg rsize treated rsize_15 $hhvar15 $arond, vce(cl groupement) 
		outreg2 using "CTa_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: yield */
/* OLS */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_1, replace): ///
	reg ryield treated $arond, vce(cl groupement)
		outreg2 using "CTb_jde.doc", ///
		title (Table 5b: Treatment effects on yield) ctitle (OLS) ///
		drop(_a*) se dec(3) label replace nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No)
	
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_2, replace): ///
	reg ryield treated $hhvar15 $arond, vce(cl groupement)
		outreg2 using "CTb_jde.doc", append ctitle (OLS) ///
		drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_3, replace): ///
	reg ryield treated ryield_15 $arond, vce(cl groupement)
		outreg2 using "CTb_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* ryield_15) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No) 

ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_ryield_4, replace): ///
	reg ryield treated ryield_15 $hhvar15 $arond, vce(cl groupement)
		outreg2 using "CTb_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Market participation */
/* OLS */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_1, replace): ///
	reg psold treated $arond, vce(cl groupement)
		outreg2 using "CTc_jde.doc", ///
		title (Table 5c: Treatment effects on market participation) ctitle (OLS) ///
		drop(_a*) se dec(3) label replace nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No)
	
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_2, replace): ///
	reg psold treated $hhvar15 $arond, vce(cl groupement)
		outreg2 using "CTc_jde.doc", append ctitle (OLS) ///
		drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_3, replace): ///
	reg psold treated psold_15 $arond, vce(cl groupement)
		outreg2 using "CTc_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* psold_15) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No) 

ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_psold_4, replace): ///
	reg psold treated psold_15 $hhvar15 $arond, vce(cl groupement)
		outreg2 using "CTc_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
/* Panel D: Income */
/* OLS */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_1, replace): ///
	reg inc_pc treated inc_pc_15 $arond, vce(cl groupement)
		outreg2 using "CTd_jde.doc", ///
		title (Table 5d: Treatment effects on income) ctitle (OLS) ///
		drop(_a*) se dec(3) label replace nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No) 

ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_2, replace): ///
	reg inc_pc treated $hhvar15 $arond, vce(cl groupement)
		outreg2 using "CTd_jde.doc", append ctitle (OLS) ///
		drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_3, replace): ///
	reg inc_pc treated inc_pc_15 $arond, vce(cl groupement)
		outreg2 using "CTd_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* inc_pc_15) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, No) 

ritest treated _b[treated], cluster(groupement) strata(depart arrond) reps(1000) seed(0) nodots saving(ri_inc_4, replace): ///
	reg inc_pc treated inc_pc_15 $hhvar15 $arond, vce(cl groupement)
		outreg2 using "CTd_jde.doc", append ctitle (ANCOVA) ///
		drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
		rice_forma assoc) se dec(3) label nocons ///
		addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

*Mean dependent variable in control group
sum rsize ryield psold inc_pc if treated == 0 & year == 2016

/* MHT */
mhtexp rsize ryield psold inc_pc, treatment(treated)

restore

/* Lee Bounds */
eststo: quietly leebounds rsize treated
eststo: quietly leebounds ryield treated
eststo: quietly leebounds psold treated
eststo: quietly leebounds inc_pc treated

esttab using Lee1.csv, cells(b(star fmt(3)) se(par fmt(3))) ///
	noobs label replace

/****************************************************************************
*																			*
* Table 6, 7, B.6, & B.7: Treatment effects of each contract characteristic	*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4

/* Panel A: Rice area */
/* OLS */
reg rsize i.typ_contract $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPa_jde.doc", ///
	title (Table 6a: Treatment effects on Rice area) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg rsize i.typ_contract  $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPa_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
/* ANCOVA */
reg rsize i.typ_contract rsize_15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPa_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg rsize i.typ_contract rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPa_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)

/* Panel B: yield */
/* OLS */
reg ryield i.typ_contract $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPb_jde.doc", ///
	title (Table 6b: Treatment effects on yield) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg ryield i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPb_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
/* ANCOVA */
reg ryield i.typ_contract ryield_15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPb_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg ryield i.typ_contract ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPb_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)

/* Panel C: Market participation */
/* OLS */
reg psold i.typ_contract $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPc_jde.doc", ///
	title (Table 6c: Treatment effects on market participation) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg psold i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPc_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
/* ANCOVA */
reg psold i.typ_contract psold_15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPc_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg psold i.typ_contract psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPc_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)

/* Panel D: Income */
/* OLS */
reg inc_pc i.typ_contract $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPd_jde.doc", ///
	title (Table 6d: Treatment effects on income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg inc_pc i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPd_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
/* ANCOVA */
reg inc_pc i.typ_contract inc_pc_15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPd_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* inc_pc_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg inc_pc i.typ_contract inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPd_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	  
test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
/* MHT */
mhtexp rsize ryield psold inc_pc, treatment(typ_contract)

/* Randomization inference */
/* Panel A: Rice area */
/* OLS */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_11, replace): ///
	reg rsize i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_rsize_12, replace): ///
	reg rsize i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_13, replace): ///
	reg rsize i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_21, replace): ///
	reg rsize i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_rsize_22, replace): ///
	reg rsize i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_23, replace): ///
	reg rsize i.treatment $hhvar15 $arond, vce(cl groupement)

/* ANCOVA */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_31, replace): ///
	reg rsize i.treatment rsize_15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_rsize_32, replace): ///
	reg rsize i.treatment rsize_15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_33, replace): ///
	reg rsize i.treatment rsize_15 $arond, vce(cl groupement)
	
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_rsize_41, replace): ///
	reg rsize i.treatment rsize_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3)  saving(ri_rsize_42, replace): ///
	reg rsize i.treatment rsize_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_rsize_43, replace): ///
	reg rsize i.treatment rsize_15 $hhvar15 $arond, vce(cl groupement)

/* Panel B: yield */
/* OLS */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_11, replace): ///
	reg ryield i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_12, replace): ///
	reg ryield i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_13, replace): ///
	reg ryield i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_21, replace): ///
	reg ryield i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_22, replace): ///
	reg ryield i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_23, replace): ///
	reg ryield i.treatment $hhvar15 $arond, vce(cl groupement)

/* ANCOVA */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_31, replace): ///
	reg ryield i.treatment ryield_15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_32, replace): ///
	reg ryield i.treatment ryield_15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_33, replace): ///
	reg ryield i.treatment ryield_15 $arond, vce(cl groupement)
	
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_ryield_41, replace): ///
	reg ryield i.treatment ryield_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_ryield_42, replace): ///
	reg ryield i.treatment ryield_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_ryield_43, replace): ///
	reg ryield i.treatment ryield_15 $hhvar15 $arond, vce(cl groupement)

/* Panel C: Market participation */
/* OLS */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_11, replace): ///
	reg psold i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_12, replace): ///
	reg psold i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_13, replace): ///
	reg psold i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_21, replace): ///
	reg psold i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_22, replace): ///
	reg psold i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_23, replace): ///
	reg psold i.treatment $hhvar15 $arond, vce(cl groupement)

/* ANCOVA */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_31, replace): ///
	reg psold i.treatment psold_15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_32, replace): ///
	reg psold i.treatment psold_15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_33, replace): ///
	reg psold i.treatment psold_15 $arond, vce(cl groupement)
	
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_psold_41, replace): ///
	reg psold i.treatment psold_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(5000) seed(0) nodots fix(1 3) saving(ri_psold_42, replace): ///
	reg psold i.treatment psold_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_psold_43, replace): ///
	reg psold i.treatment psold_15 $hhvar15 $arond, vce(cl groupement)

/* Panel D: Income */
/* OLS */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_11, replace): ///
	reg inc_pc i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_12, replace): ///
	reg inc_pc i.treatment $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_13, replace): ///
	reg inc_pc i.treatment $arond, vce(cl groupement)

ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_21, replace): ///
	reg inc_pc i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_22, replace): ///
	reg inc_pc i.treatment $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_23, replace): ///
	reg inc_pc i.treatment $hhvar15 $arond, vce(cl groupement)

/* ANCOVA */
ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_31, replace): ///
	reg inc_pc i.treatment inc_pc_15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_32, replace): ///
	reg inc_pc i.treatment inc_pc_15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_33, replace): ///
	reg inc_pc i.treatment inc_pc_15 $arond, vce(cl groupement)

ritest treatment _b[1.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(2 3) saving(ri_inc_41, replace): ///
	reg inc_pc i.treatment inc_pc_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[2.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 3) saving(ri_inc_42, replace): ///
	reg inc_pc i.treatment inc_pc_15 $hhvar15 $arond, vce(cl groupement)
	
ritest treatment _b[3.treatment], cluster(groupement) strata(depart arrond) ///
	reps(1000) seed(0) nodots fix(1 2) saving(ri_inc_43, replace): ///
	reg inc_pc i.treatment inc_pc_15 $hhvar15 $arond, vce(cl groupement)

restore

/****************************************************************************
*																			*
*			Table 8: Heterogeneity of treatment effects of contract			*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4

/* Contract Indicator */
/* Column 1: Rice area */
reg rsize i.treated##c.hhsize_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ///
	title (Table 8: Treatment effects on Rice area) ctitle (hhsize) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg rsize i.treated##c.age_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ctitle (age) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg rsize i.treated##i.sex_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ctitle (sex) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg rsize i.treated##c.rice_exp_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ctitle (exp) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg rsize i.treated##i.primeduc_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ctitle (edu) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg rsize i.treated##i.agric_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ctitle (activity) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg rsize i.treated##i.rice_forma_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ctitle (train) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg rsize i.treated##i.assoc_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_area_jde.doc", ctitle (assoc) ///
	drop(_a* rsize_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
/* Column 2: yield */
reg ryield i.treated##c.hhsize_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ///
	title (Table 8: Treatment effects on Rice yield) ctitle (hhsize) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg ryield i.treated##c.age_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ctitle (age) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg ryield i.treated##i.sex_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ctitle (sex) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg ryield i.treated##c.rice_exp_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ctitle (exp) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg ryield i.treated##i.primeduc_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ctitle (edu) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg ryield i.treated##i.agric_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ctitle (activity) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg ryield i.treated##i.rice_forma_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ctitle (train) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg ryield i.treated##i.assoc_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_yield_jde.doc", ctitle (assoc) ///
	drop(_a* ryield_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
/* Column 3: Market participation */
reg psold i.treated##c.hhsize_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ///
	title (Table 8: Treatment effects on Market participation) ctitle (hhsize) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg psold i.treated##c.age_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ctitle (age) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg psold i.treated##i.sex_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ctitle (sex) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg psold i.treated##c.rice_exp_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ctitle (exp) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg psold i.treated##i.primeduc_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ctitle (edu) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg psold i.treated##i.agric_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ctitle (activity) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg psold i.treated##i.rice_forma_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ctitle (train) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg psold i.treated##i.assoc_15 psold_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_market_jde.doc", ctitle (assoc) ///
	drop(_a* psold_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
/* Column 4: Income */
reg inc_pc i.treated##c.hhsize_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)	
	outreg2 using "h_inc_jde.doc", ///
	title (Table 8: Treatment effects on Income) ctitle (hhsize) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg inc_pc i.treated##c.age_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)	
	outreg2 using "h_inc_jde.doc", ctitle (age) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg inc_pc i.treated##i.sex_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_inc_jde.doc", ctitle (sex) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg inc_pc i.treated##c.rice_exp_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)	
	outreg2 using "h_inc_jde.doc", ctitle (exp) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)
	
reg inc_pc i.treated##i.primeduc_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)
	outreg2 using "h_inc_jde.doc", ctitle (edu) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg inc_pc i.treated##i.agric_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)	
	outreg2 using "h_inc_jde.doc", ctitle (activity) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg inc_pc i.treated##i.rice_forma_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)	
	outreg2 using "h_inc_jde.doc", ctitle (train) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg inc_pc i.treated##i.assoc_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)	
	outreg2 using "h_inc_jde.doc", ctitle (assoc) ///
	drop(_a* inc_pc_15 hhsize_15 age_15 sex_15 rice_exp_15 primeduc_15 ///
	agric_15 rice_forma_15 assoc_15) se dec(3) label append nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

/****************************************************************************
*																			*
*		Figures 4-6: Heterogeneity of treatment effects of contract			*
*																			*
****************************************************************************/

/* Contract Types */
/* Column 1: Rice area */
qui reg rsize i.typ_contract##c.hhsize_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(hhsize_15=(0(2)20))
marginsplot, title("") ytitle("Rice area (ha)") xtitle("Household size (n)") ///
	legend(pos(6) col(2) row(2)) saving(rsize_1, replace)
graph export "rsize_1.emf", as(emf) replace

qui reg rsize i.typ_contract##c.rice_exp_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(rice_exp_15=(0(5)30))
marginsplot, title("") ytitle("Rice area (ha)") xtitle("Exp. producing rice (years)") ///
	legend(pos(6) col(2) row(2)) saving(rsize_4, replace)
	graph export "rsize_4.emf", as(emf) replace

qui reg rsize i.typ_contract##i.rice_forma_15 rsize_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, over(rice_forma_15)
marginsplot, title("") ytitle("Rice area (ha)") xtitle("Farming contracts") ///
	xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(rsize_7, replace)
	graph export "rsize_7.emf", as(emf) replace

/* Column 2: yield */
qui reg ryield i.typ_contract##c.hhsize_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(hhsize_15=(0(2)20))
marginsplot, title("") ytitle("Yield (kg/ha)") xtitle("Household size (n)") ///
	legend(pos(6) col(2) row(2)) saving(ryield_1, replace)
	graph export "ryield_1.emf", as(emf) replace

qui reg ryield i.typ_contract##c.rice_exp_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(rice_exp_15=(0(5)30))
marginsplot, title("") ytitle("Yield (kg/ha)") xtitle("Exp. producing rice (years)") ///
	legend(pos(6) col(2) row(2)) saving(ryield_4, replace)
	graph export "ryield_4.emf", as(emf) replace

qui reg ryield i.typ_contract##i.rice_forma_15 ryield_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, over(rice_forma_15)
marginsplot, title("") ytitle("Yield (kg/ha)") xtitle("Farming contracts") ///
	xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(ryield_7, replace)
	graph export "ryield_7.emf", as(emf) replace

/* Column 3: market participation */
qui reg psold i.typ_contract##c.hhsize_15 psold_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(hhsize_15=(0(2)20))
marginsplot, title("") ytitle("Market participation (%)") xtitle("Household size (n)") ///
	legend(pos(6) col(2) row(2)) saving(psold_1, replace)
	graph export "psold_1.emf", as(emf) replace

qui reg psold i.typ_contract##c.rice_exp_15 psold_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(rice_exp_15=(0(5)30))
marginsplot, title("") ytitle("Market participation (%)") xtitle("Exp. producing rice (years)") ///
	legend(pos(6) col(2) row(2)) saving(psold_4, replace)
	graph export "psold_4.emf", as(emf) replace

qui reg psold i.typ_contract##i.rice_forma_15 psold_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, over(rice_forma_15)
marginsplot, title("") ytitle("Market participation (%)") xtitle("Farming contracts") ///
	xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(psold_7, replace)
	graph export "psold_7.emf", as(emf) replace

/* Column 4: income */
qui reg inc_pc i.typ_contract##c.hhsize_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(hhsize_15=(0(2)20))
marginsplot, title("") ytitle("Income per capita (US$)") xtitle("Household size (n)") ///
	legend(pos(6) col(2)) saving(inc_pc_1, replace)
	graph export "inc_pc_1.emf", as(emf) replace

qui reg inc_pc i.typ_contract##c.rice_exp_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, at(rice_exp_15=(0(5)30))
marginsplot, title("") ytitle("Income per capita (US$)") xtitle("Exp. producing rice (years)") ///
	legend(pos(6) col(2)) saving(inc_pc_4, replace)
	graph export "inc_pc_4.emf", as(emf) replace

qui reg inc_pc i.typ_contract##i.rice_forma_15 inc_pc_15 $hhvar15 $arond, vce(cl groupement)
qui margins typ_contract, over(rice_forma_15)
marginsplot, title("") ytitle("Income per capita (US$)") xtitle("Farming contracts") ///
	xlabel(0 "[C]" 1 "[T1]" 2 "[T2]" 3 "[T3]") legend(pos(6) col(2)) saving(inc_pc_7, replace)
	graph export "inc_pc_7.emf", as(emf) replace

*Combine selected graphs
gr combine rsize_1.gph ryield_1.gph psold_1.gph inc_pc_1.gph, col(2) iscale(.5) commonscheme
	graph export "hetero_hhsize.emf", as(emf) replace

gr combine rsize_4.gph ryield_4.gph psold_4.gph inc_pc_4.gph, col(2) iscale(.5) commonscheme
	graph export "hetero_exp.emf", as(emf) replace

gr combine rsize_7.gph ryield_7.gph psold_7.gph inc_pc_7.gph, col(2) iscale(.5) commonscheme
	graph export "hetero_train.emf", as(emf) replace

restore

/****************************************************************************
*																			*
*			Table B.3-B.3: Treatment effect on input use					*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4

/* Table B.3: Treated */
reg lns treated lns15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTin_jde.doc", ///
	title (Table B.3: Treatment effects on Inputs) ctitle (ln(seed)) ///
	drop(_a* lns15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

reg lnf treated lnf15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTin_jde.doc", append ctitle (ln(fert)) ///
	drop(_a* lnf15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 	

reg lnp treated lnp15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTin_jde.doc", append ctitle (ln(pest)) ///
	drop(_a* lnp15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

reg lnh treated lnh15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTin_jde.doc", append ctitle (ln(herb)) ///
	drop(_a* lnh15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 	

reg lnl treated lnl15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTin_jde.doc", append ctitle (ln(labr)) ///
	drop(_a* lnl15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 		 		

*Mean dependent variable in control group
sum lns lnf lnp lnh lnl if treated == 0 & year == 2016
	
/* Table B.4: Treatment Types */
reg lns i.typ_contract lns15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPin_jde.doc", ///
	title (Table B.4: Treatment effects on Inputs) ctitle (ln(seed)) ///
	drop(_a* lns15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)
	
reg lnf i.typ_contract lnf15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPin_jde.doc", append ctitle (ln(fert)) ///
	drop(_a* lnf15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 	

test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)

reg lnp i.typ_contract lnp15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPin_jde.doc", append ctitle (ln(pest)) ///
	drop(_a* lnp15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)

reg lnh i.typ_contract lnh15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPin_jde.doc", append ctitle (ln(herb)) ///
	drop(_a* lnh15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)

test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)	

reg lnl i.typ_contract lnl15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "iTYPin_jde.doc", append ctitle (ln(labr)) ///
	drop(_a* lnl15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 	

test (1.typ_contract=2.typ_contract) (2.typ_contract=3.typ_contract) ///
	(1.typ_contract=3.typ_contract), common cons mtest(b)	 		

restore
	
/****************************************************************************
*																			*
*		Table B.10: Treatment effect of price and quantity guarantee			*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
drop if T1C == .

/* Panel A: Rice area */
/* OLS */
reg rsize T1C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1a_jde.doc", ///
	title (Table B.10a: Treatment effects on Rice area) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
reg rsize T1C $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CT1a_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg rsize T1C rsize_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg rsize T1C rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 


/* Panel B: yield */
/* OLS */
reg ryield T1C $arond if year == 2016, vce(cl groupement)
	outreg2 using "CT1b_jde.doc", ///
	title (Table B.10b: Treatment effects on yield) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg ryield T1C $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CT1b_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg ryield T1C ryield_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
reg ryield T1C ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 


/* Panel C: Market participation */
/* OLS */
reg psold T1C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1c_jde.doc", ///
	title (Table B.10c: Treatment effects on market participation) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T1C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1c_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg psold T1C psold_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T1C psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel D: Income */
/* OLS */
reg inc_pc T1C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1d_jde.doc", ///
	title (Table B.10d: Treatment effects on income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No, ///
	Household FE, No, Year FE, No) 
	
reg inc_pc T1C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1d_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc)  se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg inc_pc T1C inc_pc_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T1C inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CT1d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc)  se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

restore

/****************************************************************************
*																			*
*		Table B.11: Treatment effect of price and extension training			*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
drop if T2C == .

/* Panel A: Rice area */
/* OLS */
reg rsize T2C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Ca_jde.doc", ///
	title (Table B.11a: Treatment effects on Rice area) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
reg rsize T2C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Ca_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg rsize T2C rsize_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Ca_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg rsize T2C rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Ca_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: yield */
/* OLS */
reg ryield T2C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cb_jde.doc", ///
	title (Table B.11b: Treatment effects on yield) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T2C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cb_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg ryield T2C ryield_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cb_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T2C ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cb_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Market participation */
/* OLS */
reg psold T2C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cc_jde.doc", ///
	title (Table B.11c: Treatment effects on market participation) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T2C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cc_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg psold T2C psold_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cc_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T2C psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cc_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel D: Income */
/* OLS */
reg inc_pc T2C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cd_jde.doc", ///
	title (Table B.11d: Treatment effects on income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T2C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cd_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc)  se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg inc_pc T2C inc_pc_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cd_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T2C inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2Cd_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc)  se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

restore

/****************************************************************************
*																			*
* 	Table B.12: Treatment effect of price, extension training, and loans	*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
drop if T3C == .

/* Panel A: Rice area */
/* OLS */
reg rsize T3C $arond if year == 2016, vce(cl groupement)
	outreg2 using "T3Ca_jde.doc", ///
	title (Table B.12a: Treatment effects on Rice area) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
reg rsize T3C $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "T3Ca_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg rsize T3C rsize_15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "T3Ca_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg rsize T3C rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "T3Ca_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: yield */
/* OLS */
reg ryield T3C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cb_jde.doc", ///
	title (Table B.12b: Treatment effects on yield) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T3C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cb_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg ryield T3C ryield_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cb_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T3C ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cb_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes)  

/* Panel C: Market participation */
/* OLS */
reg psold T3C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cc_jde.doc", ///
	title (Table B.12c: Treatment effects on market participation) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T3C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cc_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg psold T3C psold_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cc_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T3C psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cc_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
/* Panel D: Income */
/* OLS */
reg inc_pc T3C $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cd_jde.doc", ///
	title (Table B.12d: Treatment effects on income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T3C $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cd_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc)  se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg inc_pc T3C inc_pc_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cd_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T3C inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T3Cd_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc)  se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

restore
	
/****************************************************************************
*																			*
*			Table B.13: Treatment effects of extension training				*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
drop if T2T1 == .

/* Panel A: Rice area */
/* OLS */
reg rsize T2T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2a_jde.doc", ///
	title (Table B.13a: Treatment effects on Rice area) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
reg rsize T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2a_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg rsize T2T1 rsize_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg rsize T2T1 rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: yield */
/* OLS */
reg ryield T2T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2b_jde.doc", ///
	title (Table B.13b: Treatment effects on yield) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2b_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg ryield T2T1 ryield_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T2T1 ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Market participation */
/* OLS */
reg psold T2T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2c_jde.doc", ///
	title (Table B.13c: Treatment effects on market participation) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2c_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg psold T2T1 psold_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T2T1 psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel D: Income */
/* OLS */
reg inc_pc T2T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2d_jde.doc", ///
	title (Table B.13d: Treatment effects on income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T2T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2d_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg inc_pc T2T1 inc_pc_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T2T1 inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T2d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

sum rsize ryield psold inc_pc if treatment == 1 & year == 2016

restore
	
/****************************************************************************
*																			*
*				Table B.14: Treatment effects of input loans				*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
drop if T3T2 == .

/* Panel A: Rice area */
/* OLS */
reg rsize T3T2 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3a_jde.doc", ///
	title (Table B.14a: Treatment effects on Rice area) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
reg rsize T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3a_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg rsize T3T2 rsize_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg rsize T3T2 rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: Yield */
/* OLS */
reg ryield T3T2 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3b_jde.doc", ///
	title (Table B.14b: Treatment effects on yield) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3b_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc ) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg ryield T3T2 ryield_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T3T2 ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Market participation */
/* OLS */
reg psold T3T2 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3c_jde.doc", ///
	title (Table B.14c: Treatment effects on market participation) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3c_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg psold T3T2 psold_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T3T2 psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel D: Income */
/* OLS */
reg inc_pc T3T2 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3d_jde.doc", ///
	title (Table B.14d: Treatment effects on income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T3T2 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3d_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg inc_pc T3T2 inc_pc_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T3T2 inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T2T3d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

sum rsize ryield psold inc_pc if treatment == 2 & year == 2016

restore

/****************************************************************************
*																			*
*			Table B.15: Treatment effects of extension and loans			*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
drop if T3T1 == .

/* Panel A: Rice area */
/* OLS */
reg rsize T3T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3a_jde.doc", ///
	title (Table B.15a: Treatment effects on Rice area) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No)
	
reg rsize T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3a_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg rsize T3T1 rsize_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg rsize T3T1 rsize_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3a_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* rsize_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: Yield */
/* OLS */
reg ryield T3T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3b_jde.doc", ///
	title (Table B.15b: Treatment effects on yield) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3b_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg ryield T3T1 ryield_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg ryield T3T1 ryield_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3b_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* ryield_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Market participation */
/* OLS */
reg psold T3T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3c_jde.doc", ///
	title (Table B.15c: Treatment effects on market participation) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3c_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg psold T3T1 psold_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg psold T3T1 psold_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3c_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* psold_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel D: Income */
/* OLS */
reg inc_pc T3T1 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3d_jde.doc", ///
	title (Table B.15d: Treatment effects on income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T3T1 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3d_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg inc_pc T3T1 inc_pc_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg inc_pc T3T1 inc_pc_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "T1T3d_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* inc_pc_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

restore

/****************************************************************************
*																			*
*				Table B.16: Treatment effect on income						*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4

/* Panel A: Rice Income */
/* OLS */
reg lnri treated $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTri_jde.doc", ///
	title (Table B.16a: Treatment effects on Rice Income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnri treated $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTri_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnri treated lnri15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTri_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnri15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 	
	
reg lnri treated lnri15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTri_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnri15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
/* Panel B: Other Farm Income */
/* OLS */
reg lnfi treated $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTfi_jde.doc", ///
	title (Table B.16b: Treatment effects on Other Farm Income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnfi treated $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTfi_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnfi treated lnfi15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTfi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnfi15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnfi treated lnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTfi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Non-Farm Income */
/* OLS */
reg lnnfi treated $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTnfi_jde.doc", ///
	title (Table B.16c: Treatment effects on Non-Farm Income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnnfi treated $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTnfi_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnnfi treated lnnfi15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTnfi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnnfi15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 	

	reg lnnfi treated lnnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTnfi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
*Mean dependent variable in control group
sum lnri lnfi lnnfi if treated == 0 & year == 2016

restore

/****************************************************************************
*																			*
*		Table B.17: Treatment effect of each contract on income				*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
	
/* Panel A: Rice Income */
/* OLS */
reg lnri i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPri_jde.doc", ///
	title (Table B.17a: Treatment effects on Rice Income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnri i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPri_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnri i.typ_contract lnri15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPri_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* lnri15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg lnri i.typ_contract lnri15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPri_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnri15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: Other Farm Income */
/* OLS */
reg lnfi i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPfi_jde.doc", ///
	title (Table B.17b: Treatment effects on Other Farm Income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnfi i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPfi_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnfi i.typ_contract lnfi15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPfi_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* lnfi15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg lnfi i.typ_contract lnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPfi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Non-Farm Income */
/* OLS */
reg lnnfi i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPnfi_jde.doc", ///
	title (Table B.17c: Treatment effects on Non-Farm Income) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnnfi i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPnfi_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnnfi i.typ_contract lnnfi15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPnfi_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* lnnfi15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg lnnfi i.typ_contract lnnfi15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPnfi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnnfi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

restore

/****************************************************************************
*																			*
*				Table B.18: Treatment effect on food security				*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4

/* Panel A: Food Insecurity */
/* OLS */
reg hfias treated $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTe_jde.doc", ///
	title (Table B.18a: Treatment effects on HFIAS) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg hfias treated $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTe_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg hfias treated hfias_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTe_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* hfias_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg hfias treated hfias_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTe_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* hfias_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: Food Security */
/* OLS */
reg fcs treated $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTf_jde.doc", ///
	title (Table B.18b: Treatment effects on FCS) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg fcs treated $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTf_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg fcs treated fcs_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTf_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* fcs_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg fcs treated fcs_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTf_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* fcs_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 	

*Mean dependent variable in control group
sum hfias fcs if treated == 0 & year == 2016

restore

/****************************************************************************
*																			*
*	Table B.19: Treatment effect of each contract on food security			*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
	
/* Panel A: Food Insecurity */
/* OLS */
reg hfias i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPe_jde.doc", ///
	title (Table B.19a: Treatment effects on HFIAS) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg hfias i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPe_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg hfias i.typ_contract hfias_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPe_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* hfias_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg hfias i.typ_contract hfias_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPe_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* hfias_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: Food Security */
/* OLS */
reg fcs i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPf_jde.doc", ///
	title (Table B.19b: Treatment effects on HFIAS) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg fcs i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPf_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg fcs i.typ_contract fcs_15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPf_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* fcs_15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg fcs i.typ_contract fcs_15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPf_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* fcs_15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

restore

/****************************************************************************
*																			*
*				Table B.20: Treatment effect on profit						*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4

/* Panel A: Profits - low wages */
/* OLS */
reg lnpil treated $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpil_jde.doc", ///
	title (Table B.20a: Treatment effects on Profits - low wages) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnpil treated $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpil_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnpil treated lnpil15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpil_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpil15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 	
	
reg lnpil treated lnpil15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpil_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpil15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 
	
/* Panel B: Profits - high wages */
/* OLS */
reg lnpih treated $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpih_jde.doc", ///
	title (Table B.20b: Treatment effects on Profits - high wages) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnpih treated $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpih_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnpih treated lnpih15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpih_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpih15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnpih treated lnpih15 $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpih_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpih15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Profits - self-reported wages */
/* OLS */
reg lnpi treated $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpi_jde.doc", ///
	title (Table B.20c: Treatment effects on Profits - self-reported wages) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnpi treated $hhvar15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpi_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnpi treated lnpi15 $arond if year == 2016, vce(cl groupement)
	outreg2 using "CTpi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpi15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 	

	reg lnpi treated lnpi15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "CTpi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

*Mean dependent variable in control group
	sum lnpil lnpih lnpi if treated == 0 & year == 2016

restore

/****************************************************************************
*																			*
*		Table B.21: Treatment effect of each contract on profit				*
*																			*
****************************************************************************/

*Drop attritors
preserve
drop if typ_contract == 4
	
/* Panel A: Profits - low wages */
/* OLS */
reg lnpil i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPil_jde.doc", ///
	title (Table B.21a: Treatment effects on Profits - low wages) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnpil i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPil_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnpil i.typ_contract lnpil15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPil_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* lnpil15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg lnpil i.typ_contract lnpil15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPil_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpil15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel B: Profits - high wages */
/* OLS */
reg lnpih i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPih_jde.doc", ///
	title (Table B.21b: Treatment effects on Profits - high wages) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnpih i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPih_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnpih i.typ_contract lnpih15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPih_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* lnpih15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg lnpih i.typ_contract lnpih15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPih_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpih15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* Panel C: Profits - self-reported wages */
/* OLS */
reg lnpi i.typ_contract $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPi_jde.doc", ///
	title (Table B.21c: Treatment effects on Profits - self-reported wages) ctitle (OLS) ///
	drop(_a*) se dec(3) label replace nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 
	
reg lnpi i.typ_contract $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPi_jde.doc", append ctitle (OLS) ///
	drop(_a* hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

/* ANCOVA */
reg lnpi i.typ_contract lnpi15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPi_jde.doc", append ctitle (ANCOVA) ///
	drop( _a* lnpi15) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, No) 

reg lnpi i.typ_contract lnpi15 $hhvar15 $arond if year == 2016, vce(cl groupement)	
	outreg2 using "iTYPi_jde.doc", append ctitle (ANCOVA) ///
	drop(_a* lnpi15 hhsize* age* sex* rice_exp* primeduc* agric* rice_forma* assoc* ///
	rice_forma assoc) se dec(3) label nocons ///
	addtext(Arrondissement FE, Yes, Household Covariates, Yes) 

restore


/* END */